Recent changes to 2688: %e^^A returns element-by-element exponenthttps://sourceforge.net/p/maxima/bugs/2688/2014-02-19T20:04:24Z#2688 %e^^A returns element-by-element exponent2014-02-19T20:04:24Z2014-02-19T20:04:24ZRupert Swarbrickhttps://sourceforge.net/u/rswarbrick/https://sourceforge.net7148570132b79f719f17658db20f27711749eb95<div class="markdown_content"><p>I guess there are two things here. Firstly, Maxima is giving a wrong answer. That turns out to be because of some bogus logic in SIMPNCEXPT. I'm pushing a patch that fixes that problem now.</p>
<p>The second is that Maxima isn't automatically calculating e^^A when it could. Speaking as the author of most of the revamped diag.mac, it's not all that brilliant. Basically, it does linear algebra like a 1st year undergrad, calculating Jordan normal forms and passing them around all over the place.</p>
<p>(1) We'd have to think quite hard about run time costs before enabling that sort of simplification by default.</p>
<p>(2) We would probably also want to move diag.mac from share/contrib to share and/or reimplement the functionality with a less horrible API.</p></div>%e^^A returns element-by-element exponent2014-02-10T08:40:29Z2014-02-10T08:40:29ZRobert Dodierhttps://sourceforge.net/u/userid-501686/https://sourceforge.net7a550a08c2e5b5c44ab3eec4727c37cdf1ed7815<div class="markdown_content"><p>As expected, %e^A returns the element-by-element exponent (where A is a matrix). But %e^^A also returns the element-by-element exponent. I think that's a bug; it seems like "^^" should return the matrix exponent, as calculated, for example, by mat_function(exp, A).</p></div>