Wrapper module for the calculation of the regular and irregular solutions
Note
preprocessor options:
* test_runs: uncomment the line #define test_run
to prepare for test
* test_prep: uncomment the line #define test_prep
to write out the solutions to file
C a l c u l a t i o n
o f
t h e
r e g u l a r
a n d
i r r e g u l a r
s o l u t i o n s
f o r
t h e
i m p u r i t y
c o d e
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=dp) | :: | rpanbound(0:npan) | ||||
real(kind=dp) | :: | rmesh(nrmax) |
Integration matrix from left ( CS_LC^-1 in eq. 5.53) Same from right ( CS_RC^-1 in eq. 5.54) Radial mesh point |
|||
complex(kind=dp) | :: | vll(lmsize*nvec,lmsize*nvec,nrmax) |
irr. volterra sol. reg. fredholm sol. t-matrix potential term in 5.7 Bauer, PhD |
|||
complex(kind=dp) | :: | rll(lmsize2,lmsize,nrmax) |
irr. volterra sol. reg. fredholm sol. t-matrix potential term in 5.7 Bauer, PhD |
|||
complex(kind=dp) | :: | sll(lmsize2,lmsize,nrmax) |
irr. volterra sol. reg. fredholm sol. t-matrix potential term in 5.7 Bauer, PhD |
|||
complex(kind=dp) | :: | tllp(lmsize,lmsize) |
irr. volterra sol. reg. fredholm sol. t-matrix potential term in 5.7 Bauer, PhD |
|||
integer | :: | ncheb |
number of chebyshev nodes |
|||
integer | :: | npan |
number of panels |
|||
integer | :: | lmsize |
lm-components * nspin |
|||
integer | :: | lmsize2 |
lmsize * nvec |
|||
integer | :: | lbessel | ||||
integer | :: | nrmax |
total number of rad. mesh points |
|||
integer | :: | nvec |
spinor integer: nvec=1 non-rel, nvec=2 for sra and dirac |
|||
integer | :: | jlk_index(nvec*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 |
prefactor of green function (non-rel: = kappa = sqrt e, rel: including mass correction) |
|||
character(len=1) | :: | cmoderll |
These define the op(V(r)) in the eqs. above (comment in the beginning of this subroutine) cmoderll ="1" : op( )=identity for reg. solution cmoderll ="T" : op( )=transpose in L for reg. solution cmodesll: same for irregular |
|||
character(len=1) | :: | cmodesll |
These define the op(V(r)) in the eqs. above (comment in the beginning of this subroutine) cmoderll ="1" : op( )=identity for reg. solution cmoderll ="T" : op( )=transpose in L for reg. solution cmodesll: same for irregular |
|||
character(len=1) | :: | cmodetest |
These define the op(V(r)) in the eqs. above (comment in the beginning of this subroutine) cmoderll ="1" : op( )=identity for reg. solution cmoderll ="T" : op( )=transpose in L for reg. solution cmodesll: same for irregular |
|||
integer | :: | use_sratrick1 | ||||
complex(kind=dp) | :: | alphaget(lmsize,lmsize) |
H e l p e r
r o u t i n e
h e r e
o n l y
f o r
h o s t
s i n c e
m
o
d
_
r
l
l
s
l
l
u
t
i
l
s
d o e s
n o t
e x s i s t
i n
t h e
h o s t
c o d e
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
integer | :: | nmat | ||||
complex(kind=dp) | :: | mat(nmat,nmat) |