On 04/03/19 13:57 -0500, Ed Smith-Rowland wrote:
This is actually PR libstdc++/86655.
Thank you for reminding me Andre.
I remove the throw for m > l and just return 0. This is also done for
sph_legendre.
This build and tests clean on x86_64-linux.
OK?
OK for trunk, thanks.
This is actually PR libstdc++/86655.
Thank you for reminding me Andre.
I remove the throw for m > l and just return 0. This is also done for
sph_legendre.
This build and tests clean on x86_64-linux.
OK?
Ed
2018-03-04 Edward Smith-Rowland <3dw...@verizon.net>
PR libstdc++/86655
The return value specified in "8.1.2 associated Legendre polynomials"
of ISO/IEC JTC 1/SC 22/WG 21 N3060 (which is identical to the
expression in the doxygen comment of the patched function) is well-
defined for m>l: it is always zero because $ P_l(x) $ is a polynomial
of degree l.
The standard do