Changeset 2857
- Timestamp:
- 01/15/10 13:02:04 (15 years ago)
- Location:
- issm/trunk/src/c/Bamgx/objects
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/c/Bamgx/objects/ListofIntersectionTriangles.cpp
r2854 r2857 56 56 // test if the point b is in the same side 57 57 if (det(v0.i,v1.i,b)>=0) { 58 //cout << " All the edge " << A << B << endl;59 58 TriangleAdjacent edge=CloseBoundaryEdge(a,t,ba,bb); 60 59 Vertex & v0 = *edge.EdgeVertex(0), & v1 = *edge.EdgeVertex(1); … … 87 86 double ba,bb; 88 87 long int verbosity=2; 89 if (verbosity>7)90 cout << " SplitEdge: All the edge " << A << B << nbegin << det(vi,vj,b)91 << " deti= " << deti << " detj=" <<detj << endl;92 88 TriangleAdjacent edge=CloseBoundaryEdge(a,t,ba,bb); 93 89 Vertex & v0 = *edge.EdgeVertex(0), & v1 = *edge.EdgeVertex(1); … … 102 98 // we revert vi,vj because vi,vj is def in Adj triangle 103 99 if ( det(vi,vj,b)>=0) { 104 if (verbosity>7)105 cout << " SplitEdge: all AB outside " << A << B << endl;106 100 t=tbegin; 107 101 Real8 ba,bb; … … 221 215 TriangleAdjacent edge=CloseBoundaryEdge(b,t,ba,bb); 222 216 NewItem(B,Metric(ba,*edge.EdgeVertex(0),bb,*edge.EdgeVertex(1))); 223 // cout << " return " << ba << " " << bb << endl;224 // ajoute le 03 frev 1997 par F. hecht225 217 return; 226 218 } … … 271 263 /*FUNCTION ListofIntersectionTriangles::Length{{{1*/ 272 264 Real8 ListofIntersectionTriangles::Length(){ 273 // cout << " n= " << Size << ":" ;274 265 if (Size<=0){ 275 266 throw ErrorException(__FUNCT__,exprintf("Size<=0")); … … 292 283 SegI->lBegin=s; 293 284 294 for (jj=0,ii=1;ii<Size;jj=ii++) 295 { 285 for (jj=0,ii=1;ii<Size;jj=ii++) { 296 286 // seg jj,ii 297 287 x=y; … … 300 290 Mx = lIntTria[ii].m; 301 291 My = lIntTria[jj].m; 302 // Real8 &sx= lIntTria[ii].sp; // previous seg303 // Real8 &sy= lIntTria[jj].sn; // next seg304 // sx = Mx(xy);305 // sy = My(xy);306 // sxy = (Mx(xy)+ My(xy))/2.0;307 292 sxy = LengthInterpole(Mx,My,xy); 308 293 s += sxy; … … 314 299 SegI->lBegin=s; 315 300 316 // cout << ii << " " << jj << x<< y <<xy << s << lIntTria[ii].m ;317 301 } 318 302 len = s; 319 303 SegI->lEnd=s; 320 304 321 // cout << " len= " << s << endl;322 305 return s; 323 306 } … … 369 352 vertices[nbv].r = C; 370 353 vertices[nbv++].m = Metric(cx,lIntTria[ii-1].m,cy,lIntTria[ii].m); 371 if((verbosity/100%10)==2)372 cout << " -- Add point " << nbv-1 << " " << vertices[nbv-1] << " " << vertices[nbv-1].m << endl;373 374 354 } 375 355 else return nbv-nbvold; -
issm/trunk/src/c/Bamgx/objects/Triangles.cpp
r2856 r2857 2219 2219 long int verbosity=2; 2220 2220 Int4 nbtold(nbt),nbvold(nbv); 2221 if (verbosity>2) 2222 cout << " -- Triangles::NewPoints "; 2223 if (verbosity>3)cout << " nbv (in) on Boundary = " << nbv <<endl; 2221 if (verbosity>2) printf(" Triangles::NewPoints\n"); 2224 2222 Int4 i,k; 2225 2223 int j; … … 2229 2227 nbtold = nbt; 2230 2228 2231 if (KeepBackVertex && (&Bh != this) && (nbv+Bh.nbv< nbvx)) 2232 { 2229 if (KeepBackVertex && (&Bh != this) && (nbv+Bh.nbv< nbvx)){ 2233 2230 // Bh.SetVertexFieldOn(); 2234 2231 for (i=0;i<Bh.nbv;i++) … … 2242 2239 Bh.ReMakeTriangleContainingTheVertex(); 2243 2240 InsertNewPoints(nbvold,NbTSwap) ; 2244 if (verbosity>2)2245 cout << " (Nb of Points from background mesh = "2246 << nbv-nbvold << " / " << nbv1-nbvold << ")" << endl;2247 2241 } 2248 2242 else … … 2272 2266 { // for each triangle t 2273 2267 // we can change first_np_or_next_t[i] 2274 // cout << " Do the triangle " << i << " Next_t=" << next_t << endl;2275 2268 if (i<0 || i>=nbt ){ 2276 2269 throw ErrorException(__FUNCT__,exprintf("i<0 || i>=nbt")); … … 2354 2347 */ 2355 2348 NbTSwap += NbSwapf ; 2356 if (verbosity>3) cout << " " ;2357 if (verbosity>2)2358 cout << " Nb Of Vertices =" << nbv << " Nb of triangles = " << nbt-NbOutT2359 << " NbSwap final = " << NbSwapf << " Nb Total Of Swap = " << NbTSwap << endl;2360 2361 2362 2349 } 2363 2350 /*}}}1*/
Note:
See TracChangeset
for help on using the changeset viewer.