Hi everyone, I've been working on representations of operators and kets in the position and momentum bases. So far, from what I can tell, the desired behavior is to return something like DiracDelta(x_1-x_2) for a representation of |x_1> (position ket), and something like x_1*DiracDelta(x_1-x_2) for the representation of the X operator (position operator).
My question, then, is what should the representation of a full expression look like. Lets say we execute the following lines: >> x_ket = XKet() >> x_op = XOp() >> represent(x_op*x_ket, basis=XOp()) With the current behavior (well, once that desired behavior is fully implemented), it would return >> x_1*DiracDelta(x_1-x_2)**2 because it simply returns the representations of each QExpr individually, so we get two delta functions out. Is this what we actually expect the output of representing an expression like that to be? Comments are greatly appreciated! Thanks Tomo -- 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.
