From: Peter Munch Date: Mon, 6 Jun 2022 13:08:35 +0000 (+0200) Subject: VT::interpolate_boundary_values() for pyramids and wedges X-Git-Tag: v9.4.0-rc1~17^2 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=9dcf6804a3862497544456d51927e4679c2863af;p=dealii.git VT::interpolate_boundary_values() for pyramids and wedges --- 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) {