Changeset 4040

Timestamp:

06/11/10 15:57:13 (15 years ago)

Author:

jschierm

Message:

ProfileProjectOntoMesh.m: Allow profile to start and/or end externally and/or cross holes in the model.

File:

• issm/trunk/src/m/utils/Mesh/ProfileProjectOntoMesh.m (modified) (3 diffs)

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

@@ -25 +25 @@
 processed_y=[];
 for i=1:numel(indices),
-        if(((indices(i)-floor(indices(i)))~=0) & ((ceil(indices(i))-indices(i))~=0))
+        if(((indices(i)-floor(indices(i)))~=0) && ((ceil(indices(i))-indices(i))~=0))
                 processed_indices=[processed_indices;floor(indices(i))];
                 processed_x=[processed_x;x0(i)];
@@ -32 +32 @@
 end
 
-%now merge profile.x,profile.y with processed_x,processed_y, at locatoins processed_indices:
+%now merge profile.x,profile.y with processed_x,processed_y, at locations processed_indices:
 newx=profile.x;
 newy=profile.y;
@@ -49 +49 @@
 node_in_element=NodeInElement(newx,newy,md.elements,md.x,md.y,md.nodeconnectivity);
 
+% eliminate nodes that don't fall in any element
+% (profile may start and/or end externally and/or cross holes in the model)
+
+ind=find(node_in_element(:,9)>0);
+newx=newx(ind,:);
+newy=newy(ind,:);
+node_in_element=node_in_element(ind,:);
+
 mesh_profile=[newx(1:end-1) newy(1:end-1) newx(2:end) newy(2:end) zeros(length(newy(2:end)),1)];