fix typo
This commit is contained in:
gabrielkheisa
parent
436d13ab14
commit
115bcdb4c8
242
tugas.f
Normal file
242
tugas.f
Normal file
@ -0,0 +1,242 @@
|
||||
C*********************************************************************
|
||||
C
|
||||
C This benchmark test program is measuring a cpu performance
|
||||
C of floating point operation by a Poisson equation solver.
|
||||
CC
|
||||
C If you have any question, please ask me via email.
|
||||
C written by Ryutaro HIMENO, November 26, 2001.
|
||||
C Version 3.0
|
||||
C ----------------------------------------------
|
||||
C Ryutaro Himeno, Dr. of Eng.
|
||||
C Head of Computer Information Division,
|
||||
C RIKEN (The Institute of Pysical and Chemical Research)
|
||||
C Email : himeno@postman.riken.go.jp
|
||||
C ---------------------------------------------------------------
|
||||
C You can adjust the size of this benchmark code to fit your target
|
||||
C computer. In that case, please chose following sets of
|
||||
C (mimax,mjmax,mkmax):
|
||||
C small : 65,33,33
|
||||
C small : 129,65,65
|
||||
C midium: 257,129,129
|
||||
C large : 513,257,257
|
||||
C ext.large: 1025,513,513
|
||||
C This program is to measure a computer performance in MFLOPS
|
||||
C by using a kernel which appears in a linear solver of pressure
|
||||
C Poisson eq. which appears in an incompressible Navier-Stokes solver.
|
||||
C A point-Jacobi method is employed in this solver as this method can
|
||||
C be easyly vectrized and be parallelized.
|
||||
C ------------------
|
||||
C Finite-difference method, curvilinear coodinate system
|
||||
C Vectorizable and parallelizable on each grid point
|
||||
C No. of grid points : imax x jmax x kmax including boundaries
|
||||
C ------------------
|
||||
C A,B,C:coefficient matrix, wrk1: source term of Poisson equation
|
||||
C wrk2 : working area, OMEGA : relaxation parameter
|
||||
C BND:control variable for boundaries and objects ( = 0 or 1)
|
||||
C P: pressure
|
||||
C -------------------
|
||||
C -------------------
|
||||
C "use portlib" statement on the next line is for Visual fortran
|
||||
C to use UNIX libraries. Please remove it if your system is UNIX.
|
||||
C -------------------
|
||||
! use portlib
|
||||
use omp_lib
|
||||
IMPLICIT REAL*4(a-h,o-z)
|
||||
real*8 t1,t2
|
||||
C
|
||||
C PARAMETER (mimax=513,mjmax=257,mkmax=257)
|
||||
C PARAMETER (mimax=257,mjmax=129,mkmax=129)
|
||||
PARAMETER (mimax=129,mjmax=65,mkmax=65)
|
||||
C PARAMETER (mimax=65,mjmax=33,mkmax=33)
|
||||
C
|
||||
C ttarget specifys the measuring period in sec
|
||||
PARAMETER (ttarget=60.0)
|
||||
CC Arrey
|
||||
common /pres/ p(mimax,mjmax,mkmax)
|
||||
common /mtrx/ a(mimax,mjmax,mkmax,4),
|
||||
+ b(mimax,mjmax,mkmax,3),c(mimax,mjmax,mkmax,3)
|
||||
common /bound/ bnd(mimax,mjmax,mkmax)
|
||||
common /work/ wrk1(mimax,mjmax,mkmax),wrk2(mimax,mjmax,mkmax)
|
||||
CC Other constants
|
||||
common /others/ imax,jmax,kmax,omega
|
||||
C
|
||||
dimension time0(2),time1(2)
|
||||
C
|
||||
omega=0.8
|
||||
imax=mimax-1
|
||||
jmax=mjmax-1
|
||||
kmax=mkmax-1
|
||||
CC Initializing matrixes
|
||||
call initmt
|
||||
write(*,*) ' mimax=',mimax,' mjmax=',mjmax,' mkmax=',mkmax
|
||||
write(*,*) ' imax=',imax,' jmax=',jmax,' kmax=',kmax
|
||||
CC Start measuring
|
||||
C
|
||||
nn=10000
|
||||
write(*,*) ' Start rehearsal measurement process.'
|
||||
write(*,*) ' Measure the performance in 10000 times.'
|
||||
C
|
||||
cpu0=dtime(time0)
|
||||
t1 = omp_get_wtime()
|
||||
!$OMP PARALLEL DO
|
||||
C
|
||||
C Jacobi iteration
|
||||
call jacobi(nn,gosa)
|
||||
C
|
||||
cpu1= dtime(time1)
|
||||
!$OMP END PARALLEL DO
|
||||
t2 = omp_get_wtime()
|
||||
cpu = cpu1
|
||||
cpu = t2-t1
|
||||
flop=real(kmax-2)*real(jmax-2)*real(imax-2)*34.0*real(nn)
|
||||
xmflops2=flop/cpu*1.0e-6
|
||||
write(*,*) ' MFLOPS:',xmflops2,' time(s):',cpu,gosa
|
||||
C
|
||||
C end the test loop
|
||||
! nn=ifix(ttarget/(cpu/3.0))
|
||||
! write(*,*) 'Now, start the actual measurement process.'
|
||||
! write(*,*) 'The loop will be excuted in',nn,' times.'
|
||||
! write(*,*) 'This will take about one minute.'
|
||||
! write(*,*) 'Wait for a while.'
|
||||
C
|
||||
C Jacobi iteration
|
||||
! cpu0=dtime(time0)
|
||||
! call jacobi(nn,gosa)
|
||||
C
|
||||
! cpu1= dtime(time1)
|
||||
! cpu = cpu1
|
||||
! flop=real(kmax-2)*real(jmax-2)*real(imax-2)*34.0*real(nn)
|
||||
! xmflops2=flop*1.0e-6/cpu
|
||||
C
|
||||
CCC xmflops2=nflop/cpu*1.0e-6*float(nn)
|
||||
C
|
||||
! write(*,*) ' Loop executed for ',nn,' times'
|
||||
! write(*,*) ' Gosa :',gosa
|
||||
! write(*,*) ' MFLOPS:',xmflops2, ' time(s):',cpu
|
||||
! score=xmflops2/82.84
|
||||
! write(*,*) ' Score based on Pentium III 600MHz :',score
|
||||
C
|
||||
! pause
|
||||
stop
|
||||
END
|
||||
C
|
||||
C
|
||||
C**************************************************************
|
||||
subroutine initmt
|
||||
C**************************************************************
|
||||
IMPLICIT REAL*4(a-h,o-z)
|
||||
C
|
||||
C PARAMETER (mimax=513,mjmax=257,mkmax=257)
|
||||
C PARAMETER (mimax=257,mjmax=129,mkmax=129)
|
||||
PARAMETER (mimax=129,mjmax=65,mkmax=65)
|
||||
C PARAMETER (mimax=65,mjmax=33,mkmax=33)
|
||||
C
|
||||
|
||||
|
||||
CC Arrey
|
||||
common /pres/ p(mimax,mjmax,mkmax)
|
||||
common /mtrx/ a(mimax,mjmax,mkmax,4),
|
||||
+ b(mimax,mjmax,mkmax,3),c(mimax,mjmax,mkmax,3)
|
||||
common /bound/ bnd(mimax,mjmax,mkmax)
|
||||
common /work/ wrk1(mimax,mjmax,mkmax),wrk2(mimax,mjmax,mkmax)
|
||||
CC other constants
|
||||
common /others/ imax,jmax,kmax,omega
|
||||
C
|
||||
do k=1,mkmax
|
||||
do j=1,mjmax
|
||||
do i=1,mimax
|
||||
a(i,j,k,1)=0.0
|
||||
a(i,j,k,2)=0.0
|
||||
a(i,j,k,3)=0.0
|
||||
a(i,j,k,4)=0.0
|
||||
b(i,j,k,1)=0.0
|
||||
b(i,j,k,2)=0.0
|
||||
b(i,j,k,3)=0.0
|
||||
c(i,j,k,1)=0.0
|
||||
c(i,j,k,2)=0.0
|
||||
c(i,j,k,3)=0.0
|
||||
p(i,j,k) =0.0
|
||||
wrk1(i,j,k)=0.0
|
||||
bnd(i,j,k)=0.0
|
||||
enddo
|
||||
enddo
|
||||
enddo
|
||||
C
|
||||
do k=1,kmax
|
||||
do j=1,jmax
|
||||
do i=1,imax
|
||||
a(i,j,k,1)=1.0
|
||||
a(i,j,k,2)=1.0
|
||||
a(i,j,k,3)=1.0
|
||||
a(i,j,k,4)=1.0/6.0
|
||||
b(i,j,k,1)=0.0
|
||||
b(i,j,k,2)=0.0
|
||||
b(i,j,k,3)=0.0
|
||||
c(i,j,k,1)=1.0
|
||||
c(i,j,k,2)=1.0
|
||||
c(i,j,k,3)=1.0
|
||||
p(i,j,k) =float((k-1)*(k-1))/float((kmax-1)*(kmax-1))
|
||||
wrk1(i,j,k)=0.0
|
||||
bnd(i,j,k)=1.0
|
||||
enddo
|
||||
enddo
|
||||
enddo
|
||||
C
|
||||
return
|
||||
end
|
||||
C
|
||||
C*************************************************************
|
||||
subroutine jacobi(nn,gosa)
|
||||
C*************************************************************
|
||||
IMPLICIT REAL*4(a-h,o-z)
|
||||
C
|
||||
C PARAMETER (mimax=513,mjmax=257,mkmax=257)
|
||||
C PARAMETER (mimax=257,mjmax=129,mkmax=129)
|
||||
PARAMETER (mimax=129,mjmax=65,mkmax=65)
|
||||
C PARAMETER (mimax=65,mjmax=33,mkmax=33)
|
||||
C
|
||||
|
||||
CC Arrey
|
||||
common /pres/ p(mimax,mjmax,mkmax)
|
||||
common /mtrx/ a(mimax,mjmax,mkmax,4),
|
||||
+ b(mimax,mjmax,mkmax,3),c(mimax,mjmax,mkmax,3)
|
||||
common /bound/ bnd(mimax,mjmax,mkmax)
|
||||
common /work/ wrk1(mimax,mjmax,mkmax),wrk2(mimax,mjmax,mkmax)
|
||||
CC other constants
|
||||
common /others/ imax,jmax,kmax,omega
|
||||
C
|
||||
C
|
||||
DO loop=1,nn
|
||||
gosa=0.0
|
||||
DO K=2,kmax-1
|
||||
DO J=2,jmax-1
|
||||
DO I=2,imax-1
|
||||
S0=a(I,J,K,1)*p(I+1,J,K)+a(I,J,K,2)*p(I,J+1,K)
|
||||
1 +a(I,J,K,3)*p(I,J,K+1)
|
||||
2 +b(I,J,K,1)*(p(I+1,J+1,K)-p(I+1,J-1,K)
|
||||
3 -p(I-1,J+1,K)+p(I-1,J-1,K))
|
||||
4 +b(I,J,K,2)*(p(I,J+1,K+1)-p(I,J-1,K+1)
|
||||
5 -p(I,J+1,K-1)+p(I,J-1,K-1))
|
||||
6 +b(I,J,K,3)*(p(I+1,J,K+1)-p(I-1,J,K+1)
|
||||
7 -p(I+1,J,K-1)+p(I-1,J,K-1))
|
||||
8 +c(I,J,K,1)*p(I-1,J,K)+c(I,J,K,2)*p(I,J-1,K)
|
||||
9 +c(I,J,K,3)*p(I,J,K-1)+wrk1(I,J,K)
|
||||
SS=(S0*a(I,J,K,4)-p(I,J,K))*bnd(I,J,K)
|
||||
GOSA=GOSA+SS*SS
|
||||
wrk2(I,J,K)=p(I,J,K)+OMEGA *SS
|
||||
enddo
|
||||
enddo
|
||||
enddo
|
||||
C
|
||||
DO K=2,kmax-1
|
||||
DO J=2,jmax-1
|
||||
DO I=2,imax-1
|
||||
p(I,J,K)=wrk2(I,J,K)
|
||||
enddo
|
||||
enddo
|
||||
enddo
|
||||
C
|
||||
enddo
|
||||
CC End of iteration
|
||||
return
|
||||
end
|
||||
Loading…
Reference in New Issue
Block a user