[petsc-users] Advice on creating vectors defined on lower dimensional manifolds of a DMPlex

[Aldo Bonfiglioli]

Dear all,

what is the best approach for defining vectors that "sit" on the 
(vertices and/or faces) of a given stratum of the "Face Sets" of a DMPlex?

> DM Object: 3D plex 1 MPI process
>  type: plex
> 3D plex in 3 dimensions:
>  Number of 0-cells per rank: 9261
>  Number of 1-cells per rank: 59660
>  Number of 2-cells per rank: 98400
>  Number of 3-cells per rank: 48000
> Labels:
>  marker: 1 strata with value/size (1 (14402))
>  celltype: 4 strata with value/size (0 (9261), 1 (59660), 3 (98400), 6 
> (48000))
>  depth: 4 strata with value/size (0 (9261), 1 (59660), 2 (98400), 3 
> (48000))
>  Face Sets: 6 strata with value/size (1 (800), 2 (800), 3 (800), 4 
> (800), 5 (800), 6 (800))
>
These vectors are going to be used (for example) to store stresses and 
heat flux on solid surfaces.

To be more specific: suppose stratum 3 of the "Face Sets" is a solid wall.

I want to create a vector that that stores quantities computed on the 
(800) faces of that wall OR the vertices of that wall.

Thanks,

Aldo

-- 
Dr. Aldo Bonfiglioli
Associate professor of Fluid Mechanics
Dipartimento di Ingegneria
Universita' della Basilicata
V.le dell'Ateneo Lucano, 10 85100 Potenza ITALY
tel:+39.0971.205203 fax:+39.0971.205215
web:https://urldefense.us/v3/__http://docenti.unibas.it/site/home/docente.html?m=002423__;!!G_uCfscf7eWS!dx5g28NqJW34oxLLKP1Fjtp65c0KkvUjelPzjza0lBJtf6uu5ROFqpa2GTX5Cle8L7S_YjHssSDqe6szXd2PEYvYVHHq5mtW8EU$ 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20251006/ce4e612d/attachment.html>