Hej, thanks for your prompt answer! Yes, your interpretation of my poor description is correct. So then I just need to create a new FEValues within the cell loop containing the points where I want to evaluate the shape functions, right? How do I add this points?
Also if the points are in the physical quad, do I need to map them to the reference quad [-1,1]x[-1,1] before adding them, or I just add the points directly to the FEValues? El viernes, 17 de agosto de 2018, 18:18:41 (UTC+2), Daniel Arndt escribió: > > Juan, > > I would like to ask what would be the best/fastest way to evaluate > *fe_v.shape_value >> (i,x), **fe_v.shape_grad* *(i,x)*, in points that are not quadrature >> points but that are guaranteed to lie within the curved cell/quad. >> > If I understand the situation correctly, you have a FEValues object that > is initialized on a cell and now you want to evaluate it at an additional > point that is not included in the mapped quadrature points. > Is that the case? If so, I would either add that point to the existing > FEValues object or build a separate one that includes the quadrature point > (with respect to the unit cell) you want to evaluate the shape functions at. > It is sufficient to initialize this separate object for cells for which > you actually need it. > Maybe, I misunderstand what you want. Can you specify the issue you are > trying to solve? > > Best, > Daniel > > -- 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. For more options, visit https://groups.google.com/d/optout.
