Changeset 23883
- Timestamp:
- 04/22/19 11:34:41 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
TabularUnified issm/trunk-jpl/src/c/classes/Elements/Tria.cpp ¶
r23880 r23883 1319 1319 if(gl[0]*gl[1]>0){ //Nodes 0 and 1 are similar, so points must be found on segment 0-2 and 1-2 1320 1320 /*Coordinates of point 2: same as initial point 2*/ 1321 xyz_bis[2][0]= *(xyz_list+3*2+0);1322 xyz_bis[2][1]= *(xyz_list+3*2+1);1323 xyz_bis[2][2]= *(xyz_list+3*2+2);1321 xyz_bis[2][0]=xyz_list[3*2+0]; 1322 xyz_bis[2][1]=xyz_list[3*2+1]; 1323 xyz_bis[2][2]=xyz_list[3*2+2]; 1324 1324 1325 1325 /*Portion of the segments*/ … … 1328 1328 1329 1329 /*New point 1*/ 1330 xyz_bis[1][0]= *(xyz_list+3*2+0)+s1*(*(xyz_list+3*1+0)-*(xyz_list+3*2+0));1331 xyz_bis[1][1]= *(xyz_list+3*2+1)+s1*(*(xyz_list+3*1+1)-*(xyz_list+3*2+1));1332 xyz_bis[1][2]= *(xyz_list+3*2+2)+s1*(*(xyz_list+3*1+2)-*(xyz_list+3*2+2));1330 xyz_bis[1][0]=xyz_list[3*2+0]+s1*(xyz_list[3*1+0]-xyz_list[3*2+0]); 1331 xyz_bis[1][1]=xyz_list[3*2+1]+s1*(xyz_list[3*1+1]-xyz_list[3*2+1]); 1332 xyz_bis[1][2]=xyz_list[3*2+2]+s1*(xyz_list[3*1+2]-xyz_list[3*2+2]); 1333 1333 1334 1334 /*New point 0*/ 1335 xyz_bis[0][0]= *(xyz_list+3*2+0)+s2*(*(xyz_list+3*0+0)-*(xyz_list+3*2+0));1336 xyz_bis[0][1]= *(xyz_list+3*2+1)+s2*(*(xyz_list+3*0+1)-*(xyz_list+3*2+1));1337 xyz_bis[0][2]= *(xyz_list+3*2+2)+s2*(*(xyz_list+3*0+2)-*(xyz_list+3*2+2));1335 xyz_bis[0][0]=xyz_list[3*2+0]+s2*(xyz_list[3*0+0]-xyz_list[3*2+0]); 1336 xyz_bis[0][1]=xyz_list[3*2+1]+s2*(xyz_list[3*0+1]-xyz_list[3*2+1]); 1337 xyz_bis[0][2]=xyz_list[3*2+2]+s2*(xyz_list[3*0+2]-xyz_list[3*2+2]); 1338 1338 } 1339 1339 else if(gl[1]*gl[2]>0){ //Nodes 1 and 2 are similar, so points must be found on segment 0-1 and 0-2 1340 1340 /*Coordinates of point 0: same as initial point 2*/ 1341 xyz_bis[0][0]= *(xyz_list+3*0+0);1342 xyz_bis[0][1]= *(xyz_list+3*0+1);1343 xyz_bis[0][2]= *(xyz_list+3*0+2);1341 xyz_bis[0][0]=xyz_list[3*0+0]; 1342 xyz_bis[0][1]=xyz_list[3*0+1]; 1343 xyz_bis[0][2]=xyz_list[3*0+2]; 1344 1344 1345 1345 /*Portion of the segments*/ … … 1348 1348 1349 1349 /*New point 1*/ 1350 xyz_bis[1][0]= *(xyz_list+3*0+0)+s1*(*(xyz_list+3*1+0)-*(xyz_list+3*0+0));1351 xyz_bis[1][1]= *(xyz_list+3*0+1)+s1*(*(xyz_list+3*1+1)-*(xyz_list+3*0+1));1352 xyz_bis[1][2]= *(xyz_list+3*0+2)+s1*(*(xyz_list+3*1+2)-*(xyz_list+3*0+2));1350 xyz_bis[1][0]=xyz_list[3*0+0]+s1*(xyz_list[3*1+0]-xyz_list[3*0+0]); 1351 xyz_bis[1][1]=xyz_list[3*0+1]+s1*(xyz_list[3*1+1]-xyz_list[3*0+1]); 1352 xyz_bis[1][2]=xyz_list[3*0+2]+s1*(xyz_list[3*1+2]-xyz_list[3*0+2]); 1353 1353 1354 1354 /*New point 2*/ 1355 xyz_bis[2][0]= *(xyz_list+3*0+0)+s2*(*(xyz_list+3*2+0)-*(xyz_list+3*0+0));1356 xyz_bis[2][1]= *(xyz_list+3*0+1)+s2*(*(xyz_list+3*2+1)-*(xyz_list+3*0+1));1357 xyz_bis[2][2]= *(xyz_list+3*0+2)+s2*(*(xyz_list+3*2+2)-*(xyz_list+3*0+2));1355 xyz_bis[2][0]=xyz_list[3*0+0]+s2*(xyz_list[3*2+0]-xyz_list[3*0+0]); 1356 xyz_bis[2][1]=xyz_list[3*0+1]+s2*(xyz_list[3*2+1]-xyz_list[3*0+1]); 1357 xyz_bis[2][2]=xyz_list[3*0+2]+s2*(xyz_list[3*2+2]-xyz_list[3*0+2]); 1358 1358 } 1359 1359 else if(gl[0]*gl[2]>0){ //Nodes 0 and 2 are similar, so points must be found on segment 1-0 and 1-2 1360 1360 /*Coordinates of point 1: same as initial point 2*/ 1361 xyz_bis[1][0]= *(xyz_list+3*1+0);1362 xyz_bis[1][1]= *(xyz_list+3*1+1);1363 xyz_bis[1][2]= *(xyz_list+3*1+2);1361 xyz_bis[1][0]=xyz_list[3*1+0]; 1362 xyz_bis[1][1]=xyz_list[3*1+1]; 1363 xyz_bis[1][2]=xyz_list[3*1+2]; 1364 1364 1365 1365 /*Portion of the segments*/ … … 1368 1368 1369 1369 /*New point 0*/ 1370 xyz_bis[0][0]= *(xyz_list+3*1+0)+s1*(*(xyz_list+3*0+0)-*(xyz_list+3*1+0));1371 xyz_bis[0][1]= *(xyz_list+3*1+1)+s1*(*(xyz_list+3*0+1)-*(xyz_list+3*1+1));1372 xyz_bis[0][2]= *(xyz_list+3*1+2)+s1*(*(xyz_list+3*0+2)-*(xyz_list+3*1+2));1370 xyz_bis[0][0]=xyz_list[3*1+0]+s1*(xyz_list[3*0+0]-xyz_list[3*1+0]); 1371 xyz_bis[0][1]=xyz_list[3*1+1]+s1*(xyz_list[3*0+1]-xyz_list[3*1+1]); 1372 xyz_bis[0][2]=xyz_list[3*1+2]+s1*(xyz_list[3*0+2]-xyz_list[3*1+2]); 1373 1373 1374 1374 /*New point 2*/ 1375 xyz_bis[2][0]= *(xyz_list+3*1+0)+s2*(*(xyz_list+3*2+0)-*(xyz_list+3*1+0));1376 xyz_bis[2][1]= *(xyz_list+3*1+1)+s2*(*(xyz_list+3*2+1)-*(xyz_list+3*1+1));1377 xyz_bis[2][2]= *(xyz_list+3*1+2)+s2*(*(xyz_list+3*2+2)-*(xyz_list+3*1+2));1375 xyz_bis[2][0]=xyz_list[3*1+0]+s2*(xyz_list[3*2+0]-xyz_list[3*1+0]); 1376 xyz_bis[2][1]=xyz_list[3*1+1]+s2*(xyz_list[3*2+1]-xyz_list[3*1+1]); 1377 xyz_bis[2][2]=xyz_list[3*1+2]+s2*(xyz_list[3*2+2]-xyz_list[3*1+2]); 1378 1378 } 1379 1379 else _error_("case not possible");
Note:
See TracChangeset
for help on using the changeset viewer.