[Octave-bug-tracker] [bug #61246] matrix inversions give wrong results (inv, lu, mldivide)

[Michael Leitner]

Follow-up Comment #2, bug #61246 (project octave):

These are the limitations of practical implementations of floating-point
numerical linear algebra. Yes, the problem a\b would have an exact solution,
namely x, but for practical efficiency the matrix computations are done in
finite precision, and so the result is not exact, but off by some amount
(which however is very small you have to admit). For better precision, you
would have to do the computations with higher internal precision, which is
possible, but would be vastly less efficient. There is software to do that,
but octave rather has those users who value efficiency over analytical
exactness as target group.

    _______________________________________________________

Reply to this item at:

  <https://savannah.gnu.org/bugs/?61246>

_______________________________________________
  Message sent via Savannah
  https://savannah.gnu.org/