Amazingly, I have discovered a potentially serious bug in the routines
which perform edge collapse in the BDS_Mesh.

Owing to this, it is possible for a point to collapse onto a linear
boundary creating slivers (triangles with collinear points).

The problem is that the sides of the triangle, a and b, are not
recalculated after the point is moved.  Thus, area_final and area_init
are always equal!  So the test for slivers is defeated.