mod_sinwk Module

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.


Subroutines

public subroutine sinwk(f, fint, ipan, ircut)

Author
B. Drittler
License
Creative Commons License
Category
numerical-tools, KKRhost

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

Read more…

Arguments

Type IntentOptional 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