10 |
|
|
11 |
|
#include <values.h> |
12 |
|
|
13 |
+ |
#define F_TINY 1e-5 |
14 |
+ |
#define FZERO(x) ((x) < F_TINY && (x) > -F_TINY) |
15 |
+ |
#define FEQUAL(a,b) FZERO((a) - (b)) |
16 |
+ |
|
17 |
|
#define ZERO(x) ((x) < FTINY && (x) > -FTINY) |
18 |
|
#define EQUAL(a,b) ZERO((a) - (b)) |
19 |
|
|
37 |
|
typedef long BCOORD; |
38 |
|
typedef long BDIR; |
39 |
|
typedef long TINT; |
40 |
< |
#define BITS_BCOORD (BITS(long)>>1) |
41 |
< |
#define SHIFT_MAXBCOORD (BITS_BCOORD-2) |
42 |
< |
#define MAXBCOORD (1L << SHIFT_MAXBCOORD) |
40 |
> |
|
41 |
> |
#define BITS_BCOORD (BITS(long)-2) |
42 |
> |
#define SHIFT_MAXBCOORD (BITS_BCOORD-1) |
43 |
> |
#define MAXBCOORD ((1L << BITS_BCOORD)-1) |
44 |
|
#define MAXBCOORD2 (MAXBCOORD>>1) |
45 |
< |
#define MAXBDIR MAXBCOORD |
41 |
< |
#define MAXT MAXBCOORD |
42 |
< |
#define HUGET MAXLONG |
45 |
> |
#define MAXBCOORD4 (MAXBCOORD2>>1) |
46 |
|
|
47 |
|
#define M_2_3_PI PI*2/3 |
48 |
|
|
59 |
|
#define GT_ADJACENT 32 |
60 |
|
#define GT_OUT 64 |
61 |
|
|
62 |
+ |
#define IADDV3(v,a) ((v)[0] += (a)[0],(v)[1] += (a)[1],(v)[2] += (a)[2]) |
63 |
+ |
#define ISUBV3(v,a) ((v)[0] -= (a)[0],(v)[1] -= (a)[1],(v)[2] -= (a)[2]) |
64 |
+ |
#define ISCALEV3(v,a) ((v)[0] *= (a),(v)[1] *= (a),(v)[2] *= (a)) |
65 |
+ |
#define IDIVV3(v,a) ((v)[0] /= (a),(v)[1] /= (a),(v)[2] /= (a)) |
66 |
+ |
|
67 |
+ |
|
68 |
+ |
#define ADDV3(v,a,b) ((v)[0] = (a)[0]+(b)[0],(v)[1] = (a)[1]+(b)[1],\ |
69 |
+ |
(v)[2] = (a)[2]+(b)[2]) |
70 |
+ |
#define SUBV3(v,a,b) ((v)[0] = (a)[0]-(b)[0],(v)[1] = (a)[1]-(b)[1],\ |
71 |
+ |
(v)[2] = (a)[2]-(b)[2]) |
72 |
+ |
#define SUMV3(v,a,b,s) ((v)[0] = (a)[0]+(s)*(b)[0],(v)[1]=(a)[1]+(s)*(b)[1],\ |
73 |
+ |
(v)[2] = (a)[2]+(s)*(b)[2]) |
74 |
+ |
#define SCALEV3(v,a,s) ((v)[0]=(a)[0]*(s),(v)[1]=(a)[1]*(s),(v)[2]=(a)[2]*(s)) |
75 |
|
#define ZERO_VEC3(v) (ZERO(v[0]) && ZERO(v[1]) && ZERO(v[2]) ) |
76 |
< |
#define EQUAL_VEC3(a,b) (EQUAL(a[0],b[0])&&EQUAL(a[1],b[1])&&EQUAL(a[2],b[2])) |
76 |
> |
#define EQUAL_VEC3(a,b) (FEQUAL(a[0],b[0])&&FEQUAL(a[1],b[1])&&FEQUAL(a[2],b[2])) |
77 |
> |
#define OPP_EQUAL_VEC3(a,b) (EQUAL(a[0],-b[0])&&EQUAL(a[1],-b[1])&&EQUAL(a[2],-b[2])) |
78 |
> |
#define FZERO_VEC3(v) (FZERO(v[0]) && FZERO(v[1]) && FZERO(v[2]) ) |
79 |
> |
#define FEQUAL_VEC3(a,b) (FEQUAL(a[0],b[0])&&FEQUAL(a[1],b[1])&&FEQUAL(a[2],b[2])) |
80 |
|
#define NEGATE_VEC3(v) ((v)[0] *= -1.0,(v)[1] *= -1.0,(v)[2] *= -1.0) |
81 |
|
#define COPY_VEC2(v1,v2) ((v1)[0]=(v2)[0],(v1)[1]=(v2)[1]) |
82 |
|
#define DIST(a,b) (sqrt(((a)[0]-(b)[0])*((a)[0]-(b)[0]) + \ |
86 |
|
((a)[1]-(b)[1])*((a)[1]-(b)[1]) + \ |
87 |
|
((a)[2]-(b)[2])*((a)[2]-(b)[2])) |
88 |
|
|
89 |
+ |
#define SIGN(x) ((x<0)?-1:1) |
90 |
|
#define CROSS_VEC2(v1,v2) (((v1)[0]*(v2)[1]) - ((v1)[1]*(v2)[0])) |
91 |
|
#define DOT_VEC2(v1,v2) ((v1)[0]*(v2)[0] + (v1)[1]*(v2)[1]) |
92 |
|
|
93 |
< |
#define EDGE_MIDPOINT_VEC3(a,v1,v2) ((a)[0]=((v1)[0]+(v2)[0])*0.5, \ |
93 |
> |
#define EDGE_MIDPOINT(a,v1,v2) ((a)[0]=((v1)[0]+(v2)[0])*0.5, \ |
94 |
|
(a)[1]=((v1)[1]+(v2)[1])*0.5,(a)[2] = ((v1)[2]+(v2)[2])*0.5) |
95 |
|
|
96 |
|
#define MIN_VEC3(v) ((v)[0]<(v)[1]?((v)[0]<(v)[2]?(v)[0]:v[2]): \ |
137 |
|
double tri_normal(); |
138 |
|
/* double spherical_edge_normal(FVECT v0,FVECT v1,FVECT n,char norm) */ |
139 |
|
double spherical_edge_normal(); |
140 |
+ |
double point_on_sphere(); |
141 |
|
|
142 |
|
#define point_in_stri_n(n0,n1,n2,p) \ |
143 |
< |
((DOT(n0,p)<=0.0)&&(DOT(n1,p)<=0.0)&&(DOT(n2,p)<=0.0)) |
143 |
> |
((DOT(n0,p)<=FTINY)&&(DOT(n1,p)<=FTINY)&&(DOT(n2,p)<=FTINY)) |
144 |
|
|
145 |
|
#define PT_ON_PLANE(p,peq) (DOT(FP_N(peq),p)+FP_D(peq)) |
146 |
+ |
|