Home
Manual
Packages
Global Index
Keywords
Quick Reference
|
all functions - g
gadget_merge
|
returns the merge gadget snapshot with
a possible rotation rot=[anglexy,[angleyz],[anglezx]]
translation trans=[ , ,]
velocity translation vtrans=[ , ,]
velocity flip = 1/0 or void
interpreted function, defined at ./gadget.i line 1003
|
SEE ALSO:
|
|
gadget_read
|
gadget_read(fname) return a gadget structure
q= gadget_read("./snapshot_000")
for instance *q.pos holds the positions;
modified to be gadget 1.1 compatible
interpreted function, defined at ./gadget.i line 296
|
gadget_write
|
T
writes a gadget structure into the file
er= gadget_write(s,"test.dat");
interpreted function, defined at ./gadget.i line 425
|
gadgetcenter
|
recenters a given gadget snapshot on a given coordinate
interpreted function, defined at ./gadget.i line 509
|
SEE ALSO:
|
|
genRt
|
for a given 6d hypercube of 3 pts corelations
compute the residual 3 pt cor versus R and theta
EXAMPLE
//==== this can't work !!! n1 is not defined ????
interpreted function, defined at ./correl.i line 602
|
SEE ALSO:
|
|
genRt2D
|
for a given 4d hypercube of 3 pts corelations
compute the residual 3 pt cor versus R and theta
EXAMPLE
interpreted function, defined at ./correl.i line 642
|
SEE ALSO:
|
|
genYlm
|
NT
generates Ylm tables for latter use
(dimensions nn*nn).
interpreted function, defined at ./harmonic.i line 56
|
genrandfield1D
|
UMENT genrandfield(n);
generates Gaussian Random Field of dims [1,n]
which obeys the power spectrum given by fun (default k^-5/3)
EXAMPLE compute correlation function in 1D
include, "Eric/histo.i";
nm=512; ns=512;
z=genrandfield1D(ns,fun=PS); zr=float(fft(z,-1));
z2=float(fft(z*conj(z),-1));
zc=array(0.,ns); for(i=1;i<=nm;i++){z=genrandfield1D(ns,fun=PS); z2=z*conj(z);zc+=z2/float(nm);};
logxy,1,1;
kk=2*pi*(indgen(ns)-1)/float(ns);
plg,float(zc),kk; plg, PS(kk),kk,color="blue",width=4;
write, " power spectrum slope",(log(float(zc)(5:125))(dif)/log(indgen(ns)(5:125))(dif))(avg)
interpreted function, defined at ./randfield.i line 193
|
SEE
|
ALSO,
fft_indgen,,
fft.
|
genrandfield2D
|
UMENT genrandfield(n);
generates Gaussian Random Field of dims [2,n,n]
which obeys the power spectrum given by fun (default k^-5/3)
EXAMPLE
include, "Eric/histo.i";
ns=128;
z=genrandfield2D(ns); zr=float(fft(z,[-1,-1])); //pli,zr;
z2=float(fft(z*conj(z),[-1,-1])); // pli,z2;
zc=array(0.,ns,ns); for(i=1;i<=50;i++){z=genrandfield2D(ns); z2=z*conj(z);zc+=float(z2)/50.;}; //fft_pli,zc;
x1=(indgen(ns)-1)*2*pi/float(ns);
y1=(indgen(ns)-1)*2*pi/float(ns);
r1= abs(x1, y1(-,));
py= histo2(r1, px, weight=zc, average=1, interp=1,binsize=r1(dif)(min));
plg, py(2:), px(2:), color="red"; //azimuthal average
plg, PS(px(2:)),px(2:),color="blue",width=4; logxy,1,1;
write," power spectrum slope",(log(py(5:25))(dif)/log(px(5:25))(dif))(2:12:10)(avg)
interpreted function, defined at ./randfield.i line 140
|
SEE
|
ALSO,
fft_indgen,,
fft.
|
grid_read
|
grid_read(fname) return a grid structure
interpreted function, defined at ./gadget.i line 192
|
|