(JEC) 7/10/15 update lagsht_testsuite.cc

src/lagsht_testsuite.cc

@@ -26,6 +26,7 @@ struct PARAM {
   int Lmax;
   int N;
   r_8 R;
+  int Pmax; //JEC 7/10/15
   int alpha;
   string geometry;
   int ntheta;
@@ -288,6 +289,7 @@ void TestLaguerre2Bessel() {
 
   int Nmax = param.N;
   int Lmax = param.Lmax;
+  int Pmax = param.Pmax;
   r_8 Rmax = param.R;
   string geometry = param.geometry;
   int ntheta = param.ntheta;
@@ -349,7 +351,6 @@ void TestLaguerre2Bessel() {
 
   cout << "tstack 2 Start" <<endl;
   tstack_push("processing");
-  int Pmax = 64*Nmax; //To see 
   cout << "tstack 2a Start" <<endl;
   tstack_push("Init Laguerre 2 Bessel");
   Laguerre2Bessel lag2bess(sphere,lagTrans,Nmax,Pmax,Rmax);
@@ -499,6 +500,7 @@ int main(int narg, char *arg[]) {
   int N = 128; 
   r_8 R = 1.;
   int Lmax = 128;
+  int Pmax = 0; //JEC 7/10/15
 
   int test=1;
 
@@ -510,10 +512,12 @@ int main(int narg, char *arg[]) {
   while (ka<narg) {
     if (strcmp(arg[ka],"-h")==0) {
       cout << "usage: ./Objs/lagsht_testsuite -t <test number> [1]\n"
-	   << " [-n <Nmax> [5]] [-l <Lmax> [10]]\n" 
+	   << " [-n <Nmax>] [-l <Lmax>]\n" 
 	   << " [-g <geometry> Gauss|Fejer1|Healpix [Gauss]] \n"
 	   << " [-ntheta <number of theta rings> [determined by the geometry]\n     in case of Healpix gives the nside parameter]\n"
 	   << " [-nphi <number of pixel per rings> [determined by the geometry]]\n"
+	   << " [-p <Pmax for Fourier-Bessel> [Nmax] ]\n"
+	   << " [-r <Rmax> [1.]]"
 	   << endl; 
       return 0;
     }
@@ -545,17 +549,36 @@ int main(int narg, char *arg[]) {
       R = atof(arg[ka+1]);
       ka+=2;      
     }    
+    else if (strcmp(arg[ka],"-p")==0) {
+      Pmax = atof(arg[ka+1]);
+      ka+=2;      
+    }    
     else ka++;
   }//eo while
 
   param.Lmax = Lmax;
   param.N    = N;
+  if(Pmax<N){Pmax = 32*N; cout << "(Warning): Pmax<Nmax => modify the value to "<< Pmax << endl;}
+  param.Pmax = Pmax;
   param.R    = R;
   param.alpha = 2;
   param.geometry = geometry;
   param.ntheta = ntheta;
   param.nphi = nphi;
 
+  cout << "Configuration parameters are set to: "
+       << " Test number : " << test << "\n"
+       << " Lmax, Nmax, Pmax, Rmax , alpha: " 
+       << param.Lmax << ", "
+       << param.N << ", "
+       << param.Pmax << ", "
+       << param.R << ", "
+       << param.alpha << "\n"
+       << "Geometry: " << param.geometry << " Ntheta, Nphi: " 
+       << param.ntheta << ", " << param.nphi
+       << endl;
+
+
   int rc=0;
   try {