aboutsummaryrefslogtreecommitdiff
path: root/CarpetExtra/IDHydroToy
diff options
context:
space:
mode:
authorschnetter <>2003-09-20 11:47:00 +0000
committerschnetter <>2003-09-20 11:47:00 +0000
commitcae05e25c454165864950f3204cfb947e2a54703 (patch)
treefdee7ddac53724e9c4341d6757a5e11a3122e7cb /CarpetExtra/IDHydroToy
parent8fbb718b22bdbd5318f3d6d9674dfcacd8d4c250 (diff)
These files are Fortran 77, not Fortran 90.
darcs-hash:20030920114727-07bb3-a6a3356b90a00f50bf96ef24e7b1b98c157c5516.gz
Diffstat (limited to 'CarpetExtra/IDHydroToy')
-rw-r--r--CarpetExtra/IDHydroToy/src/erf.f13
-rw-r--r--CarpetExtra/IDHydroToy/src/gammln.f22
-rw-r--r--CarpetExtra/IDHydroToy/src/gammp.f16
-rw-r--r--CarpetExtra/IDHydroToy/src/gcf.f30
-rw-r--r--CarpetExtra/IDHydroToy/src/gser.f28
-rw-r--r--CarpetExtra/IDHydroToy/src/make.code.defn4
6 files changed, 2 insertions, 111 deletions
diff --git a/CarpetExtra/IDHydroToy/src/erf.f b/CarpetExtra/IDHydroToy/src/erf.f
deleted file mode 100644
index 30bc84520..000000000
--- a/CarpetExtra/IDHydroToy/src/erf.f
+++ /dev/null
@@ -1,13 +0,0 @@
- FUNCTION erf(x)
- implicit none
- REAL*8 erf,x
-CU USES gammp
- REAL*8 gammp
- if(x.lt.0)then
- erf=-gammp(.5d0,x**2)
- else
- erf=gammp(.5d0,x**2)
- endif
- return
- END
-C (C) Copr. 1986-92 Numerical Recipes Software t4.
diff --git a/CarpetExtra/IDHydroToy/src/gammln.f b/CarpetExtra/IDHydroToy/src/gammln.f
deleted file mode 100644
index b99fb6aa2..000000000
--- a/CarpetExtra/IDHydroToy/src/gammln.f
+++ /dev/null
@@ -1,22 +0,0 @@
- FUNCTION gammln(xx)
- implicit none
- REAL*8 gammln,xx
- INTEGER j
- DOUBLE PRECISION ser,stp,tmp,x,y,cof(6)
- SAVE cof,stp
- DATA cof,stp/76.18009172947146d0,-86.50532032941677d0,
- *24.01409824083091d0,-1.231739572450155d0,.1208650973866179d-2,
- *-.5395239384953d-5,2.5066282746310005d0/
- x=xx
- y=x
- tmp=x+5.5d0
- tmp=(x+0.5d0)*log(tmp)-tmp
- ser=1.000000000190015d0
- do 11 j=1,6
- y=y+1.d0
- ser=ser+cof(j)/y
-11 continue
- gammln=tmp+log(stp*ser/x)
- return
- END
-C (C) Copr. 1986-92 Numerical Recipes Software t4.
diff --git a/CarpetExtra/IDHydroToy/src/gammp.f b/CarpetExtra/IDHydroToy/src/gammp.f
deleted file mode 100644
index 99a938857..000000000
--- a/CarpetExtra/IDHydroToy/src/gammp.f
+++ /dev/null
@@ -1,16 +0,0 @@
- FUNCTION gammp(a,x)
- implicit none
- REAL*8 a,gammp,x
-CU USES gcf,gser
- REAL*8 gammcf,gamser,gln
- if(x.lt.0.or.a.le.0)pause 'bad arguments in gammp'
- if(x.lt.a+1)then
- call gser(gamser,a,x,gln)
- gammp=gamser
- else
- call gcf(gammcf,a,x,gln)
- gammp=1-gammcf
- endif
- return
- END
-C (C) Copr. 1986-92 Numerical Recipes Software t4.
diff --git a/CarpetExtra/IDHydroToy/src/gcf.f b/CarpetExtra/IDHydroToy/src/gcf.f
deleted file mode 100644
index 09fd58aa0..000000000
--- a/CarpetExtra/IDHydroToy/src/gcf.f
+++ /dev/null
@@ -1,30 +0,0 @@
- SUBROUTINE gcf(gammcf,a,x,gln)
- implicit none
- INTEGER ITMAX
- REAL*8 a,gammcf,gln,x,EPS,FPMIN
- PARAMETER (ITMAX=100,EPS=3.d-7,FPMIN=1.d-30)
-CU USES gammln
- INTEGER i
- REAL*8 an,b,c,d,del,h,gammln
- gln=gammln(a)
- b=x+1-a
- c=1/FPMIN
- d=1/b
- h=d
- do 11 i=1,ITMAX
- an=-i*(i-a)
- b=b+2
- d=an*d+b
- if(abs(d).lt.FPMIN)d=FPMIN
- c=b+an/c
- if(abs(c).lt.FPMIN)c=FPMIN
- d=1/d
- del=d*c
- h=h*del
- if(abs(del-1).lt.EPS)goto 1
-11 continue
- pause 'a too large, ITMAX too small in gcf'
-1 gammcf=exp(-x+a*log(x)-gln)*h
- return
- END
-C (C) Copr. 1986-92 Numerical Recipes Software t4.
diff --git a/CarpetExtra/IDHydroToy/src/gser.f b/CarpetExtra/IDHydroToy/src/gser.f
deleted file mode 100644
index fad2be01f..000000000
--- a/CarpetExtra/IDHydroToy/src/gser.f
+++ /dev/null
@@ -1,28 +0,0 @@
- SUBROUTINE gser(gamser,a,x,gln)
- implicit none
- INTEGER ITMAX
- REAL*8 a,gamser,gln,x,EPS
- PARAMETER (ITMAX=100,EPS=3.d-7)
-CU USES gammln
- INTEGER n
- REAL*8 ap,del,sum,gammln
- gln=gammln(a)
- if(x.le.0)then
- if(x.lt.0)pause 'x < 0 in gser'
- gamser=0
- return
- endif
- ap=a
- sum=1/a
- del=sum
- do 11 n=1,ITMAX
- ap=ap+1
- del=del*x/ap
- sum=sum+del
- if(abs(del).lt.abs(sum)*EPS)goto 1
-11 continue
- pause 'a too large, ITMAX too small in gser'
-1 gamser=sum*exp(-x+a*log(x)-gln)
- return
- END
-C (C) Copr. 1986-92 Numerical Recipes Software t4.
diff --git a/CarpetExtra/IDHydroToy/src/make.code.defn b/CarpetExtra/IDHydroToy/src/make.code.defn
index 3b82fb82f..a63bda2f5 100644
--- a/CarpetExtra/IDHydroToy/src/make.code.defn
+++ b/CarpetExtra/IDHydroToy/src/make.code.defn
@@ -1,8 +1,8 @@
# Main make.code.defn file for thorn IDHydroToy
-# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/IDHydroToy/src/make.code.defn,v 1.2 2001/03/19 21:30:37 eschnett Exp $
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/IDHydroToy/src/make.code.defn,v 1.3 2003/09/20 13:47:27 schnetter Exp $
# Source files in this directory
-SRCS = InitialData.F77 Startup.F77 erf.f gammln.f gammp.f gcf.f gser.f
+SRCS = InitialData.F77 Startup.F77 erf.f77 gammln.f77 gammp.f77 gcf.f77 gser.f77
# Subdirectories containing source files
SUBDIRS =