double ElmanNetwork::CalculateYNeuron (int i) {
double res = 0;
for (int j = 0; j < mHS; j++) {
res += mHL [j] -> getValue (i);
}
return mf2 (res);
}
complex EWSMcache::Bf_Mz2_0_mf2_mf2(const Particle f) const
{
int NumPar = 2;
double params[] = {SM.getMz(), mf(f, SM.getMz())};
if (params[1] == 0.0)
throw std::runtime_error("Error in EWSMcache::Bf_Mz_0_mf_mf()");
int ind = f.getIndex();
if (CacheCheck(Bf_Mz2_0_mf2_mf2_cache[ind], NumPar, params))
return complex(Bf_Mz2_0_mf2_mf2_cache[ind][NumPar],
Bf_Mz2_0_mf2_mf2_cache[ind][NumPar + 1], false);
else {
complex newResult = PV.Bf(SM.getMz() * SM.getMz(), 0.0, mf2(f, SM.getMz()), mf2(f, SM.getMz()));
newCacheForComplex(Bf_Mz2_0_mf2_mf2_cache[ind], NumPar, params, newResult);
return newResult;
}
}
int main(int argc, char *argv[])
{
std::string line;
std::ifstream mf("inp.log");
std::vector<CIWaveFunc*> roots = ParseCIRoots(mf);
//mf.seekg(mf.beg);
std::ifstream mf2("inp.log");
std::vector<NatOrbOcc*> occs = ParseNatOrbOcc(mf2);
// //std::cout<<roots.size()<<std::endl;
// //std::cout<<GetCIMaxmimumWeight(roots)<<std::endl;
//// for(int q=0; q<roots.size(); q++)
//// {
//// float avr = GetConfAverageWeight(roots.at(q));
//// float max = GetConfMaximumWeight(roots.at(q));
//// CIWaveFunc* root = roots.at(q);
//// float enr = root->energy - roots.at(0)->energy;
//// for(int i=0; i<root->conf.size();i++)
//// {
//// CIWaveFuncConf* conf = root->conf.at(i);
//// if(conf->conf_weight>max*0.75)
//// {
//// std::cout<<GetConfSpaceType(conf->conf_space)<<"("<<conf->conf_weight<<")";
//// }
//// }
//// std::cout<<" "<<enr*220000<<std::endl;
//// }
// for(int i=0; i<occs.size(); i++)
// AnalyzeOcc(occs,i);
AnalyzeFDStates(roots,0);
}
Obj(int x = 0) {
f();
mf();
mf2();
mf3();
if (x == 0){
std::cout << "def" << std::endl;
return; // this is the thing of this test
}
else{
std::cout << "non def" << std::endl;
}
}
complex EWSMcache::Bfp_Mz2_Mz2_mf2_mf2(const Particle f) const
{
int NumPar = 2;
double params[] = {SM.getMz(), mf(f, SM.getMz())};
int ind = f.getIndex();
if (CacheCheck(Bfp_Mz2_Mz2_mf2_mf2_cache[ind], NumPar, params))
return complex(Bfp_Mz2_Mz2_mf2_mf2_cache[ind][NumPar],
Bfp_Mz2_Mz2_mf2_mf2_cache[ind][NumPar + 1], false);
else {
complex newResult = PV.Bfp(SM.getMz() * SM.getMz(), SM.getMz() * SM.getMz(), mf2(f, SM.getMz()), mf2(f, SM.getMz()));
newCacheForComplex(Bfp_Mz2_Mz2_mf2_mf2_cache[ind], NumPar, params, newResult);
return newResult;
}
}