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

Timestamp:

09/09/11 14:58:09 (14 years ago)

Author:

seroussi

Message:

keep building mesh

File:

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

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

@@ -19 +19 @@
                 error(' ''mesh.elementconnectivity'' option does not have thge right size.');
         else
-                mesh.elementconnectivity=ElementConnectivity(md.elements,md.mesh.vertexconnectivity);
+                mesh.elementconnectivity=ElementConnectivity(md.mesh.elements,md.mesh.vertexconnectivity);
         end
 end
@@ -43 +43 @@
 
                 %get nodes of el1
-                nods1=md.elements(el1,:);
+                nods1=md.mesh.elements(el1,:);
 
                 %find the common vertices to the two elements connected to el1 (1 or 2)
-                flag=intersect(md.elements(els2(1),:),md.elements(els2(2),:));
+                flag=intersect(md.mesh.elements(els2(1),:),md.mesh.elements(els2(2),:));
 
                 %get the vertices on the boundary and build segment
@@ -53 +53 @@
 
                 %swap segment nodes if necessary
-                ord1=find(nods1(1)==md.elements(el1,:));
-                ord2=find(nods1(2)==md.elements(el1,:));
+                ord1=find(nods1(1)==md.mesh.elements(el1,:));
+                ord2=find(nods1(2)==md.mesh.elements(el1,:));
                 if ( (ord1==1 & ord2==2) | (ord1==2 & ord2==3) | (ord1==3 & ord2==1) ),
                         temp=segments(count,1);
@@ -66 +66 @@
         else
                 %get nodes of el1
-                nods1=md.elements(el1,:);
+                nods1=md.mesh.elements(el1,:);
 
                 %find the vertex  the el1 to not share with els2
-                flag=setdiff(nods1,md.elements(els2,:));
+                flag=setdiff(nods1,md.mesh.elements(els2,:));
 
                 for j=1:3,
@@ -78 +78 @@
 
                                 %swap segment nodes if necessary
-                                ord1=find(nods(1)==md.elements(el1,:));
-                                ord2=find(nods(2)==md.elements(el1,:));
+                                ord1=find(nods(1)==md.mesh.elements(el1,:));
+                                ord2=find(nods(2)==md.mesh.elements(el1,:));
                                 if ( (ord1==1 & ord2==2) | (ord1==2 & ord2==3) | (ord1==3 & ord2==1) ),
                                         temp=segments(count,1);