Address PR feedback: EnergyNormSolver -> WeightedHCurlNormSolver

Author: sebastiangrimberg

palace/linalg/hcurl.cpp

@@ -15,13 +15,12 @@
 namespace palace
 {
 
-EnergyNormSolver::EnergyNormSolver(const MaterialOperator &mat_op,
-                                   mfem::ParFiniteElementSpaceHierarchy &nd_fespaces,
-                                   mfem::ParFiniteElementSpaceHierarchy &h1_fespaces,
-                                   const std::vector<mfem::Array<int>> &nd_dbc_tdof_lists,
-                                   const std::vector<mfem::Array<int>> &h1_dbc_tdof_lists,
-                                   double tol, int max_it, int print,
-                                   int pa_order_threshold)
+WeightedHCurlNormSolver::WeightedHCurlNormSolver(
+    const MaterialOperator &mat_op, mfem::ParFiniteElementSpaceHierarchy &nd_fespaces,
+    mfem::ParFiniteElementSpaceHierarchy &h1_fespaces,
+    const std::vector<mfem::Array<int>> &nd_dbc_tdof_lists,
+    const std::vector<mfem::Array<int>> &h1_dbc_tdof_lists, double tol, int max_it,
+    int print, int pa_order_threshold)
 {
   constexpr int skip_zeros = 0;
   constexpr auto MatTypeMuInv = MaterialPropertyType::INV_PERMEABILITY;

palace/linalg/hcurl.hpp

@@ -27,7 +27,7 @@ class MaterialOperator;
 //
 // This solver implements a solver for the operator K + M in a Nedelec space.
 //
-class EnergyNormSolver
+class WeightedHCurlNormSolver
 {
 private:
   // H(curl) norm operator A = K + M and its projection Gᵀ A G.
@@ -37,12 +37,12 @@ class EnergyNormSolver
   std::unique_ptr<KspSolver> ksp;
 
 public:
-  EnergyNormSolver(const MaterialOperator &mat_op,
-                   mfem::ParFiniteElementSpaceHierarchy &nd_fespaces,
-                   mfem::ParFiniteElementSpaceHierarchy &h1_fespaces,
-                   const std::vector<mfem::Array<int>> &nd_dbc_tdof_lists,
-                   const std::vector<mfem::Array<int>> &h1_dbc_tdof_lists, double tol,
-                   int max_it, int print, int pa_order_threshold);
+  WeightedHCurlNormSolver(const MaterialOperator &mat_op,
+                          mfem::ParFiniteElementSpaceHierarchy &nd_fespaces,
+                          mfem::ParFiniteElementSpaceHierarchy &h1_fespaces,
+                          const std::vector<mfem::Array<int>> &nd_dbc_tdof_lists,
+                          const std::vector<mfem::Array<int>> &h1_dbc_tdof_lists,
+                          double tol, int max_it, int print, int pa_order_threshold);
 
   const Operator &GetOperator() { return *A; }
 

palace/models/romoperator.cpp

@@ -114,7 +114,7 @@ RomOperator::RomOperator(const IoData &iodata, SpaceOperator &spaceop) : spaceop
   if (iodata.solver.driven.adaptive_metric_aposteriori)
   {
     constexpr int curlcurl_verbose = 0;
-    kspKM = std::make_unique<EnergyNormSolver>(
+    kspKM = std::make_unique<WeightedHCurlNormSolver>(
         spaceop.GetMaterialOp(), spaceop.GetNDSpaces(), spaceop.GetH1Spaces(),
         spaceop.GetNDDbcTDofLists(), spaceop.GetH1DbcTDofLists(), iodata.solver.linear.tol,
         iodata.solver.linear.max_it, curlcurl_verbose, iodata.solver.pa_order_threshold);

palace/models/romoperator.hpp

@@ -42,7 +42,7 @@ class RomOperator
   std::unique_ptr<ComplexKspSolver> ksp;
 
   // Linear solver for inner product solves for error metric.
-  std::unique_ptr<EnergyNormSolver> kspKM;
+  std::unique_ptr<WeightedHCurlNormSolver> kspKM;
 
   // PROM matrices and vectors.
   Eigen::MatrixXcd Kr, Mr, Cr, Ar;