[Trilinos-Users] AztecOO residual vector

[Danielle Catherine Maddix]

I got the initialization working AztecOO_StatusTestResNorm restartResNorm(*Trilinos::K, *Trilinos::X,

(Epetra_Vector&) Trilinos::F[0], tol);

restartResNorm.DefineResForm(AztecOO_StatusTestResNorm::Implicit,

AztecOO_StatusTestResNorm::TwoNorm);

Solver.SetStatusTest(&restartResNorm);


I am unclear on how to get the intermediate norms though because it only returns the final one on the last call to check status.

double  GetResNormValue<https://trilinos.org/docs/r8.0/packages/aztecoo/doc/html/classAztecOO__StatusTestResNorm.html#d1c5d2fed099c119384bd7c7dc21649f> () const
        Returns the residual norm value, [$ \|r\| $] , computed in most recent call to CheckStatus.

How do we call checkStatus within the status = Solver.Iterate(numMaxIterations, tol); the inner iteration loops.
Any example of pseudocode would be helpful.

Thanks,
Danielle


________________________________
From: Duk-Soon Oh <duksoon at cims.nyu.edu>
Sent: Sunday, January 31, 2016 2:03 PM
To: Danielle Catherine Maddix
Subject: Re: [Trilinos-Users] AztecOO residual vector

Hi,

I had a similar issue previously. Please take a look at the AztecOO_StatusTest object. You may need to define your own AztecOO_StatusTest object to access the intermediate residual vectors. I hope it would be helpful.

Best,

Duk-Soon


On Sun, Jan 31, 2016 at 11:31 AM, Danielle Catherine Maddix <dcmaddix at stanford.edu<mailto:dcmaddix at stanford.edu>> wrote:

Hello,


I was wondering if it would be possible to access the residual vector at each iteration of the AztecOO linear solver.  I know that it returns the final residual.  I want to pass back to our in house finite element code the ratio of the final residual and the residual at the beginning of the last restart iteration in gmres.


Thank you,

Danielle Maddix

_______________________________________________
Trilinos-Users mailing list
Trilinos-Users at trilinos.org<mailto:Trilinos-Users at trilinos.org>
https://trilinos.org/mailman/listinfo/trilinos-users


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://trilinos.org/pipermail/trilinos-users/attachments/20160201/c7334756/attachment.html>