Discrete Jacobi had Nquad = 3N, I changed this to `max(300, 3N)` to fix a problem in the wave function project. But it's a bit hacky. Best thing to do might be to just use a Gauss quadrature rule and then it will always be exact unless the user specifies the distribution.
Discrete Jacobi had Nquad = 3N, I changed this to `max(300, 3N)` to fix a problem in the wave function project. But it's a bit hacky. Best thing to do might be to just use a Gauss quadrature rule and then it will always be exact unless the user specifies the distribution.