Changeset 14440

Timestamp:

03/11/13 14:48:00 (12 years ago)

Author:

Mathieu Morlighem

Message:

BUG: queue might be empty

File:

• issm/trunk-jpl/src/m/exp/flowlines.m (modified) (4 diffs)

issm/trunk-jpl/src/m/exp/flowlines.m

@@ -38 +38 @@
 
 %check seed points
-%tria=TriaSearch(index,x,y,x0,y0);
-tria=tsearch(x,y,index,x0,y0);
+tria=TriaSearch(index,x,y,x0,y0);
+%tria=tsearch(x,y,index,x0,y0);
 pos=find(isnan(tria));
 x0(pos)=[];
@@ -70 +70 @@
         %find current triangle
         queue=find(~done);
-        %tria=TriaSearch(index,x,y,X(queue),Y(queue));
-        tria=tsearch(x,y,index,X(queue),Y(queue));
+        tria=TriaSearch(index,x,y,X(queue),Y(queue));
+        %tria=tsearch(x,y,index,X(queue),Y(queue));
 
         %check that the point is actually inside a triangle of the mesh
@@ -112 +112 @@
 
 %same process but reverse (vel=-vel) to have a vcomplete flow line
+queue=[];
 counter=1;
 X=x0; Y=y0;
@@ -120 +121 @@
         %find current triangle
         queue=find(~done);
-        %tria=TriaSearch(index,x,y,X(queue),Y(queue));
-        tria=tsearch(x,y,index,X(queue),Y(queue));
+        tria=TriaSearch(index,x,y,X(queue),Y(queue));
+        %tria=tsearch(x,y,index,X(queue),Y(queue));
 
         %check that the point is actually inside a triangle of the mesh