N_Vector: Support PETSc (Block)Vector

Rewrite some elementwise operations to work with PETSc as well

diff --git a/include/deal.II/sundials/n_vector.templates.h b/include/deal.II/sundials/n_vector.templates.h
index 413ffb45a4df31733fa1d2044ef537742ce7c66c..a45c4683f9636371d588b56c53c79f5b8bed5f50 100644 (file)
--- a/include/deal.II/sundials/n_vector.templates.h
+++ b/include/deal.II/sundials/n_vector.templates.h
@@ -28,6 +28,8 @@
 #  include <deal.II/lac/la_parallel_block_vector.h>
 #  include <deal.II/lac/la_parallel_vector.h>
 #  include <deal.II/lac/la_vector.h>
+#  include <deal.II/lac/petsc_block_vector.h>
+#  include <deal.II/lac/petsc_vector.h>
 #  include <deal.II/lac/trilinos_parallel_block_vector.h>
 #  include <deal.II/lac/trilinos_vector.h>
 #  include <deal.II/lac/vector_memory.h>
@@ -651,9 +653,21 @@ template <typename VectorType,
 realtype
 SUNDIALS::internal::NVectorOperations::min_element(N_Vector x)
 {
-  auto *     vector    = unwrap_nvector_const<VectorType>(x);
-  const auto local_min = *std::min_element(vector->begin(), vector->end());
-  return Utilities::MPI::min(local_min, get_communicator<VectorType>(x));
+  auto *vector = unwrap_nvector_const<VectorType>(x);
+
+
+  const auto indexed_less_than = [&](const IndexSet::size_type idxa,
+                                     const IndexSet::size_type idxb) {
+    return (*vector)[idxa] < (*vector)[idxb];
+  };
+
+  auto local_elements = vector->locally_owned_elements();
+
+  const auto local_min = *std::min_element(local_elements.begin(),
+                                           local_elements.end(),
+                                           indexed_less_than);
+  return Utilities::MPI::min((*vector)[local_min],
+                             get_communicator<VectorType>(x));
 }
 
 
@@ -673,12 +687,22 @@ SUNDIALS::internal::NVectorOperations::min_element(N_Vector x)
 
   for (unsigned i = 0; i < vector->n_blocks(); ++i)
     {
+      const auto indexed_less_than = [&](const IndexSet::size_type idxa,
+                                         const IndexSet::size_type idxb) {
+        return vector->block(i)[idxa] < vector->block(i)[idxb];
+      };
+
+      auto local_elements = vector->block(i).locally_owned_elements();
+
       const auto block_local_min_element =
-        std::min_element(vector->block(i).begin(), vector->block(i).end());
+        std::min_element(local_elements.begin(),
+                         local_elements.end(),
+                         indexed_less_than);
 
       // guard against empty blocks on this processor
-      if (block_local_min_element != vector->block(i).end())
-        proc_local_min = std::min(proc_local_min, *block_local_min_element);
+      if (block_local_min_element != local_elements.end())
+        proc_local_min =
+          std::min(proc_local_min, vector->block(i)[*block_local_min_element]);
     }
 
   return Utilities::MPI::min(proc_local_min, get_communicator<VectorType>(x));
@@ -715,11 +739,12 @@ SUNDIALS::internal::NVectorOperations::elementwise_div(N_Vector x,
   AssertDimension(x_dealii->size(), z_dealii->size());
   AssertDimension(x_dealii->size(), y_dealii->size());
 
-  std::transform(x_dealii->begin(),
-                 x_dealii->end(),
-                 y_dealii->begin(),
-                 z_dealii->begin(),
-                 std::divides<>{});
+  auto x_ele = x_dealii->locally_owned_elements();
+  for (const auto idx : x_ele)
+    {
+      (*z_dealii)[idx] = (*x_dealii)[idx] / (*y_dealii)[idx];
+    }
+  z_dealii->compress(VectorOperation::insert);
 }
 
 
