Changeset 5451

Timestamp:

08/20/10 12:22:20 (15 years ago)

Author:

Mathieu Morlighem

Message:

Final fix to edges extraction

File:

• issm/trunk/src/m/classes/public/modelextract.m (modified) (1 diff)

issm/trunk/src/m/classes/public/modelextract.m

@@ -160 +160 @@
                 pos=find(md2.edges(:,4)==0);
                 md2.edges(pos,4)=NaN;
-                %Invert NaN of the third column with last column
+                %Invert NaN of the third column with last column (Also invert first two columns!!)
                 pos=find(isnan(md2.edges(:,3)));
                 md2.edges(pos,3)=md2.edges(pos,4);
                 md2.edges(pos,4)=NaN;
+                values=md2.edges(pos,2);
+                md2.edges(pos,2)=md2.edges(pos,1);
+                md2.edges(pos,1)=values;
                 %Finally remove edges that do not belong to any element
                 pos=find(isnan(md2.edges(:,3)) & isnan(md2.edges(:,4)));