[petsc-users] petsc section layout

Dharmendar Reddy dharmareddy84 at gmail.com
Thu Feb 6 16:49:07 CST 2014


Hello,
         I need to store some geometry information of a mesh in dm and
alter access it via getclosure (or a some other method you may
suggest).

For a traiangle  with vertices:  v1, v2, v3
say the circumcenter of the triangle is O
and mid point of each edges is : m1,m2,m3
The data i need to store is:
For each node a double value, which will the fraction of area of triangle:
for example: for v1 , area of quad v1,m1,O,m3

each edge has two values: edge length and distance from mid point to O

For triangle, i will have to save 9 variables, of which 6 are cell
specific and 3 are edge specific.  How should define the dof lay out
to use in DMPlexCreateSection


More information about the petsc-users mailing list