From 9dcf6804a3862497544456d51927e4679c2863af Mon Sep 17 00:00:00 2001 From: Peter Munch Date: Mon, 6 Jun 2022 15:08:35 +0200 Subject: [PATCH] VT::interpolate_boundary_values() for pyramids and wedges --- include/deal.II/numerics/vector_tools_boundary.templates.h | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/include/deal.II/numerics/vector_tools_boundary.templates.h b/include/deal.II/numerics/vector_tools_boundary.templates.h index ec71f01e86..0a5be3580f 100644 --- a/include/deal.II/numerics/vector_tools_boundary.templates.h +++ b/include/deal.II/numerics/vector_tools_boundary.templates.h @@ -283,8 +283,9 @@ namespace VectorTools // dofs on this face face_dofs.resize(fe.n_dofs_per_face(face_no)); face->get_dof_indices(face_dofs, cell->active_fe_index()); - const std::vector> &dof_locations = + std::vector> dof_locations = fe_values.get_quadrature_points(); + dof_locations.resize(fe.n_dofs_per_face(face_no)); if (fe_is_system) { -- 2.39.5