<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Thu, Mar 31, 2016 at 10:56 AM, Afanasiev  Michael <span dir="ltr"><<a href="mailto:michael.afanasiev@erdw.ethz.ch" target="_blank">michael.afanasiev@erdw.ethz.ch</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi,<br>
<br>
I’d like to define vectors on a subset of points defined by a DM (created by DMPlex). For example, if I’m modelling wave propagation in a model with both solid and fluid regions, I’d like to define ux, uy, uz vectors in the solid part, and p in the acoustic part. Basically, just a varying number of dofs per integration point. Looking through the documentation I’ve found a pretty comprehensive discussion of setting up the PetscSection object. Is an inroad here to define a custom PetscSection on a subset of elements, faces, edges, vertices, etc….? Or is there a better way to approach this problem? Or is this functionality supported at all? We’re pretty familiar with PLEX by now, just looking for ways to define vectors on subsets of the whole domain.<br></blockquote><div><br></div><div>You can certainly have 0 dofs on any mesh point, so defining a vector on part of the domain is easy.</div><div>This is enough for your case I think.</div><div><br></div><div>If you wanted to do this for many very small pieces, then in order for PetscSection to be efficient, you</div><div>would have to consecutively number the points on those pieces, which seems too hard. For this case,</div><div>I would make separate subDMs.</div><div><br></div><div>  Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Thanks,<br>
Mike.<br>
--<br>
Michael Afanasiev<br>
Ph.D. Candidate<br>
Computational Seismology<br>
Institut für Geophysik<br>
ETH Zürich<br>
<br>
Sonneggstrasse 5, NO H 39.2<br>
CH 8092 Zürich<br>
<a href="mailto:michael.afanasiev@erdw.ethz.ch">michael.afanasiev@erdw.ethz.ch</a><br>
<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature">What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div>
</div></div>