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