- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
I'm using CMKL 9.0 TT routines of staggered cosine transforms. But currently can not get results I want. The details are:
The matrix of finite difference of 1d three-point approximations Poisson equation with Neumann boundary conditions is:
1,-1,0...
-1,2,-1,...
.....
..... -1,2,-1
......,0,1,-1
Eigenvalues of this matrix is 2-2*cos((i-1)*pi/n), and the Eigenvectors are corresponding to staggered cosine transform. But I found in the example of TT routines d_trig_tforms_bvp.f90 and d_trig_tforms_slae.f90, the eigenvalues seems be given to (2.0D0*dsin(0.25D0*pi*(2*k-1)/n))**2. Am I wrong?
And manuals of different version provide different description to staggered cosine transforms. Should I followthe latest version released in march 2009 to use CMKL 9.0?
Thanks.
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
Link Copied
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
Thanks.
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
Yes, the first eigenvalue is zero. This is something concerning singularity of poisson equation with Neumann boundary condition. That means if u is a solution, u+c is also solution. c is a constant.
With best regards.
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
But I'm afraid that MKL_COSINE_TRANSFORM= 1 correspond to following matrix:
2,-2,0...
-1,2,-1,...
.....
..... -1,2,-1
......,0,2,-2
It's different from our matrix because the first and last row is twice. The eigenvalues and eigenvectors are different between two matrix.
It's so strange that I can use Numerical Recipes cosine transformer codes, cosft2, and get correct results, but cannot using TT routines. Our matrix is educed from three point finite difference of 1d poisson equation in staggered grid.
Best regards!
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
For staggered cosine transform, Forward transform of cosft2 (provided by Numerical Recipes) correspends with Backward transform provieded in MKL.
Another my question is how to do multi-column transform quickly. It seems that TT routines does not support multi-column transform.
Thanks all.
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Printer Friendly Page