8 |
|
|
9 |
|
/* Assumes included after standard.h */ |
10 |
|
|
11 |
+ |
#include <values.h> |
12 |
+ |
|
13 |
+ |
#define F_TINY 1e-3 |
14 |
+ |
|
15 |
+ |
#define FZERO(x) ((x) < F_TINY && (x) > -F_TINY) |
16 |
+ |
#define FEQUAL(a,b) FZERO((a) - (b)) |
17 |
+ |
|
18 |
|
#define ZERO(x) ((x) < FTINY && (x) > -FTINY) |
19 |
|
#define EQUAL(a,b) ZERO((a) - (b)) |
20 |
|
|
23 |
|
#define FALSE 0 |
24 |
|
#endif |
25 |
|
|
26 |
+ |
typedef struct _FPEQ { |
27 |
+ |
FVECT n; |
28 |
+ |
double d; |
29 |
+ |
char x,y,z; |
30 |
+ |
}FPEQ; |
31 |
+ |
|
32 |
+ |
#define FP_N(f) ((f).n) |
33 |
+ |
#define FP_D(f) ((f).d) |
34 |
+ |
#define FP_X(f) ((f).x) |
35 |
+ |
#define FP_Y(f) ((f).y) |
36 |
+ |
#define FP_Z(f) ((f).z) |
37 |
+ |
|
38 |
+ |
typedef long BCOORD; |
39 |
+ |
typedef long BDIR; |
40 |
+ |
typedef long TINT; |
41 |
+ |
#define BITS_BCOORD (BITS(long)>>1) |
42 |
+ |
#define SHIFT_MAXBCOORD (BITS_BCOORD-2) |
43 |
+ |
#define MAXBCOORD (1L << SHIFT_MAXBCOORD) |
44 |
+ |
#define MAXBCOORD2 (MAXBCOORD>>1) |
45 |
+ |
#define MAXBDIR MAXBCOORD |
46 |
+ |
#define MAXT MAXBCOORD |
47 |
+ |
#define HUGET MAXLONG |
48 |
+ |
|
49 |
|
#define M_2_3_PI PI*2/3 |
50 |
|
|
51 |
+ |
#ifndef INVALID |
52 |
+ |
#define INVALID -1 |
53 |
+ |
#endif |
54 |
+ |
|
55 |
|
#define GT_INVALID 0 |
56 |
|
#define GT_VERTEX 1 |
57 |
|
#define GT_EDGE 2 |
63 |
|
|
64 |
|
#define ZERO_VEC3(v) (ZERO(v[0]) && ZERO(v[1]) && ZERO(v[2]) ) |
65 |
|
#define EQUAL_VEC3(a,b) (EQUAL(a[0],b[0])&&EQUAL(a[1],b[1])&&EQUAL(a[2],b[2])) |
66 |
+ |
#define FZERO_VEC3(v) (FZERO(v[0]) && FZERO(v[1]) && FZERO(v[2]) ) |
67 |
+ |
#define FEQUAL_VEC3(a,b) (FEQUAL(a[0],b[0])&&FEQUAL(a[1],b[1])&&FEQUAL(a[2],b[2])) |
68 |
|
#define NEGATE_VEC3(v) ((v)[0] *= -1.0,(v)[1] *= -1.0,(v)[2] *= -1.0) |
69 |
|
#define COPY_VEC2(v1,v2) ((v1)[0]=(v2)[0],(v1)[1]=(v2)[1]) |
70 |
|
#define DIST(a,b) (sqrt(((a)[0]-(b)[0])*((a)[0]-(b)[0]) + \ |
71 |
|
((a)[1]-(b)[1])*((a)[1]-(b)[1]) + \ |
72 |
|
((a)[2]-(b)[2])*((a)[2]-(b)[2]))) |
73 |
+ |
#define DIST_SQ(a,b) (((a)[0]-(b)[0])*((a)[0]-(b)[0]) + \ |
74 |
+ |
((a)[1]-(b)[1])*((a)[1]-(b)[1]) + \ |
75 |
+ |
((a)[2]-(b)[2])*((a)[2]-(b)[2])) |
76 |
|
|
77 |
|
#define CROSS_VEC2(v1,v2) (((v1)[0]*(v2)[1]) - ((v1)[1]*(v2)[0])) |
78 |
|
#define DOT_VEC2(v1,v2) ((v1)[0]*(v2)[0] + (v1)[1]*(v2)[1]) |
82 |
|
|
83 |
|
#define MIN_VEC3(v) ((v)[0]<(v)[1]?((v)[0]<(v)[2]?(v)[0]:v[2]): \ |
84 |
|
(v)[1]<(v)[2]?(v)[1]:(v)[2]) |
85 |
< |
|
85 |
> |
#define MAX3(a,b,c) (((b)>(a))?((b) > (c))?(b):(c):((a)>(c))?(a):(c)) |
86 |
> |
#define MIN3(a,b,c) (((b)<(a))?((b) < (c))?(b):(c):((a)<(c))?(a):(c)) |
87 |
> |
#define MAX(a,b) (((b)>(a))?(b):(a)) |
88 |
> |
#define MIN(a,b) (((b)<(a))?(b):(a)) |
89 |
|
|
90 |
|
#define SUM_3VEC3(r,a,b,c) ((r)[0]=(a)[0]+(b)[0]+(c)[0], \ |
91 |
|
(r)[1]=(a)[1]+(b)[1]+(c)[1],(r)[2]=(a)[2]+(b)[2]+(c)[2]) |
125 |
|
/* double spherical_edge_normal(FVECT v0,FVECT v1,FVECT n,char norm) */ |
126 |
|
double spherical_edge_normal(); |
127 |
|
|
128 |
+ |
#define point_in_stri_n(n0,n1,n2,p) \ |
129 |
+ |
((DOT(n0,p)<=F_TINY)&&(DOT(n1,p)<=F_TINY)&&(DOT(n2,p)<=F_TINY)) |
130 |
|
|
131 |
+ |
#define PT_ON_PLANE(p,peq) (DOT(FP_N(peq),p)+FP_D(peq)) |