![]() |
(git:5ef3a49)
|
Calculates special integrals. More...
Functions/Subroutines | |
| subroutine, public | whittaker_c0a (wc, r, expa, erfa, alpha, l1, l2, n) |
| int(y^(2+l1+l2) * exp(-alpha*y*y),y=0..x) / x^(l2+1); wc(:) :: output r(:) :: coordinate expa(:) :: exp(-alpha*r(:)**2) erfa(:) :: erf(sqrt(alpha)*r(:)) alpha :: exponent l1, l2 :: L-quantum number n :: number of points | |
| subroutine, public | whittaker_ci (wc, r, expa, alpha, l, n) |
| int(y^(l+1) * exp(-alpha*y*y),y=x..infinity); | |
Variables | |
| real(kind=dp), parameter | epsilon = 1.e-2_dp |
Calculates special integrals.
| subroutine, public whittaker::whittaker_c0a | ( | real(kind=dp), dimension(n), intent(out) | wc, |
| real(kind=dp), dimension(n), intent(in) | r, | ||
| real(kind=dp), dimension(n), intent(in) | expa, | ||
| real(kind=dp), dimension(n), intent(in) | erfa, | ||
| real(kind=dp), intent(in) | alpha, | ||
| integer, intent(in) | l1, | ||
| integer, intent(in) | l2, | ||
| integer, intent(in) | n | ||
| ) |
int(y^(2+l1+l2) * exp(-alpha*y*y),y=0..x) / x^(l2+1); wc(:) :: output r(:) :: coordinate expa(:) :: exp(-alpha*r(:)**2) erfa(:) :: erf(sqrt(alpha)*r(:)) alpha :: exponent l1, l2 :: L-quantum number n :: number of points
| wc | ... |
| r | ... |
| expa | ... |
| erfa | ... |
| alpha | ... |
| l1 | ... |
| l2 | ... |
| n | ... |
Definition at line 51 of file whittaker.F.
| subroutine, public whittaker::whittaker_ci | ( | real(kind=dp), dimension(n), intent(out) | wc, |
| real(kind=dp), dimension(n), intent(in) | r, | ||
| real(kind=dp), dimension(n), intent(in) | expa, | ||
| real(kind=dp), intent(in) | alpha, | ||
| integer, intent(in) | l, | ||
| integer, intent(in) | n | ||
| ) |
int(y^(l+1) * exp(-alpha*y*y),y=x..infinity);
(-1 - 1/2 l~) 2
1/2 alpha GAMMA(1/2 l + 1, alpha x )
wc(:) :: output
r(:) :: coordinate
expa(:) :: exp(-alpha*r(:)**2)
alpha :: exponent
l :: L-quantum number
n :: number of points
| wc | ... |
| r | ... |
| expa | ... |
| alpha | ... |
| l | ... |
| n | ... |
Definition at line 339 of file whittaker.F.
| real(kind=dp), parameter whittaker::epsilon = 1.e-2_dp |
Definition at line 23 of file whittaker.F.