Hi,
You may want to look at :
https://www.mcs.anl.gov/petsc/petsc-current/include/petsc/private/petscfeimpl.h.html#PetscFEInterpolate_Static
It only gives the structure of what you are looking for.
In order to use it, you have to use the FE implementation of Petsc, by
looking at the examples.
I hope it helps,
Yann
Le 23/10/2018 à 09:11, Josh L a écrit :
Hi,
My FEM codes calculate the stress field on gauss point, but when I
output the stress field, I want it to on vertex. My old codes find the
gauss points that surround a vertex and calculate the average.
Is there any better way to do it with DMPlex?
The PetscSF given by DMPlexDistribute has the info how mesh points are
distributed.
Can I create a mapping from it or it can be used to map mesh point in
natural ordering to global ordering or the other way around?
Thanks,
Josh
--
___________________________
Yann JOBIC
HPC engineer
IUSTI-CNRS UMR 7343 - Polytech Marseille
Technopôle de Château Gombert
5 rue Enrico Fermi
13453 Marseille cedex 13
Tel : (33) 4 91 10 69 43
Fax : (33) 4 91 10 69 69