@@ -741,11 +766,15 @@ SUNDIALS::internal::NVectorOperations::elementwise_div(N_Vector x,
   AssertDimension(x_dealii->n_blocks(), y_dealii->n_blocks());
 
   for (unsigned i = 0; i < x_dealii->n_blocks(); ++i)
-    std::transform(x_dealii->block(i).begin(),
-                   x_dealii->block(i).end(),
-                   y_dealii->block(i).begin(),
-                   z_dealii->block(i).begin(),
-                   std::divides<>{});
+    {
+      auto x_ele = x_dealii->block(i).locally_owned_elements();
+      for (const auto idx : x_ele)
+        {
+          z_dealii->block(i)[idx] =
+            x_dealii->block(i)[idx] / y_dealii->block(i)[idx];
+        }
+    }
+  z_dealii->compress(VectorOperation::insert);
 }
 
 
@@ -760,10 +789,12 @@ SUNDIALS::internal::NVectorOperations::elementwise_inv(N_Vector x, N_Vector z)
 
   AssertDimension(x_dealii->size(), z_dealii->size());
 
-  std::transform(x_dealii->begin(),
-                 x_dealii->end(),
-                 z_dealii->begin(),
-                 [](typename VectorType::value_type v) { return 1.0 / v; });
+  auto x_ele = x_dealii->locally_owned_elements();
+  for (const auto idx : x_ele)
+    {
+      (*z_dealii)[idx] = 1.0 / (*x_dealii)[idx];
+    }
+  z_dealii->compress(VectorOperation::insert);
 }
 
 
@@ -780,10 +811,15 @@ SUNDIALS::internal::NVectorOperations::elementwise_inv(N_Vector x, N_Vector z)
   AssertDimension(x_dealii->n_blocks(), z_dealii->n_blocks());
 
   for (unsigned i = 0; i < x_dealii->n_blocks(); ++i)
-    std::transform(x_dealii->block(i).begin(),
-                   x_dealii->block(i).end(),
-                   z_dealii->block(i).begin(),
-                   [](typename VectorType::value_type v) { return 1.0 / v; });
+    {
+      auto x_ele = x_dealii->block(i).locally_owned_elements();
+      for (const auto idx : x_ele)
+        {
+          z_dealii->block(i)[idx] = 1.0 / x_dealii->block(i)[idx];
+        }
+    }
+
+  z_dealii->compress(VectorOperation::insert);
 }
 
 
@@ -798,12 +834,13 @@ SUNDIALS::internal::NVectorOperations::elementwise_abs(N_Vector x, N_Vector z)
 
   AssertDimension(x_dealii->size(), z_dealii->size());
 
-  std::transform(x_dealii->begin(),
-                 x_dealii->end(),
-                 z_dealii->begin(),
-                 [](typename VectorType::value_type v) {
-                   return std::fabs(v);
-                 });
+  auto x_ele = x_dealii->locally_owned_elements();
+  for (const auto idx : x_ele)
+    {
+      (*z_dealii)[idx] = std::fabs((*x_dealii)[idx]);
+    }
+
+  z_dealii->compress(VectorOperation::insert);
 }
 
 
@@ -820,12 +857,15 @@ SUNDIALS::internal::NVectorOperations::elementwise_abs(N_Vector x, N_Vector z)
   AssertDimension(x_dealii->n_blocks(), z_dealii->n_blocks());
 
   for (unsigned i = 0; i < x_dealii->n_blocks(); ++i)
-    std::transform(x_dealii->block(i).begin(),
-                   x_dealii->block(i).end(),
-                   z_dealii->block(i).begin(),
-                   [](typename VectorType::value_type v) {
-                     return std::fabs(v);
-                   });
+    {
+      auto x_ele = x_dealii->block(i).locally_owned_elements();
+      for (const auto idx : x_ele)
+        {
+          z_dealii->block(i)[idx] = std::fabs(x_dealii->block(i)[idx]);
+        }
+    }
+
+  z_dealii->compress(VectorOperation::insert);
 }
 
 

