On Sun, May 9, 2021 at 12:06 PM Barry Smith <bsm...@petsc.dev> wrote:
>
> Matt,
>
>
> if (centroid) for (d = 0; d < dim; ++d) centroid[d] /= (vsum*4);
> csum[d] /= (tdim+1)*vsum;
>
> Also
>
> for (d = 0; d < dim; d++) n[d] /= norm;
> etc
>
> 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?
>
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.
Thanks,
Matt
--
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
https://www.cse.buffalo.edu/~knepley/ <http://www.cse.buffalo.edu/~knepley/>
