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

[Joohwa Lee]

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

I tried to solve matrix system equations with several methods, i.e., inv,
mldivide, lu. They all give wrong results

n = 7;
a = randn(n);
x = (1:n)';
b = a * x;
r = x - a\b;

The last line can be changed to r = x - inv(a)*b;
The results should be all zeros. However, they are not zeros.
When n = 6 (or less), the code gives correct answers.

    _______________________________________________________

Reply to this item at:

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

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