<div dir="ltr">Hi all,<br><br>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?<br>
<br>Thanks,<br>Justin<br></div>