Changeset 22213
- Timestamp:
- 11/02/17 14:28:36 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/bamg/Mesh.cpp
r21931 r22213 1134 1134 1135 1135 }/*}}}*/ 1136 void 1136 void Mesh::BoundAnisotropy(double anisomax,double hminaniso) {/*{{{*/ 1137 1137 /*Original code from Frederic Hecht <hecht@ann.jussieu.fr> (BAMG v1.01, Metric.cpp/BoundAnisotropy)*/ 1138 1138 … … 2335 2335 /*Original code from Frederic Hecht <hecht@ann.jussieu.fr> (BAMG v1.01, Mesh2.cpp/FindSubDomain)*/ 2336 2336 2337 longint verbose=0;2337 int verbose=0; 2338 2338 2339 2339 if (verbose >2){ 2340 2340 if (OutSide) _printf_(" Find all external sub-domain\n"); 2341 2341 else _printf_(" Find all internal sub-domain\n"); 2342 2342 } 2343 2343 short * HeapArete = new short[nbt]; 2344 2344 Triangle ** HeapTriangle = new Triangle* [nbt]; … … 2346 2346 long k,it; 2347 2347 2348 for (int itt=0;itt<nbt;itt++)2349 triangles[itt].link=0; // par defaut pas de couleur2348 /*No color by default*/ 2349 for(int itt=0;itt<nbt;itt++) triangles[itt].link=0; 2350 2350 2351 2351 long NbSubDomTot =0; 2352 for (it=0;it<nbt;it++) {2353 if ( ! triangles[it].link ){2352 for(it=0;it<nbt;it++) { 2353 if( !triangles[it].link ){ 2354 2354 t = triangles + it; 2355 2355 NbSubDomTot++;; // new composante connexe … … 2360 2360 HeapArete[i] = 3; 2361 2361 2362 while 2363 { while ( 2362 while(i >= 0) // boucle sur la pile 2363 { while (HeapArete[i]--) // boucle sur les 3 aretes 2364 2364 { 2365 2365 int na = HeapArete[i]; … … 2445 2445 _error_("k!=nbsubdomains"); 2446 2446 } 2447 if(OutSide) 2448 { 2447 if(OutSide){ 2449 2448 // to remove all the sub domain by parity adjacents 2450 2449 // because in this case we have only the true boundary edge 2451 // so tehboundary is manifold2450 // so the boundary is manifold 2452 2451 long nbk = nbsubdomains; 2453 2452 while (nbk) … … 3084 3083 Triangle* ta=tadjj; 3085 3084 3086 //if the adjacent triangle is a boundary triangle, continu r3085 //if the adjacent triangle is a boundary triangle, continue 3087 3086 if (ta->det<0) continue; 3088 3087
Note:
See TracChangeset
for help on using the changeset viewer.