int compid_phi2, compid_phi2sus;
compid_phi2 = oPhi2.computeEasy();
- compid_phi2sus = oPhi2.computeJack(obs_phi2::phi2Sus, &(Sim->LSIZE4));
+ compid_phi2sus = oPhi2.computeJack(obs_phi2::phi2Sus, &(Sim->LSIZE2));
*out << "\t" << oPhi2.getMean(compid_phi2) << "\t" << oPhi2.getErr(compid_phi2);
*out << "\t" << oPhi2.getMean(compid_phi2sus) << "\t" << oPhi2.getErr(compid_phi2sus);
{
OM->phi2 = 0;
- for (int ix = 0; ix < Sim->LSIZE4; ix++)
+ for (int ix = 0; ix < Sim->LSIZE2; ix++)
OM->phi2 += norm( Sim->conf[ix].phi );
- OM->phi2 /= Sim->LSIZE4;
+ OM->phi2 /= Sim->LSIZE2;
}
double obs_phi2::phi2Sus(vector< vector<double> > *vals, void *para) {