[Octave-patch-tracker] [patch #9282] add minres for sparse linear systems

[Marco Caliari]

Follow-up Comment #11, patch #9282 (project octave):

You can remove the try..catch inside the iteration: if the preconditioner is
singular in the ## Initiation stage, the code will return. Cristiano used
ishermitian only in some tests, I think you can remove from the code.

I think you should change the stagnation criterion. It is almost impossible
that two consecutive iterations coincide (with ==). You should use something
like norm(resvec(k+1)-resvec(k)) <= eps * norm(resvec(k+1)).

I see you use a left preconditioner, while Matlab a split one.


    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/patch/?9282>

_______________________________________________
  Message sent via/by Savannah
  http://savannah.gnu.org/