Update piston example

fdrmrc · May 24, 2024 · c6a058b · c6a058b
1 parent 258dd8f
commit c6a058b
Showing 1 changed file with 3 additions and 5 deletions.
diff --git 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;