[Trilinos-Users] FFT in Trilinos?

[Heroux, Michael A]

Jonas,

There is not one at this point.  We have certainly discussed it, but have not had concrete application driver to help get the interface right.  If this is something you are interested in discussing, please let me know (off-list)


On 11/17/08 4:35 AM, "Jonas Thies" <J.Thies at rug.nl> wrote:

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

_______________________________________________
Trilinos-Users mailing list
Trilinos-Users at software.sandia.gov
http://software.sandia.gov/mailman/listinfo/trilinos-users


-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://software.sandia.gov/mailman/private/trilinos-users/attachments/20081117/e40ba842/attachment-0001.html