Changeset 1795

Timestamp:

08/24/09 09:18:53 (16 years ago)

Author:

seroussi

Message:

fixed averaging in 3d

File:

• issm/trunk/src/m/utils/Mesh/GetAreas.m (modified) (2 diffs)

issm/trunk/src/m/utils/Mesh/GetAreas.m

@@ -16 +16 @@
 nels=size(index,1);
 nods=length(x);
+if nargin==4, z=varargin{1}; end
 
 %some checks
@@ -46 +47 @@
 else
         %V=area(triangle)*1/3(z1+z2+z3)
-        thickness=mean(z(index(:,4:6)))-mean(z(index(:,1:3)));
-        areas=abs(0.5*((x2-x1).*(y3-y1)-(y2-y1).*(x3-x1)))*thickness;
+        thickness=mean(z(index(:,4:6)),2)-mean(z(index(:,1:3)),2);
+        areas=abs(0.5*((x2-x1).*(y3-y1)-(y2-y1).*(x3-x1))).*thickness;
 end