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