[Trilinos-Users] PETSc interface

Lucie Freret lfreret at arrow.utias.utoronto.ca
Wed Apr 30 15:52:02 MDT 2014

Hello Jed,

Thanks a lot for your answer.

I have already configured PETSc with --download-ml and run a test with  
-pc_type ml but I'm not fully satisfied with this option.

Instead of configuring PETSc with ML, I would like to use Trilinos  
with PETSc interface. The example given in EpetraExt package deals  
with a linear system using KSPsolve with ML as a preconditioner, which  
is exactly what I'm looking for. For the time being I just want to  
compile the test example, but unfortunately I'm getting some errors  
when compiling it.


Quoting "Jed Brown" <jed at jedbrown.org>:

> Lucie Freret <lfreret at arrow.utias.utoronto.ca> writes:
>> Hello All,
>> I try to apply a Trilinos preconditioner (ML) to a PETSc linear system.
>> As a first test, I consider the example given in EpetraExt package:
>> http://trilinos.sandia.gov/packages/docs/dev/packages/epetraext/doc/html/epetraext_petsc_cpp.html
>> Unfortunately, I can't compile it and get some errors (see below).
>> I have tried with both petsc-3.4.4 and petsc-3.1.3 (and trilinos 11.4)
>> It is said here
>> (http://trilinos.sandia.gov/packages/docs/dev/packages/epetraext/doc/html/epetraext_petsc_interface.html) that the Epetra_PETScAIJMatrix class was first released with Trilinos 9.0  
>> and
>> PETSc-2.3.3
>> Does anyone already tried to use Petsc interface with older versions
>> (of Trilinos and/or PETSc) ?
> I'm not sure I understand your situation.  You have an existing PETSc
> Mat that you want to solve using ML?  I suggest configuring PETSc with
> --download-ml and running with -pc_type ml.  This works with the usual
> KSPSolve; you don't need to write any special code.

This message was sent using IMP, the Internet Messaging Program.

More information about the Trilinos-Users mailing list