From c6a058bae5f5f326d84b5b2de339b0c37a802060 Mon Sep 17 00:00:00 2001
From: Marco Feder <marco.feder@sissa.it>
Date: Thu, 23 May 2024 22:19:33 +0000
Subject: [PATCH] Update piston example

---
 examples/3D_piston.cc | 8 +++-----
 1 file changed, 3 insertions(+), 5 deletions(-)

diff --git a/examples/3D_piston.cc b/examples/3D_piston.cc
index 40a5e06b..14060945 100644
--- a/examples/3D_piston.cc
+++ b/examples/3D_piston.cc
@@ -526,21 +526,19 @@ DiffusionReactionProblem<dim>::assemble_system()
                            QGauss<dim - 1>(face_quadrature_degree),
                            update_JxW_values);
 
-  DynamicSparsityPattern sparsity_pattern;
+  TrilinosWrappers::SparsityPattern sparsity_pattern;
   // double                 start_setup, stop_setup;
   // start_setup = MPI_Wtime();
   ah->distribute_agglomerated_dofs(fe_dg);
   ah->create_agglomeration_sparsity_pattern(sparsity_pattern);
+  sparsity_pattern.compress();
+  system_matrix.reinit(sparsity_pattern);
   // stop_setup = MPI_Wtime();
 
 
   locally_owned_dofs    = ah->agglo_dh.locally_owned_dofs();
   locally_relevant_dofs = DoFTools::extract_locally_relevant_dofs(ah->agglo_dh);
 
-  system_matrix.reinit(locally_owned_dofs,
-                       locally_owned_dofs,
-                       sparsity_pattern,
-                       comm);
   system_rhs.reinit(locally_owned_dofs, comm);
 
   std::unique_ptr<const RightHandSide<dim>> rhs_function;