← Back to branch summary

~ubuntu-branches/ubuntu/trusty/nwchem/trusty-proposed

« back to all changes in this revision

Viewing changes to src/tools/ga-5-1/global/testing/perform.F

  • Committer: Package Import Robot
  • Author(s): Michael Banck, Daniel Leidert, Andreas Tille, Michael Banck
  • Date: 2013-07-04 12:14:55 UTC
  • mfrom: (1.1.2)
  • Revision ID: package-import@ubuntu.com-20130704121455-5tvsx2qabor3nrui
Tags: 6.3-1
http://bugs.debian.org/696361
* New upstream release.
* Fixes anisotropic properties (Closes: #696361).
* New features include:
  + Multi-reference coupled cluster (MRCC) approaches
  + Hybrid DFT calculations with short-range HF 
  + New density-functionals including Minnesota (M08, M11) and HSE hybrid
    functionals
  + X-ray absorption spectroscopy (XAS) with TDDFT
  + Analytical gradients for the COSMO solvation model
  + Transition densities from TDDFT 
  + DFT+U and Electron-Transfer (ET) methods for plane wave calculations
  + Exploitation of space group symmetry in plane wave geometry optimizations
  + Local density of states (LDOS) collective variable added to Metadynamics
  + Various new XC functionals added for plane wave calculations, including
    hybrid and range-corrected ones
  + Electric field gradients with relativistic corrections 
  + Nudged Elastic Band optimization method
  + Updated basis sets and ECPs 

[ Daniel Leidert ]
* debian/watch: Fixed.

[ Andreas Tille ]
* debian/upstream: References

[ Michael Banck ]
* debian/upstream (Name): New field.
* debian/patches/02_makefile_flags.patch: Refreshed.
* debian/patches/06_statfs_kfreebsd.patch: Likewise.
* debian/patches/07_ga_target_force_linux.patch: Likewise.
* debian/patches/05_avoid_inline_assembler.patch: Removed, no longer needed.
* debian/patches/09_backported_6.1.1_fixes.patch: Likewise.
* debian/control (Build-Depends): Added gfortran-4.7 and gcc-4.7.
* debian/patches/10_force_gcc-4.7.patch: New patch, explicitly sets
  gfortran-4.7 and gcc-4.7, fixes test suite hang with gcc-4.8 (Closes:
  #701328, #713262).
* debian/testsuite: Added tests for COSMO analytical gradients and MRCC.
* debian/rules (MRCC_METHODS): New variable, required to enable MRCC methods.

Show diffs side-by-side

added added

removed removed

Lines of Context:
src/tools/ga-5-1/global/testing/perform.F
1
 
#if HAVE_CONFIG_H
2
 
#   include "config.fh"
3
 
#endif
4
 
c $Id: perform.F,v 1.9 2000-05-25 01:09:20 d3h325 Exp $
5
 
      program perform
6
 
c***
7
 
c***  Last modification: Fri Jan 13 12:13:27 PST 1995
8
 
c***
9
 
      implicit none
10
 
#include "mafdecls.fh"
11
 
#include "global.fh"
12
 
      integer heap
13
 
c
14
 
c***  Intitialize a message passing library
15
 
c
16
 
#include "mp3.fh"
17
 
c
18
 
c***  Intitialize the GA package
19
 
      call ga_initialize()
20
 
      if(ga_nnodes().ne.4 .and. ga_nodeid().eq.0)
21
 
     $   call ga_error('Program requires 4 GA processes',ga_nnodes())
22
 
c
23
 
c***  Initialize the MA package
24
 
      heap = 450000
25
 
      if (.not. ma_init(MT_DBL, heap,heap))
26
 
     $     call ga_error('ma init failed',2*heap)
27
 
c
28
 
      call testit()
29
 
c
30
 
      if(ga_nodeid().eq.0) print *, 'All tests successful'
31
 
c
32
 
      call ga_terminate()
33
 
c
34
 
      call MP_FINALIZE()
35
 
      end
36
 
 
37
 
 
38
 
c-----------------
39
 
 
40
 
      subroutine testit()
41
 
      implicit none
42
 
#include "mafdecls.fh"
43
 
#include "global.fh"
44
 
c
45
 
      integer n, nn, num_chunks
46
 
      parameter (n = 710, nn = n*n/4, num_chunks=12)
47
 
      double precision buf(nn)
48
 
c
49
 
      integer g_a
50
 
      integer ilo, ihi, jlo, jhi
51
 
      integer nproc, me, loop
52
 
      integer chunk(num_chunks)
53
 
      data    chunk /1,3,4,9,16,30,64,91,128,171,256,353/
54
 
c     
55
 
      nproc = ga_nnodes()
56
 
      me = ga_nodeid()
57
 
c
58
 
c***  Create global array
59
 
      if (.not. ga_create(MT_DBL, n, n, 'a', 0, 0, g_a))
60
 
     $     call ga_error(' ga_create failed ',1)
61
 
c     
62
 
      do loop=1,nn
63
 
         buf(loop) = .01d0
64
 
      enddo
65
 
      call ga_zero(g_a) 
66
 
c
67
 
      if (me .eq. 0) then
68
 
        write(*,*)' '
69
 
        print *,'> Performance of ga_get, ga_put & ga_acc n = ', n
70
 
        print *,' '
71
 
      endif
72
 
c
73
 
      
74
 
c     do loop=1,2
75
 
c
76
 
c***  local ops
77
 
c
78
 
      call ga_distribution(g_a, me, ilo, ihi, jlo, jhi)
79
 
      call TestPutGetAcc
80
 
     &     (g_a, n, chunk, num_chunks, buf, ilo, ihi, jlo, jhi)
81
 
c
82
 
c***  remote ops
83
 
c
84
 
      call ga_distribution(g_a, nproc-1, ilo, ihi, jlo, jhi)
85
 
      call TestPutGetAcc
86
 
     &     (g_a, n, chunk, num_chunks, buf, ilo, ihi, jlo, jhi)
87
 
 
88
 
c     enddo
89
 
      end
90
 
 
91
 
 
92
 
      subroutine TestPutGetAcc
93
 
     &          (g_a, n, chunk, num_chunks, buf, ilo, ihi, jlo,jhi)
94
 
      implicit none
95
 
#include "global.fh"
96
 
#include "testutil.fh"
97
 
c
98
 
      integer num_chunks, chunk(num_chunks)
99
 
      integer n, ilo, ihi, jlo,jhi,g_a
100
 
      double precision buf(*), tg, tp, ta
101
 
      double precision time_acc, time_get, time_put
102
 
c
103
 
      integer me
104
 
      integer loop, jump, count, bytes
105
 
c
106
 
      me = ga_nodeid()
107
 
      if (me .eq. 0) then
108
 
        write(6,*)' '
109
 
        write(6,'(21X,8hACCESS [,i3,1h:,i4,1h,,i3,1h:,i4,1h])') 
110
 
     &        ilo,ihi,jlo,jhi
111
 
        write(6,*)'bytes  loop         get                    put',
112
 
     &           '                 accumulate'
113
 
        call flush(6)
114
 
      endif
115
 
      call ga_sync()
116
 
c
117
 
      do loop = 1, num_chunks
118
 
        bytes = util_mdtob(1)*chunk(loop)*chunk(loop) !how much data is accessed
119
 
        jump  =  n/(60*loop) ! jump distance between consecutive patches
120
 
c
121
 
c       everybody touches own data 
122
 
        call ga_fill_patch(g_a, 1, n, 1, n , 1d0*me*loop)
123
 
        if (me .eq. 0) then
124
 
          tg= time_get(g_a,ilo,ihi,jlo,jhi,buf,chunk(loop),jump,count)
125
 
        else
126
 
          call sleep(2)
127
 
        endif
128
 
c
129
 
c       everybody touches own data
130
 
        call ga_fill_patch(g_a, 1, n, 1, n , 1d0*me*loop)
131
 
        if (me .eq. 0) then
132
 
          tp= time_put(g_a,ilo,ihi,jlo,jhi,buf,chunk(loop),jump,count)
133
 
        else
134
 
          call sleep(2)
135
 
        endif
136
 
c
137
 
c       everybody touches own data
138
 
        call ga_fill_patch(g_a, 1, n, 1, n , 1d0*me*loop)
139
 
        if (me .eq. 0) then
140
 
          ta= time_acc(g_a,ilo,ihi,jlo,jhi,buf,chunk(loop),jump,count)
141
 
        else
142
 
          call sleep(2)
143
 
        endif
144
 
c
145
 
        if (me .eq. 0) then
146
 
          write(6,77)bytes, count, tg, 1d-6*bytes/tg,
147
 
     &               tp, 1d-6*bytes/tp, ta, 1d-6*bytes/ta
148
 
          call flush(6)
149
 
        endif
150
 
      enddo
151
 
c
152
 
77    format(i6, i5, 3(d10.3, d9.2,'MB/s'))
153
 
      end
154
 
 
155
 
 
156
 
 
157
 
      double precision function
158
 
     &       time_acc(g_a, is, ie, js, je, buf, chunk, jump, count)
159
 
c
160
 
      implicit none
161
 
#include "global.fh"
162
 
#include "testutil.fh"
163
 
c
164
 
 
165
 
      integer g_a, chunk, jump, count, is, js, ie, je
166
 
c
167
 
      integer ilo, ihi, jlo, jhi
168
 
      double precision  seconds, buf 
169
 
c
170
 
      count = 0
171
 
      seconds = util_timer()
172
 
c
173
 
c       distance between consecutive patches increased by jump
174
 
c       to destroy locality of reference
175
 
        do ilo = is, ie -chunk-jump +1, chunk+jump
176
 
           ihi = ilo + chunk -1
177
 
           do jlo = js, je -chunk-jump +1, chunk+jump
178
 
              jhi = jlo + chunk -1
179
 
              count = count + 1
180
 
              call ga_acc(g_a, ilo, ihi, jlo, jhi, buf, chunk, 1d0)
181
 
           enddo
182
 
        enddo
183
 
      seconds = util_timer() - seconds
184
 
c
185
 
      time_acc = seconds/count
186
 
      end
187
 
 
188
 
 
189
 
      double precision function
190
 
     &       time_get(g_a, is, ie, js, je, buf, chunk, jump, count)
191
 
c
192
 
      implicit none
193
 
#include "global.fh"
194
 
#include "testutil.fh"
195
 
c
196
 
      integer g_a, chunk, jump, count, is, js, ie, je
197
 
c
198
 
      integer ilo, ihi, jlo, jhi
199
 
      double precision  seconds, buf
200
 
c
201
 
      count = 0
202
 
      seconds = util_timer()
203
 
c
204
 
c       distance between consecutive patches increased by jump
205
 
c       to destroy locality of reference
206
 
        do ilo = is, ie -chunk-jump +1, chunk+jump
207
 
           ihi = ilo + chunk -1
208
 
           do jlo = js, je -chunk-jump +1, chunk+jump
209
 
              jhi = jlo + chunk -1
210
 
              count = count + 1
211
 
              call ga_get(g_a, ilo, ihi, jlo, jhi, buf, chunk)
212
 
           enddo
213
 
        enddo
214
 
      seconds = util_timer() - seconds
215
 
c
216
 
      time_get = seconds/count
217
 
      end
218
 
 
219
 
 
220
 
 
221
 
      double precision function
222
 
     &       time_put(g_a, is, ie, js, je, buf, chunk, jump, count)
223
 
c
224
 
      implicit none
225
 
#include "global.fh"
226
 
#include "testutil.fh"
227
 
c
228
 
      integer g_a, chunk, jump, count, is, js, ie, je
229
 
c
230
 
      integer ilo, ihi, jlo, jhi
231
 
      double precision  seconds, buf
232
 
c
233
 
      count = 0
234
 
      seconds = util_timer()
235
 
c
236
 
c       distance between consecutive patches increased by jump
237
 
c       to destroy locality of reference
238
 
        do ilo = is, ie -chunk-jump +1, chunk+jump
239
 
           ihi = ilo + chunk -1
240
 
           do jlo = js, je -chunk-jump +1, chunk+jump
241
 
              jhi = jlo + chunk -1
242
 
              count = count + 1
243
 
              call ga_put(g_a, ilo, ihi, jlo, jhi, buf, chunk)
244
 
           enddo
245
 
        enddo
246
 
      seconds = util_timer() - seconds
247
 
c
248
 
      time_put = seconds/count
249
 
      end