diff --git a/test/unit/solver/test_cuda_lu_decomposition.cpp b/test/unit/solver/test_cuda_lu_decomposition.cpp
index d9761e1d6..4c0ea5116 100644
--- a/test/unit/solver/test_cuda_lu_decomposition.cpp
+++ b/test/unit/solver/test_cuda_lu_decomposition.cpp
@@ -68,11 +68,11 @@ void testRandomMatrix(size_t n_grids)
         for (std::size_t i_block = 0; i_block < n_grids; ++i_block)
           A[i_block][i][j] = get_double();
 
-// micm::LuDecomposition cpu_lud(A);
-// auto cpu_LU = micm::LuDecomposition::GetLUMatrices(A, 1.0e-30);
-// cpu_lud.Decompose<double, SparseMatrixPolicy>(A, cpu_LU.first, cpu_LU.second);
-// check_results<double, SparseMatrixPolicy>(
-//     A, cpu_LU.first, cpu_LU.second, [&](const double a, const double b) -> void { EXPECT_NEAR(a, b, 1.0e-5); });
+micm::LuDecomposition cpu_lud(A);
+auto cpu_LU = micm::LuDecomposition::GetLUMatrices(A, 1.0e-30);
+cpu_lud.Decompose<double, SparseMatrixPolicy>(A, cpu_LU.first, cpu_LU.second);
+check_results<double, SparseMatrixPolicy>(
+    A, cpu_LU.first, cpu_LU.second, [&](const double a, const double b) -> void { EXPECT_NEAR(a, b, 1.0e-5); });
  
   micm::CUDALuDecomposition gpu_lud(A); 
   auto gpu_LU = micm::CUDALuDecomposition::GetLUMatrices(A, 1.0e-30);