Division in file Mesh/BDS.cpp fails
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;