mptoolbox_1.3

[CdeMills]

Hello,

I started working on mptoolbox, thanks to the version I got from the
original author. So far so good, I've added implementation for matrix
inversion using Gauss-Jordan elimination and pivoting. I also added 'hilb',
starting from Octave own version, as objects may not be used to implement a
range used by 'for' loop.

Taking the inverse of hilb(mp(20)) is a bit long (15 seconds on a Intel Mac)
but at least succeeds. The discrepancy between the product of this matrix by
its inverse and the identity matrix is 5e-48. Impressive!
Also, matrix flagged as singular are really singular.

I'll start playing with equation systems where rcond is greater than 10^30.

Regards

Pascal



--
View this message in context: 
http://octave.1599824.n4.nabble.com/mptoolbox-1-3-tp4654871.html
Sent from the Octave - Maintainers mailing list archive at Nabble.com.