[Trilinos-Users] [EXTERNAL] How to implement matrix-free operation?

chenliang wang hi181904665 at msn.com
Tue Jul 17 01:56:20 MDT 2012


Thanks for your reply, you make me clear about the idea. But I am still 
confused with these interface. May I ask for a more sightly detailed 
instruction about progress or a complete tiny example?
Many thanks!
Best,
Chen-Liang Wang

On 07/17/2012 10:19 AM, Heroux, Michael A wrote:
> All Trilinos solvers perform the matrix operator using a linear operator
> interface.  If you are using Epetra, you can implement your own
> matrix-free operation by providing an adapter for Epetra_Operator.  If in
> addition you want matrix coefficient access (needed if you want to use
> Trilinos preconditioners), you can provide your own adapter for
> Epetra_RowMatrix (although I would suggest starting with
> Epetra_BasicRowMatrix).
>
> Mike
>
> On 7/16/12 9:09 PM, "chenliang wang" <hi181904665 at msn.com> wrote:
>
>> Hi,there
>> Is it any possible to implement matrix-free operation to compute pcg or
>> any other linear iterative solver? I know FLENS and MTL4 can do this,
>> but we have already write full matrix version at cluster. And we want to
>> transplant this code into pc . We still want to use trilinos . Thanks in
>> advance.
>>
>>
>> Best,
>> Chen-Liang Wang
>>
>> _______________________________________________
>> Trilinos-Users mailing list
>> Trilinos-Users at software.sandia.gov
>> http://software.sandia.gov/mailman/listinfo/trilinos-users
>
>
>





More information about the Trilinos-Users mailing list