Expose the field basis label list for pybind11; this adds previously missing functionality

Author: Martin D. Weinberg

expui/BasisFactory.H

@@ -259,6 +259,10 @@ namespace BasisClasses
       if (Naccel > 0) pseudo = currentAccel(time);
     }
 
+    //! Get the field label vector
+    std::vector<std::string> getFieldLabels(void)
+    { return getFieldLabels(coordinates); }
+
   };
 
   using BasisPtr = std::shared_ptr<Basis>;

expui/FieldBasis.H

@@ -98,9 +98,13 @@ namespace BasisClasses
     virtual std::vector<double>
     crt_eval(double x, double y, double z);
 
+    //! Get the field labels
+    std::vector<std::string> getFieldLabels(const Coord ctype)
+    { return fieldLabels; }
+  
   public:
 
-    //! Constructor from YAML node
+   //! Constructor from YAML node
     FieldBasis(const YAML::Node& conf,
 	       const std::string& name="FieldBasis") : Basis(conf, name)
     { configure(); }
@@ -158,10 +162,6 @@ namespace BasisClasses
     {
     }
 
-    //! Get the field labels
-    std::vector<std::string> getFieldLabels(const Coord ctype)
-    { return fieldLabels; }
-  
     //! Return current maximum harmonic order in expansion
     int getLmax() { return lmax; }
 

pyEXP/BasisWrappers.cc

@@ -928,8 +928,25 @@ void BasisFactoryClasses(py::module &m)
          Returns
          -------
          None
-         )");
+         )")
+    .def("getFieldLabels",
+	 [](BasisClasses::Basis& A)
+	 {
+	   return A.getFieldLabels();
+	 },
+	 R"(
+         Provide the field labels for the basis functions
 
+         Parameters
+         ----------
+         None
+
+         Returns
+         -------
+         list: str
+           list of basis function labels
+      )"
+    );
 
   py::class_<BasisClasses::BiorthBasis, std::shared_ptr<BasisClasses::BiorthBasis>, PyBiorthBasis, BasisClasses::Basis>
     (m, "BiorthBasis")
@@ -2044,11 +2061,11 @@ void BasisFactoryClasses(py::module &m)
       // orthoCheck is not in the base class and needs to have
       // different parameters depending on the basis type.  Here the
       // user can and will often need to specify a quadrature value.
-      .def("orthoCheck", [](BasisClasses::FieldBasis& A)
+    .def("orthoCheck", [](BasisClasses::FieldBasis& A)
       {
 	return A.orthoCheck();
       },
-	R"(
+      R"(
         Check orthgonality of basis functions by quadrature
 
         Inner-product matrix of orthogonal functions
@@ -2062,10 +2079,10 @@ void BasisFactoryClasses(py::module &m)
         numpy.ndarray
 	    orthogonality matrix
         )"
-	);
+      );
 
   py::class_<BasisClasses::VelocityBasis, std::shared_ptr<BasisClasses::VelocityBasis>, BasisClasses::FieldBasis>(m, "VelocityBasis")
-      .def(py::init<const std::string&>(),
+    .def(py::init<const std::string&>(),
 	 R"(
          Create a orthogonal velocity-field basis