On 6/10/21 7:10 AM, Wolfgang Bangerth wrote:

Indeed,
         std::vector<Vector<double>> velocity_qpt_vals(n_face_qpts, Vector<double>(dim));          fiv.get_fe_face_values(1).get_function_values(velocity_field, velocity_qpt_vals);

accesses the values; the issue was I hadn't explicitly instantiated the Vector<double>s with length dim, which was causing a segfault.

I forgot that earlier. Could you come up with a small program that shows the segfault? We should have caught the wrong vector size via an assertion, and I'd like to add that assertion.

Best
 W.

--
------------------------------------------------------------------------
Wolfgang Bangerth          email:                 [email protected]
                           www: http://www.math.colostate.edu/~bangerth/

--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see 
https://groups.google.com/d/forum/dealii?hl=en
--- You received this message because you are subscribed to the Google Groups "deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to [email protected].
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/d80934eb-2b05-84d6-944a-d9929bce3b09%40colostate.edu.

Reply via email to