Status: Accepted
Owner: [email protected]
Labels: Type-Defect Priority-Medium WrongResult
New issue 2423 by [email protected]: Rotation.d function incorrect
evaluation
http://code.google.com/p/sympy/issues/detail?id=2423
from sympy.physics.quantum.spin import Rotation
Rotation.d(1,0,1,pi/2)
2**(1/2)/2
Rotation.d(1,0,1,-pi/2)
2**(1/2)/2
Rotation.d(1,1,0,pi/2)
-2**(1/2)/2
Rotation.d(1,1,0,-pi/2)
-2**(1/2)/2
This is the wrong result for the small d-matrix elements, namely we should
have Rotation.d(1,0,1,pi/2) == -2**(1/2)/2 and Rotation.d(1,1,0,pi/2) ==
2**(1/2)/2.
The Rotation.d function uses Varshalovich 4.3.2 Eq 7, and as far as I can
tell, it implements the equation properly. Running this equation, in
addition to the other three differential relations given by Varshalovich,
through Mathematica gives the same result. However, disagrees with the
matrix elements given in Table 4.4, which gives these elements as
Rotation.d(1,1,0,beta)=-sin(beta)/sqrt(2) and
Rotation.d(1,0,1,beta)=sin(beta)/sqrt(2).
--
You received this message because you are subscribed to the Google Groups
"sympy-issues" group.
To post to this group, send email to [email protected].
To unsubscribe from this group, send email to
[email protected].
For more options, visit this group at
http://groups.google.com/group/sympy-issues?hl=en.
