Dear all, I am approximating a scalar function F of one single argument by discretizing the argument over the interval [lb;ub]. Then, I interpolate the nodal values utilizing VectorTools::interpolate. To generate the triangulation (a line), I call GridGenerator::hyber_cube(triangulation, lb, ub) and refine globally up to a desired level. Either FE_Q(1) or FE_Q(2) elements are used.
The issue is that there is one 'point' p in [lb;ub] for which the function must be Zero. In case of FE_Q(1) elements, I simply shifted the vertex which is closest to p, to p. For FE_Q(2) elements, the closest point to p is a vertex *or* the support point in the center of the element. I think it is not possible to 'shift' a support point as the triangulation is made up of vertices only, right? Best Simon -- 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/e4c1dc91-7e1c-4317-93ca-ef8985f1313dn%40googlegroups.com.
