Ndata = len(yData[:, 0])
del yData
yData = empty((Ndata, Nconf))
for idata in range(Ndata):
for iconf in range(Nconf):
yData[idata, iconf] = factor
for imlt in range(Nmlt):
tmpyData, tmpxData, Dummy = input_bin_data(iFname[imlt])
if (tmpxData is tmpyData is Dummy is None):
quit()
for idata in range(Ndata):
if (xData[idata] != tmpxData[idata]):
print("\nERROR: xData is differ, exit.\n")
quit()
for iconf in range(Nconf):
if (opera[imlt] == "x"):
yData[idata, iconf] *= tmpyData[idata, iconf]
elif (opera[imlt] == "/"):
yData[idata, iconf] /= tmpyData[idata, iconf]
else:
print("\nERROR: Invalid operator '%s', exit.\n" %
opera[imlt])
quit()
del tmpyData, tmpxData
output_bin_data(oFname, yData, xData)
Rcorr[0, 0, idata, ieq, iconf] = tmpRcorr[idata, iconf]
tmpLapRcorr, xDataTmp, dummy = input_bin_data(Fname_LapRcorr[ieq])
for idata in range(Ndata):
if (xData[idata] != xDataTmp[idata]):
print("\nERROR: x-data is different, exit.\n")
quit()
for iconf in range(Nconf):
LapRcorr[0, 0, idata, ieq, iconf] = tmpLapRcorr[idata, iconf]
Pot__LO = np.empty((Ndata, Nconf))
Pot_NLO = np.empty((Ndata, Nconf))
for iconf in range(Nconf):
tmpPot__LO, tmpPot_NLO = calc_potential_NLO(
Pot[:, :, :, :, iconf], Rcorr[:, :, :, :, iconf],
LapRcorr[:, :, :, :, iconf])
if (tmpPot__LO is tmpPot_NLO is None):
quit()
for idata in range(Ndata):
Pot__LO[idata, iconf] = tmpPot__LO[0, 0, idata]
Pot_NLO[idata, iconf] = tmpPot_NLO[0, 0, idata]
output_bin_data(oFname__LO, Pot__LO, xData)
output_bin_data(oFname_NLO, Pot_NLO, xData)
for idata in range(Ndata):
mean0, err0 = make_mean_err(Pot__LO[idata, :])
mean1, err1 = make_mean_err(Pot_NLO[idata, :])
print("%lf %e %e %e %e" % (xData[idata], mean0, err0, mean1, err1))
Pot__LO = np.empty((Nch, Nch, Ndata, Nconf))
Pot_NLO = np.empty((Nch, Nch, Ndata, Nconf))
for iconf in range(Nconf):
tmpPot__LO, tmpPot_NLO = calc_potential_NLO(Pot[:,:,:,:,iconf],
Rcr[:,:,:,:,iconf],
LpR[:,:,:,:,iconf])
if (tmpPot__LO is tmpPot_NLO is None):
quit()
for idata in range(Ndata):
for i in range(Nch):
for j in range(Nch):
Pot__LO[i, j, idata, iconf] = tmpPot__LO[i, j, idata]
Pot_NLO[i, j, idata, iconf] = tmpPot_NLO[i, j, idata]
output_bin_data(oFbase+"/LN-LN_CCP__LO_t%d"%(time), Pot__LO[0, 0, :, :], xData)
output_bin_data(oFbase+"/LN-LN_CCP_NLO_t%d"%(time), Pot_NLO[0, 0, :, :], xData)
output_bin_data(oFbase+"/LN-SN_CCP__LO_t%d"%(time), Pot__LO[0, 1, :, :], xData)
output_bin_data(oFbase+"/LN-SN_CCP_NLO_t%d"%(time), Pot_NLO[0, 1, :, :], xData)
output_bin_data(oFbase+"/SN-LN_CCP__LO_t%d"%(time), Pot__LO[1, 0, :, :], xData)
output_bin_data(oFbase+"/SN-LN_CCP_NLO_t%d"%(time), Pot_NLO[1, 0, :, :], xData)
output_bin_data(oFbase+"/SN-SN_CCP__LO_t%d"%(time), Pot__LO[1, 1, :, :], xData)
output_bin_data(oFbase+"/SN-SN_CCP_NLO_t%d"%(time), Pot_NLO[1, 1, :, :], xData)
for i in range(Nch):
for j in range(Nch):
for idata in range(Ndata):
mean0, err0 = make_mean_err(Pot__LO[i, j, idata, :])
mean1, err1 = make_mean_err(Pot_NLO[i, j, idata, :])
print("%lf %e %e %e %e" % (xData[idata], mean0, err0, mean1, err1))
print("\n")
for i in range(Ndata):
if (xData1[i] != xData2[i] or
xData1[i] != xData3[i] or
xData1[i] != xData4[i] or
xData2[i] != xData3[i] or
xData2[i] != xData4[i] or
xData3[i] != xData4[i]):
print("\nERROR: x-data is differ, exit.\n"); quit()
for iconf in range(Nconf):
V_NLO[i, iconf] = FitFunc(xData1[i], *Params[:, iconf]) / (yData3[i, iconf] - yData4[i, iconf])
tmpV_LO1[iconf] = yData1[i, iconf] - yData3[i, iconf] * V_NLO[i, iconf]
tmpV_LO2[iconf] = yData2[i, iconf] - yData4[i, iconf] * V_NLO[i, iconf]
# Selection:
V__LO[i, iconf] = tmpV_LO1[iconf]
#V__LO[i, iconf] = tmpV_LO2[iconf]
#V__LO[i, iconf] = (tmpV_LO1[iconf] + tmpV_LO2[iconf]) / 2.0
m__LO, e__LO = make_mean_err(V__LO[i, :])
m_NLO, e_NLO = make_mean_err(V_NLO[i, :])
# For Debug:
#m__LO, e__LO = make_mean_err(tmpV_LO1)
#m_NLO, e_NLO = make_mean_err(tmpV_LO2)
print("%lf %1.16e %1.16e %1.16e %1.16e" % (xData1[i], m__LO, e__LO, m_NLO, e_NLO))
output_bin_data(oFname__LO, V__LO, xData1)
output_bin_data(oFname_NLO, V_NLO, xData1)