diff --git a/source/sundials/n_vector.inst.in b/source/sundials/n_vector.inst.in
index 0255496f81feb0bd8198173073891b1804571411..256d649bd06de82de91424169b13caaa56de8dfa 100644 (file)
--- a/source/sundials/n_vector.inst.in
+++ b/source/sundials/n_vector.inst.in
@@ -68,3 +68,21 @@ for (V : DEAL_II_VEC_TEMPLATES)
       const TrilinosWrappers::MPI::V>;
   }
 #endif
+
+#ifdef DEAL_II_WITH_PETSC
+// PETScWrappers Vector and BlockVector
+for (V : DEAL_II_VEC_TEMPLATES)
+  {
+    template SUNDIALS::internal::NVectorView<PETScWrappers::MPI::V>
+    SUNDIALS::internal::make_nvector_view<>(PETScWrappers::MPI::V &);
+    template SUNDIALS::internal::NVectorView<const PETScWrappers::MPI::V>
+    SUNDIALS::internal::make_nvector_view<>(const PETScWrappers::MPI::V &);
+    template PETScWrappers::MPI::V
+      *SUNDIALS::internal::unwrap_nvector<PETScWrappers::MPI::V>(N_Vector);
+    template const PETScWrappers::MPI::V *
+      SUNDIALS::internal::unwrap_nvector_const<PETScWrappers::MPI::V>(N_Vector);
+
+    template class SUNDIALS::internal::NVectorView<PETScWrappers::MPI::V>;
+    template class SUNDIALS::internal::NVectorView<const PETScWrappers::MPI::V>;
+  }
+#endif
\ No newline at end of file

diff --git a/tests/sundials/n_vector.cc b/tests/sundials/n_vector.cc
index e8cc368632be3b36ec545e5ebdaa07bb87616919..89e5414129a0a3b43e25f1ba26dcbfe9d44ea0f5 100644 (file)
--- a/tests/sundials/n_vector.cc
+++ b/tests/sundials/n_vector.cc
@@ -132,6 +132,47 @@ namespace
     return vector;
   }
 
+  template <>
+  TrilinosWrappers::MPI::BlockVector
+  create_test_vector(double value)
+  {
+    const unsigned n_processes =
+      Utilities::MPI::n_mpi_processes(MPI_COMM_WORLD);
+    const unsigned block_size =
+      ((n_processes - 1) * (n_processes)) / 2 + n_processes;
+    const auto partitioning =
+      create_parallel_index_set().split_by_block({block_size, block_size});
+    TrilinosWrappers::MPI::BlockVector vector(partitioning, MPI_COMM_WORLD);
+    vector = value;
+    return vector;
+  }
+
+  template <>
+  PETScWrappers::MPI::Vector
+  create_test_vector(double value)
+  {
+    IndexSet local_dofs = create_parallel_index_set();
+
+    PETScWrappers::MPI::Vector vector(local_dofs, MPI_COMM_WORLD);
+    vector = value;
+    return vector;
+  }
+
+  template <>
+  PETScWrappers::MPI::BlockVector
+  create_test_vector(double value)
+  {
+    const unsigned n_processes =
+      Utilities::MPI::n_mpi_processes(MPI_COMM_WORLD);
+    const unsigned block_size =
+      ((n_processes - 1) * (n_processes)) / 2 + n_processes;
+    const auto partitioning =
+      create_parallel_index_set().split_by_block({block_size, block_size});
+    PETScWrappers::MPI::BlockVector vector(partitioning, MPI_COMM_WORLD);
+    vector = value;
+    return vector;
+  }
+
   template <typename VectorType>
   bool
   vector_equal(const VectorType &a, const VectorType &b)
