Merge branch 'SLboundaries' of github.com:EXP-code/EXP into SLboundaries

Author: Martin D. Weinberg

CMakeLists.txt

@@ -2,7 +2,7 @@ cmake_minimum_required(VERSION 3.25) # Needed for CUDA, MPI, and CTest features
 
 project(
   EXP
-  VERSION "7.9.0"
+  VERSION "7.9.1"
   HOMEPAGE_URL https://github.com/EXP-code/EXP
   LANGUAGES C CXX Fortran)
 

doc/exp.cfg

@@ -48,7 +48,7 @@ PROJECT_NAME           = EXP
 # could be handy for archiving the generated documentation or if some version
 # control system is used.
 
-PROJECT_NUMBER         = 7.9.0
+PROJECT_NUMBER         = 7.9.1
 
 # Using the PROJECT_BRIEF tag one can provide an optional one line description
 # for a project that appears at the top of each page and should give viewer a

expui/BiorthBasis.H

@@ -163,6 +163,10 @@ namespace BasisClasses
 
     //! Clear the particle selector callback
     void clrSelector() { ftor = nullptr; }
+
+    //! Evaluate acceleration in Cartesian coordinates in centered
+    //! coordinate system
+    virtual std::vector<double> getAccel(double x, double y, double z) = 0;
   };
 
   /**
@@ -307,6 +311,9 @@ namespace BasisClasses
       return ret;
     }
 
+    //! Evaluate acceleration in Cartesian coordinates in centered
+    //! coordinate system
+    virtual std::vector<double> getAccel(double x, double y, double z);
   };
 
   /**
@@ -542,6 +549,10 @@ namespace BasisClasses
       return ret;
     }
 
+    //! Evaluate acceleration in Cartesian coordinates in centered
+    //! coordinate system
+    std::vector<double> getAccel(double x, double y, double z);
+
   };
 
   /**
@@ -693,6 +704,10 @@ namespace BasisClasses
       return ret;
     }
 
+    //! Evaluate acceleration in Cartesian coordinates in centered
+    //! coordinate system
+    std::vector<double> getAccel(double x, double y, double z);
+
   };
 
   /**
@@ -841,10 +856,18 @@ namespace BasisClasses
     bool orthoTest()
     {
       auto [ret, worst, lworst] = orthoCompute(sl->orthoCheck());
-      // For the CTest log
-      std::cout << "---- Cylindrical::orthoTest: worst=" << worst << std::endl;
+      if (myid==0) {
+	// For the CTest log
+	std::cout << "---- Cylindrical::orthoTest: worst=" << worst
+		  << std::endl;
+      }
       return ret;
     }
+
+    //! Evaluate acceleration in Cartesian coordinates in centered
+    //! coordinate system
+    std::vector<double> getAccel(double x, double y, double z);
+
   };
 
   /**
@@ -991,6 +1014,10 @@ namespace BasisClasses
       return true;
     }
 
+    //! Evaluate acceleration in Cartesian coordinates in centered
+    //! coordinate system
+    std::vector<double> getAccel(double x, double y, double z);
+
   };
 
   /**
@@ -1113,6 +1140,10 @@ namespace BasisClasses
       return true;
     }
 
+    //! Evaluate acceleration in Cartesian coordinates in centered
+    //! coordinate system
+    std::vector<double> getAccel(double x, double y, double z);
+
   };