From 1d148f84bbf9f5ccdeee8168e4c82bb7722e82e2 Mon Sep 17 00:00:00 2001 From: Anthony Scemama Date: Sun, 21 Mar 2021 21:06:41 +0100 Subject: [PATCH] erf -> derf --- src/utils/integration.irp.f | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/src/utils/integration.irp.f b/src/utils/integration.irp.f index c907e425..974417b1 100644 --- a/src/utils/integration.irp.f +++ b/src/utils/integration.irp.f @@ -447,7 +447,7 @@ double precision function rint(n,rho) else u_inv=1.d0/dsqrt(rho) u=rho*u_inv - rint=0.5d0*u_inv*sqpi*erf(u) + rint=0.5d0*u_inv*sqpi*derf(u) endif return endif @@ -463,7 +463,7 @@ double precision function rint(n,rho) endif u=rho*u_inv two_rho_inv = 0.5d0*u_inv*u_inv - val0=0.5d0*u_inv*sqpi*erf(u) + val0=0.5d0*u_inv*sqpi*derf(u) rint=(val0-v)*two_rho_inv do k=2,n rint=(rint*dfloat(k+k-1)-v)*two_rho_inv @@ -496,7 +496,7 @@ double precision function rint_sum(n_pt_out,rho,d1) else u_inv=1.d0/dsqrt(rho) u=rho*u_inv - rint_sum=0.5d0*u_inv*sqpi*erf(u) *d1(0) + rint_sum=0.5d0*u_inv*sqpi*derf(u) *d1(0) endif do i=2,n_pt_out,2 @@ -515,7 +515,7 @@ double precision function rint_sum(n_pt_out,rho,d1) u_inv=1.d0/dsqrt(rho) u=rho*u_inv two_rho_inv = 0.5d0*u_inv*u_inv - val0=0.5d0*u_inv*sqpi*erf(u) + val0=0.5d0*u_inv*sqpi*derf(u) rint_sum=val0*d1(0) rint_tmp=(val0-v)*two_rho_inv di = 3.d0