ttim.besselnumba_old ==================== .. py:module:: ttim.besselnumba_old Attributes ---------- .. autoapisummary:: ttim.besselnumba_old.tiny ttim.besselnumba_old.c ttim.besselnumba_old.fac ttim.besselnumba_old.nrange ttim.besselnumba_old.a ttim.besselnumba_old.b ttim.besselnumba_old.fac ttim.besselnumba_old.b ttim.besselnumba_old.a ttim.besselnumba_old.gam ttim.besselnumba_old.afar ttim.besselnumba_old.fac ttim.besselnumba_old.bot ttim.besselnumba_old.fac ttim.besselnumba_old.psi ttim.besselnumba_old.psi ttim.besselnumba_old.a1 ttim.besselnumba_old.b1 ttim.besselnumba_old.twologhalf ttim.besselnumba_old.wg ttim.besselnumba_old.xg ttim.besselnumba_old.exprange ttim.besselnumba_old.anew ttim.besselnumba_old.bnew Functions --------- .. autoapisummary:: ttim.besselnumba_old.besselk0near ttim.besselnumba_old.besselk0cheb ttim.besselnumba_old.besselk0 ttim.besselnumba_old.bessells_int ttim.besselnumba_old.bessells_gauss ttim.besselnumba_old.bessellsuni ttim.besselnumba_old.bessellsuniv ttim.besselnumba_old.circle_line_intersection ttim.besselnumba_old.bessellsv2 ttim.besselnumba_old.find_d1d2 ttim.besselnumba_old.bessells ttim.besselnumba_old.bessells_gauss_ho ttim.besselnumba_old.bessells_gauss_ho_d1d2 ttim.besselnumba_old.isinside ttim.besselnumba_old.bessellsqxqyv2 ttim.besselnumba_old.bessellsqxqy ttim.besselnumba_old.bessells_gauss_ho_qxqy_d1d2 ttim.besselnumba_old.lapld_int_ho_d1d2 ttim.besselnumba_old.lapld_int_ho ttim.besselnumba_old.bessells_gauss_ho_qxqy ttim.besselnumba_old.besselk1cheb ttim.besselnumba_old.besselk1 ttim.besselnumba_old.besselk1near ttim.besselnumba_old.besselldv2 ttim.besselnumba_old.besselld ttim.besselnumba_old.besselld_gauss_ho_d1d2 ttim.besselnumba_old.besselld_gauss_ho ttim.besselnumba_old.besselldqxqyv2 ttim.besselnumba_old.besselldqxqy ttim.besselnumba_old.besselld_gauss_ho_qxqy_d1d2 ttim.besselnumba_old.besselld_gauss_ho_qxqy ttim.besselnumba_old.besselldpart ttim.besselnumba_old.lapld_int_ho_wdis_d1d2 ttim.besselnumba_old.lapld_int_ho_wdis Module Contents --------------- .. py:data:: tiny :value: 1e-10 .. py:data:: c .. py:data:: fac :value: 1.0 .. py:data:: nrange .. py:data:: a .. py:data:: b .. py:data:: fac .. py:data:: b .. py:data:: a .. py:data:: gam .. py:data:: afar .. py:data:: fac :value: 1.0 .. py:data:: bot .. py:data:: fac .. py:data:: psi .. py:data:: psi .. py:data:: a1 .. py:data:: b1 .. py:data:: twologhalf .. py:data:: wg .. py:data:: xg .. py:function:: besselk0near(z, Nt) besselk0near. implicit none complex(kind=8), intent(in) :: z integer, intent(in) :: Nt complex(kind=8) :: omega complex(kind=8) :: rsq, log1, term integer :: n .. !! processed by numpydoc !! .. py:function:: besselk0cheb(z, Nt) besselk0cheb. implicit none complex(kind=8), intent(in) :: z integer, intent(in) :: Nt complex(kind=8) :: omega integer :: n, n2, ts real(kind=8) :: a, b, c, A3, u complex(kind=8) :: A1, A2, cn, cnp1, cnp2, cnp3 complex(kind=8) :: z1, z2, S, T .. !! processed by numpydoc !! .. py:function:: besselk0(x, y, lab) besselk0. implicit none real(kind=8), intent(in) :: x,y complex(kind=8), intent(in) :: lab complex(kind=8) :: z, omega real(kind=8) :: cond .. !! processed by numpydoc !! .. py:data:: exprange .. py:data:: anew .. py:data:: bnew .. py:function:: bessells_int(x, y, z1, z2, lab) bessells_int. implicit none real(kind=8), intent(in) :: x,y complex(kind=8), intent(in) :: z1,z2,lab real(kind=8) :: biglab, biga, L, ang, tol complex(kind=8) :: zeta, zetabar, omega, log1, log2, term1, term2, d1minzeta, d2minzeta complex(kind=8), dimension(0:20) :: zminzbar, anew, bnew, exprange complex(kind=8), dimension(0:20,0:20) :: gamnew, gam2 complex(kind=8), dimension(0:40) :: alpha, beta, alpha2 integer :: n .. !! processed by numpydoc !! .. py:function:: bessells_gauss(x, y, z1, z2, lab) bessells_gauss. implicit none real(kind=8), intent(in) :: x,y complex(kind=8), intent(in) :: z1,z2 complex(kind=8), intent(in) :: lab complex(kind=8) :: omega integer :: n real(kind=8) :: L, x0 complex(kind=8) :: bigz, biglab .. !! processed by numpydoc !! .. py:function:: bessellsuni(x, y, z1, z2, lab) Bessellsuni. # Uniform strength implicit none real(kind=8), intent(in) :: x,y complex(kind=8), intent(in) :: z1,z2 complex(kind=8), intent(in) :: lab complex(kind=8) :: omega integer :: Nls, n real(kind=8) :: Lnear, L complex(kind=8) :: z, delz, za, zb .. !! processed by numpydoc !! .. py:function:: bessellsuniv(x, y, z1, z2, lab, rzero) Bessellsuniv. # Uniform strength implicit none real(kind=8), intent(in) :: x,y complex(kind=8), intent(in) :: z1,z2 integer, intent(in) :: nlab complex(kind=8), dimension(nlab), intent(in) :: lab complex(kind=8), dimension(nlab), intent(inout) :: omega integer :: n .. !! processed by numpydoc !! .. py:function:: circle_line_intersection(z1, z2, zc, R) circle_line_intersection. implicit none complex(kind=8), intent(in) :: z1, z2, zc real(kind=8), intent(in) :: R real(kind=8), intent(inout) :: xouta, youta, xoutb, youtb integer, intent(inout) :: N real(kind=8) :: Lover2, d, xa, xb complex(kind=8) :: bigz, za, zb .. !! processed by numpydoc !! .. py:function:: bessellsv2(x, y, z1, z2, lab, order, R) bessellsv2. implicit none integer, intent(in) :: order real(kind=8), intent(in) :: x,y,R complex(kind=8), intent(in) :: z1,z2 integer, intent(in) :: nlab real(kind=8) :: d1, d2 complex(kind=8), dimension(nlab), intent(in) :: lab complex(kind=8), dimension(order+1,nlab) :: omega integer :: n, nterms .. !! processed by numpydoc !! .. py:function:: find_d1d2(z1, z2, zc, R) find_d1d2. implicit none complex(kind=8), intent(in) :: z1, z2, zc real(kind=8), intent(in) :: R real(kind=8), intent(inout) :: d1, d2 real(kind=8) :: Lover2, d, xa, xb complex(kind=8) :: bigz .. !! processed by numpydoc !! .. py:function:: bessells(x, y, z1, z2, lab, order, d1in, d2in) Bessells. implicit none integer, intent(in) :: order real(kind=8), intent(in) :: x,y,d1in,d2in complex(kind=8), intent(in) :: z1,z2 complex(kind=8), intent(in) :: lab complex(kind=8), dimension(0:order) :: omega integer :: Nls, n real(kind=8) :: Lnear, L, d1, d2, delta complex(kind=8) :: z, delz, za, zb .. !! processed by numpydoc !! .. py:function:: bessells_gauss_ho(x, y, z1, z2, lab, order) bessells_gauss_ho. implicit none integer, intent(in) :: order real(kind=8), intent(in) :: x,y complex(kind=8), intent(in) :: z1,z2 complex(kind=8), intent(in) :: lab complex(kind=8), dimension(0:order) :: omega integer :: n, p real(kind=8) :: L, x0 complex(kind=8) :: bigz, biglab complex(kind=8), dimension(8) :: k0 .. !! processed by numpydoc !! .. py:function:: bessells_gauss_ho_d1d2(x, y, z1, z2, lab, order, d1, d2) Returns integral from d1 to d2 along real axis. While strength is still Delta^order from -1 to +1. implicit none integer, intent(in) :: order real(kind=8), intent(in) :: x,y,d1,d2 complex(kind=8), intent(in) :: z1,z2,lab complex(kind=8), dimension(0:order) :: omega, omegac integer :: n, m real(kind=8) :: xp, yp, dc, fac complex(kind=8) :: z1p,z2p,bigz1,bigz2 .. !! processed by numpydoc !! .. py:function:: isinside(z1, z2, zc, R) Checks whether point zc is within oval with 'radius' R from line element. implicit none complex(kind=8), intent(in) :: z1, z2, zc real(kind=8), intent(in) :: R integer :: irv real(kind=8) :: Lover2, d, xa, xb complex(kind=8) :: bigz .. !! processed by numpydoc !! .. py:function:: bessellsqxqyv2(x, y, z1, z2, lab, order, R) bessellsqxqyv2. implicit none integer, intent(in) :: order real(kind=8), intent(in) :: x,y,R complex(kind=8), intent(in) :: z1,z2 integer, intent(in) :: nlab real(kind=8) :: d1, d2 complex(kind=8), dimension(nlab), intent(in) :: lab complex(kind=8), dimension(2*(order+1),nlab) :: qxqy complex(kind=8), dimension(0:2*order+1) :: qxqylab integer :: n, nterms, nhalf .. !! processed by numpydoc !! .. py:function:: bessellsqxqy(x, y, z1, z2, lab, order, d1in, d2in) Bessellsqxqy. implicit none integer, intent(in) :: order real(kind=8), intent(in) :: x,y,d1in,d2in complex(kind=8), intent(in) :: z1,z2 complex(kind=8), intent(in) :: lab complex(kind=8), dimension(0:2*order+1) :: qxqy integer :: Nls, n real(kind=8) :: Lnear, L, d1, d2, delta complex(kind=8) :: z, delz, za, zb .. !! processed by numpydoc !! .. py:function:: bessells_gauss_ho_qxqy_d1d2(x, y, z1, z2, lab, order, d1, d2) Returns integral from d1 to d2 along real axis. While strength is still Delta^order from -1 to +1. implicit none integer, intent(in) :: order real(kind=8), intent(in) :: x,y,d1,d2 complex(kind=8), intent(in) :: z1,z2,lab complex(kind=8), dimension(0:2*order+1) :: qxqy, qxqyc integer :: n, m real(kind=8) :: xp, yp, dc, fac complex(kind=8) :: z1p,z2p,bigz1,bigz2 .. !! processed by numpydoc !! .. py:function:: lapld_int_ho_d1d2(x, y, z1, z2, order, d1, d2) lapld_int_ho_d1d2. Near field only Returns integral from d1 to d2 along real axis while strength is still Delta^order from -1 to +1 implicit none integer, intent(in) :: order real(kind=8), intent(in) :: x,y,d1,d2 complex(kind=8), intent(in) :: z1,z2 complex(kind=8), dimension(0:order) :: omega, omegac integer :: n, m real(kind=8) :: xp, yp, dc, fac complex(kind=8) :: z1p,z2p,bigz1,bigz2 .. !! processed by numpydoc !! .. py:function:: lapld_int_ho(x, y, z1, z2, order) lapld_int_ho. ! Near field only implicit none integer, intent(in) :: order real(kind=8), intent(in) :: x,y complex(kind=8), intent(in) :: z1,z2 complex(kind=8), dimension(0:order) :: omega, qm integer :: m, n real(kind=8) :: L complex(kind=8) :: z, zplus1, zmin1 .. !! processed by numpydoc !! .. py:function:: bessells_gauss_ho_qxqy(x, y, z1, z2, lab, order) bessells_gauss_ho_qxqy. implicit none integer, intent(in) :: order real(kind=8), intent(in) :: x,y complex(kind=8), intent(in) :: z1,z2 complex(kind=8), intent(in) :: lab complex(kind=8), dimension(0:2*order+1) :: qxqy integer :: n, p real(kind=8) :: L, bigy, angz complex(kind=8) :: bigz, biglab real(kind=8), dimension(8) :: r, xmind complex(kind=8), dimension(8) :: k1 complex(kind=8), dimension(0:order) :: qx,qy .. !! processed by numpydoc !! .. py:function:: besselk1cheb(z, Nt) besselk1cheb. implicit none complex(kind=8), intent(in) :: z integer, intent(in) :: Nt complex(kind=8) :: omega integer :: n, n2, ts real(kind=8) :: a, b, c, A3, u complex(kind=8) :: A1, A2, cn, cnp1, cnp2, cnp3 complex(kind=8) :: z1, z2, S, T .. !! processed by numpydoc !! .. py:function:: besselk1(x, y, lab) besselk1. implicit none real(kind=8), intent(in) :: x,y complex(kind=8), intent(in) :: lab complex(kind=8) :: z, omega real(kind=8) :: cond .. !! processed by numpydoc !! .. py:function:: besselk1near(z, Nt) besselk1near. implicit none complex(kind=8), intent(in) :: z integer, intent(in) :: Nt complex(kind=8) :: omega complex(kind=8) :: zsq, log1, term integer :: n .. !! processed by numpydoc !! .. py:function:: besselldv2(x, y, z1, z2, lab, order, R) besselldv2. implicit none integer, intent(in) :: order real(kind=8), intent(in) :: x,y,R complex(kind=8), intent(in) :: z1,z2 integer, intent(in) :: nlab real(kind=8) :: d1, d2 complex(kind=8), dimension(nlab), intent(in) :: lab complex(kind=8), dimension(order+1,nlab) :: omega integer :: n, nterms .. !! processed by numpydoc !! .. py:function:: besselld(x, y, z1, z2, lab, order, d1in, d2in) Besselld. implicit none integer, intent(in) :: order real(kind=8), intent(in) :: x,y,d1in,d2in complex(kind=8), intent(in) :: z1,z2 complex(kind=8), intent(in) :: lab complex(kind=8), dimension(0:order) :: omega integer :: Nls, n real(kind=8) :: Lnear, L, d1, d2, delta complex(kind=8) :: z, delz, za, zb .. !! processed by numpydoc !! .. py:function:: besselld_gauss_ho_d1d2(x, y, z1, z2, lab, order, d1, d2) besselld_gauss_ho_d1d2. # Returns integral from d1 to d2 along real axis while strength is still # Delta^order from -1 to +1 implicit none integer, intent(in) :: order real(kind=8), intent(in) :: x,y,d1,d2 complex(kind=8), intent(in) :: z1,z2,lab complex(kind=8), dimension(0:order) :: omega, omegac integer :: n, m real(kind=8) :: xp, yp, dc, fac complex(kind=8) :: z1p,z2p,bigz1,bigz2 .. !! processed by numpydoc !! .. py:function:: besselld_gauss_ho(x, y, z1, z2, lab, order) besselld_gauss_ho. implicit none integer, intent(in) :: order real(kind=8), intent(in) :: x,y complex(kind=8), intent(in) :: z1,z2 complex(kind=8), intent(in) :: lab complex(kind=8), dimension(0:order) :: omega integer :: n, p real(kind=8) :: L, x0, r complex(kind=8) :: bigz, biglab complex(kind=8), dimension(8) :: k1overr .. !! processed by numpydoc !! .. py:function:: besselldqxqyv2(x, y, z1, z2, lab, order, R) besselldqxqyv2. implicit none integer, intent(in) :: order real(kind=8), intent(in) :: x,y,R complex(kind=8), intent(in) :: z1,z2 integer, intent(in) :: nlab real(kind=8) :: d1, d2 complex(kind=8), dimension(nlab), intent(in) :: lab complex(kind=8), dimension(2*(order+1),nlab) :: qxqy complex(kind=8), dimension(0:2*order+1) :: qxqylab integer :: n, nterms, nhalf .. !! processed by numpydoc !! .. py:function:: besselldqxqy(x, y, z1, z2, lab, order, d1in, d2in) Besselldqxqy. implicit none integer, intent(in) :: order real(kind=8), intent(in) :: x,y,d1in,d2in complex(kind=8), intent(in) :: z1,z2 complex(kind=8), intent(in) :: lab complex(kind=8), dimension(0:2*order+1) :: qxqy integer :: Nls, n real(kind=8) :: Lnear, L, d1, d2, delta complex(kind=8) :: z, delz, za, zb .. !! processed by numpydoc !! .. py:function:: besselld_gauss_ho_qxqy_d1d2(x, y, z1, z2, lab, order, d1, d2) Returns integral from d1 to d2 along real axis. While strength is still Delta^order from -1 to +1. implicit none integer, intent(in) :: order real(kind=8), intent(in) :: x,y,d1,d2 complex(kind=8), intent(in) :: z1,z2,lab complex(kind=8), dimension(0:2*order+1) :: qxqy, qxqyc integer :: n, m real(kind=8) :: xp, yp, dc, fac complex(kind=8) :: z1p,z2p,bigz1,bigz2 .. !! processed by numpydoc !! .. py:function:: besselld_gauss_ho_qxqy(x, y, z1, z2, lab, order) besselld_gauss_ho_qxqy. implicit none integer, intent(in) :: order real(kind=8), intent(in) :: x,y complex(kind=8), intent(in) :: z1,z2 complex(kind=8), intent(in) :: lab complex(kind=8), dimension(0:2*order+1) :: qxqy integer :: n, p real(kind=8) :: L, bigy, angz complex(kind=8) :: bigz, biglab real(kind=8), dimension(8) :: r, xmind complex(kind=8), dimension(8) :: k0,k1 complex(kind=8), dimension(0:order) :: qx,qy .. !! processed by numpydoc !! .. py:function:: besselldpart(x, y, z1, z2, lab, order, d1, d2) Besselldpart. implicit none integer, intent(in) :: order real(kind=8), intent(in) :: x,y,d1,d2 complex(kind=8), intent(in) :: z1,z2,lab complex(kind=8), dimension(0:order) :: omega real(kind=8) :: biglab, biga, L, ang, tol, bigy complex(kind=8) :: zeta, zetabar, log1, log2, term1, term2, d1minzeta, d2minzeta, bigz complex(kind=8) :: cm, biglabcomplex complex(kind=8), dimension(0:20) :: zminzbar, anew, bnew, exprange complex(kind=8), dimension(0:20,0:20) :: gamnew, gam2 complex(kind=8), dimension(0:40) :: alpha, beta, alpha2 complex(kind=8), dimension(0:50) :: alphanew, betanew, alphanew2 ! Order fixed to 10 integer :: m, n, p .. !! processed by numpydoc !! .. py:function:: lapld_int_ho_wdis_d1d2(x, y, z1, z2, order, d1, d2) lapld_int_ho_wdis_d1d2. # Near field only # Returns integral from d1 to d2 along real axis while strength is still # Delta^order from -1 to +1 implicit none integer, intent(in) :: order real(kind=8), intent(in) :: x,y,d1,d2 complex(kind=8), intent(in) :: z1,z2 complex(kind=8), dimension(0:order) :: wdis, wdisc integer :: n, m real(kind=8) :: xp, yp, dc, fac complex(kind=8) :: z1p,z2p,bigz1,bigz2 .. !! processed by numpydoc !! .. py:function:: lapld_int_ho_wdis(x, y, z1, z2, order) lapld_int_ho_wdis. # Near field only implicit none integer, intent(in) :: order real(kind=8), intent(in) :: x,y complex(kind=8), intent(in) :: z1,z2 complex(kind=8), dimension(0:order) :: wdis complex(kind=8), dimension(0:10) :: qm # Max order is 10 integer :: m, n complex(kind=8) :: z, zplus1, zmin1, term1, term2, zterm .. !! processed by numpydoc !!