diff options
author | schnetter <> | 2003-09-20 11:47:00 +0000 |
---|---|---|
committer | schnetter <> | 2003-09-20 11:47:00 +0000 |
commit | cae05e25c454165864950f3204cfb947e2a54703 (patch) | |
tree | fdee7ddac53724e9c4341d6757a5e11a3122e7cb /CarpetExtra/IDHydroToy | |
parent | 8fbb718b22bdbd5318f3d6d9674dfcacd8d4c250 (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.f | 13 | ||||
-rw-r--r-- | CarpetExtra/IDHydroToy/src/gammln.f | 22 | ||||
-rw-r--r-- | CarpetExtra/IDHydroToy/src/gammp.f | 16 | ||||
-rw-r--r-- | CarpetExtra/IDHydroToy/src/gcf.f | 30 | ||||
-rw-r--r-- | CarpetExtra/IDHydroToy/src/gser.f | 28 | ||||
-rw-r--r-- | CarpetExtra/IDHydroToy/src/make.code.defn | 4 |
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 = |