@@ -575,8 +616,6 @@ test_min_element()
   auto nv_b = make_nvector_view(vector_b);
 
   auto result = N_VMin(nv_a);
-  std::cout << "result: " << result << std::endl;
-  vector_a.print(std::cout);
   Assert(std::fabs(result - 2.0) < 1e-12, NVectorTestError());
 
   result = N_VMin(nv_b);
@@ -654,4 +693,15 @@ main(int argc, char **argv)
   run_all_tests<LinearAlgebra::distributed::BlockVector<double>>(
     "LinearAlgebra::distributed::BlockVector<double>");
   run_all_tests<TrilinosWrappers::MPI::Vector>("TrilinosWrappers::MPI::Vector");
+  run_all_tests<TrilinosWrappers::MPI::BlockVector>(
+    "TrilinosWrappers::MPI::BlockVector");
+  run_all_tests<PETScWrappers::MPI::Vector>("PETScWrappers::MPI::Vector");
+  run_all_tests<PETScWrappers::MPI::BlockVector>(
+    "PETScWrappers::MPI::BlockVector");
+
+  // although the memory would be cleared in ~MPI_InitFinalize it needs to be
+  // done manually to satisfy the PETSc memory check inside ~MPILogInitAll,
+  // which is invoked first
+  GrowingVectorMemory<PETScWrappers::MPI::Vector>::release_unused_memory();
+  GrowingVectorMemory<PETScWrappers::MPI::BlockVector>::release_unused_memory();
 }

diff --git a/tests/sundials/n_vector.with_sundials.geq.5.0.0.with_mpi=true.with_trilinos=true.mpirun=2.output b/tests/sundials/n_vector.with_sundials.geq.5.0.0.with_mpi=true.mpirun=2.with_trilinos=true.with_petsc=true.output
similarity index 62%
rename from tests/sundials/n_vector.with_sundials.geq.5.0.0.with_mpi=true.with_trilinos=true.mpirun=2.output
rename to tests/sundials/n_vector.with_sundials.geq.5.0.0.with_mpi=true.mpirun=2.with_trilinos=true.with_petsc=true.output
index 391c5897897ada83faf7daca6d07f8bc33b9d810..137c7ab77ce6ab160833fda91150b57909b09edd 100644 (file)
--- a/tests/sundials/n_vector.with_sundials.geq.5.0.0.with_mpi=true.with_trilinos=true.mpirun=2.output
+++ b/tests/sundials/n_vector.with_sundials.geq.5.0.0.with_mpi=true.mpirun=2.with_trilinos=true.with_petsc=true.output
@@ -89,6 +89,60 @@ DEAL:0:TrilinosWrappers::MPI::Vector::test_weighted_rms_norm OK
 DEAL:0:TrilinosWrappers::MPI::Vector::test_max_norm OK
 DEAL:0:TrilinosWrappers::MPI::Vector::test_min_element OK
 DEAL:0:TrilinosWrappers::MPI::Vector::test_scale OK
