Changeset 8070

Timestamp:

Nov 16, 2010, 7:35:01 AM (14 years ago)

Author:

steve

Message:

Added test of zero area triangle in check_integrity

File:

• trunk/anuga_core/source/anuga/abstract_2d_finite_volumes/neighbour_mesh.py (modified) (1 diff)

trunk/anuga_core/source/anuga/abstract_2d_finite_volumes/neighbour_mesh.py

@@ -635 +635 @@
             assert abs((area - ref)/area) < epsilon, msg
 
+            msg = 'Triangle (%f,%f), (%f,%f), (%f, %f)' % (x0,y0,x1,y1,x2,y2)
+            msg += ' is degenerate:  area == %f' % self.areas[i]
+            assert area > 0.0, msg
+
             # Check that points are arranged in counter clock-wise order
             v0 = [x1-x0, y1-y0]