69 |
|
f->area = normalize(f->norm); |
70 |
|
if (f->area == 0.0) { |
71 |
|
objerror(o, WARNING, "zero area"); /* used to be fatal */ |
72 |
< |
f->const = 0.0; |
72 |
> |
f->offset = 0.0; |
73 |
|
f->ax = 0; |
74 |
|
return(f); |
75 |
|
} |
76 |
|
f->area *= 0.5; |
77 |
< |
/* compute constant */ |
77 |
> |
/* compute offset */ |
78 |
|
badvert = 0; |
79 |
< |
f->const = DOT(f->norm, VERTEX(f,0)); |
79 |
> |
f->offset = DOT(f->norm, VERTEX(f,0)); |
80 |
|
for (i = 1; i < f->nv; i++) { |
81 |
|
d1 = DOT(f->norm, VERTEX(f,i)); |
82 |
< |
badvert += fabs(d1 - f->const/i) > VERTEPS; |
83 |
< |
f->const += d1; |
82 |
> |
badvert += fabs(d1 - f->offset/i) > VERTEPS; |
83 |
> |
f->offset += d1; |
84 |
|
} |
85 |
< |
f->const /= f->nv; |
85 |
> |
f->offset /= f->nv; |
86 |
|
if (badvert) |
87 |
|
objerror(o, WARNING, "non-planar vertex"); |
88 |
|
/* find axis */ |
90 |
|
if (fabs(f->norm[2]) > fabs(f->norm[f->ax])) |
91 |
|
f->ax = 2; |
92 |
|
|
93 |
< |
(FACE *)o->os = f; /* save face */ |
93 |
> |
o->os = (char *)f; /* save face */ |
94 |
|
return(f); |
95 |
|
} |
96 |
|
|