|
|
@ -8,9 +8,9 @@
|
|
|
|
double precision HNC,EPS
|
|
|
|
double precision HNC,EPS
|
|
|
|
double precision PRB, BOUND
|
|
|
|
double precision PRB, BOUND
|
|
|
|
double precision, parameter :: infinity = 37.0d0
|
|
|
|
double precision, parameter :: infinity = 37.0d0
|
|
|
|
Cf2py integer, intent(hide), depend(rho) :: N = len(rho)
|
|
|
|
Cf2py integer, intent(hide), depend(rho) :: N = len(rho)
|
|
|
|
Cf2py depend(N) a
|
|
|
|
Cf2py depend(N) a
|
|
|
|
Cf2py depend(N) b
|
|
|
|
Cf2py depend(N) b
|
|
|
|
Cf2py integer, optional :: NDF = 0
|
|
|
|
Cf2py integer, optional :: NDF = 0
|
|
|
|
Cf2py double precision, optional :: abseps = 0.001
|
|
|
|
Cf2py double precision, optional :: abseps = 0.001
|
|
|
|
Cf2py double precision, optional :: HNC = 0.24
|
|
|
|
Cf2py double precision, optional :: HNC = 0.24
|
|
|
@ -56,7 +56,7 @@ C D(Ndim) = D(K)
|
|
|
|
CALL MVSTUD(NDF,B,A,RHO,ABSEPS,Ndim,INFIN,D,IERC,HNC,
|
|
|
|
CALL MVSTUD(NDF,B,A,RHO,ABSEPS,Ndim,INFIN,D,IERC,HNC,
|
|
|
|
& PRB,BOUND,IFAULT)
|
|
|
|
& PRB,BOUND,IFAULT)
|
|
|
|
|
|
|
|
|
|
|
|
C CALL MVNPRD(A, B, BPD, EPS, N, INF, IERC, HINC, PROB, BOUND,
|
|
|
|
C CALL MVNPRD(A, B, BPD, EPS, N, INF, IERC, HINC, PROB, BOUND,
|
|
|
|
C * IFAULT)
|
|
|
|
C * IFAULT)
|
|
|
|
return
|
|
|
|
return
|
|
|
|
end subroutine prbnormtndpc
|
|
|
|
end subroutine prbnormtndpc
|
|
|
@ -73,9 +73,9 @@ C * IFAULT)
|
|
|
|
double precision,intent(out) :: abserr,prb
|
|
|
|
double precision,intent(out) :: abserr,prb
|
|
|
|
integer, intent(out) :: IFT
|
|
|
|
integer, intent(out) :: IFT
|
|
|
|
|
|
|
|
|
|
|
|
Cf2py integer, intent(hide), depend(rho) :: N = len(rho)
|
|
|
|
Cf2py integer, intent(hide), depend(rho) :: N = len(rho)
|
|
|
|
Cf2py depend(N) a
|
|
|
|
Cf2py depend(N) a
|
|
|
|
Cf2py depend(N) b
|
|
|
|
Cf2py depend(N) b
|
|
|
|
Cf2py double precision, optional :: abseps = 0.001
|
|
|
|
Cf2py double precision, optional :: abseps = 0.001
|
|
|
|
Cf2py double precision, optional :: releps = 0.001
|
|
|
|
Cf2py double precision, optional :: releps = 0.001
|
|
|
|
Cf2py logical, optional :: useBreakPoints =1
|
|
|
|
Cf2py logical, optional :: useBreakPoints =1
|
|
|
|