Re: octave can not solve x=A\B, but freemat can

[tomdylan]

Thank you very much for all the hints, especially the techniques manual!

I reviewed my mixed FE code and found out that the number of elements which
is calculated from the number of nodes was 48.5. I rounded that off to the
next integer ("floor").  And this caused in octave some in building up the
coefficient matrix confusion that led to a kind of singular system.

Apparently freemat did not bother about that "wrong" number of elements and
yielded correct results!

I changed my discretization (uneven number of nodes) and octave gave my
accurate results too!

Without having found out, what exactly happened in octave through the
"wrong" number of elements I am very lucky, to have a "running" octave
system. As I assumed octave selects the appropriate matrix factorization
depending on the properties of the matrix itself.

Hope you too have a spring sun shining through your window!




--
View this message in context: 
http://octave.1599824.n4.nabble.com/octave-can-not-solve-x-A-B-but-freemat-can-tp4663428p4663442.html
Sent from the Octave - General mailing list archive at Nabble.com.