Add gmgpolar namespace

include/ConfigParser/config_parser.h

@@ -13,6 +13,9 @@
 #include "../../include/GMGPolar/gmgpolar.h"
 #include "test_selection.h"
 
+namespace gmgpolar
+{
+
 class ConfigParser
 {
 public:
@@ -118,3 +121,4 @@ class ConfigParser
     void selectTestCase(GeometryType geometry_type, ProblemType problem_type, AlphaCoeff alpha_type,
                         BetaCoeff beta_type, double Rmax, double kappa_eps, double delta_e, double alpha_jump);
 };
+} // namespace gmgpolar

include/ConfigParser/test_selection.h

@@ -4,8 +4,13 @@
 
 #include "../../include/GMGPolar/test_cases.h"
 
+namespace gmgpolar
+{
+
 using DomainGeometryVariant = std::variant<CircularGeometry, ShafranovGeometry, CzarnyGeometry, CulhamGeometry>;
 
 using DensityProfileCoefficientsVariant =
     std::variant<PoissonCoefficients, SonnendruckerCoefficients, SonnendruckerGyroCoefficients, ZoniCoefficients,
                  ZoniGyroCoefficients, ZoniShiftedCoefficients, ZoniShiftedGyroCoefficients>;
+
+} // namespace gmgpolar

include/Definitions/equals.h

@@ -4,9 +4,14 @@
 #include <cmath>
 #include <algorithm>
 
+namespace gmgpolar
+{
+
 template <typename T>
 inline std::enable_if_t<std::is_floating_point_v<T>, bool> equals(T lhs, T rhs)
 {
     return std::abs(lhs - rhs) <=
            (1e3 * std::numeric_limits<T>::epsilon()) * std::max(1.0, std::max(std::abs(lhs), std::abs(rhs)));
-}
+}
+
+} // namespace gmgpolar

include/DirectSolver/DirectSolver-COO-MUMPS-Give/directSolverGive.h

@@ -4,6 +4,9 @@
 
 #ifdef GMGPOLAR_USE_MUMPS
 
