Changeset 1800

Timestamp:

Sep 7, 2005, 4:29:49 PM (19 years ago)

Author:

ole

Message:

Error message in check_integrity

File:

• inundation/pyvolution/mesh.py (modified) (1 diff)

inundation/pyvolution/mesh.py

@@ -536 +536 @@
 
             for u, v in [ (v0, normal2), (v1, normal0), (v2, normal1) ]:
-                assert u[0]*v[0] + u[1]*v[1] < epsilon
+
+                #Normalise
+                from math import sqrt
+                l_u = sqrt(u[0]*u[0] + u[1]*u[1])
+                l_v = sqrt(v[0]*v[0] + v[1]*v[1])                
+
+                x = (u[0]*v[0] + u[1]*v[1])/l_u/l_v #Inner product
+                
+                msg = 'Normal vector (%f,%f) is not perpendicular to' %tuple(v)
+                msg += ' edge (%f,%f) in triangle %d.' %(tuple(u) + (i,))
+                msg += ' Inner product is %e.' %x                
+                assert x < epsilon, msg