+DEAL:0:TrilinosWrappers::MPI::BlockVector::test_nvector_view_unwrap OK
+DEAL:0:TrilinosWrappers::MPI::BlockVector::test_get_vector_id OK
+DEAL:0:TrilinosWrappers::MPI::BlockVector::test_clone OK
+DEAL:0:TrilinosWrappers::MPI::BlockVector::test_destroy OK
+DEAL:0:TrilinosWrappers::MPI::BlockVector::test_get_communicator OK
+DEAL:0:TrilinosWrappers::MPI::BlockVector::test_length OK
+DEAL:0:TrilinosWrappers::MPI::BlockVector::test_linear_sum OK
+DEAL:0:TrilinosWrappers::MPI::BlockVector::test_dot_product OK
+DEAL:0:TrilinosWrappers::MPI::BlockVector::test_set_constant OK
+DEAL:0:TrilinosWrappers::MPI::BlockVector::test_add_constant OK
+DEAL:0:TrilinosWrappers::MPI::BlockVector::test_elementwise_product OK
+DEAL:0:TrilinosWrappers::MPI::BlockVector::test_elementwise_div OK
+DEAL:0:TrilinosWrappers::MPI::BlockVector::test_elementwise_inv OK
+DEAL:0:TrilinosWrappers::MPI::BlockVector::test_elementwise_abs OK
+DEAL:0:TrilinosWrappers::MPI::BlockVector::test_weighted_rms_norm OK
+DEAL:0:TrilinosWrappers::MPI::BlockVector::test_max_norm OK
+DEAL:0:TrilinosWrappers::MPI::BlockVector::test_min_element OK
+DEAL:0:TrilinosWrappers::MPI::BlockVector::test_scale OK
+DEAL:0:PETScWrappers::MPI::Vector::test_nvector_view_unwrap OK
+DEAL:0:PETScWrappers::MPI::Vector::test_get_vector_id OK
+DEAL:0:PETScWrappers::MPI::Vector::test_clone OK
+DEAL:0:PETScWrappers::MPI::Vector::test_destroy OK
+DEAL:0:PETScWrappers::MPI::Vector::test_get_communicator OK
+DEAL:0:PETScWrappers::MPI::Vector::test_length OK
+DEAL:0:PETScWrappers::MPI::Vector::test_linear_sum OK
+DEAL:0:PETScWrappers::MPI::Vector::test_dot_product OK
+DEAL:0:PETScWrappers::MPI::Vector::test_set_constant OK
+DEAL:0:PETScWrappers::MPI::Vector::test_add_constant OK
+DEAL:0:PETScWrappers::MPI::Vector::test_elementwise_product OK
+DEAL:0:PETScWrappers::MPI::Vector::test_elementwise_div OK
+DEAL:0:PETScWrappers::MPI::Vector::test_elementwise_inv OK
+DEAL:0:PETScWrappers::MPI::Vector::test_elementwise_abs OK
+DEAL:0:PETScWrappers::MPI::Vector::test_weighted_rms_norm OK
+DEAL:0:PETScWrappers::MPI::Vector::test_max_norm OK
+DEAL:0:PETScWrappers::MPI::Vector::test_min_element OK
+DEAL:0:PETScWrappers::MPI::Vector::test_scale OK
+DEAL:0:PETScWrappers::MPI::BlockVector::test_nvector_view_unwrap OK
+DEAL:0:PETScWrappers::MPI::BlockVector::test_get_vector_id OK
+DEAL:0:PETScWrappers::MPI::BlockVector::test_clone OK
+DEAL:0:PETScWrappers::MPI::BlockVector::test_destroy OK
+DEAL:0:PETScWrappers::MPI::BlockVector::test_get_communicator OK
+DEAL:0:PETScWrappers::MPI::BlockVector::test_length OK
+DEAL:0:PETScWrappers::MPI::BlockVector::test_linear_sum OK
+DEAL:0:PETScWrappers::MPI::BlockVector::test_dot_product OK
+DEAL:0:PETScWrappers::MPI::BlockVector::test_set_constant OK
+DEAL:0:PETScWrappers::MPI::BlockVector::test_add_constant OK
+DEAL:0:PETScWrappers::MPI::BlockVector::test_elementwise_product OK
+DEAL:0:PETScWrappers::MPI::BlockVector::test_elementwise_div OK
+DEAL:0:PETScWrappers::MPI::BlockVector::test_elementwise_inv OK
+DEAL:0:PETScWrappers::MPI::BlockVector::test_elementwise_abs OK
+DEAL:0:PETScWrappers::MPI::BlockVector::test_weighted_rms_norm OK
+DEAL:0:PETScWrappers::MPI::BlockVector::test_max_norm OK
+DEAL:0:PETScWrappers::MPI::BlockVector::test_min_element OK
+DEAL:0:PETScWrappers::MPI::BlockVector::test_scale OK
 
 DEAL:1:Vector<double>::test_nvector_view_unwrap OK
 DEAL:1:Vector<double>::test_get_vector_id OK