+namespace gmgpolar
+{
+
 template <class LevelCacheType>
 class DirectSolver_COO_MUMPS_Give : public DirectSolver<LevelCacheType>
 {
@@ -84,4 +87,5 @@ class DirectSolver_COO_MUMPS_Give : public DirectSolver<LevelCacheType>
     #include "directSolverGive.inl"
     #include "matrixStencil.inl"
 
+} // namespace gmgpolar
 #endif

include/DirectSolver/DirectSolver-COO-MUMPS-Take/directSolverTake.h

@@ -4,6 +4,9 @@
 
 #ifdef GMGPOLAR_USE_MUMPS
 
+namespace gmgpolar
+{
+
 template <class LevelCacheType>
 class DirectSolver_COO_MUMPS_Take : public DirectSolver<LevelCacheType>
 {
@@ -83,4 +86,5 @@ class DirectSolver_COO_MUMPS_Take : public DirectSolver<LevelCacheType>
     #include "directSolverTake.inl"
     #include "matrixStencil.inl"
 
+} // namespace gmgpolar
 #endif

include/DirectSolver/DirectSolver-CSR-LU-Give/directSolverGiveCustomLU.h

@@ -2,6 +2,9 @@
 
 #include "../directSolver.h"
 
+namespace gmgpolar
+{
+
 template <class LevelCacheType>
 class DirectSolver_CSR_LU_Give : public DirectSolver<LevelCacheType>
 {
@@ -64,3 +67,4 @@ class DirectSolver_CSR_LU_Give : public DirectSolver<LevelCacheType>
 #include "buildSolverMatrix.inl"
 #include "directSolverGive.inl"
 #include "matrixStencil.inl"
+} // namespace gmgpolar

include/DirectSolver/DirectSolver-CSR-LU-Take/directSolverTakeCustomLU.h

@@ -2,6 +2,9 @@
 
 #include "../directSolver.h"
 
+namespace gmgpolar
+{
+
 template <class LevelCacheType>
 class DirectSolver_CSR_LU_Take : public DirectSolver<LevelCacheType>
 {
@@ -63,3 +66,4 @@ class DirectSolver_CSR_LU_Take : public DirectSolver<LevelCacheType>
 #include "buildSolverMatrix.inl"
 #include "directSolverTake.inl"
 #include "matrixStencil.inl"
+} // namespace gmgpolar

include/DirectSolver/directSolver.h

@@ -19,6 +19,9 @@
     #include "mpi.h"
 #endif
 
+namespace gmgpolar
+{
+
 template <class LevelCacheType>
 class DirectSolver
 {
@@ -43,3 +46,4 @@ class DirectSolver
     const bool DirBC_Interior_;
     const int num_omp_threads_;
 };
+} // namespace gmgpolar

include/ExtrapolatedSmoother/ExtrapolatedSmootherGive/extrapolatedSmootherGive.h

@@ -2,6 +2,9 @@
 
 #include "../extrapolatedSmoother.h"
 
+namespace gmgpolar
+{
+
 // The extrapolated smoother implements the coupled circle-radial smoothing with coarse node elimination.
 // Coarse nodes are excluded from the smoothing iteration by:
 // 1. Setting their diagonal entries in A_sc to 1.0 (identity mapping)
@@ -196,3 +199,4 @@ class ExtrapolatedSmootherGive : public ExtrapolatedSmoother<LevelCacheType>
 #include "applyAscOrtho.inl"
 #include "solveAscSystem.inl"
 #include "initializeMumps.inl"
+} // namespace gmgpolar

include/ExtrapolatedSmoother/ExtrapolatedSmootherTake/extrapolatedSmootherTake.h

@@ -2,6 +2,9 @@
 
 #include "../extrapolatedSmoother.h"
 
+namespace gmgpolar
+{
+
 // The extrapolated smoother implements the coupled circle-radial smoothing with coarse node elimination.
 // Coarse nodes are excluded from the smoothing iteration by:
 // 1. Setting their diagonal entries in A_sc to 1.0 (identity mapping)
@@ -200,3 +203,4 @@ class ExtrapolatedSmootherTake : public ExtrapolatedSmoother<LevelCacheType>
 #include "buildTridiagonalAsc.inl"
 #include "applyAscOrtho.inl"
 #include "solveAscSystem.inl"
+} // namespace gmgpolar

include/ExtrapolatedSmoother/extrapolatedSmoother.h

@@ -20,6 +20,9 @@
     #include "mpi.h"
 #endif
 
+namespace gmgpolar
+{
+
 template <class LevelCacheType>
 class ExtrapolatedSmoother
 {
@@ -42,3 +45,4 @@ class ExtrapolatedSmoother
     const bool DirBC_Interior_;
     const int num_omp_threads_;
 };
+} // namespace gmgpolar

include/GMGPolar/gmgpolar.h

@@ -12,6 +12,9 @@
 
 #include "igmgpolar.h"
 
+namespace gmgpolar
+{
+
 template <concepts::DomainGeometry DomainGeometry, concepts::DensityProfileCoefficients DensityProfileCoefficients>
 class GMGPolar : public IGMGPolar
 {
@@ -210,3 +213,4 @@ class GMGPolar : public IGMGPolar
 #include "MultigridMethods/extrapolated_multigrid_V_Cycle.h"
 #include "MultigridMethods/multigrid_F_Cycle.h"
 #include "MultigridMethods/multigrid_W_Cycle.h"
+} // namespace gmgpolar

include/GMGPolar/igmgpolar.h

@@ -18,6 +18,9 @@
 #include "../PolarGrid/polargrid.h"
 #include "../Definitions/global_definitions.h"
 
+namespace gmgpolar
+{
+
 class IGMGPolar
 {
 public:
@@ -222,3 +225,4 @@ class IGMGPolar
     double t_avg_MGC_residual_;
     double t_avg_MGC_directSolver_;
 };
+} // namespace gmgpolar

include/InputFunctions/BoundaryConditions/cartesianR2_Boundary_CircularGeometry.h

@@ -4,6 +4,9 @@
 
 #include "../boundaryConditions.h"
 
+namespace gmgpolar
+{
+
 class CartesianR2_Boundary_CircularGeometry
 {
 public:
@@ -17,3 +20,4 @@ class CartesianR2_Boundary_CircularGeometry
 };
 
 static_assert(concepts::BoundaryConditions<CartesianR2_Boundary_CircularGeometry>);
+} // namespace gmgpolar

include/InputFunctions/BoundaryConditions/cartesianR2_Boundary_CzarnyGeometry.h

@@ -4,6 +4,9 @@
 
 #include "../boundaryConditions.h"
 
+namespace gmgpolar
+{
+
 class CartesianR2_Boundary_CzarnyGeometry
 {
 public:
@@ -24,3 +27,4 @@ class CartesianR2_Boundary_CzarnyGeometry
 };
 
 static_assert(concepts::BoundaryConditions<CartesianR2_Boundary_CzarnyGeometry>);
+} // namespace gmgpolar

include/InputFunctions/BoundaryConditions/cartesianR2_Boundary_ShafranovGeometry.h

@@ -4,6 +4,9 @@
 
 #include "../boundaryConditions.h"
 
+namespace gmgpolar
+{
+
 class CartesianR2_Boundary_ShafranovGeometry
 {
 public:
@@ -19,3 +22,4 @@ class CartesianR2_Boundary_ShafranovGeometry
 };
 
 static_assert(concepts::BoundaryConditions<CartesianR2_Boundary_ShafranovGeometry>);
+} // namespace gmgpolar

include/InputFunctions/BoundaryConditions/cartesianR6_Boundary_CircularGeometry.h

@@ -4,6 +4,9 @@
 
 #include "../boundaryConditions.h"
 
+namespace gmgpolar
+{
+
 class CartesianR6_Boundary_CircularGeometry
 {
 public:
@@ -17,3 +20,4 @@ class CartesianR6_Boundary_CircularGeometry
 };
 
 static_assert(concepts::BoundaryConditions<CartesianR6_Boundary_CircularGeometry>);
+} // namespace gmgpolar

include/InputFunctions/BoundaryConditions/cartesianR6_Boundary_CzarnyGeometry.h

@@ -4,6 +4,9 @@
 
 #include "../boundaryConditions.h"
 
+namespace gmgpolar
+{
+
 class CartesianR6_Boundary_CzarnyGeometry
 {
 public:
@@ -24,3 +27,4 @@ class CartesianR6_Boundary_CzarnyGeometry
 };
 
 static_assert(concepts::BoundaryConditions<CartesianR6_Boundary_CzarnyGeometry>);
+} // namespace gmgpolar

include/InputFunctions/BoundaryConditions/cartesianR6_Boundary_ShafranovGeometry.h

@@ -4,6 +4,9 @@
 
 #include "../boundaryConditions.h"
 
+namespace gmgpolar
+{
+
 class CartesianR6_Boundary_ShafranovGeometry
 {
 public:
@@ -19,3 +22,4 @@ class CartesianR6_Boundary_ShafranovGeometry
 };
 
 static_assert(concepts::BoundaryConditions<CartesianR6_Boundary_ShafranovGeometry>);
+} // namespace gmgpolar

include/InputFunctions/BoundaryConditions/polarR6_Boundary_CircularGeometry.h

@@ -4,6 +4,9 @@
 
 #include "../boundaryConditions.h"
 
+namespace gmgpolar
+{
+
 class PolarR6_Boundary_CircularGeometry
 {
 public:
@@ -17,3 +20,4 @@ class PolarR6_Boundary_CircularGeometry
 };
 
 static_assert(concepts::BoundaryConditions<PolarR6_Boundary_CircularGeometry>);
+} // namespace gmgpolar

include/InputFunctions/BoundaryConditions/polarR6_Boundary_CulhamGeometry.h

@@ -4,6 +4,9 @@
 
 #include "../boundaryConditions.h"
 
+namespace gmgpolar
+{
+
 class PolarR6_Boundary_CulhamGeometry
 {
 public:
@@ -17,3 +20,4 @@ class PolarR6_Boundary_CulhamGeometry
 };
 
 static_assert(concepts::BoundaryConditions<PolarR6_Boundary_CulhamGeometry>);
+} // namespace gmgpolar

include/InputFunctions/BoundaryConditions/polarR6_Boundary_CzarnyGeometry.h

@@ -4,6 +4,9 @@
 
 #include "../boundaryConditions.h"
 
+namespace gmgpolar
+{
+
 class PolarR6_Boundary_CzarnyGeometry
 {
 public:
@@ -23,3 +26,4 @@ class PolarR6_Boundary_CzarnyGeometry
 };
 
 static_assert(concepts::BoundaryConditions<PolarR6_Boundary_CzarnyGeometry>);
+} // namespace gmgpolar

include/InputFunctions/BoundaryConditions/polarR6_Boundary_ShafranovGeometry.h

@@ -4,6 +4,9 @@
 
 #include "../boundaryConditions.h"
 
+namespace gmgpolar
+{
+
 class PolarR6_Boundary_ShafranovGeometry
 {
 public:
@@ -19,3 +22,4 @@ class PolarR6_Boundary_ShafranovGeometry
 };
 
 static_assert(concepts::BoundaryConditions<PolarR6_Boundary_ShafranovGeometry>);
+} // namespace gmgpolar

include/InputFunctions/BoundaryConditions/refined_Boundary_CircularGeometry.h

@@ -4,6 +4,9 @@
 
 #include "../boundaryConditions.h"
 
+namespace gmgpolar
+{
+
 class Refined_Boundary_CircularGeometry
 {
 public:
@@ -17,3 +20,4 @@ class Refined_Boundary_CircularGeometry
 };
 
 static_assert(concepts::BoundaryConditions<Refined_Boundary_CircularGeometry>);
+} // namespace gmgpolar

include/InputFunctions/BoundaryConditions/refined_Boundary_CulhamGeometry.h

@@ -4,6 +4,9 @@
 
 #include "../boundaryConditions.h"
 
+namespace gmgpolar
+{
+
 class Refined_Boundary_CulhamGeometry
 {
 public:
@@ -17,3 +20,4 @@ class Refined_Boundary_CulhamGeometry
 };
 
 static_assert(concepts::BoundaryConditions<Refined_Boundary_CulhamGeometry>);
+} // namespace gmgpolar

include/InputFunctions/BoundaryConditions/refined_Boundary_CzarnyGeometry.h

@@ -4,6 +4,9 @@
 
 #include "../boundaryConditions.h"
 
+namespace gmgpolar
+{
+
 class Refined_Boundary_CzarnyGeometry
 {
 public:
@@ -23,3 +26,4 @@ class Refined_Boundary_CzarnyGeometry
 };
 
 static_assert(concepts::BoundaryConditions<Refined_Boundary_CzarnyGeometry>);
+} // namespace gmgpolar