implemented: Mesh.get_area()

leuraph · Jan 16, 2024 · db66d61 · db66d61
1 parent a4efc50
commit db66d61
Showing 1 changed file with 6 additions and 2 deletions.
diff --git a/p1afempy/mesh.py b/p1afempy/mesh.py
@@ -33,8 +33,12 @@ def get_area(self) -> np.ndarray:
         """
         returns the area of each element
         """
-        # TODO implement
-        return np.zeros(self.elements.shape[0])
+        c1 = self.coordinates[self.elements[:, 0], :]
+        d21 = self.coordinates[self.elements[:, 1], :] - c1
+        d31 = self.coordinates[self.elements[:, 2], :] - c1
+
+        # vector of element areas 4*|T|
+        return 0.5 * (d21[:, 0]*d31[:, 1] - d21[:, 1] * d31[:, 0])
 
 
 class BoundaryCondition: