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,
)