@@ -180,4 +234,58 @@ DEAL:1:TrilinosWrappers::MPI::Vector::test_weighted_rms_norm OK
 DEAL:1:TrilinosWrappers::MPI::Vector::test_max_norm OK
 DEAL:1:TrilinosWrappers::MPI::Vector::test_min_element OK
 DEAL:1:TrilinosWrappers::MPI::Vector::test_scale OK
+DEAL:1:TrilinosWrappers::MPI::BlockVector::test_nvector_view_unwrap OK
+DEAL:1:TrilinosWrappers::MPI::BlockVector::test_get_vector_id OK
+DEAL:1:TrilinosWrappers::MPI::BlockVector::test_clone OK
+DEAL:1:TrilinosWrappers::MPI::BlockVector::test_destroy OK
+DEAL:1:TrilinosWrappers::MPI::BlockVector::test_get_communicator OK
+DEAL:1:TrilinosWrappers::MPI::BlockVector::test_length OK
+DEAL:1:TrilinosWrappers::MPI::BlockVector::test_linear_sum OK
+DEAL:1:TrilinosWrappers::MPI::BlockVector::test_dot_product OK
+DEAL:1:TrilinosWrappers::MPI::BlockVector::test_set_constant OK
+DEAL:1:TrilinosWrappers::MPI::BlockVector::test_add_constant OK
+DEAL:1:TrilinosWrappers::MPI::BlockVector::test_elementwise_product OK
+DEAL:1:TrilinosWrappers::MPI::BlockVector::test_elementwise_div OK
+DEAL:1:TrilinosWrappers::MPI::BlockVector::test_elementwise_inv OK
+DEAL:1:TrilinosWrappers::MPI::BlockVector::test_elementwise_abs OK
+DEAL:1:TrilinosWrappers::MPI::BlockVector::test_weighted_rms_norm OK
+DEAL:1:TrilinosWrappers::MPI::BlockVector::test_max_norm OK
+DEAL:1:TrilinosWrappers::MPI::BlockVector::test_min_element OK
+DEAL:1:TrilinosWrappers::MPI::BlockVector::test_scale OK
+DEAL:1:PETScWrappers::MPI::Vector::test_nvector_view_unwrap OK
+DEAL:1:PETScWrappers::MPI::Vector::test_get_vector_id OK
+DEAL:1:PETScWrappers::MPI::Vector::test_clone OK
+DEAL:1:PETScWrappers::MPI::Vector::test_destroy OK
+DEAL:1:PETScWrappers::MPI::Vector::test_get_communicator OK
+DEAL:1:PETScWrappers::MPI::Vector::test_length OK
+DEAL:1:PETScWrappers::MPI::Vector::test_linear_sum OK
+DEAL:1:PETScWrappers::MPI::Vector::test_dot_product OK
+DEAL:1:PETScWrappers::MPI::Vector::test_set_constant OK
+DEAL:1:PETScWrappers::MPI::Vector::test_add_constant OK
+DEAL:1:PETScWrappers::MPI::Vector::test_elementwise_product OK
+DEAL:1:PETScWrappers::MPI::Vector::test_elementwise_div OK
+DEAL:1:PETScWrappers::MPI::Vector::test_elementwise_inv OK
+DEAL:1:PETScWrappers::MPI::Vector::test_elementwise_abs OK
+DEAL:1:PETScWrappers::MPI::Vector::test_weighted_rms_norm OK
+DEAL:1:PETScWrappers::MPI::Vector::test_max_norm OK
+DEAL:1:PETScWrappers::MPI::Vector::test_min_element OK
+DEAL:1:PETScWrappers::MPI::Vector::test_scale OK
+DEAL:1:PETScWrappers::MPI::BlockVector::test_nvector_view_unwrap OK
+DEAL:1:PETScWrappers::MPI::BlockVector::test_get_vector_id OK
+DEAL:1:PETScWrappers::MPI::BlockVector::test_clone OK
+DEAL:1:PETScWrappers::MPI::BlockVector::test_destroy OK
+DEAL:1:PETScWrappers::MPI::BlockVector::test_get_communicator OK
+DEAL:1:PETScWrappers::MPI::BlockVector::test_length OK
+DEAL:1:PETScWrappers::MPI::BlockVector::test_linear_sum OK
+DEAL:1:PETScWrappers::MPI::BlockVector::test_dot_product OK
+DEAL:1:PETScWrappers::MPI::BlockVector::test_set_constant OK
+DEAL:1:PETScWrappers::MPI::BlockVector::test_add_constant OK
+DEAL:1:PETScWrappers::MPI::BlockVector::test_elementwise_product OK
+DEAL:1:PETScWrappers::MPI::BlockVector::test_elementwise_div OK
+DEAL:1:PETScWrappers::MPI::BlockVector::test_elementwise_inv OK
+DEAL:1:PETScWrappers::MPI::BlockVector::test_elementwise_abs OK
+DEAL:1:PETScWrappers::MPI::BlockVector::test_weighted_rms_norm OK
+DEAL:1:PETScWrappers::MPI::BlockVector::test_max_norm OK
+DEAL:1:PETScWrappers::MPI::BlockVector::test_min_element OK
+DEAL:1:PETScWrappers::MPI::BlockVector::test_scale OK
 

