|
|
Up |
|
|
  |
Author: FatemehFatemeh Date: Jul 5, 2008 02:23
Dear all;
I'm going to use rlft3.f (3D Fast Fourier transform subroutine in
Numerical Recipes book )
I have constructed a 107*107*107 mesh for data.
I should call rlft3(data,speg,nn1,nn2,nn3,isign).
-----------------------------------------------------------------------------
SUBROUTINE rlft3(data,speq,nn1,nn2,nn3,isign)
INTEGER isign,nn1,nn2,nn3
COMPLEX data(nn1/2,nn2,nn3),speq(nn2,nn3)
-----------------------------------------------------------------------------
1- I don't understand why for data is written nn1/2 ? I have a
107*107*107 mesh!!!
2- Is it true to put nn1=nn2=nn3=107
3- I don't know what should I put instead of speq ?
I'm anxiously looking forward your reply and guidlines.
With my best regards,
Fatemeh
|
| |
|
| | 1 Comment |
|
  |
|
|
  |
Author: DamianDamian Date: Jul 5, 2008 15:18
On Jul 5, 2:23 am, Fatemeh gmail.com> wrote:
> Dear all;
>
> I'm going to use rlft3.f (3D Fast Fourier transform subroutine in
> Numerical Recipes book )
> I have constructed a 107*107*107 mesh for data.
> I should call rlft3(data,speg,nn1,nn2,nn3,isign).
> -----------------------------------------------------------------------------
> SUBROUTINE rlft3(data,speq,nn1,nn2,nn3,isign)
> INTEGER isign,nn1,nn2,nn3
> COMPLEX data(nn1/2,nn2,nn3),speq(nn2,nn3)
> -----------------------------------------------------------------------------
> 1- I don't understand why for data is written nn1/2 ? I have a
> 107*107*107 mesh!!!
> 2- Is it true to put nn1=nn2=nn3=107
> 3- I don't know what should I put instead of speq ?
>
> I'm anxiously looking forward your reply and guidlines.
> With my best regards,
> Fatemeh ...
|
| Show full article (1.79Kb) |
|
| | no comments |
|
|