Miscellanous numerical methods, mostly calculus


Subroutine Description The routine at work ... Author, Creation date
angint evaluates the integral over the solid angle of a product of three spherical harmonics via clebsch gordon coefficients: tangint.f, tangint.txt Heiko Appel, 5.18.99
arbder finds the derivative of a function on a non-uniformly spaced set of points, using two-point formulas missing! Kieron Burke, 9.15.99
arbint finds the integral of a function on a non-uniformly spaced set of points, using two-point formulas missing! Kieron Burke, 9.15.99
clebschg constructs the clebsch gordan coefficient for the special case tclebschg.f, tclebschg.txt Heiko Appel, 5.18.99
derfx DERFX returns the value of the error function over x missing! Heiko Appel, 5.18.99, updated by Kieron, 8.15.99
derfxp yields the derivative d/dx erfx(x) missing! Heiko Appel, 5.18.99
dgintgrl1 integrates a function from 0 to x using two successive grids: (y=1/(1+r/r0) to scale into [1,r0/(r0+x)] and a uniform grid to scale to [0,1] tdgintgrl1.f, tdgintgrl1.txt Kieron Burke
dgintgrl2 integrates a function from x to Infinity using two successive grids: (y=1/(1+r/r0) to scale into [1,r0/(r0+x)] and a uniform grid to scale to [0,1] tdgintgrl2.f, tdgintgrl2.txt Kieron Burke
dspline is a copy of Numerical recipes spline routine, in real*8. It sets up coefficients for a spline, in array Y2. It need only be called once. missing! Adapted from Numerical Recipes, Kieron Burke, 9.15.99
dsplint is a copy of NUM REC routine splint, in real*8. It gives the value of the splined function at point x. DSPLINE must be called beforehand. missing! Adapted from Numerical Recipes, Kieron Burke, 9.15.99
eigenf constructs the eigenfunctions for a basis-set expansion on the gridpoints rho with coefficents from the corresponding eigenvectors missing! Heiko Appel, 11.23.99
extrplam takes a function defined for various values of lambda, and interpolates the last value from the two before missing! Kieron, 9.22.99
odff ODFF does an open-ended differentiation: Just like stddff, but not using the y=0 point. This is helpful for r integrations, to avoid r=infinity. missing! Kieron, 8.14.99
rdff differentiates a function of r from 0 to infinity on n points using y=1/(1+r/r0), using 6-point formulas missing! updated by Kieron, 8.14.99
rntgr2 simply integrates a function from 0 to infinity on n+1 uniformly spaced points of y=1/(1+r/r0), producing the integral at all points inbetween. missing! Kieron Burke, 1.29.99
rntgrl simply integrates a function from 0 to infinity on n+1 uniformly spaced points of y=1/(1+r/r0). trntgrl.f, trntgrl.txt Kieron Burke, 1.29.99
sntgr2 simply integrates a function from 0 to 1 on n+1 uniformly spaced points, producing the integral at all points inbetween. missing! Adapted from Numerical Recipes, 4.1.14, Kieron Burke, 1.29.99
sntgrl simply integrates a function from 0 to 1 on n+1 uniformly spaced points; n must be at least 8. tsntgrl.f, tsntgrl.txt Adapted from Numerical Recipes, 4.1.14, Kieron Burke, 1.29.99
stddff simply differentiates a function from 0 to 1 on n+1 uniformly spaced points, using 6-point formulas tstddff.f, tstddff.txt Kieron Burke, 11.22.99
stdntgr2 simply integrates a function from 0 to 1 on n+1 uniformly spaced points, producing the integral at all points inbetween. missing! ?
stdntgrl simply integrates a function from 0 to 1 on n+1 uniformly spaced points. missing! ?

[ last page | main page ]



Last Modified: January 14, 2000
Heiko Appel, email: appel@physics.rutgers.edu