Commit ef93b61d authored by Jean-Eric Campagne's avatar Jean-Eric Campagne
Browse files

(JEC) 25/4/15 polish

parent 2c2d72af
......@@ -36,42 +36,7 @@ LaguerreTransform::LaguerreTransform(int N, r_8 R, int alpha) :
int alphaFact = 1;
for(int i=1;i<=alpha_; i++) alphaFact *= i;
// //verification
// r_16 sumR = 0.;
// r_16 sumW = 0.;
// r_16 tau = nodes_[N_-1]/708.0; //double number 708 ~ ln(minimum positive double)
// for (int i=0;i<N_;i++){
// sumR += nodes_[i];
// sumW += weights_[i]*exp(-nodes_[i]/tau);
// }
// r_16 sumRTh = (r_16)(N_*(N_+alpha_));
// r_16 sumWTh = (r_16)(alphaFact)*pow(tau,(r_16)(alpha_+1));
// #if DEBUG >= 1
// cout << "Sum roots = " << sumR << " -> diff with theory: " << sumR - sumRTh << endl;
// cout << "Sum weights = " << sumW << " -> diff with theory: " << sumW - sumWTh << endl;
// #endif
// if (fabs(sumR-sumRTh)>(0.000001*sumRTh)) {
// cout << " !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!\n"
// << "WARNING!!!! verify LaguerrePolyQuad sum of nodes not accurate: "
// << " current value " << sumR << " theoretical value : " << sumRTh
// << " !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!\n"
// << endl;
// }
// if (fabs(sumW-sumWTh)>(0.000001*sumWTh)) {
// cout << " !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!\n"
// << "WARNING!!!! verify LaguerrePolyQuad sum of weights not accurate: "
// << " current value " << sumW << " theoretical value : " << sumWTh
// << " !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!\n"
// << endl;
// }
alphaFact_ = sqrt((r_8)alphaFact); //sqrt(alpha!)
}//Ctor
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment