From: Wolfgang Bangerth <bangerth@colostate.edu>
Date: Sun, 2 Jul 2023 16:18:05 +0000 (-0600)
Subject: Get it closer with hanging nodes.
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=f23b7974a9f0766a71690122b9522b0c9bcd772e;p=dealii.git

Get it closer with hanging nodes.
---

diff --git a/examples/step-86/step-86.cc b/examples/step-86/step-86.cc
index 6981ae53d4..4ab5db30ef 100644
--- a/examples/step-86/step-86.cc
+++ b/examples/step-86/step-86.cc
@@ -199,12 +199,10 @@ namespace Step86
     laplace_matrix.reinit(dof_handler.locally_owned_dofs(), dsp, MPI_COMM_SELF);
     system_matrix.reinit(dof_handler.locally_owned_dofs(), dsp, MPI_COMM_SELF);
 
-    MatrixCreator::create_mass_matrix(dof_handler,
-                                      QGauss<dim>(fe.degree + 1),
-                                      mass_matrix);
-    MatrixCreator::create_laplace_matrix(dof_handler,
-                                         QGauss<dim>(fe.degree + 1),
-                                         laplace_matrix);
+    MatrixCreator::create_mass_matrix(
+      dof_handler, QGauss<dim>(fe.degree + 1), mass_matrix, {}, constraints);
+    MatrixCreator::create_laplace_matrix(
+      dof_handler, QGauss<dim>(fe.degree + 1), laplace_matrix, {}, constraints);
 
     solution.reinit(dof_handler.locally_owned_dofs(), MPI_COMM_SELF);
     old_solution.reinit(dof_handler.locally_owned_dofs(), MPI_COMM_SELF);
@@ -245,9 +243,9 @@ namespace Step86
     data_out.set_flags(DataOutBase::VtkFlags(time, timestep_number));
 
     const std::string filename =
-      "solution-" + Utilities::int_to_string(timestep_number, 3) + ".vtk";
+      "solution-" + Utilities::int_to_string(timestep_number, 3) + ".vtu";
     std::ofstream output(filename);
-    data_out.write_vtk(output);
+    data_out.write_vtu(output);
   }
 
 
@@ -362,7 +360,7 @@ namespace Step86
         system_matrix.copy_from(mass_matrix);
         system_matrix.add(theta * time_step, laplace_matrix);
 
-        //        constraints.condense(system_matrix, system_rhs);
+        constraints.condense(system_rhs);
 
         {
           BoundaryValues<dim> boundary_values_function;