I am debugging a program using the function
*'get_generalized_support_points()' *( where**has_support_points()=0,
while has_generalized_support_points()=1*)*. My FE system is defined as
*'FESystem<3> fe(FE_Nedelec<3>(0), 2);'*, therefore, each active cell
has 12*2 dofs. So I would also expect 'get_generalized_support_points()' can
return the support points with a vector of size 24 (of course the value will
repeat once). However, it only has 12 valid Point<3> values, the other 12 are
zero or some crazy number.
My question is, is this reasonable, or there is sth wrong with my understand
on this.
Longying & Jochen:
the latter. You assume that the array you are reading from has 12*2 entries,
but the array really only has 12 entries and as a consequence when you output
24 elements, the latter half is accessing invalid memory.
Generalized support points are the ones at which you need to know the values
of a function to compute some kind of interpolant or projection. It is enough
to know a function (which would have 2 components) at 12 points to determine
the 24 coefficients of the interpolant because in your specific case, the two
elements that describe the two components have support points at the same
location.
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: bange...@colostate.edu
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 dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/edfc2483-14b1-4aad-efdc-baaa27abfcf8%40colostate.edu.
