Hi, Daniel:
Thank you for your help. But this seems not to be convenient. I am thinking
to run my code like the following (psuedocode)
for (i=0; i != number_of_nodes; ++i){
Point<dim> pi = get_point_from_global_index(i);
Quadrature quad(pi);
FEValues<dim> fe_values (fe, quad, |update_gradients);
... ...//do some computations with fe_values.shape_grad(i, q)
}
In this way, I can easily run over all the triangulation nodes without
worrying about worrying about recomputation at some nodes if I instead run
over each cell.
But the difficullty here that I can't find a function which can do the same
job as "get_point_from_global_index(i)", namely give me the point with
global node index i of the triangulation.
Is there any function in dealii with this functionality?
Great thanks!
Best,
David
On Sunday, May 22, 2016 at 6:11:47 PM UTC+8, Daniel Arndt wrote:
>
> David,
>
> assuming that you want to access the gradients at the support points of
> the element you are using, you can simply create a
> appropriate Quadrature object to initialize FEValues[1]. In this case the
> quadrature points correspond to the local DoF indices.
>
> Apart from that, you can use DataPostprocessor[2] for postprocessing. Its
> use is illustrated in step-29, step-32, and step-33
>
>
--
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].
For more options, visit https://groups.google.com/d/optout.
