wigner_argument_mpi
Graham Dennis and Joe Hope
Simulation of the Wigner function for an anharmonic oscillator with the initial state being a coherent state.
real Uint_hbar_on16;
/* derived constants */ Uint_hbar_on16 = Uint_hbar/16.0;
t
W
W = 2.0/M_PI * exp(-2.0*(y*y + (x-alpha_0)*(x-alpha_0)));
damping
if (sqrt(x*x + y*y) > _max_x-width) damping = 0.0; else damping = 1.0;
50
Lx Ly Lxxx Lxxy Lxyy Lyyy
Lx = i*kx; Ly = i*ky; Lxxx = -i*kx*kx*kx; Lxxy = -i*kx*kx*ky; Lxyy = -i*kx*ky*ky; Lyyy = -i*ky*ky*ky;
main
dampConstants
real rotation = omega + Uint_hbar*(-1.0 + x*x + y*y);

dW_dt = damping * ( rotation * (x*Ly[W] - y*Lx[W]) - Uint_hbar_on16*( x*(Lxxy[W] + Lyyy[W]) - y*(Lxyy[W] + Lxxx[W]) ) );
WR WI
main
_SAMPLE_COMPLEX(W);