Make locally inlined functions static (#172)

julianlitz · web-flow · commit 7748cf59e693 · 2026-02-10T15:39:07.000+01:00
diff --git a/src/DirectSolver/DirectSolver-COO-MUMPS-Give/buildSolverMatrix.cpp b/src/DirectSolver/DirectSolver-COO-MUMPS-Give/buildSolverMatrix.cpp
@@ -4,7 +4,8 @@
 
 #ifdef GMGPOLAR_USE_MUMPS
 
-inline void updateMatrixElement(SparseMatrixCOO<double>& matrix, int ptr, int offset, int row, int col, double val)
+static inline void updateMatrixElement(SparseMatrixCOO<double>& matrix, int ptr, int offset, int row, int col,
+                                       double val)
 {
     matrix.row_index(ptr + offset) = row;
     matrix.col_index(ptr + offset) = col;
diff --git a/src/DirectSolver/DirectSolver-COO-MUMPS-Take/buildSolverMatrix.cpp b/src/DirectSolver/DirectSolver-COO-MUMPS-Take/buildSolverMatrix.cpp
@@ -2,7 +2,8 @@
 
 #ifdef GMGPOLAR_USE_MUMPS
 
-inline void updateMatrixElement(SparseMatrixCOO<double>& matrix, int ptr, int offset, int row, int col, double val)
+static inline void updateMatrixElement(SparseMatrixCOO<double>& matrix, int ptr, int offset, int row, int col,
+                                       double val)
 {
     matrix.row_index(ptr + offset) = row;
     matrix.col_index(ptr + offset) = col;
diff --git a/src/DirectSolver/DirectSolver-CSR-LU-Give/buildSolverMatrix.cpp b/src/DirectSolver/DirectSolver-CSR-LU-Give/buildSolverMatrix.cpp
@@ -2,7 +2,7 @@
 
 #include "../../../include/common/geometry_helper.h"
 
-inline void updateMatrixElement(SparseMatrixCSR<double>& matrix, int offset, int row, int col, double val)
+static inline void updateMatrixElement(SparseMatrixCSR<double>& matrix, int offset, int row, int col, double val)
 {
     matrix.row_nz_index(row, offset) = col;
     matrix.row_nz_entry(row, offset) += val;
diff --git a/src/DirectSolver/DirectSolver-CSR-LU-Take/buildSolverMatrix.cpp b/src/DirectSolver/DirectSolver-CSR-LU-Take/buildSolverMatrix.cpp
@@ -1,6 +1,6 @@
 #include "../../../include/DirectSolver/DirectSolver-CSR-LU-Take/directSolverTakeCustomLU.h"
 
-inline void updateMatrixElement(SparseMatrixCSR<double>& matrix, int offset, int row, int col, double val)
+static inline void updateMatrixElement(SparseMatrixCSR<double>& matrix, int offset, int row, int col, double val)
 {
     matrix.row_nz_index(row, offset) = col;
     matrix.row_nz_entry(row, offset) = val;
diff --git a/src/ExtrapolatedSmoother/ExtrapolatedSmootherGive/buildAscMatrices.cpp b/src/ExtrapolatedSmoother/ExtrapolatedSmootherGive/buildAscMatrices.cpp
@@ -3,7 +3,8 @@
 #include "../../../include/common/geometry_helper.h"
 
 /* Tridiagonal matrices */
-inline void updateTridiagonalElement(SymmetricTridiagonalSolver<double>& matrix, int row, int column, double value)
+static inline void updateTridiagonalElement(SymmetricTridiagonalSolver<double>& matrix, int row, int column,
+                                            double value)
 {
     if (row == column)
         matrix.main_diagonal(row) += value;
@@ -14,23 +15,23 @@ inline void updateTridiagonalElement(SymmetricTridiagonalSolver<double>& matrix,
 }
 
 /* Diagonal matrices */
-inline void updateDiagonalElement(DiagonalSolver<double>& matrix, int row, int column, double value)
+static inline void updateDiagonalElement(DiagonalSolver<double>& matrix, int row, int column, double value)
 {
     matrix.diagonal(row) += value;
 }
 
 /* Inner Boundary COO/CSR matrix */
 #ifdef GMGPOLAR_USE_MUMPS
-inline void updateCOOCSRMatrixElement(SparseMatrixCOO<double>& matrix, int ptr, int offset, int row, int col,
-                                      double val)
+static inline void updateCOOCSRMatrixElement(SparseMatrixCOO<double>& matrix, int ptr, int offset, int row, int col,
+                                             double val)
 {
     matrix.row_index(ptr + offset) = row;
     matrix.col_index(ptr + offset) = col;
     matrix.value(ptr + offset) += val;
 }
 #else
-inline void updateCOOCSRMatrixElement(SparseMatrixCSR<double>& matrix, int ptr, int offset, int row, int col,
-                                      double val)
+static inline void updateCOOCSRMatrixElement(SparseMatrixCSR<double>& matrix, int ptr, int offset, int row, int col,
+                                             double val)
 {
     matrix.row_nz_index(row, offset) = col;
     matrix.row_nz_entry(row, offset) += val;
diff --git a/src/ExtrapolatedSmoother/ExtrapolatedSmootherGive/smootherSolver.cpp b/src/ExtrapolatedSmoother/ExtrapolatedSmootherGive/smootherSolver.cpp
@@ -2,10 +2,10 @@
 
 #include "../../../include/common/geometry_helper.h"
 
-inline void nodeApplyAscOrthoCircleGive(int i_r, int i_theta, const PolarGrid& grid, bool DirBC_Interior,
-                                        SmootherColor smoother_color, ConstVector<double>& x, ConstVector<double>& rhs,
-                                        Vector<double>& result, double arr, double att, double art, double detDF,
-                                        double coeff_beta)
+static inline void nodeApplyAscOrthoCircleGive(int i_r, int i_theta, const PolarGrid& grid, bool DirBC_Interior,
+                                               SmootherColor smoother_color, ConstVector<double>& x,
+                                               ConstVector<double>& rhs, Vector<double>& result, double arr, double att,
+                                               double art, double detDF, double coeff_beta)
 {
     assert(i_r >= 0 && i_r <= grid.numberSmootherCircles());
 
@@ -394,10 +394,10 @@ inline void nodeApplyAscOrthoCircleGive(int i_r, int i_theta, const PolarGrid& g
     }
 }
 
-inline void nodeApplyAscOrthoRadialGive(int i_r, int i_theta, const PolarGrid& grid, bool DirBC_Interior,
-                                        SmootherColor smoother_color, ConstVector<double>& x, ConstVector<double>& rhs,
-                                        Vector<double>& result, double arr, double att, double art, double detDF,
-                                        double coeff_beta)
+static inline void nodeApplyAscOrthoRadialGive(int i_r, int i_theta, const PolarGrid& grid, bool DirBC_Interior,
+                                               SmootherColor smoother_color, ConstVector<double>& x,
+                                               ConstVector<double>& rhs, Vector<double>& result, double arr, double att,
+                                               double art, double detDF, double coeff_beta)
 {
     assert(i_r >= grid.numberSmootherCircles() - 1 && i_r < grid.nr());
 
diff --git a/src/ExtrapolatedSmoother/ExtrapolatedSmootherTake/buildAscMatrices.cpp b/src/ExtrapolatedSmoother/ExtrapolatedSmootherTake/buildAscMatrices.cpp
@@ -1,7 +1,8 @@
 #include "../../../include/ExtrapolatedSmoother/ExtrapolatedSmootherTake/extrapolatedSmootherTake.h"
 
 /* Tridiagonal matrices */
-inline void updateTridiagonalElement(SymmetricTridiagonalSolver<double>& matrix, int row, int column, double value)
+static inline void updateTridiagonalElement(SymmetricTridiagonalSolver<double>& matrix, int row, int column,
+                                            double value)
 {
     if (row == column)
         matrix.main_diagonal(row) = value;
@@ -12,23 +13,23 @@ inline void updateTridiagonalElement(SymmetricTridiagonalSolver<double>& matrix,
 }
 
 /* Diagonal matrices */
-inline void updateDiagonalElement(DiagonalSolver<double>& matrix, int row, int column, double value)
+static inline void updateDiagonalElement(DiagonalSolver<double>& matrix, int row, int column, double value)
 {
     matrix.diagonal(row) = value;
 }
 
 /* Inner Boundary COO/CSR matrix */
 #ifdef GMGPOLAR_USE_MUMPS
-inline void updateCOOCSRMatrixElement(SparseMatrixCOO<double>& matrix, int ptr, int offset, int row, int col,
-                                      double val)
+static inline void updateCOOCSRMatrixElement(SparseMatrixCOO<double>& matrix, int ptr, int offset, int row, int col,
+                                             double val)
 {
     matrix.row_index(ptr + offset) = row;
     matrix.col_index(ptr + offset) = col;
     matrix.value(ptr + offset)     = val;
 }
 #else
-inline void updateCOOCSRMatrixElement(SparseMatrixCSR<double>& matrix, int ptr, int offset, int row, int col,
-                                      double val)
+static inline void updateCOOCSRMatrixElement(SparseMatrixCSR<double>& matrix, int ptr, int offset, int row, int col,
+                                             double val)
 {
     matrix.row_nz_index(row, offset) = col;
     matrix.row_nz_entry(row, offset) = val;
diff --git a/src/ExtrapolatedSmoother/ExtrapolatedSmootherTake/smootherSolver.cpp b/src/ExtrapolatedSmoother/ExtrapolatedSmootherTake/smootherSolver.cpp
@@ -1,10 +1,11 @@
 #include "../../../include/ExtrapolatedSmoother/ExtrapolatedSmootherTake/extrapolatedSmootherTake.h"
 
-inline void nodeApplyAscOrthoCircleTake(int i_r, int i_theta, const PolarGrid& grid, bool DirBC_Interior,
-                                        SmootherColor smoother_color, ConstVector<double>& x, ConstVector<double>& rhs,
-                                        Vector<double>& result, ConstVector<double>& arr, ConstVector<double>& att,
-                                        ConstVector<double>& art, ConstVector<double>& detDF,
-                                        ConstVector<double>& coeff_beta)
+static inline void nodeApplyAscOrthoCircleTake(int i_r, int i_theta, const PolarGrid& grid, bool DirBC_Interior,
+                                               SmootherColor smoother_color, ConstVector<double>& x,
+                                               ConstVector<double>& rhs, Vector<double>& result,
+                                               ConstVector<double>& arr, ConstVector<double>& att,
+                                               ConstVector<double>& art, ConstVector<double>& detDF,
+                                               ConstVector<double>& coeff_beta)
 {
     assert(i_r >= 0 && i_r <= grid.numberSmootherCircles());
 
@@ -182,11 +183,12 @@ inline void nodeApplyAscOrthoCircleTake(int i_r, int i_theta, const PolarGrid& g
     }
 }
 
-inline void nodeApplyAscOrthoRadialTake(int i_r, int i_theta, const PolarGrid& grid, bool DirBC_Interior,
-                                        SmootherColor smoother_color, ConstVector<double>& x, ConstVector<double>& rhs,
-                                        Vector<double>& result, ConstVector<double>& arr,
-                                        const ConstVector<double>& att, ConstVector<double>& art,
-                                        const ConstVector<double>& detDF, ConstVector<double>& coeff_beta)
+static inline void nodeApplyAscOrthoRadialTake(int i_r, int i_theta, const PolarGrid& grid, bool DirBC_Interior,
+                                               SmootherColor smoother_color, ConstVector<double>& x,
+                                               ConstVector<double>& rhs, Vector<double>& result,
+                                               ConstVector<double>& arr, const ConstVector<double>& att,
+                                               ConstVector<double>& art, const ConstVector<double>& detDF,
+                                               ConstVector<double>& coeff_beta)
 {
     assert(i_r >= grid.numberSmootherCircles() - 1 && i_r < grid.nr());
 
diff --git a/src/Interpolation/extrapolated_prolongation.cpp b/src/Interpolation/extrapolated_prolongation.cpp
@@ -46,9 +46,9 @@
  *
  */
 
-inline void fineNodeExtrapolatedProlongation(int i_r, int i_theta, int i_r_coarse, int i_theta_coarse,
-                                             const PolarGrid& coarse_grid, const PolarGrid& fine_grid,
-                                             Vector<double>& fine_result, ConstVector<double>& coarse_values)
+static inline void fineNodeExtrapolatedProlongation(int i_r, int i_theta, int i_r_coarse, int i_theta_coarse,
+                                                    const PolarGrid& coarse_grid, const PolarGrid& fine_grid,
+                                                    Vector<double>& fine_result, ConstVector<double>& coarse_values)
 {
     if (i_r & 1) {
         if (i_theta & 1) { /* (odd, odd) -> node in center of coarse cell */
diff --git a/src/Interpolation/extrapolated_restriction.cpp b/src/Interpolation/extrapolated_restriction.cpp
@@ -21,9 +21,9 @@
  *  - Radial direction: check domain boundaries
  */
 
-inline void coarseNodeExtrapolatedRestriction(int i_r_coarse, int i_theta_coarse, const PolarGrid& fine_grid,
-                                              const PolarGrid& coarse_grid, Vector<double>& coarse_result,
-                                              ConstVector<double>& fine_values)
+static inline void coarseNodeExtrapolatedRestriction(int i_r_coarse, int i_theta_coarse, const PolarGrid& fine_grid,
+                                                     const PolarGrid& coarse_grid, Vector<double>& coarse_result,
+                                                     ConstVector<double>& fine_values)
 {
     int i_r     = i_r_coarse * 2;
     int i_theta = i_theta_coarse * 2;
diff --git a/src/Interpolation/fmg_interpolation.cpp b/src/Interpolation/fmg_interpolation.cpp
@@ -53,9 +53,9 @@
  * 
  */
 
-inline void fineNodeFMGInterpolation(int i_r, int i_theta, int i_r_coarse, int i_theta_coarse,
-                                     const PolarGrid& coarse_grid, const PolarGrid& fine_grid,
-                                     Vector<double>& fine_result, ConstVector<double>& coarse_values)
+static inline void fineNodeFMGInterpolation(int i_r, int i_theta, int i_r_coarse, int i_theta_coarse,
+                                            const PolarGrid& coarse_grid, const PolarGrid& fine_grid,
+                                            Vector<double>& fine_result, ConstVector<double>& coarse_values)
 {
     /* Case 1: On the boundary */
     if (i_r == 0 || i_r == fine_grid.nr() - 1) {
diff --git a/src/Interpolation/injection.cpp b/src/Interpolation/injection.cpp
@@ -2,9 +2,9 @@
 
 /* Remark: This injection is not scaled. */
 
-inline void coarseNodeInjection(int i_r_coarse, int i_theta_coarse, const PolarGrid& fine_grid,
-                                const PolarGrid& coarse_grid, Vector<double>& coarse_result,
-                                ConstVector<double>& fine_values)
+static inline void coarseNodeInjection(int i_r_coarse, int i_theta_coarse, const PolarGrid& fine_grid,
+                                       const PolarGrid& coarse_grid, Vector<double>& coarse_result,
+                                       ConstVector<double>& fine_values)
 {
     int i_r     = i_r_coarse * 2;
     int i_theta = i_theta_coarse * 2;
diff --git a/src/Interpolation/prolongation.cpp b/src/Interpolation/prolongation.cpp
@@ -51,9 +51,9 @@
  *  - k1, k2 in angular direction
  */
 
-inline void fineNodeProlongation(int i_r, int i_theta, int i_r_coarse, int i_theta_coarse, const PolarGrid& coarse_grid,
-                                 const PolarGrid& fine_grid, Vector<double>& fine_result,
-                                 ConstVector<double>& coarse_values)
+static inline void fineNodeProlongation(int i_r, int i_theta, int i_r_coarse, int i_theta_coarse,
+                                        const PolarGrid& coarse_grid, const PolarGrid& fine_grid,
+                                        Vector<double>& fine_result, ConstVector<double>& coarse_values)
 {
     if (i_r & 1) {
         if (i_theta & 1) { /* (odd, odd) -> fine node in center of coarse cell */
diff --git a/src/Interpolation/restriction.cpp b/src/Interpolation/restriction.cpp
@@ -22,9 +22,9 @@
  *  - Radial direction: check domain boundaries (inner/outer radius)
  */
 
-inline void coarseNodeRestriction(int i_r_coarse, int i_theta_coarse, const PolarGrid& fine_grid,
-                                  const PolarGrid& coarse_grid, Vector<double>& coarse_result,
-                                  ConstVector<double>& fine_values)
+static inline void coarseNodeRestriction(int i_r_coarse, int i_theta_coarse, const PolarGrid& fine_grid,
+                                         const PolarGrid& coarse_grid, Vector<double>& coarse_result,
+                                         ConstVector<double>& fine_values)
 {
     int i_r     = i_r_coarse * 2;
     int i_theta = i_theta_coarse * 2;
diff --git a/src/Residual/ResidualGive/applyAGive.cpp b/src/Residual/ResidualGive/applyAGive.cpp
@@ -1,8 +1,8 @@
 #include "../../../include/Residual/ResidualGive/residualGive.h"
 
-inline void node_apply_a_give(int i_r, int i_theta, double r, double theta, const PolarGrid& grid,
-                              const LevelCache& level_cache, bool DirBC_Interior, Vector<double>& result,
-                              ConstVector<double>& x)
+static inline void node_apply_a_give(int i_r, int i_theta, double r, double theta, const PolarGrid& grid,
+                                     const LevelCache& level_cache, bool DirBC_Interior, Vector<double>& result,
+                                     ConstVector<double>& x)
 {
     /* ---------------------------------------- */
     /* Compute or retrieve stencil coefficients */
diff --git a/src/Residual/ResidualTake/applyResidualTake.cpp b/src/Residual/ResidualTake/applyResidualTake.cpp
@@ -1,9 +1,10 @@
 #include "../../../include/Residual/ResidualTake/residualTake.h"
 
-inline void node_apply_residual_take(int i_r, int i_theta, const PolarGrid& grid, bool DirBC_Interior,
-                                     Vector<double>& result, ConstVector<double>& rhs, ConstVector<double>& x,
-                                     ConstVector<double>& arr, ConstVector<double>& att, ConstVector<double>& art,
-                                     ConstVector<double>& detDF, ConstVector<double>& coeff_beta)
+static inline void node_apply_residual_take(int i_r, int i_theta, const PolarGrid& grid, bool DirBC_Interior,
+                                            Vector<double>& result, ConstVector<double>& rhs, ConstVector<double>& x,
+                                            ConstVector<double>& arr, ConstVector<double>& att,
+                                            ConstVector<double>& art, ConstVector<double>& detDF,
+                                            ConstVector<double>& coeff_beta)
 {
     /* -------------------- */
     /* Node in the interior */
diff --git a/src/Smoother/SmootherGive/buildMatrix.cpp b/src/Smoother/SmootherGive/buildMatrix.cpp
@@ -3,7 +3,7 @@
 #include "../../../include/common/geometry_helper.h"
 
 /* Tridiagonal matrices */
-inline void updateMatrixElement(SymmetricTridiagonalSolver<double>& matrix, int row, int column, double value)
+static inline void updateMatrixElement(SymmetricTridiagonalSolver<double>& matrix, int row, int column, double value)
 {
     if (row == column)
         matrix.main_diagonal(row) += value;
@@ -15,16 +15,16 @@ inline void updateMatrixElement(SymmetricTridiagonalSolver<double>& matrix, int
 
 /* Inner Boundary COO/CSR matrix */
 #ifdef GMGPOLAR_USE_MUMPS
-inline void updateCOOCSRMatrixElement(SparseMatrixCOO<double>& matrix, int ptr, int offset, int row, int col,
-                                      double val)
+static inline void updateCOOCSRMatrixElement(SparseMatrixCOO<double>& matrix, int ptr, int offset, int row, int col,
+                                             double val)
 {
     matrix.row_index(ptr + offset) = row;
     matrix.col_index(ptr + offset) = col;
     matrix.value(ptr + offset) += val;
 }
 #else
-inline void updateCOOCSRMatrixElement(SparseMatrixCSR<double>& matrix, int ptr, int offset, int row, int col,
-                                      double val)
+static inline void updateCOOCSRMatrixElement(SparseMatrixCSR<double>& matrix, int ptr, int offset, int row, int col,
+                                             double val)
 {
     matrix.row_nz_index(row, offset) = col;
     matrix.row_nz_entry(row, offset) += val;
diff --git a/src/Smoother/SmootherGive/smootherSolver.cpp b/src/Smoother/SmootherGive/smootherSolver.cpp
@@ -2,10 +2,10 @@
 
 #include "../../../include/common/geometry_helper.h"
 
-inline void nodeApplyAscOrthoCircleGive(int i_r, int i_theta, const PolarGrid& grid, bool DirBC_Interior,
-                                        SmootherColor smoother_color, ConstVector<double>& x, ConstVector<double>& rhs,
-                                        Vector<double>& result, double arr, double att, double art, double detDF,
-                                        double coeff_beta)
+static inline void nodeApplyAscOrthoCircleGive(int i_r, int i_theta, const PolarGrid& grid, bool DirBC_Interior,
+                                               SmootherColor smoother_color, ConstVector<double>& x,
+                                               ConstVector<double>& rhs, Vector<double>& result, double arr, double att,
+                                               double art, double detDF, double coeff_beta)
 {
     assert(i_r >= 0 && i_r <= grid.numberSmootherCircles());
 
@@ -186,10 +186,10 @@ inline void nodeApplyAscOrthoCircleGive(int i_r, int i_theta, const PolarGrid& g
     }
 }
 
-inline void nodeApplyAscOrthoRadialGive(int i_r, int i_theta, const PolarGrid& grid, bool DirBC_Interior,
-                                        SmootherColor smoother_color, ConstVector<double>& x, ConstVector<double>& rhs,
-                                        Vector<double>& result, double arr, double att, double art, double detDF,
-                                        double coeff_beta)
+static inline void nodeApplyAscOrthoRadialGive(int i_r, int i_theta, const PolarGrid& grid, bool DirBC_Interior,
+                                               SmootherColor smoother_color, ConstVector<double>& x,
+                                               ConstVector<double>& rhs, Vector<double>& result, double arr, double att,
+                                               double art, double detDF, double coeff_beta)
 {
     assert(i_r >= grid.numberSmootherCircles() - 1 && i_r < grid.nr());
 
diff --git a/src/Smoother/SmootherTake/buildMatrix.cpp b/src/Smoother/SmootherTake/buildMatrix.cpp
diff --git a/src/Smoother/SmootherTake/smootherSolver.cpp b/src/Smoother/SmootherTake/smootherSolver.cpp

Original file line number	Diff line number	Diff line change
`@@ -4,7 +4,8 @@`
`4`	`4`
`5`	`5`	`#ifdef GMGPOLAR_USE_MUMPS`
`6`	`6`
`7`		`-inline void updateMatrixElement(SparseMatrixCOO<double>& matrix, int ptr, int offset, int row, int col, double val)`
	`7`	`+static inline void updateMatrixElement(SparseMatrixCOO<double>& matrix, int ptr, int offset, int row, int col,`
	`8`	`+ double val)`
`8`	`9`	`{`
`9`	`10`	`matrix.row_index(ptr + offset) = row;`
`10`	`11`	`matrix.col_index(ptr + offset) = col;`
Original file line number	Diff line number	Diff line change
`@@ -2,7 +2,8 @@`
`2`	`2`
`3`	`3`	`#ifdef GMGPOLAR_USE_MUMPS`
`4`	`4`
`5`		`-inline void updateMatrixElement(SparseMatrixCOO<double>& matrix, int ptr, int offset, int row, int col, double val)`
	`5`	`+static inline void updateMatrixElement(SparseMatrixCOO<double>& matrix, int ptr, int offset, int row, int col,`
	`6`	`+ double val)`
`6`	`7`	`{`
`7`	`8`	`matrix.row_index(ptr + offset) = row;`
`8`	`9`	`matrix.col_index(ptr + offset) = col;`
Original file line number	Diff line number	Diff line change
`@@ -2,7 +2,7 @@`
`2`	`2`
`3`	`3`	`#include "../../../include/common/geometry_helper.h"`
`4`	`4`
`5`		`-inline void updateMatrixElement(SparseMatrixCSR<double>& matrix, int offset, int row, int col, double val)`
	`5`	`+static inline void updateMatrixElement(SparseMatrixCSR<double>& matrix, int offset, int row, int col, double val)`
`6`	`6`	`{`
`7`	`7`	`matrix.row_nz_index(row, offset) = col;`
`8`	`8`	`matrix.row_nz_entry(row, offset) += val;`
Original file line number	Diff line number	Diff line change
`@@ -1,6 +1,6 @@`
`1`	`1`	`#include "../../../include/DirectSolver/DirectSolver-CSR-LU-Take/directSolverTakeCustomLU.h"`
`2`	`2`
`3`		`-inline void updateMatrixElement(SparseMatrixCSR<double>& matrix, int offset, int row, int col, double val)`
	`3`	`+static inline void updateMatrixElement(SparseMatrixCSR<double>& matrix, int offset, int row, int col, double val)`
`4`	`4`	`{`
`5`	`5`	`matrix.row_nz_index(row, offset) = col;`
`6`	`6`	`matrix.row_nz_entry(row, offset) = val;`
Original file line number	Diff line number	Diff line change
`@@ -3,7 +3,8 @@`
`3`	`3`	`#include "../../../include/common/geometry_helper.h"`
`4`	`4`
`5`	`5`	`/* Tridiagonal matrices */`
`6`		`-inline void updateTridiagonalElement(SymmetricTridiagonalSolver<double>& matrix, int row, int column, double value)`
	`6`	`+static inline void updateTridiagonalElement(SymmetricTridiagonalSolver<double>& matrix, int row, int column,`
	`7`	`+ double value)`
`7`	`8`	`{`
`8`	`9`	`if (row == column)`
`9`	`10`	`matrix.main_diagonal(row) += value;`
`@@ -14,23 +15,23 @@ inline void updateTridiagonalElement(SymmetricTridiagonalSolver<double>& matrix,`
`14`	`15`	`}`
`15`	`16`
`16`	`17`	`/* Diagonal matrices */`
`17`		`-inline void updateDiagonalElement(DiagonalSolver<double>& matrix, int row, int column, double value)`
	`18`	`+static inline void updateDiagonalElement(DiagonalSolver<double>& matrix, int row, int column, double value)`
`18`	`19`	`{`
`19`	`20`	`matrix.diagonal(row) += value;`
`20`	`21`	`}`
`21`	`22`
`22`	`23`	`/* Inner Boundary COO/CSR matrix */`
`23`	`24`	`#ifdef GMGPOLAR_USE_MUMPS`
`24`		`-inline void updateCOOCSRMatrixElement(SparseMatrixCOO<double>& matrix, int ptr, int offset, int row, int col,`
`25`		`- double val)`
	`25`	`+static inline void updateCOOCSRMatrixElement(SparseMatrixCOO<double>& matrix, int ptr, int offset, int row, int col,`
	`26`	`+ double val)`
`26`	`27`	`{`
`27`	`28`	`matrix.row_index(ptr + offset) = row;`
`28`	`29`	`matrix.col_index(ptr + offset) = col;`
`29`	`30`	`matrix.value(ptr + offset) += val;`
`30`	`31`	`}`
`31`	`32`	`#else`
`32`		`-inline void updateCOOCSRMatrixElement(SparseMatrixCSR<double>& matrix, int ptr, int offset, int row, int col,`
`33`		`- double val)`
	`33`	`+static inline void updateCOOCSRMatrixElement(SparseMatrixCSR<double>& matrix, int ptr, int offset, int row, int col,`
	`34`	`+ double val)`
`34`	`35`	`{`
`35`	`36`	`matrix.row_nz_index(row, offset) = col;`
`36`	`37`	`matrix.row_nz_entry(row, offset) += val;`
Original file line number	Diff line number	Diff line change
`@@ -2,10 +2,10 @@`
`2`	`2`
`3`	`3`	`#include "../../../include/common/geometry_helper.h"`
`4`	`4`
`5`		`-inline void nodeApplyAscOrthoCircleGive(int i_r, int i_theta, const PolarGrid& grid, bool DirBC_Interior,`
`6`		`- SmootherColor smoother_color, ConstVector<double>& x, ConstVector<double>& rhs,`
`7`		`- Vector<double>& result, double arr, double att, double art, double detDF,`
`8`		`- double coeff_beta)`
	`5`	`+static inline void nodeApplyAscOrthoCircleGive(int i_r, int i_theta, const PolarGrid& grid, bool DirBC_Interior,`
	`6`	`+ SmootherColor smoother_color, ConstVector<double>& x,`
	`7`	`+ ConstVector<double>& rhs, Vector<double>& result, double arr, double att,`
	`8`	`+ double art, double detDF, double coeff_beta)`
`9`	`9`	`{`
`10`	`10`	`assert(i_r >= 0 && i_r <= grid.numberSmootherCircles());`
`11`	`11`
`@@ -394,10 +394,10 @@ inline void nodeApplyAscOrthoCircleGive(int i_r, int i_theta, const PolarGrid& g`
`394`	`394`	`}`
`395`	`395`	`}`
`396`	`396`
`397`		`-inline void nodeApplyAscOrthoRadialGive(int i_r, int i_theta, const PolarGrid& grid, bool DirBC_Interior,`
`398`		`- SmootherColor smoother_color, ConstVector<double>& x, ConstVector<double>& rhs,`
`399`		`- Vector<double>& result, double arr, double att, double art, double detDF,`
`400`		`- double coeff_beta)`
	`397`	`+static inline void nodeApplyAscOrthoRadialGive(int i_r, int i_theta, const PolarGrid& grid, bool DirBC_Interior,`
	`398`	`+ SmootherColor smoother_color, ConstVector<double>& x,`
	`399`	`+ ConstVector<double>& rhs, Vector<double>& result, double arr, double att,`
	`400`	`+ double art, double detDF, double coeff_beta)`
`401`	`401`	`{`
`402`	`402`	`assert(i_r >= grid.numberSmootherCircles() - 1 && i_r < grid.nr());`
`403`	`403`
Original file line number	Diff line number	Diff line change
`@@ -1,7 +1,8 @@`
`1`	`1`	`#include "../../../include/ExtrapolatedSmoother/ExtrapolatedSmootherTake/extrapolatedSmootherTake.h"`
`2`	`2`
`3`	`3`	`/* Tridiagonal matrices */`
`4`		`-inline void updateTridiagonalElement(SymmetricTridiagonalSolver<double>& matrix, int row, int column, double value)`
	`4`	`+static inline void updateTridiagonalElement(SymmetricTridiagonalSolver<double>& matrix, int row, int column,`
	`5`	`+ double value)`
`5`	`6`	`{`
`6`	`7`	`if (row == column)`
`7`	`8`	`matrix.main_diagonal(row) = value;`
`@@ -12,23 +13,23 @@ inline void updateTridiagonalElement(SymmetricTridiagonalSolver<double>& matrix,`
`12`	`13`	`}`
`13`	`14`
`14`	`15`	`/* Diagonal matrices */`
`15`		`-inline void updateDiagonalElement(DiagonalSolver<double>& matrix, int row, int column, double value)`
	`16`	`+static inline void updateDiagonalElement(DiagonalSolver<double>& matrix, int row, int column, double value)`
`16`	`17`	`{`
`17`	`18`	`matrix.diagonal(row) = value;`
`18`	`19`	`}`
`19`	`20`
`20`	`21`	`/* Inner Boundary COO/CSR matrix */`
`21`	`22`	`#ifdef GMGPOLAR_USE_MUMPS`
`22`		`-inline void updateCOOCSRMatrixElement(SparseMatrixCOO<double>& matrix, int ptr, int offset, int row, int col,`
`23`		`- double val)`
	`23`	`+static inline void updateCOOCSRMatrixElement(SparseMatrixCOO<double>& matrix, int ptr, int offset, int row, int col,`
	`24`	`+ double val)`
`24`	`25`	`{`
`25`	`26`	`matrix.row_index(ptr + offset) = row;`
`26`	`27`	`matrix.col_index(ptr + offset) = col;`
`27`	`28`	`matrix.value(ptr + offset) = val;`
`28`	`29`	`}`
`29`	`30`	`#else`
`30`		`-inline void updateCOOCSRMatrixElement(SparseMatrixCSR<double>& matrix, int ptr, int offset, int row, int col,`
`31`		`- double val)`
	`31`	`+static inline void updateCOOCSRMatrixElement(SparseMatrixCSR<double>& matrix, int ptr, int offset, int row, int col,`
	`32`	`+ double val)`
`32`	`33`	`{`
`33`	`34`	`matrix.row_nz_index(row, offset) = col;`
`34`	`35`	`matrix.row_nz_entry(row, offset) = val;`
Original file line number	Diff line number	Diff line change
`@@ -46,9 +46,9 @@`
`46`	`46`	`*`
`47`	`47`	`*/`
`48`	`48`
`49`		`-inline void fineNodeExtrapolatedProlongation(int i_r, int i_theta, int i_r_coarse, int i_theta_coarse,`
`50`		`- const PolarGrid& coarse_grid, const PolarGrid& fine_grid,`
`51`		`- Vector<double>& fine_result, ConstVector<double>& coarse_values)`
	`49`	`+static inline void fineNodeExtrapolatedProlongation(int i_r, int i_theta, int i_r_coarse, int i_theta_coarse,`
	`50`	`+ const PolarGrid& coarse_grid, const PolarGrid& fine_grid,`
	`51`	`+ Vector<double>& fine_result, ConstVector<double>& coarse_values)`
`52`	`52`	`{`
`53`	`53`	`if (i_r & 1) {`
`54`	`54`	`if (i_theta & 1) { /* (odd, odd) -> node in center of coarse cell */`