From: Peter Munch <peterrmuench@gmail.com>
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<Point<spacedim>> &dof_locations =
+                      std::vector<Point<spacedim>> dof_locations =
                         fe_values.get_quadrature_points();
+                      dof_locations.resize(fe.n_dofs_per_face(face_no));
 
                       if (fe_is_system)
                         {