On 4/28/06, Barton Willis <willisb at unk.edu> wrote:
> Yes, Maxima has at least three packages for computing
> the matrix exponential.
Well, I wonder if we can clean up this situation a little bit.
(1) It is likely that two of the the three packages for
computing the matrix exponential are redundant
with the third. It would be best to review the packages
and move any that seem redundant to the archive.
If one package is left, we should make it autoload.
(2) SIMPNCEXPT is not careful about simplifying ^^ .
aa ^^ bb simplifies to aa ^ bb when aa is a scalar or constant,
(e.g. %e ^^ matrix => %e ^ matrix) but that's not generally correct.
We can debate whether SIMPNCEXPT ought to go to the
expense of executing the matrix exponential code,
but if not that, then SIMPNCEXPT should at least
return a valid result.
For what it's worth,
Robert Dodier