Skip to content

Faster and more accurate Gauss-Legendre nodes and weights

Martin Reinecke requested to merge gl_tests into ducc0

This algorithm was pointed out to me by Sigurd Naess.

For degrees <=100 I'm still using the old one (however now working with long double instead of double internally). This is also quick enough and avoids storing huge coefficient tables.

Merge request reports