[petsc-users] DMPlex, output value on gauss point to vertex

[Yann Jobic]

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