On Thu, Mar 1, 2012 at 8:05 AM, Ronan Lamy <[email protected]> wrote: > Le jeudi 01 mars 2012 à 08:20 -0600, Matthew Rocklin a écrit : >> I've guessed that reading the code. Would it be hard for me to >> implement the missing MatPow function ? >> >> It should be easy to implement such a class. I would suggest the name >> MatExponential. We already have a MatPow for the case A**b (A a >> matrix, b a scalar). I would suggest looking at MatPow and MatMul as >> examples. You'll also have to change the __rpow__ operator in MatExpr >> (look at the surrounding functions to see what to do). This should >> allow for the syntax e**A but not exp(A). > > x**A is a confusing notation and is very rarely used with x != E. We > should only have exp(A).
I don't see any reason to not support it. Maybe we could automatically convert it to exp(log(x)*A) in this case, to make it clear, though personally I don't even see why we would have to do that. If no one uses x != E, then it won't ever be used. Aaron Meurer > > > > -- > 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.
