Changeset 19949

Timestamp:

01/13/16 10:15:44 (9 years ago)

Author:

Eric.Larour

Message:

CHG: allow for optional removal of orphans

File:

• issm/trunk-jpl/src/m/mesh/triangle.m (modified) (1 diff)

issm/trunk-jpl/src/m/mesh/triangle.m

@@ -41 +41 @@
 
 %check that all the created nodes belong to at least one element
-orphan=find(~ismember([1:length(x)],sort(unique(elements(:)))));
-for i=1:length(orphan),
-        disp('WARNING: removing orphans');
-        %get rid of the orphan node i
-        %update x and y
-        x=[x(1:orphan(i)-(i-1)-1); x(orphan(i)-(i-1)+1:end)];
-        y=[y(1:orphan(i)-(i-1)-1); y(orphan(i)-(i-1)+1:end)];
-        %update elements
-        pos=find(elements>orphan(i)-(i-1));
-        elements(pos)=elements(pos)-1;
-        %update segments
-        pos1=find(segments(:,1)>orphan(i)-(i-1));
-        pos2=find(segments(:,2)>orphan(i)-(i-1));
-        segments(pos1,1)=segments(pos1,1)-1;
-        segments(pos2,2)=segments(pos2,2)-1;
+removeorphans=1;
+if removeorphans,
+        orphan=find(~ismember([1:length(x)],sort(unique(elements(:)))));
+        for i=1:length(orphan),
+                disp('WARNING: removing orphans');
+                %get rid of the orphan node i
+                %update x and y
+                x=[x(1:orphan(i)-(i-1)-1); x(orphan(i)-(i-1)+1:end)];
+                y=[y(1:orphan(i)-(i-1)-1); y(orphan(i)-(i-1)+1:end)];
+                %update elements
+                pos=find(elements>orphan(i)-(i-1));
+                elements(pos)=elements(pos)-1;
+                %update segments
+                pos1=find(segments(:,1)>orphan(i)-(i-1));
+                pos2=find(segments(:,2)>orphan(i)-(i-1));
+                segments(pos1,1)=segments(pos1,1)-1;
+                segments(pos2,2)=segments(pos2,2)-1;
+        end
 end