Actual source code: ex1f.F90
petsc-3.14.6 2021-03-30
1: !
2: !
3: ! Description: Demonstrates how users can augment the PETSc profiling by
4: ! inserting their own event logging.
5: !
6: !/*T
7: ! Concepts: PetscLog^user-defined event profiling (basic example);
8: ! Concepts: PetscLog^activating/deactivating events for profiling (basic example);
9: ! Processors: n
10: !T*/
11: ! -----------------------------------------------------------------------
13: program SchoolDay
14: #include <petsc/finclude/petscsys.h>
15: #include <petsc/finclude/petsclog.h>
16: use petscsys
17: implicit none
19: !====================================================================
20: ! Local Variables
22: ! Settings:
23: integer, parameter :: verbose=0 ! 0: silent, >=1 : increasing amount of debugging output
24: integer, parameter :: msgLen = 30 ! number of reals which is sent with MPI_Isend
25: PetscReal, parameter :: second=0.1; ! time is sped up by a factor 10
27: ! Codes
28: integer, parameter :: BOY=1, GIRL=2, TEACHER=0
29: PetscMPIInt, parameter :: tagMsg = 1200;
31: ! Timers
32: PetscLogEvent :: Morning, Afternoon
33: PetscLogEvent :: PlayBall, SkipRope
34: PetscLogEvent :: TidyClass
35: PetscLogEvent :: Lessons, CorrectHomework
36: PetscClassId classid
38: ! Petsc-stuff
39: PetscErrorCode :: ierr
41: ! MPI-stuff
42: PetscMPIInt :: rank, size
43: PetscReal, allocatable :: message(:,:)
44: integer :: item, maxItem
45: integer4 :: status(MPI_STATUS_SIZE)
46: PetscMPIInt req
48: ! Own stuff
49: integer4 :: role ! is this process a BOY, a GIRL or a TEACHER?
50: integer4 :: i, j
51: integer4,parameter :: one=1
52: !====================================================================
53: ! Initializations
54: call PetscInitialize(PETSC_NULL_CHARACTER,ierr)
55: if (ierr .ne. 0) then
56: print*,'Unable to initialize PETSc'
57: stop
58: endif
59: call MPI_Comm_size(PETSC_COMM_WORLD, size,ierr)
60: call MPI_Comm_rank(PETSC_COMM_WORLD, rank,ierr)
62: if (rank==0) then
63: role = TEACHER
64: else if (rank<0.4*size) then
65: role = GIRL
66: else
67: role = BOY
68: end if
70: allocate(message(msgLen,msglen))
71: do i = 1,msgLen
72: do j = 1,msgLen
73: message(i,j) = 10.0*j + i*1.0/(rank+one)
74: end do
75: end do
76: !
77: !====================================================================
78: ! Create new user-defined events
79: classid = 0
80: call PetscLogEventRegister('Morning', classid, Morning, ierr)
81: call PetscLogEventRegister('Afternoon', classid, Afternoon, ierr)
82: call PetscLogEventRegister('Play Ball', classid, PlayBall, ierr)
83: call PetscLogEventRegister('Skip Rope', classid, SkipRope, ierr)
84: call PetscLogEventRegister('Tidy Classroom', classid, TidyClass, ierr)
85: call PetscLogEventRegister('Lessons', classid, Lessons, ierr)
86: call PetscLogEventRegister('Correct Homework',classid,CorrectHomework, &
87: & ierr)
88: if (verbose>=1) then
89: print '(a,i0,a)','[',rank,'] SchoolDay events have been defined'
90: endif
92: !====================================================================
93: ! Go through the school day
94: call PetscLogEventBegin(Morning,ierr)
96: call PetscLogFlops(190000d0,ierr)
97: call PetscSleep(0.5*second,ierr)
99: call PetscLogEventBegin(Lessons,ierr)
100: call PetscLogFlops(23000d0,ierr)
101: call PetscSleep(1*second, ierr)
102: if (size>1) then
103: call MPI_Isend( message, msgLen, MPI_DOUBLE_PRECISION, &
104: & mod(rank+1,size), &
105: & tagMsg+rank, PETSC_COMM_WORLD, req, ierr)
106: call MPI_Recv( message, msgLen, MPI_DOUBLE_PRECISION, &
107: & mod(rank-1+size,size), &
108: & tagMsg+mod(rank-1+size,size), PETSC_COMM_WORLD, &
109: & status, ierr)
110: call MPI_Wait(req,MPI_STATUS_IGNORE,ierr)
111: end if
112: call PetscLogEventEnd(Lessons,ierr)
114: if (role==TEACHER) then
115: call PetscLogEventBegin(TidyClass,ierr)
116: call PetscLogFlops(600000d0,ierr)
117: call PetscSleep(0.6*second, ierr)
118: call PetscLogEventBegin(CorrectHomework,ierr)
119: call PetscLogFlops(234700d0,ierr)
120: call PetscSleep(0.4*second, ierr)
121: call PetscLogEventEnd(CorrectHomework,ierr)
122: call PetscLogEventEnd(TidyClass,ierr)
123: else if (role==BOY) then
124: call PetscLogEventBegin(SkipRope,ierr)
125: call PetscSleep(0.8*second, ierr)
126: call PetscLogEventEnd(SkipRope,ierr)
127: else
128: call PetscLogEventBegin(PlayBall,ierr)
129: call PetscSleep(0.9*second, ierr)
130: call PetscLogEventEnd(PlayBall,ierr)
131: end if
133: call PetscLogEventBegin(Lessons,ierr)
134: call PetscLogFlops(120000d0,ierr)
135: call PetscSleep(0.7*second, ierr)
136: call PetscLogEventEnd(Lessons,ierr)
138: call PetscLogEventEnd(Morning,ierr)
140: call PetscLogEventBegin(Afternoon,ierr)
142: item = rank*(3-rank)
143: call MPI_Allreduce(item, maxItem, 1, MPI_INTEGER, MPI_MAX, &
144: & PETSC_COMM_WORLD, ierr)
146: item = rank*(10-rank)
147: call MPI_Allreduce(item, maxItem, 1, MPI_INTEGER, MPI_MAX, &
148: & PETSC_COMM_WORLD, ierr)
150: call PetscLogFlops(58988d0,ierr)
151: call PetscSleep(0.6*second,ierr)
153: call PetscLogEventBegin(Lessons,ierr)
154: call PetscLogFlops(123456d0,ierr)
155: call PetscSleep(1*second, ierr)
156: call PetscLogEventEnd(Lessons,ierr)
158: if (role==TEACHER) then
159: call PetscLogEventBegin(TidyClass,ierr)
160: call PetscLogFlops(17800d0,ierr)
161: call PetscSleep(1.1*second, ierr)
162: call PetscLogEventBegin(Lessons,ierr)
163: call PetscLogFlops(72344d0,ierr)
164: call PetscSleep(0.5*second, ierr)
165: call PetscLogEventEnd(Lessons,ierr)
166: call PetscLogEventEnd(TidyClass,ierr)
167: else if (role==GIRL) then
168: call PetscLogEventBegin(SkipRope,ierr)
169: call PetscSleep(0.7*second, ierr)
170: call PetscLogEventEnd(SkipRope,ierr)
171: else
172: call PetscLogEventBegin(PlayBall,ierr)
173: call PetscSleep(0.8*second, ierr)
174: call PetscLogEventEnd(PlayBall,ierr)
175: end if
177: call PetscLogEventBegin(Lessons,ierr)
178: call PetscLogFlops(72344d0,ierr)
179: call PetscSleep(0.5*second, ierr)
180: call PetscLogEventEnd(Lessons,ierr)
182: call PetscLogEventEnd(Afternoon,ierr)
184: if (.false.) then
185: continue
186: else if (role==TEACHER) then
187: call PetscLogEventBegin(TidyClass,ierr)
188: call PetscLogFlops(612300d0,ierr)
189: call PetscSleep(1.1*second, ierr)
190: call PetscLogEventEnd(TidyClass,ierr)
191: call PetscLogEventBegin(CorrectHomework,ierr)
192: call PetscLogFlops(234700d0,ierr)
193: call PetscSleep(1.1*second, ierr)
194: call PetscLogEventEnd(CorrectHomework,ierr)
195: else
196: call PetscLogEventBegin(SkipRope,ierr)
197: call PetscSleep(0.7*second, ierr)
198: call PetscLogEventEnd(SkipRope,ierr)
199: call PetscLogEventBegin(PlayBall,ierr)
200: call PetscSleep(0.8*second, ierr)
201: call PetscLogEventEnd(PlayBall,ierr)
202: end if
204: call PetscLogEventBegin(Lessons,ierr)
205: call PetscLogFlops(120000d0,ierr)
206: call PetscSleep(0.7*second, ierr)
207: call PetscLogEventEnd(Lessons,ierr)
209: call PetscSleep(0.25*second,ierr)
211: call PetscLogEventBegin(Morning,ierr)
213: call PetscLogFlops(190000d0,ierr)
214: call PetscSleep(0.5*second,ierr)
216: call PetscLogEventBegin(Lessons,ierr)
217: call PetscLogFlops(23000d0,ierr)
218: call PetscSleep(1*second, ierr)
219: if (size>1) then
220: call MPI_Isend( message, msgLen, MPI_DOUBLE_PRECISION, &
221: & mod(rank+1,size), &
222: & tagMsg+rank, PETSC_COMM_WORLD, req, ierr)
223: call MPI_Recv( message, msgLen, MPI_DOUBLE_PRECISION, &
224: & mod(rank-1+size,size), &
225: & tagMsg+mod(rank-1+size,size), PETSC_COMM_WORLD, &
226: & status, ierr)
227: call MPI_Wait(req,MPI_STATUS_IGNORE,ierr)
228: end if
229: call PetscLogEventEnd(Lessons,ierr)
231: if (role==TEACHER) then
232: call PetscLogEventBegin(TidyClass,ierr)
233: call PetscLogFlops(600000d0,ierr)
234: call PetscSleep(1.2*second, ierr)
235: call PetscLogEventEnd(TidyClass,ierr)
236: else if (role==BOY) then
237: call PetscLogEventBegin(SkipRope,ierr)
238: call PetscSleep(0.8*second, ierr)
239: call PetscLogEventEnd(SkipRope,ierr)
240: else
241: call PetscLogEventBegin(PlayBall,ierr)
242: call PetscSleep(0.9*second, ierr)
243: call PetscLogEventEnd(PlayBall,ierr)
244: end if
246: call PetscLogEventBegin(Lessons,ierr)
247: call PetscLogFlops(120000d0,ierr)
248: call PetscSleep(0.7*second, ierr)
249: call PetscLogEventEnd(Lessons,ierr)
251: call PetscLogEventEnd(Morning,ierr)
253: deallocate(message)
255: call PetscFinalize(ierr)
257: end program SchoolDay
259: !/*TEST
260: !
261: ! testset:
262: ! args: -log_view ascii:filename.txt
263: ! output_file: output/ex1f.out
264: ! test:
265: ! suffix: 1
266: ! nsize: 1
267: ! test:
268: ! suffix: 2
269: ! nsize: 2
270: ! test:
271: ! suffix: 3
272: ! nsize: 3
273: !
274: ! testset:
275: ! suffix: detail
276: ! args: -log_view ascii:filename.txt:ascii_info_detail
277: ! output_file: output/ex1f.out
278: ! test:
279: ! suffix: 1
280: ! nsize: 1
281: ! test:
282: ! suffix: 2
283: ! nsize: 2
284: ! test:
285: ! suffix: 3
286: ! nsize: 3
287: !
288: ! testset:
289: ! suffix: xml
290: ! args: -log_view ascii:filename.xml:ascii_xml
291: ! output_file: output/ex1f.out
292: ! test:
293: ! suffix: 1
294: ! nsize: 1
295: ! test:
296: ! suffix: 2
297: ! nsize: 2
298: ! test:
299: ! suffix: 3
300: ! nsize: 3
301: !
302: !TEST*/