<div dir="ltr">I think I sent my previous question (on Dec 28th) to the wrong place <br>(<a href="mailto:petsc-users-request@mcs.anl.gov">petsc-users-request@mcs.anl.gov</a>).<br><div><br>To repeat,<br><br>I am having bit of a difficulty in understanding the introduction of</div><div>constraints in DMPlex. From a quick study of the User Manual I gather <br>that it is easiest done using DMPlexSetAnchors ? The description of this <br>routine says that there is an anchorIS that specifies the anchor points (rows in the </div><div>matrix). This is okay and easily understood.</div><div><br></div><div>There is also an anchorSection which is described as a map from constraint points </div><div>(columns ?)  to the anchor points listed in the anchorIS. Should this not be a map between <br></div><div>solution indices (i.e. indices appearing in the vectors and matrices) ?<br></div><div><br></div><div>For example I am completely unable to set up a simple constraint matrix for the following (say):<br></div><div><br>Point 1, Field A, B<br></div><div>Point 2-10 Field A<br></div><div>At point 1, Field B depends on Field A at points 1-10 <br><br></div><div>When I set it up it appears to create a matrix where field A depends on field A values at points 1-10.<br></div><div><br></div><div>How does the mapping work in this case ? Will the DMPlexSetAnchors() routine work <br>for this simple scenario ? <br> </div><div>If not, is the only recourse to create the constraint matrix oneself </div><div>using DMSetDefaultConstraints ?<br> <br></div><div>Also documentation for DMSetDefaultConstraints is incomplete.<br></div><div>The function accepts three arguments (dm, section and Mat) but<br></div><div>what the section is is not described at all. <br></div><div><br></div><div>I don't know if my question makes any sense. If it does not then it is</div><div>only a reflection of my utter confusion regarding the routine DMPlexSetAnchors :-(</div><div><br></div><div>Regards,<br></div>Rochan</div>