[Trilinos-Users] FFT in Trilinos?

Jonas Thies J.Thies at rug.nl
Mon Nov 17 03:35:07 MST 2008


Hi all,

Is there an FFT somewhere in Trilinos, or an interface to fftw?

thanks,
Jonas

Heroux, Michael A wrote:
> Great.  I will make sure that the constructor of this class checks for a 
> zero pointer in the future, so this kind of mistake can be caught.
> 
> Mike
> 
> 
> On 11/11/08 10:35 AM, "andrea3.villa at mail.polimi.it" 
> <andrea3.villa at mail.polimi.it> wrote:
> 
>     Ok I've fixed the problem, I've made 2 mistakes:
> 
>     1) The CreateTranspose function call the NEW operator on the pointer
>     Rcut. Therefore the line:
>     Epetra_CrsMatrix* Rcut = new Epetra_CrsMatrix(Copy,*mapXp,1);
>     is useless and creates a memory leak.
> 
>     2) The destructor of the class Epetra_RowMatrixTransposer DELETEs the
>     memory pointed by Rcut. This creates a conflict with:
>       delete Rcut
> 
>     Thank you very much.
>     Andrea.
> 
> 
> 
> 
> 
>     Citando "Heroux, Michael A" <maherou at sandia.gov>:
> 
>     >  Hi Andrea,
>     >
>     >  Have you checked the return code on the call to CreateTranspose?
>     >
>     >  Also, you should not need the "Rcut->FillComplete();" call since
>     >  this is done for you when CreateTranspose is called.
>     >
>     >  Finally, unless you know that Tcut has exactly 1 nonzero per row,
>     >  you should not use a "1" for the third argument to the Rcut
>     >  constructor.  Instead use a "0".
>     >
>     >  If you cannot figure out the error and have a small test case, I am
>     >  happy to look at it.
>     >
>     >  Mike
>     >
>     >
>     >  On 11/8/08 4:01 AM, "andrea3.villa at mail.polimi.it"
>     >  <andrea3.villa at mail.polimi.it> wrote:
>     >
>     >  Good morning.
>     >  I have a problem using Epetra_RowMatrixTransposer. I'm running a test
>     >  code like:
>     >
>     >  Epetra_CrsMatrix* Rcut = new Epetra_CrsMatrix(Copy,*mapXp,1);
>     >  Epetra_RowMatrixTransposer Traspositore(Tcut);
>     >  Traspositore.CreateTranspose(false,Rcut);
>     >  Rcut->FillComplete();
>     >  delete Rcut;
>     >
>     >  Where Tcut is an already filled FECrsMatrix. I get several
>     >  segmentation fault problems. I've also used valgrind and it defines
>     >  the last code line as: "Invalid free() / delete / delete[]". Did I use
>     >  the Epetra_RowMatrixTransposer properly?
>     >
>     >  Thank you very much.
>     >  Andrea.
>     >
>     >
>     >
>     >  _______________________________________________
>     >  Trilinos-Users mailing list
>     >  Trilinos-Users at software.sandia.gov
>     >  http://software.sandia.gov/mailman/listinfo/trilinos-users
>     >
>     >
>     >
> 
> 
> 
> 
> 
> 
>     _______________________________________________
>     Trilinos-Users mailing list
>     Trilinos-Users at software.sandia.gov
>     http://software.sandia.gov/mailman/listinfo/trilinos-users
> 
> 
> 
> ------------------------------------------------------------------------
> 
> _______________________________________________
> 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