ttim.besselnumba ================ .. py:module:: ttim.besselnumba Attributes ---------- .. autoapisummary:: ttim.besselnumba.tiny ttim.besselnumba.c ttim.besselnumba.fac ttim.besselnumba.nrange ttim.besselnumba.a ttim.besselnumba.b ttim.besselnumba.fac ttim.besselnumba.b ttim.besselnumba.a ttim.besselnumba.gam ttim.besselnumba.binom ttim.besselnumba.afar ttim.besselnumba.fac ttim.besselnumba.bot ttim.besselnumba.fac ttim.besselnumba.psi ttim.besselnumba.psi ttim.besselnumba.a1 ttim.besselnumba.b1 ttim.besselnumba.twologhalf ttim.besselnumba.wg ttim.besselnumba.xg ttim.besselnumba.exprange ttim.besselnumba.anew ttim.besselnumba.bnew Functions --------- .. autoapisummary:: ttim.besselnumba.besselk0near ttim.besselnumba.besselk0cheb ttim.besselnumba.besselk0 ttim.besselnumba.bessells_int ttim.besselnumba.bessells_gauss ttim.besselnumba.circle_line_intersection ttim.besselnumba.bessellsv2 ttim.besselnumba.find_d1d2 ttim.besselnumba.bessells ttim.besselnumba.bessells_gauss_ho ttim.besselnumba.bessells_gauss_ho_d1d2 ttim.besselnumba.isinside ttim.besselnumba.bessellsqxqyv2 ttim.besselnumba.bessellsqxqy ttim.besselnumba.bessells_gauss_ho_qxqy_d1d2 ttim.besselnumba.lapls_int_ho ttim.besselnumba.lapls_int_ho_wdis ttim.besselnumba.lapld_int_ho_d1d2 ttim.besselnumba.lapld_int_ho ttim.besselnumba.bessells_gauss_ho_qxqy ttim.besselnumba.besselk1cheb ttim.besselnumba.besselk1 ttim.besselnumba.besselk1near ttim.besselnumba.besselldv2 ttim.besselnumba.besselld ttim.besselnumba.besselld_gauss_ho_d1d2 ttim.besselnumba.besselld_gauss_ho ttim.besselnumba.besselldqxqyv2 ttim.besselnumba.besselldqxqy ttim.besselnumba.besselld_gauss_ho_qxqy_d1d2 ttim.besselnumba.besselld_gauss_ho_qxqy ttim.besselnumba.besselldpart ttim.besselnumba.lapld_int_ho_wdis_d1d2 ttim.besselnumba.lapld_int_ho_wdis ttim.besselnumba.Fp ttim.besselnumba.bessells_int_ho ttim.besselnumba.bessells_int_ho_qxqy ttim.besselnumba.besselld_int_ho ttim.besselnumba.besselld_int_ho_qxqy ttim.besselnumba.potbeslsv ttim.besselnumba.disbeslsv ttim.besselnumba.potbesldv ttim.besselnumba.disbesldv 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:: binom .. 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:: 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:: lapls_int_ho(x, y, z1, z2, order) lapls_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:: lapls_int_ho_wdis(x, y, z1, z2, order) Note this is W andReturns Qx - iQy. .. !! 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 !! .. py:function:: Fp(x, y, z1, z2, biga, order, d1, d2, a, b, nt) .. py:function:: bessells_int_ho(x, y, z1, z2, lab, order, d1, d2, nt=20) Docs. To come here .. !! processed by numpydoc !! .. py:function:: bessells_int_ho_qxqy(x, y, z1, z2, lab, order, d1, d2) Docs. To come here .. !! processed by numpydoc !! .. py:function:: besselld_int_ho(x, y, z1, z2, lab, order, d1, d2) Docs. To come here .. !! processed by numpydoc !! .. py:function:: besselld_int_ho_qxqy(x, y, z1, z2, lab, order, d1, d2) Docs. To come here .. !! processed by numpydoc !! .. py:function:: potbeslsv(x, y, z1, z2, lab, order, ilap, naq, R=8) Potential of line-sink for use in timml. .. !! processed by numpydoc !! .. py:function:: disbeslsv(x, y, z1, z2, lab, order, ilap, naq, R=8) Disvec of line-sink for use in timml. .. !! processed by numpydoc !! .. py:function:: potbesldv(x, y, z1, z2, lab, order, ilap, naq, R=8) Potential of line-doublet for use in timml. .. !! processed by numpydoc !! .. py:function:: disbesldv(x, y, z1, z2, lab, order, ilap, naq, R=8) Disvec of line-doublet for use in timml. .. !! processed by numpydoc !!