724: Establish common solver math layer (#727)

* 676: Decouple dense math utilities from FE

* 724: Establish common solver math layer

* 724: Align math namespace with component

* 724: Simplify core exception references

* 724: Update math documentation groups

Author: bearsontheloose

Code/Source/solver/CMakeLists.txt

@@ -177,15 +177,15 @@ if(ENABLE_COVERAGE)
 endif()
 
 add_subdirectory(Core)
+add_subdirectory(common/math)
 add_subdirectory(FE/Common)
-add_subdirectory(FE/Math)
 add_subdirectory(FE/Quadrature)
 add_subdirectory(FE/Basis)
 
 set(SVMP_ARCHITECTURE_TARGETS
   svmp_core
   svmp_fe_common
-  svmp_fe_math
+  svmp_math
   svmp_fe_basis
   svmp_fe_quadrature
 )

Code/Source/solver/FE/Basis/BasisFunction.h

@@ -6,8 +6,8 @@
 
 #include "BasisExceptions.h"
 #include "BasisTraits.h"
-#include "Math/Matrix.h"
-#include "Math/Vector.h"
+#include "math/Matrix.h"
+#include "math/Vector.h"
 #include "Types.h"
 
 #include <cstddef>

Code/Source/solver/FE/Basis/CMakeLists.txt

@@ -19,6 +19,6 @@ target_include_directories(svmp_fe_basis PUBLIC
 
 target_link_libraries(svmp_fe_basis PUBLIC
   svmp_fe_common
-  svmp_fe_math
+  svmp_math
   svmp_fe_quadrature
 )

Code/Source/solver/FE/Basis/NodeOrderingConventions.h

@@ -5,7 +5,7 @@
 #define SVMP_FE_BASIS_NODEORDERINGCONVENTIONS_H
 
 #include "BasisTraits.h"
-#include "Math/Vector.h"
+#include "math/Vector.h"
 #include "Types.h"
 
 #include <array>

Code/Source/solver/FE/Basis/SerendipityBasis.cpp

@@ -3,7 +3,8 @@
 
 #include "SerendipityBasis.h"
 #include "NodeOrderingConventions.h"
-#include "Math/DenseLinearAlgebra.h"
+#include "Core/Exception.h"
+#include "math/DenseLinearAlgebra.h"
 
 #include <algorithm>
 #include <array>
@@ -195,21 +196,21 @@ std::vector<double> build_inverse_vandermonde(
     // too.
     const double norm_v = matrix_norm_inf(vandermonde, n);
 
-    // invert_dense_matrix raises a generic FEException if the Vandermonde is
+    // invert_dense_matrix raises a generic CoreException if the Vandermonde is
     // exactly singular (a rank-deficient pivot). For a serendipity family that
     // means the node set is not unisolvent at this order -- a construction failure
     // in basis terms -- so translate it to BasisConstructionException, presenting
     // the singular and the ill-conditioned cases (below) as one catchable type in
     // one vocabulary. The matrix was just built n-by-n from n modes, so a
-    // size-mismatch FEException cannot originate here; rank deficiency is the only
-    // FEException this call can raise. (Defensive: the supported node sets are
+    // size-mismatch CoreException cannot originate here; rank deficiency is the only
+    // CoreException this call can raise. (Defensive: the supported node sets are
     // provably unisolvent, so this branch is not reachable for the shipped orders.)
     std::vector<double> inverse;
     try {
         inverse = math::invert_dense_matrix(
             std::move(vandermonde), n,
             "SerendipityBasis interpolation matrix for " + label);
-    } catch (const FEException&) {
+    } catch (const CoreException&) {
         svmp::raise<BasisConstructionException>("SerendipityBasis: " + label +
                 " interpolation matrix is singular; the serendipity node set is not "
                 "unisolvent at the requested order");

Code/Source/solver/FE/FE.h

@@ -11,7 +11,7 @@
  * This header intentionally contains no declarations. It gives Doxygen a
  * header-based home for the top-level FE group; submodule groups attach to it
  * from their own headers, including FE_Basis (Basis/BasisFunction.h),
- * FE_Common (Common/Types.h), FE_Math (Math/Vector.h), and FE_Quadrature
+ * FE_Common (Common/Types.h), shared Math (math/Vector.h), and FE_Quadrature
  * (Quadrature/QuadratureRule.h).
  */
 

Code/Source/solver/FE/Math/CMakeLists.txt

@@ -18,5 +18,5 @@ target_include_directories(svmp_fe_quadrature PUBLIC
 
 target_link_libraries(svmp_fe_quadrature PUBLIC
   svmp_fe_common
-  svmp_fe_math
+  svmp_math
 )

Code/Source/solver/FE/Quadrature/CMakeLists.txt

@@ -40,7 +40,7 @@
  */
 
 #include "FE/Common/Types.h"
-#include "FE/Math/Vector.h"
+#include "math/Vector.h"
 
 #include <cstddef>
 #include <span>

Code/Source/solver/FE/Quadrature/QuadratureRule.h

@@ -0,0 +1,22 @@
+file(GLOB_RECURSE SVMP_MATH_FILES CONFIGURE_DEPENDS
+  "${CMAKE_CURRENT_SOURCE_DIR}/include/*.h"
+  "${CMAKE_CURRENT_SOURCE_DIR}/include/*.hpp"
+  "${CMAKE_CURRENT_SOURCE_DIR}/include/*.inl"
+  "${CMAKE_CURRENT_SOURCE_DIR}/src/*.cpp"
+)
+
+add_library(svmp_math OBJECT
+  ${SVMP_MATH_FILES}
+)
+
+target_compile_features(svmp_math PUBLIC cxx_std_20)
+
+target_include_directories(svmp_math PUBLIC
+  ${CMAKE_CURRENT_SOURCE_DIR}/include
+  ${SV_SOURCE_DIR}/ThirdParty/eigen/include
+  ${SV_SOURCE_DIR}/ThirdParty/eigen/include/eigen3
+)
+
+target_link_libraries(svmp_math PUBLIC
+  svmp_core
+)