1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
|
#if HAVE_CONFIG_H
# include "config.fh"
#endif
program main
implicit double precision (a-h,o-z)
c
c FORTRAN program to test message passing routines
c
c LOG is the FORTRAN unit number for standard output.
c
parameter (LOG = 6)
parameter (MAXLEN = 262144 )
#include "msgtypesf.h"
dimension buf(MAXLEN)
integer ibuf(MAXLEN)
character*80 fname
integer dtype
c
c Always the first thing to do is call pbeginf
c
call pbeginf
call setdbg(0)
call evon
c
c who am i and how many processes
c
nproc = nnodes()
me = nodeid()
c
c now try broadcasting messages from all nodes to every other node
c send each process my id as a message
c
call evbgin('Hello test')
itype = 1 + MSGINT
do 10 iproc = 0,nproc-1
itest = me
call brdcst(itype, itest, mitob(1),iproc)
if (iproc.ne.me) then
write(LOG,1) me, itest
1 format(' me=',i3,', itest=',i3)
endif
10 continue
call evend('Hello test')
call evbgin('Counter test')
c
c now try using the shared counter
c
mproc = nproc
do 20 i = 1,10
write(LOG,*) ' process ',me,' got ',nxtval(mproc)
20 continue
junk = nxtval(-mproc)
call evend('Counter test')
c
c now time sending a message round a ring
c
if (nproc.gt.1) then
call evbgin('Ring test')
itype = 3
left = mod(me + nproc - 1, nproc)
iright = mod(me + 1, nproc)
c
lenbuf = 1
30 if (me .eq. 0) then
start = tcgtime()
call snd(itype, buf, lenbuf, left, 1)
call rcv(itype, buf, lenbuf, lenmes, iright, node, 1)
used = tcgtime() - start
if (used.gt.0d0) then
rate = 1.0d-6 * dble(nproc * lenbuf) / used
else
rate = 0.0d0
endif
write(LOG,31) lenbuf, used, rate
else
call rcv(itype, buf, lenbuf, lenmes, iright, node, 1)
call snd(itype, buf, lenbuf, left, 1)
endif
lenbuf = lenbuf * 2
if (lenbuf .le. mdtob(MAXLEN)) goto 30
31 format(' len=',i7,' bytes, used=',f8.2,' cs, rate=',f10.6,' Mb/s')
call evend('Ring test')
endif
c
c global sums
c
do i=1,MAXLEN
ibuf(i) = i*me
buf(i) = dble(ibuf(i))
enddo
dtype=1+MSGDBL
call igop(itype, ibuf, MAXLEN, "+")
call dgop(dtype, buf, MAXLEN, "+")
do i=1,MAXLEN
iresult = i*nproc*(nproc-1)/2
if (ibuf(i).ne.iresult.or.buf(i).ne.dble(iresult))
. call error('TestGlobals: global sum failed', i)
enddo
if (me.eq.0) write(LOG,*) 'global sums OK'
c
c
c Check that everyone can open, write, read and close
c a binary FORTRAN file
c
call pfname('junk',fname)
open(9,file=fname,form='unformatted',status='unknown',
& err=1000)
write(9,err=1001) buf
rewind 9
read(9,err=1002) buf
close(9,status='delete')
call event('Read file OK')
c
if (me.eq.0) call stats
c
c Always the last thing to do is call pend
c
call pend
c
c check that everyone makes it thru after pend .. NODEID
c is not actually guaranteed to work outside of pbegin/pend
c section ... it may return junk. All you should do is exit
c is some FORTRAN supported fashion
c
write(LOG,32) nodeid()
32 format(' Process ',i4,' after pend')
stop
c
c error returns for FORTRAN I/O
c
1000 call error('failed to open fortran binary file',-1)
1001 call error('failed to write fortran binary file',-1)
1002 call error('failed to read fortran binary file',-1)
c
end
subroutine pfname(name, fname)
character*(*) name, fname
c
c construct a unique filename by appending the process
c number after the stub name
c i.e. <fname> = <name>.<mynode>
c
c find last non-blank character in name
c
do 10 i = len(name),1,-1
if (name(i:i).ne.' ') goto 20
10 continue
call error('pfname: name is all blanks!',i)
c
c check that have room to store result and then write result
c
20 if (i+4.gt.len(fname))
& call error('pfname: fname too short for name.id',len(fname))
fname = name
write(fname(i+1:i+4),1) nodeid()
1 format('.',i3.3)
c
end
subroutine error(s,i)
parameter (LOG = 6)
character*(*) s
integer i
c
write(LOG,1) s,i
1 format(//
$ ' !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!'/
$ 1x,a,1x,i8/
$ ' !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!'/)
c $ 1x,a,1x,i8/
c
call parerr(i)
c
end
|
