(JEC) 19/10/15 polish

src/laguerreTransform.cc

@@ -32,7 +32,6 @@ namespace LagSHT {
   LaguerreFuncQuad lag(N_,alpha_);
   lag.QuadWeightNodes(nodes_,weights_);
 
-  //JEC 30/9/15 using alphaFact & alphaFact_ with different types is missleading
   int tmp = 1;
   for(int i=1;i<=alpha_; i++) tmp *= i;
 
@@ -55,8 +54,7 @@ namespace LagSHT {
 }//Ctor
 
 int LaguerreTransform::R2Index(r_8 r) const {
-  using namespace std;
-  //JEC 22/9/15  r_8 rscaled = r*nodes_[N_-1]/R_;
+  //  using namespace std;
   r_8 rscaled = r/tau_;
   vector<r_8> dist(N_);
   transform(nodes_.begin(),nodes_.end(),dist.begin(), distance(rscaled));
@@ -118,7 +116,7 @@ void LaguerreTransform::MultiAnalysis(const vector< complex<r_8> >& fi,
 
 }//MultiAnalysis
 
-  //JEC 9/6/15 real case
+  //JEC 9/6/15 real case: not optimized
 void LaguerreTransform::MultiAnalysisR(const vector<r_8>& fi, 
 				      vector<r_8>& fn, int stride) {