Add biorthogonal check to all reset locations

EXP-code · Oct 1, 2023 · 4979ca4 · 4979ca4
1 parent 13bebf1
commit 4979ca4
Show file tree

Hide file tree

Showing 2 changed files with 18 additions and 4 deletions.
diff --git a/src/Cylinder.cc b/src/Cylinder.cc
@@ -214,10 +214,6 @@ Cylinder::Cylinder(Component* c0, const YAML::Node& conf, MixtureBasis *m) :
   ortho = std::make_shared<CylEXP>
     (nmaxfid, lmaxfid, mmax, nmax, acyl, hcyl, ncylodd, cachename);
 
-  // Test for basis consistency
-  //
-  orthoTest(ortho->orthoCheck(), "Cylinder", "m");
-
   // Set azimuthal harmonic order restriction?
   //
   if (mlim>=0)  ortho->set_mlim(mlim);
@@ -379,6 +375,13 @@ Cylinder::Cylinder(Component* c0, const YAML::Node& conf, MixtureBasis *m) :
   }
 #endif
 
+  // Test for basis consistency
+  //
+  std::cout << "---- ";
+  orthoTest(ortho->orthoCheck(), "Cylinder", "m");
+
+  // Initialize internal variables
+  //
   ncompcyl = 0;
 
   pos.resize(nthrds);

diff --git a/src/Sphere.cc b/src/Sphere.cc
@@ -86,6 +86,11 @@ Sphere::Sphere(Component* c0, const YAML::Node& conf, MixtureBasis* m) :
   }
 
 
+  // Test for basis consistency
+  //
+  std::cout << "---- ";
+  orthoTest(ortho->orthoCheck(std::max<int>(nmax*5, 100)), "Sphere", "l");
+
   setup();
 }
 
@@ -322,6 +327,7 @@ void Sphere::make_model_bin()
 
   // Test for basis consistency
   //
+  std::cout << "---- ";
   orthoTest(ortho->orthoCheck(std::max<int>(nmax*5, 100)), "Sphere", "l");
 
   // Update time trigger
@@ -464,6 +470,11 @@ void Sphere::make_model_plummer()
   std::string cachename = outdir  + cache_file + "." + runtag;
   ortho = std::make_shared<SLGridSph>(mod, Lmax, nmax, numr, Rmin, Rmax, false, 1, 1.0, cachename);
 
+  // Test for basis consistency
+  //
+  std::cout << "---- ";
+  orthoTest(ortho->orthoCheck(std::max<int>(nmax*5, 100)), "Sphere", "l");
+
   // Update time trigger
   //
   tnext = tnow + dtime;