fix cuda mismatch template

Author: boulderdaze

include/micm/GPU.hpp

@@ -20,6 +20,7 @@
 #include <micm/cuda/util/cuda_param.hpp>
 #include <micm/cuda/util/cuda_sparse_matrix.hpp>
 #include <micm/cuda/util/cuda_util.cuh>
+#include <micm/constraint/constraint_set.hpp>
 #include <micm/solver/solver_builder.hpp>
 
 namespace micm
@@ -30,7 +31,7 @@ namespace micm
   using GpuState = CudaState<CudaDenseMatrixVector, CudaSparseMatrixVector, CudaLuDecompositionMozartInPlace>;
 
   using CudaRosenbrockVectorType = typename CudaRosenbrockSolverParameters::
-      template SolverType<CudaProcessSet<CudaDenseMatrixVector, CudaSparseMatrixVector>, CudaLinearSolverInPlace<CudaSparseMatrixVector>>;
+      template SolverType<CudaProcessSet<CudaDenseMatrixVector, CudaSparseMatrixVector>, CudaLinearSolverInPlace<CudaSparseMatrixVector>, ConstraintSet<CudaDenseMatrixVector, CudaSparseMatrixVector>>;
   using CudaRosenbrock = Solver<CudaRosenbrockVectorType, GpuState>;
 
   using GpuRosenbrockThreeStageBuilder = CudaSolverBuilderInPlace<CudaRosenbrockSolverParameters>;

include/micm/cuda/solver/cuda_solver_parameters.hpp

@@ -12,8 +12,8 @@ namespace micm
   /// @brief Parameters for the CUDA Rosenbrock solver
   struct CudaRosenbrockSolverParameters : public RosenbrockSolverParameters
   {
-    template<class RatesPolicy, class LinearSolverPolicy>
-    using SolverType = CudaRosenbrockSolver<RatesPolicy, LinearSolverPolicy>;
+    template<class RatesPolicy, class LinearSolverPolicy, class ConstraintSetPolicy>
+    using SolverType = CudaRosenbrockSolver<RatesPolicy, LinearSolverPolicy, ConstraintSetPolicy>;
 
     /// @brief Constructor from base class
     /// @param base