This subroutine does an inwards integration of llmax functions f with an extended 3-point-simpson :
The starting value for this integration at ist - 1 is determined by a 4 point lagrangian integration, coefficients given by m. abramowitz and i.a. stegun, handbook of mathematical functions, nbs applied mathematics series 55 (1968)
Warning
in case of radial integration : the weights drdi have to be multiplied before calling this subroutine. B. Drittler Mar. 1989
Modified for functions with kinks - at each kink the integration is restarted
Warning
It is supposed that is less than imt! B. Drittler july 1989
I n w a r d
i n t e g r a t i o n
o f
l l m a x
f u n c t i o n s
w i t h
e x t e n d e d
p o i n t
s i m p s o n
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
double complex | :: | F(LMMSQD,IRMIND:IRMD) | ||||
double complex | :: | FINT(LMMSQD,IRMIND:IRMD) | ||||
integer | :: | LMMSQD | ||||
integer | :: | IRMIND | ||||
integer | :: | IRMD | ||||
integer | :: | IPAN | ||||
integer | :: | IRCUT(0:IPAN) |