Changeset 5238
- Timestamp:
- 08/13/10 10:55:03 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/c/shared/TriMesh/TriMeshUtils.cpp
r3775 r5238 291 291 *we can only rely on the position (x,y) of the rift grids to create a segment:*/ 292 292 for (k=0;k<3;k++){ 293 if ((x[(int)*(index+el1*3+k)-1]==x[(int)*(segments+3*j+0)-1]) & (y[(int)*(index+el1*3+k)-1]==y[(int)*(segments+3*j+0)-1])){293 if ((x[(int)*(index+el1*3+k)-1]==x[(int)*(segments+3*j+0)-1]) && (y[(int)*(index+el1*3+k)-1]==y[(int)*(segments+3*j+0)-1])){ 294 294 *(segments+3*j+0)=*(index+el1*3+k); 295 295 break; … … 297 297 } 298 298 for (k=0;k<3;k++){ 299 if ((x[(int)*(index+el1*3+k)-1]==x[(int)*(segments+3*j+1)-1]) & (y[(int)*(index+el1*3+k)-1]==y[(int)*(segments+3*j+1)-1])){299 if ((x[(int)*(index+el1*3+k)-1]==x[(int)*(segments+3*j+1)-1]) && (y[(int)*(index+el1*3+k)-1]==y[(int)*(segments+3*j+1)-1])){ 300 300 *(segments+3*j+1)=*(index+el1*3+k); 301 301 break; … … 306 306 *(segmentmarkerlist+(nsegs+i))=*(segmentmarkerlist+j); 307 307 for (k=0;k<3;k++){ 308 if ((x[(int)*(index+el2*3+k)-1]==x[(int)*(segments+3*j+0)-1]) & (y[(int)*(index+el2*3+k)-1]==y[(int)*(segments+3*j+0)-1])){308 if ((x[(int)*(index+el2*3+k)-1]==x[(int)*(segments+3*j+0)-1]) && (y[(int)*(index+el2*3+k)-1]==y[(int)*(segments+3*j+0)-1])){ 309 309 *(segments+3*(nsegs+i)+0)=*(index+el2*3+k); 310 310 break; … … 312 312 } 313 313 for (k=0;k<3;k++){ 314 if ((x[(int)*(index+el2*3+k)-1]==x[(int)*(segments+3*j+1)-1]) & (y[(int)*(index+el2*3+k)-1]==y[(int)*(segments+3*j+1)-1])){314 if ((x[(int)*(index+el2*3+k)-1]==x[(int)*(segments+3*j+1)-1]) && (y[(int)*(index+el2*3+k)-1]==y[(int)*(segments+3*j+1)-1])){ 315 315 *(segments+3*(nsegs+i)+1)=*(index+el2*3+k); 316 316 break; … … 322 322 /*Let's update segments[j][:] using element el2 and the corresponding rift segment: */ 323 323 for (k=0;k<3;k++){ 324 if ((x[(int)*(index+el2*3+k)-1]==x[(int)*(segments+3*j+0)-1]) & (y[(int)*(index+el2*3+k)-1]==y[(int)*(segments+3*j+0)-1])){324 if ((x[(int)*(index+el2*3+k)-1]==x[(int)*(segments+3*j+0)-1]) && (y[(int)*(index+el2*3+k)-1]==y[(int)*(segments+3*j+0)-1])){ 325 325 *(segments+3*j+0)=*(index+el2*3+k); 326 326 break; … … 328 328 } 329 329 for (k=0;k<3;k++){ 330 if ((x[(int)*(index+el2*3+k)-1]==x[(int)*(segments+3*j+1)-1]) & (y[(int)*(index+el2*3+k)-1]==y[(int)*(segments+3*j+1)-1])){330 if ((x[(int)*(index+el2*3+k)-1]==x[(int)*(segments+3*j+1)-1]) && (y[(int)*(index+el2*3+k)-1]==y[(int)*(segments+3*j+1)-1])){ 331 331 *(segments+3*j+1)=*(index+el2*3+k); 332 332 break; … … 337 337 *(segmentmarkerlist+(nsegs+i))=*(segmentmarkerlist+j); 338 338 for (k=0;k<3;k++){ 339 if ((x[(int)*(index+el1*3+k)-1]==x[(int)*(segments+3*j+0)-1]) & (y[(int)*(index+el1*3+k)-1]==y[(int)*(segments+3*j+0)-1])){339 if ((x[(int)*(index+el1*3+k)-1]==x[(int)*(segments+3*j+0)-1]) && (y[(int)*(index+el1*3+k)-1]==y[(int)*(segments+3*j+0)-1])){ 340 340 *(segments+3*(nsegs+i)+0)=*(index+el1*3+k); 341 341 break; … … 343 343 } 344 344 for (k=0;k<3;k++){ 345 if ((x[(int)*(index+el1*3+k)-1]==x[(int)*(segments+3*j+1)-1]) & (y[(int)*(index+el1*3+k)-1]==y[(int)*(segments+3*j+1)-1])){345 if ((x[(int)*(index+el1*3+k)-1]==x[(int)*(segments+3*j+1)-1]) && (y[(int)*(index+el1*3+k)-1]==y[(int)*(segments+3*j+1)-1])){ 346 346 *(segments+3*(nsegs+i)+1)=*(index+el1*3+k); 347 347 break;
Note:
See TracChangeset
for help on using the changeset viewer.