GPUE  v1.0
GPU Gross-Pitaevskii Equation numerical solver for Bose-Einstein condensates
observables.py File Reference

Go to the source code of this file.

Namespaces

 observables
 

Functions

def observables.kinertrum (Psi, dx, i, quOn)
 Kinetic energy spectrum = kinertrum. More...
 
def observables.kinertrum_loop (dataName, initValue, finalValue, incr)
 
def observables.dens_struct_fact (dataName, initValue, finalValue, incr)
 
def observables.energy_total (dataName, initValue, finalValue, increment)
 
def observables.energy_kinetic (dataName, initValue, finalValue, increment)
 
def observables.energy_potential (dataName, initValue, finalValue, increment)
 
def observables.ang_mom (dataName, initValue, finalValue, incr, ev_type, imgdpi)
 
def observables.expec_val_monopole (dataName, initValue, finalValue, incr)
 
def observables.expec_val_quadrupole (dataName, initValue, finalValue, incr)
 
def observables.expec_val_ (quant_name, quantity, dataName, initValue, finalValue, incr)
 

Variables

float observables.HBAR = 1.05457148e-34
 
float observables.PI = 3.141592653589793
 
 observables.prec
 
 observables.c = ConfigParser.ConfigParser()
 
 observables.xDim = int(c.getfloat('Params','xDim'))
 
 observables.yDim = int(c.getfloat('Params','yDim'))
 
 observables.gndMaxVal = int(c.getfloat('Params','gsteps'))
 
 observables.evMaxVal = int(c.getfloat('Params','esteps'))
 
 observables.incr = int(c.getfloat('Params','print_out'))
 
tuple observables.dx = (c.getfloat('Params','dx'))
 
tuple observables.dy = (c.getfloat('Params','dy'))
 
tuple observables.dkx = (c.getfloat('Params','dpx'))
 
tuple observables.dky = (c.getfloat('Params','dpy'))
 
tuple observables.dt = (c.getfloat('Params','dt'))
 
tuple observables.xMax = (c.getfloat('Params','xMax'))
 
tuple observables.yMax = (c.getfloat('Params','yMax'))
 
tuple observables.omegaZ = (c.getfloat('Params','omegaZ'))
 
tuple observables.mass = (c.getfloat('Params','Mass'))
 
tuple observables.omega = (c.getfloat('Params','omega'))
 
tuple observables.omegaX = (c.getfloat('Params','omegaX'))
 
 observables.num_vort = int(c.getfloat('Params','Num_vort'))
 
 observables.N = int(c.getfloat('Params','atoms'))
 
 observables.data = numpy.ndarray(shape=(xDim,yDim))
 
 observables.x = np.asarray(open('x_0').read().splitlines(),dtype='f8')
 
 observables.y = np.asarray(open('y_0').read().splitlines(),dtype='f8')
 
 observables.kx = np.reshape( np.array( [np.linspace( 0, (xDim/2-1)*dkx, xDim/2), np.linspace( (-xDim/2-1)*dkx, -dkx, xDim/2)]), (xDim,1) )
 
 observables.ky = np.reshape( np.array( [np.linspace( 0, (yDim/2-1)*dky, yDim/2), np.linspace( (-yDim/2-1)*dky, -dky, yDim/2)]), (yDim,1) )
 
 observables.kxm
 
 observables.kym
 
 observables.km_mag = np.sqrt( kxm**2 + kym**2 )
 
 observables.k_mag = np.sqrt( kx**2 + ky**2 )
 
 observables.kMax = max(max(k_mag))
 
float observables.hbar = 1.05457e-34
 
float observables.m = 1.4431607e-25
 
 observables.V = np.array(open('V_0').read().splitlines(),dtype='f8')
 
 observables.K = np.array(open('K_0').read().splitlines(),dtype='f8')
 
 observables.xPy = np.array(open('xPy_0').read().splitlines(),dtype='f8')
 
 observables.yPx = np.array(open('yPx_0').read().splitlines(),dtype='f8')
 
tuple observables.g = (0.5*N)*4.0*HBAR*HBAR*PI*(4.67e-9/mass)*np.sqrt(mass*omegaZ/(2.0*PI*HBAR))