Changeset 9728 for issm/trunk/src/m/model/mesh/findsegments.m

Timestamp:

09/09/11 13:40:05 (14 years ago)

Author:

seroussi

Message:

working on mesh class

File:

• issm/trunk/src/m/model/mesh/findsegments.m (modified) (3 diffs)

issm/trunk/src/m/model/mesh/findsegments.m

@@ -3 +3 @@
 %
 %   Optional inputs:
-%      'elementconnectivity'
+%      'mesh.elementconnectivity'
 %
 %   Usage:
@@ -12 +12 @@
 
 %Get connectivity
-elementconnectivity=getfieldvalue(options,'elementconnectivity',md.elementconnectivity);
+mesh.elementconnectivity=getfieldvalue(options,'mesh.elementconnectivity',md.mesh.elementconnectivity);
 
 %Now, build the connectivity tables for this mesh if not correclty done
-if size(md.elementconnectivity,1)~=md.mesh.numberofelements,
-        if exist(options,'elementconnectivity'),
-                error(' ''elementconnectivity'' option does not have thge right size.');
+if size(md.mesh.elementconnectivity,1)~=md.mesh.numberofelements,
+        if exist(options,'mesh.elementconnectivity'),
+                error(' ''mesh.elementconnectivity'' option does not have thge right size.');
         else
-                elementconnectivity=ElementConnectivity(md.elements,md.nodeconnectivity);
+                mesh.elementconnectivity=ElementConnectivity(md.elements,md.mesh.nodeconnectivity);
         end
 end
 
 %Recreate the segments
-elementonboundary=double(elementconnectivity(:,3)==0);
+elementonboundary=double(mesh.elementconnectivity(:,3)==0);
 pos=find(elementonboundary);
 num_segments=length(pos);
@@ -37 +37 @@
 
         %get elements connected to el1
-        els2=elementconnectivity(el1,find(elementconnectivity(el1,:)));
+        els2=mesh.elementconnectivity(el1,find(mesh.elementconnectivity(el1,:)));
 
         %el1 is connected to 2 other elements