1 |
/* Copyright (c) 1998 Silicon Graphics, Inc. */ |
2 |
|
3 |
/* SCCSid "$SunId$ SGI" */ |
4 |
|
5 |
/* |
6 |
* sm_geom.h |
7 |
*/ |
8 |
|
9 |
/* Assumes included after standard.h */ |
10 |
|
11 |
#define ZERO(x) ((x) < FTINY && (x) > -FTINY) |
12 |
#define EQUAL(a,b) ZERO((a) - (b)) |
13 |
|
14 |
#ifndef TRUE |
15 |
#define TRUE 1 |
16 |
#define FALSE 0 |
17 |
#endif |
18 |
|
19 |
#define M_2_3_PI PI*2/3 |
20 |
|
21 |
#define INVALID -1 |
22 |
|
23 |
#define GT_INVALID 0 |
24 |
#define GT_VERTEX 1 |
25 |
#define GT_EDGE 2 |
26 |
#define GT_FACE 4 |
27 |
#define GT_INTERIOR 8 |
28 |
#define GT_INTERSECT 16 |
29 |
#define GT_ADJACENT 32 |
30 |
#define GT_OUT 64 |
31 |
|
32 |
#define ZERO_VEC3(v) (ZERO(v[0]) && ZERO(v[1]) && ZERO(v[2]) ) |
33 |
#define EQUAL_VEC3(a,b) (EQUAL(a[0],b[0])&&EQUAL(a[1],b[1])&&EQUAL(a[2],b[2])) |
34 |
#define NEGATE_VEC3(v) ((v)[0] *= -1.0,(v)[1] *= -1.0,(v)[2] *= -1.0) |
35 |
#define COPY_VEC2(v1,v2) ((v1)[0]=(v2)[0],(v1)[1]=(v2)[1]) |
36 |
#define DIST(a,b) (sqrt(((a)[0]-(b)[0])*((a)[0]-(b)[0]) + \ |
37 |
((a)[1]-(b)[1])*((a)[1]-(b)[1]) + \ |
38 |
((a)[2]-(b)[2])*((a)[2]-(b)[2]))) |
39 |
#define DIST_SQ(a,b) (((a)[0]-(b)[0])*((a)[0]-(b)[0]) + \ |
40 |
((a)[1]-(b)[1])*((a)[1]-(b)[1]) + \ |
41 |
((a)[2]-(b)[2])*((a)[2]-(b)[2])) |
42 |
|
43 |
#define CROSS_VEC2(v1,v2) (((v1)[0]*(v2)[1]) - ((v1)[1]*(v2)[0])) |
44 |
#define DOT_VEC2(v1,v2) ((v1)[0]*(v2)[0] + (v1)[1]*(v2)[1]) |
45 |
|
46 |
#define EDGE_MIDPOINT_VEC3(a,v1,v2) ((a)[0]=((v1)[0]+(v2)[0])*0.5, \ |
47 |
(a)[1]=((v1)[1]+(v2)[1])*0.5,(a)[2] = ((v1)[2]+(v2)[2])*0.5) |
48 |
|
49 |
#define MIN_VEC3(v) ((v)[0]<(v)[1]?((v)[0]<(v)[2]?(v)[0]:v[2]): \ |
50 |
(v)[1]<(v)[2]?(v)[1]:(v)[2]) |
51 |
|
52 |
|
53 |
#define SUM_3VEC3(r,a,b,c) ((r)[0]=(a)[0]+(b)[0]+(c)[0], \ |
54 |
(r)[1]=(a)[1]+(b)[1]+(c)[1],(r)[2]=(a)[2]+(b)[2]+(c)[2]) |
55 |
|
56 |
#define NTH_BIT(n,i) ((n) & (1<<(i))) |
57 |
#define SET_NTH_BIT(n,i) ((n) |= (1<<(i))) |
58 |
|
59 |
|
60 |
/* int convex_angle(FVECT v0,FVECT v1,FVECT v2) */ |
61 |
/* void triangle_centroid(FVECT v0,FVECT v1,FVECT v2,FVECT c) */ |
62 |
/* void triangle_plane_equation(FVECT v0,FVECT v1,FVECT v2,FVECT n,double *nd, |
63 |
char norm) */ |
64 |
/* int vec3_equal(FVECT v1,v2) */ |
65 |
/* int point_relative_to_plane(FVECT p,FVECT n, double nd) */ |
66 |
/* int point_in_circle(FVECT p,FVECT p0,FVECT p1) */ |
67 |
/* int intersect_line_plane(FVECT r,FVECT p1,FVECT p2,float *plane) */ |
68 |
/* int point_in_cone(FVECT p,FVECT p1,FVECT p2,FVECT p3,FVECT p4) */ |
69 |
/* void point_on_sphere(FVECT ps,FVECT p,FVECT c) */ |
70 |
/* int test_point_against_spherical_tri(FVECT v0,FVECT v1,FVECT v2,FVECT p, |
71 |
FVECT n,char *nset,char *which,char sides[3]) */ |
72 |
/* int test_single_point_against_spherical_tri(FVECT v0,FVECT v1,FVECT v2, |
73 |
FVECT p,char *which )*/ |
74 |
/* int test_vertices_for_tri_inclusion(FVECT tri[3],FVECT pts[3],char *nset, |
75 |
FVECT n[3],FVECT avg,char pt_sides[3][3]); */ |
76 |
/* void set_sidedness_tests(FVECT tri[3],FVECT pts[3],char test[3], |
77 |
char sides[3][3],char nset,FVECT n[3]) |
78 |
*/ |
79 |
/* int cs_spherical_edge_edge_test(FVECT n[2][3],int i,int j,FVECT avg[2]) */ |
80 |
/* int spherical_tri_tri_intersect(FVECT a1,FVECT a2,FVECT a3, |
81 |
FVECT b1,FVECT b2,FVECT b3) */ |
82 |
|
83 |
/* void calculate_view_frustum(FVECT vp,hv,vv,double horiz,vert,near,far, |
84 |
FVECT fnear[4],FVECT ffar[4]) |
85 |
*/ |
86 |
/* double triangle_normal_Newell(FVECT v0,FVECT v1,FVECT v2,FVECT n,char n)*/ |
87 |
double tri_normal(); |
88 |
/* double spherical_edge_normal(FVECT v0,FVECT v1,FVECT n,char norm) */ |
89 |
double spherical_edge_normal(); |
90 |
|
91 |
|