qi-hardware/cae-tools
1
0
Fork 0
mirror of git://projects.qi-hardware.com/cae-tools.git synced 2025-01-03 20:20:15 +02:00
Code

slicer/slice.c (add): detect and report inclined facets

Browse Source
This commit is contained in:
Werner Almesberger 2015-01-08 08:24:18 -03:00
parent c984021251
commit a30024dce1
1 changed files with 11 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
slicer/slice.c
View File

@ -71,17 +71,26 @@ static void mark_z(float z)
static void add(const struct v *a, const struct v *b, const struct v *c) static void add(const struct v *a, const struct v *b, const struct v *c)
{ {
float cross;
struct z *z; struct z *z;
struct line *line; struct line *line;
/* @@@ should check that C is above AB */
if (eq(a->x, b->x) && eq(a->y, b->y)) { if (eq(a->x, b->x) && eq(a->y, b->y)) {
fprintf(stderr, "zero point\n"); fprintf(stderr, "zero point\n");
return; return;
} }
if (c->z < a->z) if (c->z < a->z)
return; return;
cross = (b->x - a->x) * (c->y - a->y) - (b->y - a->y) * (c->x - a->x);
if (!eq(cross, 0)) {
fprintf(stderr,
"inclined facet\n\t%f %f %f\n\t%f %f %f\n\t%f %f %f\n",
a->x, a->y, a->z, b->x, b->y, b->z, c->x, c->y, c->z);
exit(1);
}
z = get_z(a->z); z = get_z(a->z);
line = alloc_type(struct line); line = alloc_type(struct line);
line->ax = a->x; line->ax = a->x;