groundstate_gaussian
Graham Dennis
Calculate the ground state of a negative gaussian potential.
t
printf("Hello world from a named filter!\n");
V1
real Vtrap = -exp(-y*y);

V1 = -i*Vtrap;

phi

// This will be automatically normalised later phi = exp(-y*y);
Ncalc
wavefunction
// Calculate the current normalisation of the wave function. Ncalc = mod2(phi);
toplevelfilter(); normalisy();
50 50 1
wavefunction normalisation
// Correct normalisation of the wavefunction phi *= sqrt(1.0/Ncalc);
T
T = -0.5*ky*ky;
wavefunction
potential
dphi_dt = T[phi] - (i*V1)*phi;
wavefunction potential
dens
wavefunction
dens = mod2(phi);
N
normalisation
N = Ncalc;
g1R g1I g2
wavefunction
complex g1 = conj(phi(y => y)) * phi(y => yp); _SAMPLE_COMPLEX(g1); g2 = mod2(phi(y => y)) * mod2(phi(y => yp));