From: Wolfgang Bangerth 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(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;