| 1 |
< |
/* Copyright (c) 1998 Silicon Graphics, Inc. */ |
| 2 |
< |
|
| 3 |
< |
/* SCCSid "$SunId$ SGI" */ |
| 4 |
< |
|
| 1 |
> |
/* RCSid: $Id$ */ |
| 2 |
|
/* |
| 3 |
|
* sm_geom.h |
| 4 |
|
*/ |
| 7 |
|
|
| 8 |
|
#include <values.h> |
| 9 |
|
|
| 10 |
+ |
#ifdef SMLFLT |
| 11 |
+ |
#define EQUALITY_EPS 1e-6 |
| 12 |
+ |
#else |
| 13 |
+ |
#define EQUALITY_EPS 1e-10 |
| 14 |
+ |
#endif |
| 15 |
+ |
|
| 16 |
|
#define F_TINY 1e-10 |
| 17 |
|
#define FZERO(x) ((x) < F_TINY && (x) > -F_TINY) |
| 18 |
|
#define FEQUAL(a,b) FZERO((a) - (b)) |
| 39 |
|
|
| 40 |
|
typedef long BCOORD; |
| 41 |
|
typedef long BDIR; |
| 39 |
– |
typedef long TINT; |
| 42 |
|
|
| 43 |
|
#define BITS_BCOORD (BITS(long)-2) |
| 44 |
|
#define SHIFT_MAXBCOORD (BITS_BCOORD-1) |
| 51 |
|
#ifndef INVALID |
| 52 |
|
#define INVALID -1 |
| 53 |
|
#endif |
| 52 |
– |
|
| 53 |
– |
#define GT_INVALID 0 |
| 54 |
– |
#define GT_VERTEX 1 |
| 55 |
– |
#define GT_EDGE 2 |
| 56 |
– |
#define GT_FACE 4 |
| 57 |
– |
#define GT_INTERIOR 8 |
| 58 |
– |
#define GT_INTERSECT 16 |
| 59 |
– |
#define GT_ADJACENT 32 |
| 60 |
– |
#define GT_OUT 64 |
| 61 |
– |
|
| 54 |
|
#define IADDV3(v,a) ((v)[0] += (a)[0],(v)[1] += (a)[1],(v)[2] += (a)[2]) |
| 55 |
|
#define ISUBV3(v,a) ((v)[0] -= (a)[0],(v)[1] -= (a)[1],(v)[2] -= (a)[2]) |
| 56 |
|
#define ISCALEV3(v,a) ((v)[0] *= (a),(v)[1] *= (a),(v)[2] *= (a)) |
| 88 |
|
#define MIN_VEC3(v) ((v)[0]<(v)[1]?((v)[0]<(v)[2]?(v)[0]:v[2]): \ |
| 89 |
|
(v)[1]<(v)[2]?(v)[1]:(v)[2]) |
| 90 |
|
#define MAX3(a,b,c) (((b)>(a))?((b) > (c))?(b):(c):((a)>(c))?(a):(c)) |
| 91 |
< |
#define MIN3(a,b,c) (((b)<(a))?((b) < (c))?(b):(c):((a)<(c))?(a):(c)) |
| 91 |
> |
#define MIN3(a,b,c) (((b)<(a))?((b) < (c))?(b):(c):((a)<(c))?(a):(c)) |
| 92 |
> |
|
| 93 |
|
#define MAX(a,b) (((b)>(a))?(b):(a)) |
| 94 |
|
#define MIN(a,b) (((b)<(a))?(b):(a)) |
| 95 |
|
|
| 99 |
|
#define NTH_BIT(n,i) ((n) & (1<<(i))) |
| 100 |
|
#define SET_NTH_BIT(n,i) ((n) |= (1<<(i))) |
| 101 |
|
|
| 102 |
+ |
#define PT_ON_PLANE(p,peq) (DOT(FP_N(peq),p)+FP_D(peq)) |
| 103 |
|
|
| 104 |
< |
/* int convex_angle(FVECT v0,FVECT v1,FVECT v2) */ |
| 105 |
< |
/* void triangle_centroid(FVECT v0,FVECT v1,FVECT v2,FVECT c) */ |
| 106 |
< |
/* void triangle_plane_equation(FVECT v0,FVECT v1,FVECT v2,FVECT n,double *nd, |
| 107 |
< |
char norm) */ |
| 108 |
< |
/* int vec3_equal(FVECT v1,v2) */ |
| 109 |
< |
/* int point_relative_to_plane(FVECT p,FVECT n, double nd) */ |
| 110 |
< |
/* int point_in_circle(FVECT p,FVECT p0,FVECT p1) */ |
| 111 |
< |
/* int intersect_line_plane(FVECT r,FVECT p1,FVECT p2,float *plane) */ |
| 112 |
< |
/* int point_in_cone(FVECT p,FVECT p1,FVECT p2,FVECT p3,FVECT p4) */ |
| 113 |
< |
/* void point_on_sphere(FVECT ps,FVECT p,FVECT c) */ |
| 114 |
< |
/* int test_point_against_spherical_tri(FVECT v0,FVECT v1,FVECT v2,FVECT p, |
| 115 |
< |
FVECT n,char *nset,char *which,char sides[3]) */ |
| 116 |
< |
/* int test_single_point_against_spherical_tri(FVECT v0,FVECT v1,FVECT v2, |
| 123 |
< |
FVECT p,char *which )*/ |
| 124 |
< |
/* int test_vertices_for_tri_inclusion(FVECT tri[3],FVECT pts[3],char *nset, |
| 125 |
< |
FVECT n[3],FVECT avg,char pt_sides[3][3]); */ |
| 126 |
< |
/* void set_sidedness_tests(FVECT tri[3],FVECT pts[3],char test[3], |
| 127 |
< |
char sides[3][3],char nset,FVECT n[3]) |
| 128 |
< |
*/ |
| 129 |
< |
/* int cs_spherical_edge_edge_test(FVECT n[2][3],int i,int j,FVECT avg[2]) */ |
| 130 |
< |
/* int spherical_tri_tri_intersect(FVECT a1,FVECT a2,FVECT a3, |
| 131 |
< |
FVECT b1,FVECT b2,FVECT b3) */ |
| 132 |
< |
|
| 133 |
< |
/* void calculate_view_frustum(FVECT vp,hv,vv,double horiz,vert,near,far, |
| 134 |
< |
FVECT fnear[4],FVECT ffar[4]) |
| 135 |
< |
*/ |
| 136 |
< |
/* double triangle_normal_Newell(FVECT v0,FVECT v1,FVECT v2,FVECT n,char n)*/ |
| 104 |
> |
/* FUNCTIONS: |
| 105 |
> |
int point_in_cone(FVECT p,a,b,c) |
| 106 |
> |
void triangle_centroid(FVECT v0,v1,v2,c) |
| 107 |
> |
double tri_normal(FVECT v0,v1,v2,n,int norm) |
| 108 |
> |
void tri_plane_equation(FVECT v0,v1,v2,FPEQ *peqptr,int norm) |
| 109 |
> |
int intersect_ray_plane(FVECT orig,dir,FPEQ peq,double *pd,FVECT r) |
| 110 |
> |
double point_on_sphere(FVECT ps,p,c) |
| 111 |
> |
int point_in_stri(FVECT v0,v1,v2,p) |
| 112 |
> |
int ray_intersect_tri(FVECT orig,dir,v0,v1,v2,pt) |
| 113 |
> |
void calculate_view_frustum(FVECT vp,hv,vv,double horiz,vert,near,far, |
| 114 |
> |
FVECT fnear[4],ffar[4]) |
| 115 |
> |
void bary2d(double x1,y1,x2,y2,x3,y3,px,py,coord) |
| 116 |
> |
*/ |
| 117 |
|
double tri_normal(); |
| 138 |
– |
/* double spherical_edge_normal(FVECT v0,FVECT v1,FVECT n,char norm) */ |
| 139 |
– |
double spherical_edge_normal(); |
| 118 |
|
double point_on_sphere(); |
| 119 |
|
|
| 142 |
– |
#define point_in_stri_n(n0,n1,n2,p) \ |
| 143 |
– |
((DOT(n0,p)<=FTINY)&&(DOT(n1,p)<=FTINY)&&(DOT(n2,p)<=FTINY)) |
| 120 |
|
|
| 145 |
– |
#define PT_ON_PLANE(p,peq) (DOT(FP_N(peq),p)+FP_D(peq)) |
| 121 |
|
|
