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
|
C -*- Mode: Fortran; -*-
C
C (C) 2012 by Argonne National Laboratory.
C See COPYRIGHT in top-level directory.
C
c------------------------------------------------------------------------------
c
c Check for correct source, tag, count, and data in test message.
c
c------------------------------------------------------------------------------
subroutine msg_check( recv_buf, source, tag, count, status, n,
* name, errs )
implicit none
include 'mpif.h'
integer n, errs
real recv_buf(n)
integer source, tag, count, rank, status(MPI_STATUS_SIZE)
character*(*) name
logical foundError
integer ierr, recv_src, recv_tag, recv_count
foundError = .false.
recv_src = status(MPI_SOURCE)
recv_tag = status(MPI_TAG)
call MPI_Comm_rank( MPI_COMM_WORLD, rank, ierr )
call MPI_Get_count(status, MPI_REAL, recv_count, ierr)
if (recv_src .ne. source) then
print *, '[', rank, '] Unexpected source:', recv_src,
* ' in ', name
errs = errs + 1
foundError = .true.
end if
if (recv_tag .ne. tag) then
print *, '[', rank, '] Unexpected tag:', recv_tag, ' in ', name
errs = errs + 1
foundError = .true.
end if
if (recv_count .ne. count) then
print *, '[', rank, '] Unexpected count:', recv_count,
* ' in ', name
errs = errs + 1
foundError = .true.
end if
call verify_test_data(recv_buf, count, n, name, errs )
end
c------------------------------------------------------------------------------
c
c Check that requests have been set to null
c
c------------------------------------------------------------------------------
subroutine rq_check( requests, n, msg )
include 'mpif.h'
integer n, requests(n)
character*(*) msg
integer i
c
do 10 i=1, n
if (requests(i) .ne. MPI_REQUEST_NULL) then
print *, 'Nonnull request in ', msg
endif
10 continue
c
end
c------------------------------------------------------------------------------
c
c Initialize test data buffer with integral sequence.
c
c------------------------------------------------------------------------------
subroutine init_test_data(buf,n)
integer n
real buf(n)
integer i
do 10 i = 1, n
buf(i) = REAL(i)
10 continue
end
c------------------------------------------------------------------------------
c
c Clear test data buffer
c
c------------------------------------------------------------------------------
subroutine clear_test_data(buf, n)
integer n
real buf(n)
integer i
do 10 i = 1, n
buf(i) = 0.
10 continue
end
c------------------------------------------------------------------------------
c
c Verify test data buffer
c
c------------------------------------------------------------------------------
subroutine verify_test_data( buf, count, n, name, errs )
implicit none
include 'mpif.h'
integer n, errs
real buf(n)
character *(*) name
integer count, ierr, i
C
do 10 i = 1, count
if (buf(i) .ne. REAL(i)) then
print 100, buf(i), i, count, name
errs = errs + 1
endif
10 continue
C
do 20 i = count + 1, n
if (buf(i) .ne. 0.) then
print 100, buf(i), i, n, name
errs = errs + 1
endif
20 continue
C
100 format('Invalid data', f6.1, ' at ', i4, ' of ', i4, ' in ', a)
C
end
C
C This routine is used to prevent the compiler from deallocating the
C array "a", which may happen in some of the tests (see the text in
C the MPI standard about why this may be a problem in valid Fortran
C codes). Without this, for example, tests fail with the Cray ftn
C compiler.
C
subroutine dummyRef( a, n, ie )
integer n, ie
real a(n)
C This condition will never be true, but the compile won't know that
if (ie .eq. -1) then
print *, a(n)
endif
return
end
|
