[Octave-bug-tracker] [bug #53869] Residue returning an inaccurate result

[Marco Caliari]

Follow-up Comment #6, bug #53869 (project octave):

What about a left preconditioning like


## Solve for the residues.
d = max (abs (A), [], 2);
r = (diag (d) \ A) \ (B ./ d);


It seems to fix the problem, for me. Of course, it would be better to
construct A and B already in a normalized way, but I do not know the details
of the algorithm.

    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/bugs/?53869>

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