almF3D
|
NT
field3D=almF3D(alm,nn)
returns the 3D recomposed field (dimensions nn*nn*nn) using alm.
field is of the form theta along absiss phi along ordinate
r along alt ordinate
EXAMPLE
nn=10;
#include "Dom/Density_utils.i"
x=span(-1,1,nn)(,-:1:nn)(,,-:1:nn);
y=span(-1,1,nn)(-:1:nn,)(,,-:1:nn);
z=span(-1,1,nn)(-:1:nn,)(-:1:nn,,);
u=1/sqrt((x-1.1)^2+(y-1.1)^2+0.1^2);
v=Projec_Halo(u,[0.,0,0],[x,y,z],1.,sz=nn*2);
alm=Falm(v,nn);
u1=almF3D(alm,nn,xyz=1,rmax=1).re;
w=where(x^2+y^2+z^2<1);
pl,(u1 -u)(w)(*);
pl,(u1)(w)(*),color=-6;
interpreted function, defined at ./harmonic.i line 143
|