Calculation of the local regular solutions
VJLR
, VNL
, VHLR
, which appear in the integrands TAU(K,IPAN)
is used instead of TAU(K,IPAN)**2
, which directly gives RLL(r)
and SLL(r)
multiplied with r
. TAU
is the radial mesh.YR
, ZR
, YI
, ZI
because of the conventions used
by Gonzalez et al, Journal of Computational Physics 134, 134-149 (1997)
a factor is included in the source terms this factor is removed by
the definition of ZSLC1SUM
given below
i.e. prepare terms , appearing in 5.11, 5.12.Then determine the local solutions solve the equations and and and i.e., solve system , see eq. 5.68.
C a l c u l a t i o n
o f
t h e
l o c a l
r e g u l a r
s o l u t i o n s
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
complex(kind=dp) | :: | vll(lmsize*nvec,lmsize*nvec,nrmax) | ||||
real(kind=dp) | :: | tau(0:ncheb) | ||||
real(kind=dp) | :: | drpan2 | ||||
real(kind=dp) | :: | cslc1(0:ncheb,0:ncheb) | ||||
real(kind=dp) | :: | slc1sum(0:ncheb) | ||||
complex(kind=dp) | :: | mrnvy(lmsize,lmsize) | ||||
complex(kind=dp) | :: | mrnvz(lmsize,lmsize) | ||||
complex(kind=dp) | :: | mrjvy(lmsize,lmsize) | ||||
complex(kind=dp) | :: | mrjvz(lmsize,lmsize) | ||||
complex(kind=dp) | :: | yrf(lmsize2,lmsize,0:ncheb) | ||||
complex(kind=dp) | :: | zrf(lmsize2,lmsize,0:ncheb) | ||||
integer | :: | ncheb | ||||
integer | :: | ipan | ||||
integer | :: | lmsize | ||||
integer | :: | lmsize2 | ||||
integer | :: | nrmax | ||||
integer | :: | nvec | ||||
integer | :: | jlk_index(2*lmsize) | ||||
complex(kind=dp) | :: | hlk(lbessel,nrmax) | ||||
complex(kind=dp) | :: | jlk(lbessel,nrmax) | ||||
complex(kind=dp) | :: | hlk2(lbessel,nrmax) | ||||
complex(kind=dp) | :: | jlk2(lbessel,nrmax) | ||||
complex(kind=dp) | :: | gmatprefactor | ||||
character(len=1) | :: | cmoderll | ||||
integer | :: | lbessel | ||||
integer | :: | use_sratrick1 |