Fix iterator to pointer issue in ublas_space.h (#12761)

KratosMultiphysics · Oct 17, 2024 · 574bd67 · 574bd67
1 parent d18b2c4
commit 574bd67
Showing 1 changed file with 11 additions and 11 deletions.
diff --git a/kratos/spaces/ublas_space.h b/kratos/spaces/ublas_space.h
@@ -710,8 +710,8 @@ class UblasSpace
     {
         const std::size_t system_size = rA.size1();
 
-        const double* Avalues = rA.value_data().begin();
-        const std::size_t* Arow_indices = rA.index1_data().begin();
+        const auto& Avalues = rA.value_data();
+        const auto& Arow_indices = rA.index1_data();
 
         // Define  zero value tolerance
         const double zero_tolerance = std::numeric_limits<double>::epsilon();
@@ -778,9 +778,9 @@ class UblasSpace
      */
     static double GetDiagonalNorm(const MatrixType& rA)
     {
-        const double* Avalues = rA.value_data().begin();
-        const std::size_t* Arow_indices = rA.index1_data().begin();
-        const std::size_t* Acol_indices = rA.index2_data().begin();
+        const auto& Avalues = rA.value_data();
+        const auto& Arow_indices = rA.index1_data();
+        const auto& Acol_indices = rA.index2_data();
 
         const double diagonal_norm = IndexPartition<std::size_t>(Size1(rA)).for_each<SumReduction<double>>([&](std::size_t Index){
             const std::size_t col_begin = Arow_indices[Index];
@@ -813,9 +813,9 @@ class UblasSpace
      */
     static double GetMaxDiagonal(const MatrixType& rA)
     {
-        const double* Avalues = rA.value_data().begin();
-        const std::size_t* Arow_indices = rA.index1_data().begin();
-        const std::size_t* Acol_indices = rA.index2_data().begin();
+        const auto& Avalues = rA.value_data();
+        const auto& Arow_indices = rA.index1_data();
+        const auto& Acol_indices = rA.index2_data();
 
         return IndexPartition<std::size_t>(Size1(rA)).for_each<MaxReduction<double>>([&](std::size_t Index){
             const std::size_t col_begin = Arow_indices[Index];
@@ -836,9 +836,9 @@ class UblasSpace
      */
     static double GetMinDiagonal(const MatrixType& rA)
     {
-        const double* Avalues = rA.value_data().begin();
-        const std::size_t* Arow_indices = rA.index1_data().begin();
-        const std::size_t* Acol_indices = rA.index2_data().begin();
+        const auto& Avalues = rA.value_data();
+        const auto& Arow_indices = rA.index1_data();
+        const auto& Acol_indices = rA.index2_data();
 
         return IndexPartition<std::size_t>(Size1(rA)).for_each<MinReduction<double>>([&](std::size_t Index){
             const std::size_t col_begin = Arow_indices[Index];