- Timestamp:
- 08/11/10 09:01:16 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/c/objects/Bamg/ListofIntersectionTriangles.cpp
r5130 r5143 254 254 ilast=ifirst; 255 255 double ba,bb; 256 TriangleAdjacentedge=CloseBoundaryEdge(a,t,ba,bb);256 AdjacentTriangle edge=CloseBoundaryEdge(a,t,ba,bb); 257 257 BamgVertex & v0 = *edge.EdgeVertex(0), & v1 = *edge.EdgeVertex(1); 258 258 NewItem(A,Metric(ba,v0,bb,v1)); … … 260 260 // test if the point b is in the same side 261 261 if (det(v0.i,v1.i,b)>=0) { 262 TriangleAdjacentedge=CloseBoundaryEdge(a,t,ba,bb);262 AdjacentTriangle edge=CloseBoundaryEdge(a,t,ba,bb); 263 263 BamgVertex & v0 = *edge.EdgeVertex(0), & v1 = *edge.EdgeVertex(1); 264 264 NewItem(A,Metric(ba,v0,bb,v1)); … … 284 284 else if (detj<=0) // go to j direction on gamma 285 285 ocut = NextEdge[ocut]; 286 TriangleAdjacenttadj =t->Adj(ocut);286 AdjacentTriangle tadj =t->Adj(ocut); 287 287 t = tadj; 288 288 iedge= tadj; … … 290 290 double ba,bb; 291 291 long int verbose=2; 292 TriangleAdjacentedge=CloseBoundaryEdge(a,t,ba,bb);292 AdjacentTriangle edge=CloseBoundaryEdge(a,t,ba,bb); 293 293 BamgVertex & v0 = *edge.EdgeVertex(0), & v1 = *edge.EdgeVertex(1); 294 294 NewItem(A,Metric(ba,v0,bb,v1)); … … 304 304 t=tbegin; 305 305 double ba,bb; 306 TriangleAdjacentedge=CloseBoundaryEdge(b,t,ba,bb);306 AdjacentTriangle edge=CloseBoundaryEdge(b,t,ba,bb); 307 307 NewItem(B,Metric(ba,*edge.EdgeVertex(0),bb,*edge.EdgeVertex(1))); 308 308 return; … … 412 412 ilast=NewItem(t, ba[0],ba[1],ba[2]); 413 413 414 TriangleAdjacentta =t->Adj(ocut);414 AdjacentTriangle ta =t->Adj(ocut); 415 415 t = ta; 416 416 iedge= ta; 417 417 if (t->det <= 0) { 418 418 double ba,bb; 419 TriangleAdjacentedge=CloseBoundaryEdge(b,t,ba,bb);419 AdjacentTriangle edge=CloseBoundaryEdge(b,t,ba,bb); 420 420 NewItem(B,Metric(ba,*edge.EdgeVertex(0),bb,*edge.EdgeVertex(1))); 421 421 return;
Note:
See TracChangeset
for help on using the changeset viewer.