def objF(x): global NofFuncCalls F = model.averaging(random.randint(0,5275500),int(round(x[0])),int(round(x[1])), x[2],x[3],x[4],x[5],x[6],x[7],x[8],x[9],x[10]) #print command #print F NofFuncCalls = NofFuncCalls + 1 if math.fmod(NofFuncCalls,5)==0.0: print NofFuncCalls fop = open("curres",'a') fop.write(str(NofFuncCalls)+ " Chi="+ str(F) + " ") fop.write(str(int(round(x[0])))+" ") fop.write(str(int(round(x[1])))+" ") fop.write(str(round(x[2],2))+" ") fop.write(str(round(x[3],2))+" ") fop.write(str(round(x[4],2))+" ") fop.write(str(round(x[5],2))+" ") fop.write(str(round(x[6],2))+" ") fop.write(str(round(x[7],2))+" ") fop.write(str(round(x[8],2))+" ") fop.write(str(round(x[9],2))+" ") fop.write(str(round(x[10],2))+"\n") fop.close() print NofFuncCalls,F.strip("\n"),round(x[0]),round(x[1]),round(x[2],2),round(x[3],2),round(x[4],2),round(x[5],2),round(x[6],2),round(x[7],2),round(x[8],2),round(x[9],2),round(x[10],2) return float(F)
def Func(n_mean,n_sigma,I_mean,I_sigma,mean_Px,sigma_Px,mean_Py,sigma_Py,pcx,pcy,Ic): global NofFuncCalls F = model.averaging(random.randint(0,5275500),n_mean,n_sigma,I_mean,I_sigma,mean_Px,sigma_Px,mean_Py,sigma_Py,pcx,pcy,Ic) #print command #print F NofFuncCalls = NofFuncCalls + 1 if math.fmod(NofFuncCalls,1000)==0.0: print NofFuncCalls fop = open("curres",'a') fop.write(str(NofFuncCalls)+"\n") fop.close() return float(F)