[Trilinos-Users] [EXTERNAL] Epetra_Operator as Stratimikos preconditioner

[Nico Schlömer]

Those appear to apply to the Epetra interface. I'm trying to get the Thyra
interface to work.

When browsing through Albany, I find the comment

>> Analog of EpetraExt::ModelEvaluator::Preconditioner does not exist
>> in Thyra yet! So problem will run for now with no
>> preconditioner...

in a couple of spots (e.g., [1]). This makes me think there's is still some
work in Piro/LOCA/NOX to do to support this.

Cheers,
Nico


[1]
https://github.com/gahansen/Albany/blob/master/src/QCAD/QCADT_CoupledPoissonSchrodinger.cpp#L636


On Wed, Jun 3, 2015 at 5:09 PM Salinger, Andrew <agsalin at sandia.gov> wrote:





Nico,


Here’s some code that you can use as a working use case. Look in:
packages/nox/src-epetra/NOX_Epetra_LinearSystem_Stratimikos.C


If you have an operator that approximates the inverse of the matrix, then
search for (precMatrixSource == UserDefined_) (See Line 518).
If you have a matrix that approximates the true matrix that you want
to use to generate a preconditioner using ILU/ML, then look for
precMatrixSource == SeparateMatrix (see line 484).


This class is often called from packages/piro/src/Piro_Epetra_NOXSolver.cpp
if you want to see how it is used.
Hope this helps,
Andy




From: Nico Schlömer <nico.schloemer at gmail.com>
Date: Wednesday, June 3, 2015 at 7:24 AM
To: "trilinos-users at trilinos.org" <trilinos-users at trilinos.org>
Subject: [EXTERNAL] [Trilinos-Users] Epetra_Operator as Stratimikos
preconditioner






Hi,


I'd like to use an Epetra_Operator (which I could wrap as
Thyra::LinearOpBase) as a preconditioner for a
Thyra::LinearOpWithSolveFactoryBase.
How can I do that?


Cheers,
Nico

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