Changeset 5391 for issm/trunk

Timestamp:

08/19/10 08:20:35 (15 years ago)

Author:

Eric.Larour

Message:

Simplified TriaSearch

File:

• issm/trunk/src/c/modules/TriaSearchx/TriaSearchx.cpp (modified) (2 diffs)

issm/trunk/src/c/modules/TriaSearchx/TriaSearchx.cpp

@@ -33 +33 @@
 
         // read background mesh 
-        if (verbose) printf("Reading mesh\n");
         Mesh Th(index,x,y,nods,nel); 
         Th.CreateSingleVertexToTriangleConnectivity();
@@ -48 +47 @@
 
                 // internal point 
-                if (tb.det>0){ 
-                        //Area coordinate
-                        areacoord[0]= (double) dete[0]/ tb.det;
-                        areacoord[1]= (double) dete[1] / tb.det;
-                        areacoord[2]= (double) dete[2] / tb.det;
-                        //3 vertices of the triangle
-                        i0=Th.GetId(tb[0]);
-                        i1=Th.GetId(tb[1]);
-                        i2=Th.GetId(tb[2]);
-                        //triangle number
-                        tria[i]=(double)Th.GetId(tb);
-                }
+                if (Th.GetId(tb)<nel)tria[i]=(double)Th.GetId(tb);
                 //external point
                 else tria[i]=NAN;
         }
 
+
         /*Assign output pointers:*/
         *ptria=tria;