I have been using octave for a couple of years now for my university
courses as a replacement for That Other Software.
Last year we had the chance to take a closer look at the Moler-Van Loan
paper  when we studied non-linear dynamics. The exponential of a
matrix was also the basic tool we used to solve the Automatic Control
That said, I would like to contribute to the gnu-octave project by
improving the matrix functions as cited in the ideas page.
I would like to know what is the (estimated) scope of the improvement to
be made and see whether my ''some mathematics'' level is up to the
In reply to this post by Jordi Gutiérrez Hermoso-2
al., 2015.eko marren 16a 21:28(e)an, Jordi Gutiérrez Hermoso igorleak
> Can you understand any of the papers in this website?
> http://www.maths.manchester.ac.uk/~higham/NAMF/ >
> All the links are dead, but perhaps you can still find the articles
I wanted to go deep into at least one of the articles to further
responding. The one I have chosen is the one named "Computing the Matrix
Cosine" by N.J. Higham and M.I. Smith.
After compiling octave, I have written a draft of a would be cosm
function. It seems to compare well to the actual funm in octave and
This has brought up to a couple of questions about testing matrix
funcions: is there any way to systematically test a (matrix) function?
The second one is about Padé approximations in octave, right now I am
using sympy (sympy.org) for that purpose.
Either way, my plan for the gsoc is to follow the articles and try to
effectively implement and test all of them.
Finally, I provide you with the m file and new (working) links to the