From 5df79c53745fde5d6c3340a2979b1429cd5892c1 Mon Sep 17 00:00:00 2001 From: Henrik Rydberg Date: Sat, 8 Oct 2011 20:30:28 +0200 Subject: Initial import of htcd system 1.0 Signed-off-by: Henrik Rydberg --- src/labat/gcor.f | 15 +++++++++++++++ 1 file changed, 15 insertions(+) create mode 100644 src/labat/gcor.f (limited to 'src/labat/gcor.f') diff --git a/src/labat/gcor.f b/src/labat/gcor.f new file mode 100644 index 0000000..afdc739 --- /dev/null +++ b/src/labat/gcor.f @@ -0,0 +1,15 @@ + subroutine gcor(a,a1,b1,b2,b3,b4,p,rs,gg,ggrs) +c called by subroutines ldaec and ldauc + implicit double precision (a-h,o-z) + p1 = p + 1.d0 + q0 = -2.d0*a*(1.d0+a1*rs) + rs12 = dsqrt(rs) + rs32 = rs12**3 + rsp = rs**p + q1 = 2.d0*a*(b1*rs12+b2*rs+b3*rs32+b4*rs*rsp) + q2 = dlog(1.d0+1.d0/q1) + gg = q0*q2 + q3 = a*(b1/rs12+2.d0*b2+3.d0*b3*rs12+2.d0*b4*p1*rsp) + ggrs = -2.d0*a*a1*q2-q0*q3/(q1**2+q1) + return + end -- cgit v1.2.3