3DPrinting/PrusaSlicer-NonPlainar
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Bugfix: overflow causing full object facets reversal when STL file has large coordinates. #1666

Browse Source
This commit is contained in:
Alessandro Ranellucci 2014-01-02 18:50:54 +01:00
parent 6afc01c3b9
commit 713bdd8055
1 changed files with 10 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

17
xs/src/admesh/util.c
View File

@ -376,19 +376,22 @@ void stl_calculate_volume(stl_file *stl)
static float get_area(stl_facet *facet) static float get_area(stl_facet *facet)
{ {
float cross[3][3]; double cross[3][3];
float sum[3]; float sum[3];
float n[3]; float n[3];
float area; float area;
int i; int i;
// cast to double before calculating cross product because large coordinates
// can result in overflowing product
// (bad area is responsible for bad volume and bad facets reversal)
for(i = 0; i < 3; i++){ for(i = 0; i < 3; i++){
cross[i][0]=((facet->vertex[i].y * facet->vertex[(i + 1) % 3].z) - cross[i][0]=(((double)facet->vertex[i].y * (double)facet->vertex[(i + 1) % 3].z) -
(facet->vertex[i].z * facet->vertex[(i + 1) % 3].y)); ((double)facet->vertex[i].z * (double)facet->vertex[(i + 1) % 3].y));
cross[i][1]=((facet->vertex[i].z * facet->vertex[(i + 1) % 3].x) - cross[i][1]=(((double)facet->vertex[i].z * (double)facet->vertex[(i + 1) % 3].x) -
(facet->vertex[i].x * facet->vertex[(i + 1) % 3].z)); ((double)facet->vertex[i].x * (double)facet->vertex[(i + 1) % 3].z));
cross[i][2]=((facet->vertex[i].x * facet->vertex[(i + 1) % 3].y) - cross[i][2]=(((double)facet->vertex[i].x * (double)facet->vertex[(i + 1) % 3].y) -
(facet->vertex[i].y * facet->vertex[(i + 1) % 3].x)); ((double)facet->vertex[i].y * (double)facet->vertex[(i + 1) % 3].x));
} }
sum[0] = cross[0][0] + cross[1][0] + cross[2][0]; sum[0] = cross[0][0] + cross[1][0] + cross[2][0];