Changeset 1062

Timestamp:

Mar 11, 2005, 12:14:59 PM (20 years ago)

Author:

duncan

Message:

added a hack to util to stop acos(1.0) and acos(-1.0) causing an error. confused as to why it was causing an error.

File:

• inundation/ga/storm_surge/pyvolution/util.py (modified) (1 diff)

inundation/ga/storm_surge/pyvolution/util.py

@@ -20 +20 @@
         theta = acos(v1)
     except ValueError, e:
-        print 'Angle acos(%s) failed: %s' %(str(v1), e)         
-        raise ValueError, e
-
+        print 'WARNING Angle acos(%s) failed: %s' %(str(v1), e)
+
+        #FIXME, hack to avoid acos(1.0) Value error
+        # why is it happening?
+        # is it catching something we should avoid?
+        s = 1e-6
+        if (v1+s >  1.0) and (v1-s < 1.0) :
+            theta = 0.0
+        elif (v1+s >  -1.0) and (v1-s < -1.0):
+            theta = 3.1415926535897931
+        print "WARNING, setting theta to ", theta   
     if v2 < 0:
         #Quadrant 3 or 4