<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Wed, Nov 20, 2013 at 7:22 PM, Geoffrey Irving <span dir="ltr"><<a href="mailto:irving@naml.us" target="_blank">irving@naml.us</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Was there any particular reason for making DMPlexMarkBoundaryFaces a<br>
quadratic time algorithm?  I realize it's hard to write a subquadratic<br>
time algorithm on top of DMLabelSetValue; maybe PETSc needs some basic<br>
integer hash tables?<br></blockquote><div><br></div><div>Why is it quadratic time? I just looked again, and it seems to be linear time to me.</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">

Note: boundaries are not necessarily sublinear in the case of strongly<br>
adapted meshes.<br>
<span class="HOEnZb"><font color="#888888"><br>
Geoffrey<br>
</font></span></blockquote></div><br><br clear="all"><div><br></div>-- <br>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>