Ross Barnett's
Coulomb and Bessel functions
Code
Simple driver program
Full driver program
Bessel functions
Complete package:
Details:
READMECOUL : comments on the test files.
last changed 7/7/95 (ARB)
These are FORTRAN 77 files for COUL90.FOR, SBESJY.FOR and RICBES.FOR Coulomb
functions, spherical Bessel functions and RiccatiBessel functions.
(The '90' refers to a year NOT to FORTRAN 90)
The subroutine COUL90 calculates Coulomb and spherical/cylindrical
Bessels.
The other two subroutines are special cases.
Simple version:
The driver program COUL90_S.FOR calls the test file COULSMPL.IN and
creates the output file COULSMPL.OUT which should be identical to
COULSMPL.REF
The 1st line of input is text.
Each line then contains : eta, x, xm, lrange, kfn in format 3F10.3 ,
2I5
xm is the minimum Lvalue desired (usually zero)
lrange is the number of additional Lvalues desired
kfn = 0 gives Coulomb functions
kfn = 1 gives spherical Bessel functions
kfn = 2 gives cylindrical Bessel functions
Where appropriate, calls to SBESJY and RICBES are also made.
Full test version:
The driver program COUL90_T.FOR calls the test file COULTEST.IN and
creates the output file COULTEST.OUT which should be identical to
COULTEST.REF
The 4letter key words are COUL SBES BESS and BESI and signify the
functions to be calculated, BESS being cylindrical Bessels and BESI being
fractional order cylindrical Bessels.
NX is the number of xvalues for the given orders.
NL is the number of orders for the given xvalues.
The 1st line of input is text.
Each 3 [or 4] lines then contains :
CHFN, NL, NX, eta, xm in format 1X, A4, 7X, 2I3, 2F10.2
CHFN a 4letter key word as above
NL integer Lvalues in free format
NX real xvalues in free format
[text to annotate the output] if NX is negative
xm is the minimum Lvalue desired (usually zero)
lrange is the number of additional Lvalues desired
kfn = 0 gives Coulomb functions
kfn = 1 gives spherical Bessel functions
kfn = 2 gives cylindrical Bessel functions
Where appropriate, calls to SBESJY and RICBES are also made.
The complete description of the programs is to be found in Chapter 9 of
Computational Atomic Physics ed. K Bartschat, which is called
The Calculation of Spherical Bessel and Coulomb
Functions by A R Barnett.
