From f23b7974a9f0766a71690122b9522b0c9bcd772e Mon Sep 17 00:00:00 2001 From: Wolfgang Bangerth Date: Sun, 2 Jul 2023 10:18:05 -0600 Subject: [PATCH] Get it closer with hanging nodes. --- examples/step-86/step-86.cc | 16 +++++++--------- 1 file changed, 7 insertions(+), 9 deletions(-) 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(fe.degree + 1), - mass_matrix); - MatrixCreator::create_laplace_matrix(dof_handler, - QGauss(fe.degree + 1), - laplace_matrix); + MatrixCreator::create_mass_matrix( + dof_handler, QGauss(fe.degree + 1), mass_matrix, {}, constraints); + MatrixCreator::create_laplace_matrix( + dof_handler, QGauss(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 boundary_values_function; -- 2.39.5