def lqmn(m,n,z):
"""Associated Legendre functions of the second kind, Qmn(z) and its
derivative, ``Qmn'(z)`` of order m and degree n. Returns two
arrays of size ``(m+1, n+1)`` containing ``Qmn(z)`` and ``Qmn'(z)`` for
all orders from ``0..m`` and degrees from ``0..n``.
z can be complex.
"""
if not isscalar(m) or (m<0):
raise ValueError("m must be a non-negative integer.")
if not isscalar(n) or (n<0):
raise ValueError("n must be a non-negative integer.")
if not isscalar(z):
raise ValueError("z must be scalar.")
m = int(m)
n = int(n)
# Ensure neither m nor n == 0
mm = max(1,m)
nn = max(1,n)
if iscomplex(z):
q,qd = specfun.clqmn(mm,nn,z)
else:
q,qd = specfun.lqmn(mm,nn,z)
return q[:(m+1),:(n+1)],qd[:(m+1),:(n+1)]
def lqmn(m, n, z):
"""Associated Legendre functions of the second kind, Qmn(z) and its
derivative, Qmn'(z) of order m and degree n. Returns two
arrays of size (m+1,n+1) containing Qmn(z) and Qmn'(z) for
all orders from 0..m and degrees from 0..n.
z can be complex.
"""
if not isscalar(m) or (m < 0):
raise ValueError, "m must be a non-negative integer."
if not isscalar(n) or (n < 0):
raise ValueError, "n must be a non-negative integer."
if not isscalar(z):
raise ValueError, "z must be scalar."
m = int(m)
n = int(n)
# Ensure neither m nor n == 0
mm = max(1, m)
nn = max(1, n)
if iscomplex(z):
q, qd = specfun.clqmn(mm, nn, z)
else:
q, qd = specfun.lqmn(mm, nn, z)
return q[:(m + 1), :(n + 1)], qd[:(m + 1), :(n + 1)]
def lqmn(m,n,z):
"""Associated Legendre functions of the second kind, Qmn(z) and its
derivative, Qmn'(z) of order m and degree n. Returns two
arrays of size (m+1,n+1) containing Qmn(z) and Qmn'(z) for
all orders from 0..m and degrees from 0..n.
z can be complex.
"""
if not isscalar(m) or (m<0):
raise ValueError, "m must be a non-negative integer."
if not isscalar(n) or (n<0):
raise ValueError, "n must be a non-negative integer."
if not isscalar(z):
raise ValueError, "z must be scalar."
m = int(m)
n = int(n)
if (m*n == 0):
mm = max(1,m)
nn = max(1,n)
if iscomplex(z):
q,qd = specfun.clqmn(mm,nn,z)
else:
q,qd = specfun.lqmn(mm,nn,z)
return q[:(m+1),:(n+1)],qd[:(m+1),:(n+1)]