[petsc-dev] should plexgeometry.c check for divide by zero?
knepley at gmail.com
Sun May 9 12:12:08 CDT 2021
On Sun, May 9, 2021 at 12:06 PM Barry Smith <bsmith at petsc.dev> wrote:
> if (centroid) for (d = 0; d < dim; ++d) centroid[d] /= (vsum*4);
> csum[d] /= (tdim+1)*vsum;
> for (d = 0; d < dim; d++) n[d] /= norm;
> 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
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the petsc-dev