diff --git a/tests/sundials/n_vector.with_sundials.geq.5.0.0.output b/tests/sundials/n_vector.with_sundials.geq.5.0.0.with_trilinos=true.with_petsc=true.output
similarity index 62%
rename from tests/sundials/n_vector.with_sundials.geq.5.0.0.output
rename to tests/sundials/n_vector.with_sundials.geq.5.0.0.with_trilinos=true.with_petsc=true.output
index b9285624988886a35ae653e1a6f26b52334b8088..7cb9f49ebdc80d98eef3f67bdb91c276f70eb312 100644 (file)
--- a/tests/sundials/n_vector.with_sundials.geq.5.0.0.output
+++ b/tests/sundials/n_vector.with_sundials.geq.5.0.0.with_trilinos=true.with_petsc=true.output
@@ -89,3 +89,57 @@ DEAL:0:TrilinosWrappers::MPI::Vector::test_weighted_rms_norm OK
 DEAL:0:TrilinosWrappers::MPI::Vector::test_max_norm OK
 DEAL:0:TrilinosWrappers::MPI::Vector::test_min_element OK
 DEAL:0:TrilinosWrappers::MPI::Vector::test_scale OK
+DEAL:0:TrilinosWrappers::MPI::BlockVector::test_nvector_view_unwrap OK
+DEAL:0:TrilinosWrappers::MPI::BlockVector::test_get_vector_id OK
+DEAL:0:TrilinosWrappers::MPI::BlockVector::test_clone OK
+DEAL:0:TrilinosWrappers::MPI::BlockVector::test_destroy OK
+DEAL:0:TrilinosWrappers::MPI::BlockVector::test_get_communicator OK
+DEAL:0:TrilinosWrappers::MPI::BlockVector::test_length OK
+DEAL:0:TrilinosWrappers::MPI::BlockVector::test_linear_sum OK
+DEAL:0:TrilinosWrappers::MPI::BlockVector::test_dot_product OK
+DEAL:0:TrilinosWrappers::MPI::BlockVector::test_set_constant OK
+DEAL:0:TrilinosWrappers::MPI::BlockVector::test_add_constant OK
+DEAL:0:TrilinosWrappers::MPI::BlockVector::test_elementwise_product OK
+DEAL:0:TrilinosWrappers::MPI::BlockVector::test_elementwise_div OK
+DEAL:0:TrilinosWrappers::MPI::BlockVector::test_elementwise_inv OK
+DEAL:0:TrilinosWrappers::MPI::BlockVector::test_elementwise_abs OK
+DEAL:0:TrilinosWrappers::MPI::BlockVector::test_weighted_rms_norm OK
+DEAL:0:TrilinosWrappers::MPI::BlockVector::test_max_norm OK
+DEAL:0:TrilinosWrappers::MPI::BlockVector::test_min_element OK
+DEAL:0:TrilinosWrappers::MPI::BlockVector::test_scale OK
+DEAL:0:PETScWrappers::MPI::Vector::test_nvector_view_unwrap OK
+DEAL:0:PETScWrappers::MPI::Vector::test_get_vector_id OK
+DEAL:0:PETScWrappers::MPI::Vector::test_clone OK
+DEAL:0:PETScWrappers::MPI::Vector::test_destroy OK
+DEAL:0:PETScWrappers::MPI::Vector::test_get_communicator OK
+DEAL:0:PETScWrappers::MPI::Vector::test_length OK
+DEAL:0:PETScWrappers::MPI::Vector::test_linear_sum OK
+DEAL:0:PETScWrappers::MPI::Vector::test_dot_product OK
+DEAL:0:PETScWrappers::MPI::Vector::test_set_constant OK
+DEAL:0:PETScWrappers::MPI::Vector::test_add_constant OK
+DEAL:0:PETScWrappers::MPI::Vector::test_elementwise_product OK
+DEAL:0:PETScWrappers::MPI::Vector::test_elementwise_div OK
+DEAL:0:PETScWrappers::MPI::Vector::test_elementwise_inv OK
+DEAL:0:PETScWrappers::MPI::Vector::test_elementwise_abs OK
+DEAL:0:PETScWrappers::MPI::Vector::test_weighted_rms_norm OK
+DEAL:0:PETScWrappers::MPI::Vector::test_max_norm OK
+DEAL:0:PETScWrappers::MPI::Vector::test_min_element OK
+DEAL:0:PETScWrappers::MPI::Vector::test_scale OK
+DEAL:0:PETScWrappers::MPI::BlockVector::test_nvector_view_unwrap OK
+DEAL:0:PETScWrappers::MPI::BlockVector::test_get_vector_id OK
+DEAL:0:PETScWrappers::MPI::BlockVector::test_clone OK
+DEAL:0:PETScWrappers::MPI::BlockVector::test_destroy OK
+DEAL:0:PETScWrappers::MPI::BlockVector::test_get_communicator OK
+DEAL:0:PETScWrappers::MPI::BlockVector::test_length OK
+DEAL:0:PETScWrappers::MPI::BlockVector::test_linear_sum OK
+DEAL:0:PETScWrappers::MPI::BlockVector::test_dot_product OK
+DEAL:0:PETScWrappers::MPI::BlockVector::test_set_constant OK
+DEAL:0:PETScWrappers::MPI::BlockVector::test_add_constant OK
+DEAL:0:PETScWrappers::MPI::BlockVector::test_elementwise_product OK
+DEAL:0:PETScWrappers::MPI::BlockVector::test_elementwise_div OK
+DEAL:0:PETScWrappers::MPI::BlockVector::test_elementwise_inv OK
+DEAL:0:PETScWrappers::MPI::BlockVector::test_elementwise_abs OK
+DEAL:0:PETScWrappers::MPI::BlockVector::test_weighted_rms_norm OK
+DEAL:0:PETScWrappers::MPI::BlockVector::test_max_norm OK
+DEAL:0:PETScWrappers::MPI::BlockVector::test_min_element OK
+DEAL:0:PETScWrappers::MPI::BlockVector::test_scale OK