def zern_rho_coefs(n, m):
"""Compute coefficients of radial part of Zernike (n, m).
"""
from galsim.utilities import nCr
kmax = (n - abs(m)) // 2
A = np.zeros(n + 1)
for k in range(kmax + 1):
A[n - 2 * k] = (-1)**k * nCr(n - k, k) * nCr(n - 2 * k, kmax - k)
return A
def zern_rho_coefs(n, m):
"""Compute coefficients of radial part of Zernike (n, m).
"""
from galsim.utilities import nCr
kmax = (n-abs(m)) // 2
A = np.zeros(n+1)
for k in range(kmax+1):
A[n-2*k] = (-1)**k * nCr(n-k, k) * nCr(n-2*k, kmax-k)
return A