Homec4science
Diffusion COSOlver (improved_sugama) Browse

cosolver/basis_transformation/tests/test_coca_t4performance.f90
cd1286e2d410improved_sugama

Last Change
Raw File

test_coca_t4performance.f90
View Options

program test_coca_t4performance
use FMZM
use COCA, only: cp, T4
use iso_fortran_env, only: REAL64
implicit none
integer, parameter :: NUM_LOOPS = 1000
integer(kind=REAL64):: i_loop, count_start, count_rate, &
& count_stop
real(cp) :: rntm
type(FM) :: rslt
write(*,fmt="(a,I5,a)") "Testing T4[5,5,5,5] over ", &
& 10*NUM_LOOPS, " iterations."
call system_clock( count=count_start, count_rate=count_rate )
do i_loop = 1, 10*NUM_LOOPS
rslt = T4( 5, 5, 5, 5 )
end do
call system_clock( count=count_stop )
rntm = real(count_stop-count_start) / &
& real( count_rate * 10 * NUM_LOOPS )
write(*,fmt="(a,E10.4,a,F9.3,a)") "Execution time: ", rntm, &
& " sec. (", 0.0017/rntm ,"x speedup)"
write(*,fmt="(a,I5,a)") "Testing T4[10,10,10,10] over ", &
& NUM_LOOPS, " iterations."
call system_clock( count=count_start, count_rate=count_rate )
do i_loop = 1, NUM_LOOPS
rslt = T4( 10, 10, 10, 10 )
end do
call system_clock( count=count_stop )
rntm = real(count_stop-count_start) / &
& real( count_rate * NUM_LOOPS )
write(*,fmt="(a,E10.4,a,F9.3,a)") "Execution time: ", rntm, &
& " sec. (", 0.0126/rntm ,"x speedup)"
write(*,fmt="(a,I5,a)") "Testing T4[15,15,15,15] over ", &
& NUM_LOOPS, " iterations."
call system_clock( count=count_start, count_rate=count_rate )
do i_loop = 1, NUM_LOOPS
rslt = T4( 15, 15, 15, 15 )
end do
call system_clock( count=count_stop )
rntm = real(count_stop-count_start) / &
& real( count_rate * NUM_LOOPS )
write(*,fmt="(a,E10.4,a,F9.3,a)") "Execution time: ", rntm, &
& " sec. (", 0.037/rntm ,"x speedup)"
end program test_coca_t4performance
c4science · Help