Changeset 14297

Timestamp:

02/01/13 10:04:41 (12 years ago)

Author:

Mathieu Morlighem

Message:

use tsearch instead of TriaSearch

File:

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

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

@@ -38 +38 @@
 
 %check seed points
-tria=TriaSearch(index,x,y,x0,y0);
+%tria=TriaSearch(index,x,y,x0,y0);
+tria=tsearch(x,y,index,x0,y0);
 pos=find(isnan(tria));
 x0(pos)=[];
@@ -69 +70 @@
         %find current triangle
         queue=find(~done);
-        tria=TriaSearch(index,x,y,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
@@ -118 +120 @@
         %find current triangle
         queue=find(~done);
-        tria=TriaSearch(index,x,y,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