diff options
Diffstat (limited to 'apps/plugins/sdl/progs/quake/FixedPointMath.h')
-rw-r--r-- | apps/plugins/sdl/progs/quake/FixedPointMath.h | 246 |
1 files changed, 246 insertions, 0 deletions
diff --git a/apps/plugins/sdl/progs/quake/FixedPointMath.h b/apps/plugins/sdl/progs/quake/FixedPointMath.h new file mode 100644 index 0000000000..112b1a103b --- /dev/null +++ b/apps/plugins/sdl/progs/quake/FixedPointMath.h | |||
@@ -0,0 +1,246 @@ | |||
1 | #ifndef _FIXEDPOINTMATH_H | ||
2 | #define _FIXEDPOINTMATH_H | ||
3 | //Fixed point math routines (16.16) | ||
4 | //Dan East | ||
5 | //01-24-2001 | ||
6 | |||
7 | #define FPM_PI 205887L | ||
8 | #define FPM_2PI 411775L | ||
9 | #define FPM_E 178144L | ||
10 | #define FPM_ROOT2 74804L | ||
11 | #define FPM_ROOT3 113512L | ||
12 | #define FPM_GOLDEN 106039L | ||
13 | #define FPM_MAX 0x7fff0000 | ||
14 | |||
15 | |||
16 | typedef long fixedpoint_t; | ||
17 | typedef long fixedpoint8_24_t; | ||
18 | |||
19 | |||
20 | fixedpoint_t fpm_FromFloat(double f); | ||
21 | float fpm_ToFloat(fixedpoint_t fxp); | ||
22 | |||
23 | fixedpoint_t fpm_FromLong(long l); | ||
24 | long fpm_ToLong(fixedpoint_t fxp); | ||
25 | |||
26 | fixedpoint_t fpm_Add(fixedpoint_t fxp1, fixedpoint_t fxp2); | ||
27 | fixedpoint_t fpm_Add3(fixedpoint_t fxp1, fixedpoint_t fxp2, fixedpoint_t fxp3); | ||
28 | |||
29 | fixedpoint_t fpm_Sub(fixedpoint_t fxp1, fixedpoint_t fxp2); | ||
30 | |||
31 | fixedpoint_t fpm_Mul(fixedpoint_t fxp1, fixedpoint_t fxp2); | ||
32 | |||
33 | fixedpoint_t fpm_Div(fixedpoint_t fxp1, fixedpoint_t fxp2); | ||
34 | fixedpoint_t fpm_DivInt(fixedpoint_t fxp1, long l); | ||
35 | fixedpoint_t fpm_Abs(fixedpoint_t fxp); | ||
36 | |||
37 | fixedpoint_t fpm_Ceil(fixedpoint_t fxp); | ||
38 | fixedpoint_t fpm_Floor(fixedpoint_t fxp); | ||
39 | |||
40 | fixedpoint_t fpm_Sqrt(fixedpoint_t fxp); | ||
41 | fixedpoint_t fpm_Sqr(fixedpoint_t fxp); | ||
42 | |||
43 | fixedpoint_t fpm_Inv(fixedpoint_t fxp); | ||
44 | |||
45 | fixedpoint_t fpm_Sin(fixedpoint_t fxp); | ||
46 | fixedpoint_t fpm_Cos(fixedpoint_t fxp); | ||
47 | fixedpoint_t fpm_Tan(fixedpoint_t fxp); | ||
48 | fixedpoint_t fpm_ATan(fixedpoint_t fxp); | ||
49 | //These take degrees | ||
50 | fixedpoint_t fpm_SinDeg(fixedpoint_t fxp); | ||
51 | fixedpoint_t fpm_CosDeg(fixedpoint_t fxp); | ||
52 | fixedpoint_t fpm_TanDeg(fixedpoint_t fxp); | ||
53 | fixedpoint_t fpm_ATanDeg(fixedpoint_t fxp); | ||
54 | |||
55 | |||
56 | #define FPM_FROMFLOAT(f) fpm_FromFloat(f) | ||
57 | #define FPM_FROMFLOATC(f) ((long)((f) * 65536.0 )) //Constant version | ||
58 | #define FPM_TOFLOAT(fxp) fpm_ToFloat(fxp) | ||
59 | |||
60 | #define FPM_FROMLONG(l) fpm_FromLong(l) | ||
61 | #define FPM_FROMLONGC(l) ((l)<<16) //Constant version | ||
62 | #define FPM_TOLONG(l) fpm_ToLong(l) | ||
63 | |||
64 | #define FPM_ADD(f1, f2) fpm_Add(f1, f2) | ||
65 | #define FPM_ADD3(f1, f2, f3) fpm_Add3(f1, f2, f3) | ||
66 | #define FPM_INC(f1) ((f1)=FPM_ADD(f1, FPM_FROMLONG(1))) | ||
67 | |||
68 | #define FPM_SUB(f1, f2) fpm_Sub(f1, f2) | ||
69 | #define FPM_DEC(f1) ((f1)=FPM_SUB(f1, FPM_FROMLONG(1))) | ||
70 | |||
71 | #define FPM_MUL(f1, f2) fpm_Mul(f1, f2) | ||
72 | |||
73 | #define FPM_DIV(n, d) fpm_Div(n, d) | ||
74 | #define FPM_DIVINT(n, d) fpm_DivInt(n, d) | ||
75 | |||
76 | #define FPM_ABS(fxp) fpm_Abs(fxp) | ||
77 | |||
78 | #define FPM_CEIL(fxp) fpm_Ceil(fxp) | ||
79 | #define FPM_FLOOR(fxp) fpm_Floor(fxp) | ||
80 | |||
81 | #define FPM_SQRT(fxp) fpm_Sqrt(fxp) | ||
82 | #define FPM_SQR(fxp) fpm_Sqr(fxp) | ||
83 | |||
84 | #define FPM_INV(fxp) fpm_Inv(fxp) | ||
85 | |||
86 | //These take radians | ||
87 | #define FPM_SIN(r) fpm_Sin(r) | ||
88 | #define FPM_COS(r) fpm_Cos(r) | ||
89 | #define FPM_TAN(r) fpm_Tan(r) | ||
90 | #define FPM_ATAN(r) fpm_ATan(r) | ||
91 | |||
92 | //These take degrees | ||
93 | #define FPM_SIN_DEG(d) fpm_SinDeg(d) | ||
94 | #define FPM_COS_DEG(d) fpm_CosDeg(d) | ||
95 | #define FPM_TAN_DEG(d) fpm_TanDeg(d) | ||
96 | #define FPM_ATAN_DEG(d) fpm_ATanDeg(d) | ||
97 | |||
98 | fixedpoint8_24_t fpm_FromFloat(double f); | ||
99 | float fpm_ToFloat8_24(fixedpoint8_24_t fxp); | ||
100 | |||
101 | fixedpoint8_24_t fpm_FromLong8_24(long l); | ||
102 | long fpm_ToLong8_24(fixedpoint8_24_t fxp); | ||
103 | |||
104 | fixedpoint8_24_t fpm_FromFixedPoint(fixedpoint_t fxp); | ||
105 | fixedpoint_t fpm_ToFixedPoint(fixedpoint8_24_t fxp); | ||
106 | |||
107 | fixedpoint8_24_t fpm_Add8_24(fixedpoint8_24_t fxp1, fixedpoint8_24_t fxp2); | ||
108 | fixedpoint8_24_t fpm_Add38_24(fixedpoint8_24_t fxp1, fixedpoint8_24_t fxp2, fixedpoint8_24_t fxp3); | ||
109 | |||
110 | fixedpoint8_24_t fpm_Sub8_24(fixedpoint8_24_t fxp1, fixedpoint8_24_t fxp2); | ||
111 | |||
112 | fixedpoint8_24_t fpm_Mul8_24(fixedpoint8_24_t fxp1, fixedpoint8_24_t fxp2); | ||
113 | fixedpoint_t fpm_MulMixed8_24(fixedpoint8_24_t fxp1, fixedpoint_t fxp2); | ||
114 | |||
115 | fixedpoint8_24_t fpm_Div8_24(fixedpoint8_24_t fxp1, fixedpoint8_24_t fxp2); | ||
116 | fixedpoint8_24_t fpm_DivInt8_24(fixedpoint8_24_t fxp1, long l); | ||
117 | fixedpoint8_24_t fpm_DivInt64_8_24(fixedpoint8_24_t fxp1, long long l); | ||
118 | |||
119 | fixedpoint8_24_t fpm_Abs8_24(fixedpoint8_24_t fxp); | ||
120 | |||
121 | fixedpoint8_24_t fpm_Ceil8_24(fixedpoint8_24_t fxp); | ||
122 | fixedpoint8_24_t fpm_Floor8_24(fixedpoint8_24_t fxp); | ||
123 | |||
124 | fixedpoint8_24_t fpm_Sqrt8_24(fixedpoint8_24_t fxp); | ||
125 | fixedpoint8_24_t fpm_Sqr8_24(fixedpoint8_24_t fxp); | ||
126 | |||
127 | fixedpoint8_24_t fpm_Inv8_24(fixedpoint8_24_t fxp); | ||
128 | |||
129 | |||
130 | #define FPM_FROMFLOAT8_24(f) fpm_FromFloat8_24(f) | ||
131 | #define FPM_FROMFLOATC8_24(f) ((long)((f) * 16777216.0 )) //Constant version | ||
132 | #define FPM_TOFLOAT8_24(fxp) fpm_ToFloat8_24(fxp) | ||
133 | |||
134 | #define FPM_FROMLONG8_24(l) fpm_FromLong8_24(l) | ||
135 | #define FPM_FROMLONGC8_24(l) ((l)<<24) //Constant version | ||
136 | #define FPM_TOLONG8_24(l) fpm_ToLong8_24(l) | ||
137 | |||
138 | |||
139 | /* | ||
140 | extern __int64 FPM_TMPVAR_INT64; | ||
141 | |||
142 | #define FPM_FROMFLOAT(f) ((long)((f) * 65536.0 )) //+0.5 | ||
143 | #define FPM_TOFLOAT(fxp) (((float)(fxp)) / 65536.0) | ||
144 | |||
145 | #define FPM_FROMLONG(l) ((l)<<16) | ||
146 | #define FPM_TOLONG(l) ((l)<0?(-(long)((l)^0xffffffff)>>16):(((l)>>16)&0x0000ffff)) | ||
147 | |||
148 | #define FPM_ADD(f1, f2) ((f1)+(f2)) | ||
149 | #define FPM_ADD3(f1, f2, f3) ((f1)+(f2)+(f3)) | ||
150 | #define FPM_INC(f1) ((f1)=FPM_ADD(f1, FPM_FROMLONG(1))) | ||
151 | |||
152 | #define FPM_SUB(f1, f2) ((f1)-(f2)) | ||
153 | #define FPM_DEC(f1) ((f1)=FPM_SUB(f1, FPM_FROMLONG(1))) | ||
154 | |||
155 | #define FPM_MUL(f1, f2) (((f1)>>8)*((f2)>>8)) | ||
156 | //#define FPM_MUL(f1, f2) ((fixedpoint_t)((FPM_TMPVAR_INT64=(f1))*(f2))>>16) | ||
157 | //#define FPM_MUL(f1, f2) (((f1)*(f2))>>16) | ||
158 | //TODO: This needs to be done without copying to another var | ||
159 | #define FPM_DIV(n, d) ((fixedpoint_t)(((FPM_TMPVAR_INT64=(n))<<16)/d)) | ||
160 | #define FPM_DIVINT(n, d) ((fixedpoint_t)((n)/(d))) | ||
161 | //#define FPM_DIV(n, d) ((long)(((__int64)n)<<16)/(d)) | ||
162 | |||
163 | #define FPM_ABS(fxp) (abs(fxp)) | ||
164 | |||
165 | //TODO: could be more effecient | ||
166 | #define FPM_CEIL(fxp) ((fxp)&0x0000ffff?((fxp)<=0?((fxp)&0xffff0000):(((fxp)&0xffff0000)+FPM_FROMLONG(1))):(fxp)) | ||
167 | #define FPM_FLOOR(fxp) ((fxp)&0x0000ffff?((fxp)<0?(((fxp)&0xffff0000)-FPM_FROMLONG(1)):((fxp)&0xffff0000)):(fxp)) | ||
168 | |||
169 | //TODO: Implement sqrt mathematically instead of converting to float and back | ||
170 | #define FPM_SQRT(fxp) (FPM_FROMFLOAT(sqrt(FPM_TOFLOAT(fxp)))) | ||
171 | #define FPM_SQR(fxp) (FPM_MUL(fxp,fxp)>>16) | ||
172 | |||
173 | #define FPM_INV(fxp) (FPM_DIV(0x10000, fxp)) | ||
174 | //TODO: Calc trig functions (or lookup) instead of converting to float and back | ||
175 | //These take radians | ||
176 | #define FPM_SIN(r) (FPM_FROMFLOAT(sin(FPM_TOFLOAT(r)))) | ||
177 | #define FPM_COS(r) (FPM_FROMFLOAT(cos(FPM_TOFLOAT(r)))) | ||
178 | #define FPM_TAN(r) (FPM_FROMFLOAT(tan(FPM_TOFLOAT(r)))) | ||
179 | #define FPM_ATAN(r) (FPM_FROMFLOAT(atan(FPM_TOFLOAT(r)))) | ||
180 | |||
181 | //These take degrees | ||
182 | #define FPM_SIN_DEG(d) (FPM_SIN(FPM_DIV(FPM_MUL(d,FPM_PI),0xB40000))) //0xB40000 = 180.0 | ||
183 | #define FPM_COS_DEG(d) (FPM_COS(FPM_DIV(FPM_MUL(d,FPM_PI),0xB40000))) | ||
184 | #define FPM_TAN_DEG(d) (FPM_TAN(FPM_DIV(FPM_MUL(d,FPM_PI),0xB40000))) | ||
185 | #define FPM_ATAN_DEG(d) (FPM_ATAN(FPM_DIV(FPM_MUL(d,FPM_PI),0xB40000))) | ||
186 | |||
187 | */ | ||
188 | /* | ||
189 | #define FPM_PI 3.14 | ||
190 | #define FPM_2PI (3.14*2) | ||
191 | #define FPM_E 178144L | ||
192 | #define FPM_ROOT2 74804L | ||
193 | #define FPM_ROOT3 113512L | ||
194 | #define FPM_GOLDEN 106039L | ||
195 | |||
196 | |||
197 | typedef float fixedpoint_t; | ||
198 | //This variable must be declared in one of the implementation files. | ||
199 | extern __int64 FPM_TMPVAR_INT64; | ||
200 | |||
201 | #define FPM_FROMFLOAT(f) (f) | ||
202 | #define FPM_TOFLOAT(fxp) (fxp) | ||
203 | |||
204 | #define FPM_FROMLONG(l) ((float)l) | ||
205 | #define FPM_TOLONG(fxp) ((long)fxp) | ||
206 | |||
207 | #define FPM_ADD(f1, f2) ((f1)+(f2)) | ||
208 | #define FPM_ADD3(f1, f2, f3) ((f1)+(f2)+(f3)) | ||
209 | #define FPM_INC(f1) ((f1)++) | ||
210 | |||
211 | #define FPM_SUB(f1, f2) ((f1)-(f2)) | ||
212 | #define FPM_DEC(f1) ((f1)--) | ||
213 | |||
214 | #define FPM_MUL(f1, f2) ((f1)*(f2)) | ||
215 | //#define FPM_MUL(f1, f2) ((fixedpoint_t)((FPM_TMPVAR_INT64=(f1))*(f2))>>16) | ||
216 | //#define FPM_MUL(f1, f2) (((f1)*(f2))>>16) | ||
217 | //TODO: This needs to be done without copying to another var | ||
218 | #define FPM_DIV(n, d) ((n)/(d)) | ||
219 | //#define FPM_DIV(n, d) ((long)(((__int64)n)<<16)/(d)) | ||
220 | |||
221 | #define FPM_ABS(fxp) (abs(fxp)) | ||
222 | |||
223 | //TODO: Implement ceil mathematically instead of converting to float and back | ||
224 | #define FPM_CEIL(fxp) (ceil(fxp)) | ||
225 | #define FPM_FLOOR(fxp) (floor(fxp)) | ||
226 | |||
227 | //TODO: Implement sqrt mathematically instead of converting to float and back | ||
228 | #define FPM_SQRT(fxp) (sqrt(fxp)) | ||
229 | #define FPM_SQR(fxp) ((fxp)*(fxp)) | ||
230 | |||
231 | #define FPM_INV(fxp) (1/(fxp)) | ||
232 | //TODO: Calc trig functions (or lookup) instead of converting to float and back | ||
233 | //These take radians | ||
234 | #define FPM_SIN(r) (sin(r)) | ||
235 | #define FPM_COS(r) (cos(r)) | ||
236 | #define FPM_TAN(r) (tan(r)) | ||
237 | #define FPM_ATAN(r) (atan(r)) | ||
238 | |||
239 | //These take degrees | ||
240 | #define FPM_SIN_DEG(d) FPM_SIN(((d)*FPM_PI)/180.0) //0xB40000 = 180.0 | ||
241 | #define FPM_COS_DEG(d) FPM_COS(((d)*FPM_PI)/180.0) | ||
242 | #define FPM_TAN_DEG(d) FPM_TAN(((d)*FPM_PI)/180.0) | ||
243 | #define FPM_ATAN_DEG(d) FPM_ATAN(((d)*FPM_PI)/180.0) | ||
244 | */ | ||
245 | |||
246 | #endif //_FIXEDPOINTMATH_H | ||