if sigma == 0 or sigma == None: valfunc = lambda x,y,z: -x/y/z else: valfunc = lambda x,y,z: -x/y*z/(z**2+((2*math.pi/24*sigmaterm)**2)/3) elif op == "O012": valfunc = lambda x,y,z: -x/y/((2*math.pi/24*sigmaterm)**2) elif op == "O0ii": valfunc = lambda x,y,z: -x/y*2/((2*math.pi/24*sigmaterm)**2) elif op == "O01": valfunc = lambda x,y,z: -x/y*np.sqrt(2)*2/((2*math.pi/24*sigmaterm)) elif op == "O12": valfunc = lambda x,y,z: -x/y*z*2*2/((2*math.pi/24*sigmaterm)**2) print runr.getopcf(op,sigma)+","+runr.getopcf("x",sigma)+","+destdir+","+str(latT)+","+str(numBstraps)+","+str(doasym) a.fitplatval(runr.getopcf(op,sigma),runr.getopcf("x",sigma),destdir,latT,numBstraps,valfunc,doasym = doasym) elif task == "fitplatval3": doasym = 0 if taskr.istrue("asym"): doasym = 1 sigmaterm = runr.getsigma(sigma) if op == "O44": valfunc = lambda x,y,z,w: -x/y*z/(z**2+((2*math.pi/24*sigmaterm)**2)/3) elif op == "O012": valfunc = lambda x,y,z,w: -x/y/(z**2-w**2) elif op == "O0ii": valfunc = lambda x,y,z,w: -x/y*2/((2*math.pi/24*sigmaterm)**2) print runr.getopcf(op,sigma)+","+runr.getopcf("x",sigma)+","+destdir+","+str(latT)+","+str(numBstraps)+","+str(doasym)
import sys import numpy as np import os import math import os, inspect currentdir = os.path.dirname(os.path.abspath(inspect.getfile(inspect.currentframe()))) parentdir = os.path.dirname(currentdir) sys.path.insert(0, parentdir) import a sigma = 5 print "sigma: " + str(sigma) averx2val = lambda x, y, z: -x / y * 2 * (np.sqrt(3)) ** 3 / ((2 * math.pi / 24 * (sigma) / 5) ** 3) if sigma == 0: averx2val = lambda x, y, z: -x / y a.fitplatval( "/hiskp2/oehm/correlators/momentum3/A80/O0123/sigma5-8/O0123_sigma5-8", "/hiskp2/oehm/correlators/momentum3/A80/x/sigma" + str(sigma) + "/x" + str(sigma), "/hiskp2/oehm/correlators/momentum3/A80/O0123/analysis/sigma5-8/", 48, 1500, averx2val, )