Hi, The branch of the ticket https://trac.sagemath.org/ticket/33117 has been merged in Sage 9.6.beta5, so in Sage 9.6 spherical harmonics will agree with those of SymPy, SciPy, Mathematica and Wikipedia, and will have correct derivatives. There remains the issue of simplifying some sqrt(sin(theta)^2) terms which appear for odd orders m. This is now https://trac.sagemath.org/ticket/33501.
Eric. Le mercredi 5 janvier 2022 à 09:14:23 UTC+1, Eric Gourgoulhon a écrit : > Le mercredi 5 janvier 2022 à 08:27:56 UTC+1, Eric Gourgoulhon a écrit : > >> >> Actually, the difference between the two results is essentially due to a >> different convention in the Condon-Shortley phase >> (cf. >> https://en.wikipedia.org/wiki/Spherical_harmonics#Condon%E2%80%93Shortley_phase >> ), >> which makes Sage's spherical harmonics Y_l^m differ from Wikipedia and >> Mathematica ones by a factor (-1)^m. >> The other difference in the above example is a lack of simplification of >> sqrt(sin(theta)^2). >> >> I would vote for including the Condon-Shortley phase in Sage's spherical >> harmonics, since this is standard in quantum mechanics and this would make >> Sage agree with Wikipedia and Mathematica. >> > > I've opened > https://trac.sagemath.org/ticket/33117 > for this. > > In doing so, I've noticed that current Sage's spherical harmonics > disagree with SymPy as well. > I've also found a very serious bug in the computation of derivatives of > spherical harmonics (see the ticket for details). This has not been seen > earlier probably because before https://trac.sagemath.org/ticket/25034 > (merged in Sage 9.3), spherical harmonics were basically not usable in > Sage. > > Eric. > -- You received this message because you are subscribed to the Google Groups "sage-devel" group. To unsubscribe from this group and stop receiving emails from it, send an email to sage-devel+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/sage-devel/3034012d-c5a8-4626-99a1-9944fc264234n%40googlegroups.com.
