You cannot select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
70 lines
2.6 KiB
FortranFixed
70 lines
2.6 KiB
FortranFixed
15 years ago
|
PROGRAM TST_RIND
|
||
|
|
USE RINDMOD
|
||
|
|
*
|
||
|
|
* Test program for rindmod2007
|
||
|
|
* gfortran -W -Wall -pedantic-errors -fbounds-check -Werror -o test_rindmod.exe intmodule.f jacobmod.f swapmod.f fimod.f rindmod.f test_rindmod.f
|
||
|
|
|
||
|
|
*
|
||
|
|
DOUBLE PRECISION ABSEPS, RELEPS
|
||
|
|
INTEGER N, Nc,Nt,Nx, I, J, K, MAXPTS
|
||
|
|
PARAMETER ( N = 5, Nc=0, Nt=5, Nx=1, MAXPTS = 5000*N*N*N )
|
||
|
|
PARAMETER ( ABSEPS = 0.00005, RELEPS = 0 )
|
||
|
|
DOUBLE PRECISION CORREL(N,N), LOW(N), UP(N), BLOW(Nx,N),BUP(Nx,N)
|
||
|
|
DOUBLE PRECISION Ex(N), Vals(Nx),ERR(Nx), TERR(Nx),Xc(Nc,Nx)
|
||
|
|
INTEGER INFIN(N), INDI(N+1)
|
||
|
|
* Chen Problem
|
||
|
|
DATA ( UP(I), I=1,N) /.0, 1.5198, 1.7817, 1.4755, 1.5949/
|
||
|
|
DATA (LOW(I), I=1,N) /.0, .0 , -1.7817, -1.4755, -1.5949/
|
||
|
|
DATA (INFIN(I), I=1,N)/ 1, 2 , 1 , 1 , 0 /
|
||
|
|
INDI = (/(I,I=0,N)/)
|
||
|
|
CORREL(:,:) = 0.0d0
|
||
|
|
CORREL(1,2) = -0.707107
|
||
|
|
CORREL(2,3) = 1.d0
|
||
|
|
CORREL(2:3,4:5) = 0.5d0
|
||
|
|
CORREL(4,5) = 0.5d0
|
||
|
|
Ex(:) = 0.d0
|
||
|
|
DO I = 1, N
|
||
|
|
CORREL(I,I) = 1.0d0
|
||
|
|
DO j = I+1, N
|
||
|
|
CORREL(I,J) = 0.3
|
||
|
|
CORREL(J,I) = CORREL(I,J)
|
||
|
|
ENDDO
|
||
|
|
ENDDO
|
||
|
|
CALL setConstants(METHOD=3,ABSEPS=abseps,RELEPS=releps,
|
||
|
|
! MAXPTS=maxpts)
|
||
|
|
PRINT '('' Test of MVNDST'')'
|
||
|
|
PRINT '(12X, ''Requested Accuracy '',F8.5)', MAX(ABSEPS,RELEPS)
|
||
|
|
PRINT '('' Number of Dimensions is '',I2)', N
|
||
|
|
PRINT '('' Maximum # of Function Values is '',I7)', MAXPTS
|
||
|
|
*
|
||
|
|
DO K = 1, 3
|
||
|
|
PRINT '(/'' I Limits'')'
|
||
|
|
PRINT '(4X,''Lower Upper Lower Left of Correlation Matrix'')'
|
||
|
|
|
||
|
|
DO I = 1, N
|
||
|
|
IF ( INFIN(I) .LT. 0 ) THEN
|
||
|
|
PRINT '(I2, '' -infin infin '', 7F9.5)',
|
||
|
|
& I, ( CORREL(I,J), J = 1,I )
|
||
|
|
ELSE IF ( INFIN(I) .EQ. 0 ) THEN
|
||
|
|
PRINT '(I2, '' -infin'', F7.4, 1X, 7F9.5)',
|
||
|
|
& I, UP(I), ( CORREL(I,J), J = 1,I)
|
||
|
|
ELSE IF ( INFIN(I) .EQ. 1 ) THEN
|
||
|
|
PRINT '(I2, F7.4, '' infin '', 7F9.5)',
|
||
|
|
& I, LOW(I), ( CORREL(I,J), J = 1,I)
|
||
|
|
ELSE
|
||
|
|
PRINT '(I2, 2F7.4, 1X, 7F9.5)',
|
||
|
|
& I, LOW(I), UP(I), ( CORREL(I,J), J = 1,I )
|
||
|
|
ENDIF
|
||
|
|
END DO
|
||
|
|
BLow(1,:) = LOW
|
||
|
|
Bup(1,:) = UP
|
||
|
|
CALL RIND(VALS,ERR,TERR,CORREL,Ex,Xc,Nt,indI,BLow,Bup,INFIN)
|
||
|
|
|
||
|
|
|
||
|
|
PRINT '('' Results for: RINDD'')'
|
||
|
|
PRINT '('' Value : '', F12.6 )', VALS(1)
|
||
|
|
PRINT '('' Error Estimate : '', '' ('', F8.6'')'' )', ERR
|
||
|
|
INFIN(1) = INFIN(1) - 1
|
||
|
|
END DO
|
||
|
|
END
|