(JEC) 24/11/15 clean the default debuging

src/lagsht_testsuite.cc

@@ -2,6 +2,7 @@
 #include <iostream>
 #include <iomanip>      // std::setprecision
 #include <fstream>
+#include <sstream>
 #include <string.h>
 #include <limits>       // std::numeric_limits
 #include <string>
@@ -18,7 +19,7 @@ using namespace std;
 
 using namespace LagSHT;
 
-#define DEBUG 2
+#define DEBUG 0
 
 //-------- Parameters set in the main and  used in the different test functions
 struct PARAM {
@@ -236,11 +237,11 @@ void MultiSphericalLaguerreTransform() {
 #endif
   vector<r_8> fijk(NpTot);
   sphlagtrans.Synthesis(flmnOrig,fijk);
-  //#if DEBUG >= 2
+#if DEBUG >= 2
   for (int i=0; i<NpTot; i++){
     cout << "fijk("<<i<<"): " << fijk[i] << endl;
   }
-  //#endif
+#endif
   
   tstack_pop("processing part Synthesis");
   tstack_push("processing part Analysis");
@@ -572,8 +573,8 @@ void TestPixelization() {
 
   int state = 1234567 + 8912 ; //random seed
   
-  
-  std::ofstream ofs ("pixels.txt", std::ofstream::out);
+  stringstream ss; ss << "pixels-" << geometry << "-L" << Lmax << ".txt";
+  std::ofstream ofs (ss.str(), std::ofstream::out);
 
   r_8 dTheta = M_PI/ntheta/2;
   r_8 dPhi   = 2*M_PI/nphi/2;
@@ -711,13 +712,6 @@ int main(int narg, char *arg[]) {
        << param.alpha << "\n"
        << "Geometry: " << param.geometry << " Ntheta, Nphi: " 
        << param.ntheta << ", " << param.nphi<<"\n";
-  if(recomputeJlnp){
-    cout << "rebuild Jlnp: \n Clenshaw Path<"<<param.clenshawDir<<">\n and jlnp Path<"
-	 << param.jlnpDir << ">"
-	 << endl;
-  } else {
-    cout << "Jlnp path <" << param.jlnpDir << ">" << endl;
-  }
 
   int rc=0;
   try {
@@ -757,6 +751,13 @@ int main(int narg, char *arg[]) {
 
     case 5:
       {
+	if(recomputeJlnp){
+	  cout << "rebuild Jlnp: \n Clenshaw Path<"<<param.clenshawDir<<">\n and jlnp Path<"
+	       << param.jlnpDir << ">"
+	       << endl;
+	} else {
+	  cout << "Jlnp path <" << param.jlnpDir << ">" << endl;
+	}
 	if(Pmax<N){ param.Pmax = 16*N; cout << "(Warning): Pmax<Nmax => modify the value to "<<  param.Pmax << endl;}
 	tstack_push("TestLaguerre2Bessel");
 	TestLaguerre2Bessel();

src/laguerre2bessel.cc

@@ -10,7 +10,7 @@
 #include <sstream>
 
 
-#define DEBUG 1
+#define DEBUG 0
 
 namespace LagSHT {