[petsc-users] Identifying mesh points within DMPlex FEM

[Justin Chang]

Hi all,

In the discretization schemes for the DMPlex FEM examples (ex12 and ex62),
it lets us specify equations as pointwise functions. But what if I wanted
to assign "regions" to different points within the domain? In one of my
problems, I want the diffusivity tensor and/or forcing function to differ
for various layers/regions within the computational domain. I would imagine
calling something like DMPlexGetLabelValue() to identify a region
associated with a given mesh point but is it possible to do something like
this within the f0/f1/g3/etc functions?

Thanks,
Justin
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20140825/1b8b19e5/attachment.html>