[Trilinos-Users] [EXTERNAL] Re: MueLu + Belos + Epetra Runtime Error

[Prokopenko, Andrey V.]

Hi Alicia and Shieli,

Sorry for late reply. MueLuOperator for Epetra is applied through 
ApplyInverse(). MueLuOperator for Tpetra is applied throught apply().

We have an example in trilinos/packages/muelu/test/scaling/Driver.cpp. 
It works for both Epetra and Tpetra with Belos. However, we always use 
right preconditioner.

-Andrey

On 02/20/2017 09:05 AM, Klinvex, Alicia Marie wrote:
>
> Hello MueLu developers,
>
> Does the MueLu preconditioner get applied with Apply or ApplyInverse?  
> Should Shilei wrap the preconditioner in an Epetra_InvOperator before 
> passing it to Belos (so that Apply and ApplyInverse get swapped)?
>
> Thanks,
>
> Alicia
>
> *From:*Trilinos-Users [mailto:trilinos-users-bounces at trilinos.org] *On 
> Behalf Of *Lin, Shilei
> *Sent:* Saturday, February 18, 2017 5:35 PM
> *To:* trilinos-users at trilinos.org
> *Subject:* [EXTERNAL] Re: [Trilinos-Users] MueLu + Belos + Epetra 
> Runtime Error
>
> Here is the Code.
>
> ------------------------------------------------------------------------
>
> *From:*Lin, Shilei
> *Sent:* Saturday, February 18, 2017 4:33:45 PM
> *To:* trilinos-users at trilinos.org <mailto:trilinos-users at trilinos.org>
> *Subject:* MueLu + Belos + Epetra Runtime Error
>
> Dear developers / Users,
>
>
>     I am try to modify an example code from the MueLu packages. I have
>     attahed my source code. The code will compile, however, it won't
>     work when I set the Preconditioner to the Belos linear problem.
>
>     In line 152: //Problem->setLeftPrec(M);
>
>     This is the error message from runtime:
>
>     terminate called after throwing an instance of
>     'Belos::EpetraOpFailure'
>
>     what():
>      /home/shilei/Documents/Source/Trilinos/packages/belos/epetra/src/BelosEpetraAdapter.cpp:816:
>
>     Throw number = 1
>
>     Throw test that evaluated to true: info != 0
>
>     Belos::OperatorTraits::Apply (Epetra specialization): Calling the
>     Epetra_Operator object's Apply() method failed, returning a
>     nonzero error code of -1.
>
>     When I disabled preconditioner, the programm actually works.
>     Also, I can not find an example for using Belos as solver and
>     MueLu as preconditioner with Epetra.
>
> Thanks!
> Shilei Lin
>
>
>
> _______________________________________________
> Trilinos-Users mailing list
> 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/20170223/c803a778/attachment.html>