rlft3.f (3D Fast Fourier transform subroutine)???
  Home FAQ Contact Sign in
comp.lang.fortran only
 
Advanced search
POPULAR GROUPS

more...

comp.lang.fortran Profile…
 Up
rlft3.f (3D Fast Fourier transform subroutine)???         


Author: Fatemeh
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
Re: rlft3.f (3D Fast Fourier transform subroutine)???         


Author: Damian
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