Post Reply 
Precision problem with matrix inversion
07-17-2015, 08:18 AM
Post: #21
RE: Precision problem with matrix inversion
FYI, I have eventually checked from inside CAS.
b:=400e3*[1,1,1,1]; linsolve(a,b)
returns the correct answer, as well as inv(a)*b. The CAS does not support the b/a notation because the CAS parser translate / to inv, here to b*inv(a).
If you do linear algebra, I would recomment to be inside CAS, because except for a few operations (matrix * and inverse, perhaps a few more, Cyrille could confirm), HOME matrices are converted to CAS matrices and giac algorithms are called. The constraints for Cyrille inside HOME (perhaps for learning apps like the linear system solver with HOME approx data but expecting exact results at the end) are not the same as for giac inside the CAS where the user is expected to understand the underlying maths, rounding to 0 for example might explain wrong results. Since RPN users are probably mastering the maths, the conclusion is perhaps that it would be a good idea to support RPN inside CAS.
Find all posts by this user
Quote this message in a reply
Post Reply 

Messages In This Thread
RE: Precision problem with matrix inversion - parisse - 07-17-2015 08:18 AM

User(s) browsing this thread: 1 Guest(s)