Division in file Mesh/BDS.cpp fails
I use gmsh
for an application on the hawk
HPC cluster in Stuttgart. gmsh
is connected to my application using the API. Sometimes I end up with a problem in the code Mesh/BDS.cpp:
double fact = sqrt((dx * dx + dy * dy + dz * dz) / (du * du + dv * dv));
In some cases within the kernel
array exactly point p
or a point very close to p
is included. In those cases the division evaluates to nan
and the meshing procedure fails.
A possible dirty hack is to insert
if ( isnan( fact ) ) fact = 0.0;
to check fact
.