test

NCAR · Oct 22, 2023 · a685796 · a685796
1 parent 865aebc
commit a685796
Show file tree

Hide file tree

Showing 3 changed files with 4 additions and 4 deletions.
diff --git a/include/micm/solver/cuda_linear_solver.hpp b/include/micm/solver/cuda_linear_solver.hpp
@@ -8,12 +8,12 @@
 #include <micm/solver/cuda_linear_solver.cuh>
 namespace micm{
     template<typename T, template<class> class SparseMatrixPolicy, class LuDecompositionPolicy = CudaLuDecomposition>
-    class CudaLinearSolver: public LinearSolver<T, SparseMatrixPolicy, LuDecomposition> {
+    class CudaLinearSolver : public LinearSolver<T, SparseMatrixPolicy, CudaLuDecomposition> {
     public:
         //constructor
         CudaLinearSolver(){};
 
-        CudaLinearSolver(const SparseMatrixPolicy<T>& matrix, T initial_value): LinearSolver<T, SparseMatrixPolicy, LuDecomposition> (matrix, initial_value){};
+        CudaLinearSolver(const SparseMatrixPolicy<T>& matrix, T initial_value): LinearSolver<T, SparseMatrixPolicy, CudaLuDecomposition> (matrix, initial_value){};
 
         template<template<class> class MatrixPolicy> 
         requires(VectorizableDense<MatrixPolicy<T>> || VectorizableSparse<SparseMatrixPolicy<T>>)

diff --git a/include/micm/solver/linear_solver.hpp b/include/micm/solver/linear_solver.hpp
@@ -24,7 +24,7 @@ namespace micm
   template<typename T, template<class> class SparseMatrixPolicy, class LuDecompositionPolicy = LuDecomposition>
   class LinearSolver
   {
-   public:
+   protected:
     // Parameters needed to calculate L (U x) = b
     //
     // The calculation is split into calculation of L y = b where y = U x:

diff --git a/test/unit/solver/test_cuda_linear_solver.cpp b/test/unit/solver/test_cuda_linear_solver.cpp
@@ -7,9 +7,9 @@
 #include <micm/util/sparse_matrix_vector_ordering.hpp>
 #include <micm/util/vector_matrix.hpp>
 #include <micm/solver/linear_solver.hpp>
+#include <micm/solver/lu_decomposition.hpp>
 #include <micm/solver/cuda_linear_solver.hpp>
 #include <micm/solver/cuda_lu_decomposition.hpp>
-#include <micm/solver/lu_decomposition.hpp>
 #include "test_linear_solver_policy.hpp"
 
 template<class T>