Commit 057ebccb authored by Daniel Wortmann's avatar Daniel Wortmann

Added some code to extend calculation of sph. Bessel fct to lmax==0

parent 1fa7a489
......@@ -28,7 +28,14 @@ c********************************************************************
! .. Intrinsic Functions ..
INTRINSIC abs,cos,sin
IF (lmax.LT.1) CALL juDFT_error("sphbes1",calledby="sphbes")
IF (lmax==0) THEN
IF (x.GE.small) THEN
fj(0) = sin(x)/x
ELSE
fj(0) = 1 - x*x/6.* (1.-x*x/20.* (1.-x*x/42.))
END IF
RETURN
ENDIF
IF (x.LT.zero) CALL juDFT_error("sphbes2",calledby="sphbes")
xx = x*x
IF (x.GE.small) THEN
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment