Changeset 1191

Timestamp:

Apr 4, 2005, 11:05:16 AM (20 years ago)

Author:

prow

Message:

adding function 'set_to_incribed_circle' as another way to use other timestep.

File:

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

inundation/ga/storm_surge/pyvolution/mesh.py

@@ -151 +151 @@
                %(self.coordinates.shape[0], len(self), len(self.boundary))
 
+
+    def set_to_incribed_circle(self):
+        #FIXME phase out eventually
+        from math import sqrt
+        N = self.number_of_elements
+        V = self.vertex_coordinates
+        for i in range(N):
+            x0 = V[i, 0]; y0 = V[i, 1]
+            x1 = V[i, 2]; y1 = V[i, 3]
+            x2 = V[i, 4]; y2 = V[i, 5]
+            a = sqrt((x0-x1)**2+(y0-y1)**2)
+            b = sqrt((x1-x2)**2+(y1-y2)**2)
+            c = sqrt((x2-x0)**2+(y2-y0)**2)
+
+            self.radii[i]=self.areas[i]/(2*(a+b+c)) 
 
     def build_neighbour_structure(self):