testStats_prg Program

  1. testStats.f90
  2. testStats_prg
program~~teststats_prg~~UsesGraph program~teststats_prg testStats_prg module~stats_mod stats_mod module~stats_mod->program~teststats_prg module~kinds_mod kinds_mod module~kinds_mod->program~teststats_prg module~kinds_mod->module~stats_mod
Help


Test program for stats_mod

Calls

program~~teststats_prg~~CallsGraph program~teststats_prg testStats_prg proc~testrandomnormal testRandomNormal program~teststats_prg->proc~testrandomnormal proc~testrandomuniform testRandomUniform program~teststats_prg->proc~testrandomuniform proc~testsetrandomseed testSetRandomSeed program~teststats_prg->proc~testsetrandomseed proc~teststdev testStDev program~teststats_prg->proc~teststdev proc~testmean testMean program~teststats_prg->proc~testmean interface~randomnormal randomNormal proc~testrandomnormal->interface~randomnormal proc~setrandomseed setRandomSeed proc~testrandomnormal->proc~setrandomseed proc~mean mean proc~testrandomnormal->proc~mean proc~stdev stDev proc~testrandomnormal->proc~stdev proc~testrandomuniform->proc~setrandomseed proc~testrandomuniform->proc~mean proc~testrandomuniform->proc~stdev interface~randomuniform randomUniform proc~testrandomuniform->interface~randomuniform proc~testsetrandomseed->proc~setrandomseed proc~teststdev->proc~stdev proc~testmean->proc~mean
Help

Subroutines

Source Code

testStats_prg

Subroutines

subroutine testSetRandomSeed()

Test setRandomSeed to verify operation

Arguments

None

subroutine testRandomUniform()

Test randomUniform to verify approximate distribution properties

Arguments

None

subroutine testRandomNormal()

Test randomNormal to verify approximate distribution properties

Arguments

None

subroutine testMean()

Test mean to verify operation

Arguments

None

subroutine testStDev()

Test stDev to verify operation

Arguments

None

Source Code

program testStats_prg
	!! Test program for stats_mod
	use kinds_mod
	use stats_mod
	implicit none
	
	call testSetRandomSeed
	call testRandomUniform
	call testRandomNormal
	
	call testMean
	call testStDev
	
contains

	subroutine testSetRandomSeed
		!! Test setRandomSeed to verify operation
		logical,dimension(2)::results
		
		integer,parameter::N = 10
		real(wp),dimension(N)::x,y,z
		
		call setRandomSeed(1)
		call random_number(x)
		
		call setRandomSeed(2)
		call random_number(y)
		
		call setRandomSeed(1)
		call random_number(z)
		
		results(1) = any(x==y)
		results(2) = .not.all(x==z)
		
		if( any(results) ) error stop "Failed setRandomSeed check"
	end subroutine testSetRandomSeed

	subroutine testRandomUniform
		!! Test randomUniform to verify approximate distribution properties
		logical,dimension(4)::results
		
		integer,parameter::N = 1000000
		real(wp),dimension(N)::x
		
		call setRandomSeed(1)
		x = randomUniform(N)
		
		results(1) = abs(mean(x))>5.0E-3_wp
		results(2) = abs(stDev(x)-sqrt(3.0_wp)**(-1))>5.0E-3_wp
		results(3) = any(x> 1.0_wp+epsilon(1.0_wp))
		results(4) = any(x<-1.0_wp-epsilon(1.0_wp))
		
		if( any(results) ) error stop "Failed randomUniform check"
	end subroutine testRandomUniform

	subroutine testRandomNormal
		!! Test randomNormal to verify approximate distribution properties
		logical,dimension(4)::results
		
		integer,parameter::N = 1000000
		real(wp),dimension(N)::x
		
		call setRandomSeed(1)
		x = randomNormal(N)
		
		results(1) = abs(mean(x))>5.0E-3_wp
		results(2) = abs(stdev(x)-1.0_wp)>5.0E-3_wp
		results(3) = any(x> 6.0_wp+12.0_wp*epsilon(1.0_wp))
		results(4) = any(x<-6.0_wp-12.0_wp*epsilon(1.0_wp))
		
		if( any(results) ) error stop "Failed randomNormal check"
	end subroutine testRandomNormal

	subroutine testMean
		!! Test mean to verify operation
		logical,dimension(1)::results
		
		results(1) = abs(mean([0.0_wp,1.0_wp,2.0_wp,3.0_wp,4.0_wp])-2.0_wp)<2.0_wp**4*epsilon(1.0_wp)
		
		if( .not.all(results) ) error stop "Failed mean check"
	end subroutine testMean

	subroutine testStDev
		!! Test stDev to verify operation
		logical,dimension(1)::results
		
		results(1) = abs(stDev([-1.0_wp,0.0_wp,1.0_wp])-1.0_wp)<2.0_wp**4*epsilon(1.0_wp)
		
		if( .not.all(results) ) error stop "Failed mean check"
	end subroutine testStDev

end program testStats_prg

All Programs