PythonWaveFunctionGaussianIons

From PIMC++

def WaveFunction(self,R):
#we could make this more general by using loops but since we call it a lot and 
#python does bounds checking it would be slow down the code
    return math.exp(-alpha*(dot(R[0],R[0]) + dot(R[1],R[1]))