To integrate the function at each kink the integration is restarted the starting value for this integration 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). the weights have to be multiplied before calling this subroutine.
T h i s
s u b r o u t i n e
d o e s
a n
o u t w a r d s
i n t e g r a t i o n
o f
a
f u n c t i o n
w i t h
k i n k s
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=dp), | intent(in), | dimension (*) | :: | f | ||
real(kind=dp), | intent(out), | dimension (*) | :: | fint | ||
integer, | intent(in) | :: | ipan |
Number of panels in non-MT-region |
||
integer, | intent(in), | dimension (0:ipand) | :: | ircut |
R points of panel |