[Trilinos-Users] Interfacing fftw and trilinos

[Veltz Romain]

Mike,

I am sorry because I was not clear enough.

My basic goal is to apply a FFT to an Epetra_Vector in parallel.

Romain.


On Nov 28, 2012, at 1:57 PM, Heroux, Michael A wrote:

> Romain,
> 
> If I understand your situation correctly, you will want to write an
> adapter for the pure virtual class Epetra_Operator, which is used by all
> Trilinos linear solvers to apply the linear operator.  A reasonable
> starting point for such a class is in
> Trilinos/packages/AztecOO/src/AztecOO_Operator.[h|cpp].
> 
> Depending on the details of what you are doing you may need to reverse the
> roles of Apply() and ApplyInverse() from what AztecOO_Operator does. You
> will implement Apply() or ApplyInverse() using the FFT library and you
> will indicate that the other is not implemented.  AztecOO_Operator
> implements ApplyInverse() and specifies that Apply() is not implemented.
> In fact you should be able to use AztecOO_Operator::Apply() basically
> as-is for one of your own methods.
> 
> I hope this helps.
> 
> Mike
> 
> On 11/28/12 2:46 PM, "Veltz Romain" <romain.veltz at inria.fr> wrote:
> 
>> Hello,
>> 
>> 
>> I would like to interface a Fast Fourier Transform library with Epetra.
>> More precisely, I am using NOX and LOCA to solve integral equations which
>> involve the computations of Fourier Transform (in Matrix-Free formalism)
>> and sparse Matrix-Vector products at the same time.
>> 
>> There is an old fftw code for serial 1d FFT in Thyra but it does not fit
>> my needs. 
>> 
>> Could you give me some advice me please,
>> 
>> Thank you for your help.
>> 
>> Romain
>> 
>> 
>> 
>> _______________________________________________
>> Trilinos-Users mailing list
>> Trilinos-Users at software.sandia.gov
>> http://software.sandia.gov/mailman/listinfo/trilinos-users
>> 
> 
>