Great job! So lets finish the pull request with all this, write tests and get it in. If somebody finds a mistake in sympy, we can fix it later. I think that sympy is right too.
Sent from mobile phone On Jun 20, 2011 10:49 AM, "Sean Vig" <[email protected]> wrote: > On Mon, Jun 20, 2011 at 01:54, Ondrej Certik <[email protected] >wrote: > >> Can you post here your session showing the sympy (wrong) results? It >> will make it much easier for me to reproduce what you got and see what >> is wrong. >> >> Thanks >> Ondrej > > > Here's what I did >>>> from sympy.physics.quantum.spin import * >>>> JzKet(1,1).rewrite("Jx") # This agrees with the pdf > |1,-1>/2 + 2**(1/2)*|1,0>/2 + |1,1>/2 >>>> JzKet(1,0).rewrite("Jx") # This does not agree with the pdf > 2**(1/2)*|1,-1>/2 - 2**(1/2)*|1,1>/2 >>>> JzKet(1,-1).rewrite("Jx") # This agrees with the pdf > |1,-1>/2 - 2**(1/2)*|1,0>/2 + |1,1>/2 > > That, said, I think that pdf isn't right. I checked Shankar and what it has > about the rotation operator. From that I worked out the same results as > sympy returns, which I attached (note that J in that should be S, and thus > not have the hbar). > > Sean > > -- > You received this message because you are subscribed to the Google Groups "sympy" 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?hl=en. > -- You received this message because you are subscribed to the Google Groups "sympy" 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?hl=en.
