<div dir="ltr"><div dir="ltr">On Sun, May 9, 2021 at 12:06 PM Barry Smith <<a href="mailto:bsmith@petsc.dev">bsmith@petsc.dev</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><br>
  Matt,<br>
<br>
<br>
  if (centroid) for (d = 0; d < dim; ++d) centroid[d] /= (vsum*4);<br>
   csum[d] /= (tdim+1)*vsum;<br>
<br>
  Also <br>
<br>
  for (d = 0; d < dim; d++) n[d] /= norm;<br>
  etc<br>
<br>
  Would adding (some/all) tests detect errors in the geometry information earlier then later strange behavior that happens because of the infinities in the arrays?<br></blockquote></div><div><br></div><div>We could perhaps benefit from earlier checks. Here however all points would have to be the same in order for vsum == 0. which is a generacy I have never seen.</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>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><br></div><div><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>