Changeset 25559
- Timestamp:
- 09/11/20 09:44:41 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/modules/Trianglex/Trianglex.cpp
r25554 r25559 42 42 /*Create initial triangulation to call triangulate(). First number of points:*/ 43 43 in.numberofpoints=0; 44 for ( Object* & object : domain->objects){45 contour=(Contour<IssmPDouble>*) object;44 for (i=0;i<domain->Size();i++){ 45 contour=(Contour<IssmPDouble>*)domain->GetObjectByOffset(i); 46 46 in.numberofpoints+=contour->nods-1; 47 47 } 48 for ( Object* & object : rifts->objects){49 contour=(Contour<IssmPDouble>*) object;48 for (i=0;i<rifts->Size();i++){ 49 contour=(Contour<IssmPDouble>*)rifts->GetObjectByOffset(i); 50 50 in.numberofpoints+=contour->nods; 51 51 } … … 58 58 59 59 counter=0; 60 for ( Object* & object : domain->objects){61 contour=(Contour<IssmPDouble>*) object;60 for (i=0;i<domain->Size();i++){ 61 contour=(Contour<IssmPDouble>*)domain->GetObjectByOffset(i); 62 62 for (j=0;j<contour->nods-1;j++){ 63 63 in.pointlist[2*counter+0]=contour->x[j]; … … 66 66 } 67 67 } 68 for ( Object* & object : rifts->objects){69 contour=(Contour<IssmPDouble>*) object;68 for (i=0;i<rifts->Size();i++){ 69 contour=(Contour<IssmPDouble>*)rifts->GetObjectByOffset(i); 70 70 for (j=0;j<contour->nods;j++){ 71 71 in.pointlist[2*counter+0]=contour->x[j]; … … 85 85 /*Build segments. First figure out number of segments: holes and closed outlines have as many segments as vertices: */ 86 86 in.numberofsegments=0; 87 for ( Object* & object : domain->objects){88 contour=(Contour<IssmPDouble>*) object;87 for (i=0;i<domain->Size();i++){ 88 contour=(Contour<IssmPDouble>*)domain->GetObjectByOffset(i); 89 89 in.numberofsegments+=contour->nods-1; 90 90 } 91 for (Object* & object : rifts->objects){92 contour=(Contour<IssmPDouble>*) object;91 for(i=0;i<rifts->Size();i++){ 92 contour=(Contour<IssmPDouble>*)rifts->GetObjectByOffset(i); 93 93 /*for rifts, we have one less segment as we have vertices*/ 94 94 in.numberofsegments+=contour->nods-1; … … 99 99 counter=0; 100 100 backcounter=0; 101 for ( Object* & object : domain->objects){102 contour=(Contour<IssmPDouble>*) object;101 for (i=0;i<domain->Size();i++){ 102 contour=(Contour<IssmPDouble>*)domain->GetObjectByOffset(i); 103 103 for (j=0;j<contour->nods-2;j++){ 104 104 in.segmentlist[2*counter+0]=counter; … … 115 115 } 116 116 counter2=counter; 117 for ( Object* & object : rifts->objects){118 contour=(Contour<IssmPDouble>*) object;117 for (i=0;i<rifts->Size();i++){ 118 contour=(Contour<IssmPDouble>*)rifts->GetObjectByOffset(i); 119 119 for (j=0;j<(contour->nods-1);j++){ 120 120 in.segmentlist[2*counter2+0]=counter;
Note:
See TracChangeset
for help on using the changeset viewer.