DESC
DESC copied to clipboard
Quadrature grid has too many nodes for given L, M
for n in range(3,13):
basis = ZernikePolynomial(n,n)
grid1 = QuadratureGrid(n,n,0)
grid2 = ConcentricGrid(n,n,0)
print(f"{n} {grid1.num_nodes/basis.num_modes:.2f} {grid2.num_nodes/basis.num_modes:.2f}")
3 2.80 1.00
4 3.00 1.00
5 3.14 1.00
6 3.25 1.00
7 3.33 1.00
8 3.40 1.00
9 3.45 1.00
10 3.50 1.00
11 3.54 1.00
12 3.57 1.00
We expect the concentric grid to have as many nodes as the basis has modes (which is True), and the quadrature grid to have double that. However in practice the quadrature grid seems to have 3x-3.5x as many nodes.