1 /*
2 * Copyright 2007 Vijay Kiran Kamuju
3 * Copyright 2007 David Adam
4 *
5 * This library is free software; you can redistribute it and/or
6 * modify it under the terms of the GNU Lesser General Public
7 * License as published by the Free Software Foundation; either
8 * version 2.1 of the License, or (at your option) any later version.
9 *
10 * This library is distributed in the hope that it will be useful,
11 * but WITHOUT ANY WARRANTY; without even the implied warranty of
12 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
13 * Lesser General Public License for more details.
14 *
15 * You should have received a copy of the GNU Lesser General Public
16 * License along with this library; if not, write to the Free Software
17 * Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301, USA
18 */
19
20 #include <math.h>
21
22 #include "d3drmdef.h"
23
24 #include "wine/test.h"
25
26 #define PI (4.0f*atanf(1.0f))
27 #define admit_error 0.000001f
28
29 #define expect_mat( expectedmat, gotmat)\
30 { \
31 int i,j; \
32 BOOL equal = TRUE; \
33 for (i=0; i<4; i++)\
34 {\
35 for (j=0; j<4; j++)\
36 {\
37 if (fabs(expectedmat[i][j]-gotmat[i][j])>admit_error)\
38 {\
39 equal = FALSE;\
40 }\
41 }\
42 }\
43 ok(equal, "Expected matrix=\n(%f,%f,%f,%f\n %f,%f,%f,%f\n %f,%f,%f,%f\n %f,%f,%f,%f\n)\n\n" \
44 "Got matrix=\n(%f,%f,%f,%f\n %f,%f,%f,%f\n %f,%f,%f,%f\n %f,%f,%f,%f)\n", \
45 expectedmat[0][0],expectedmat[0][1],expectedmat[0][2],expectedmat[0][3], \
46 expectedmat[1][0],expectedmat[1][1],expectedmat[1][2],expectedmat[1][3], \
47 expectedmat[2][0],expectedmat[2][1],expectedmat[2][2],expectedmat[2][3], \
48 expectedmat[3][0],expectedmat[3][1],expectedmat[3][2],expectedmat[3][3], \
49 gotmat[0][0],gotmat[0][1],gotmat[0][2],gotmat[0][3], \
50 gotmat[1][0],gotmat[1][1],gotmat[1][2],gotmat[1][3], \
51 gotmat[2][0],gotmat[2][1],gotmat[2][2],gotmat[2][3], \
52 gotmat[3][0],gotmat[3][1],gotmat[3][2],gotmat[3][3] ); \
53 }
54
55 #define expect_quat(expectedquat,gotquat) \
56 ok( (fabs(U1(expectedquat.v).x-U1(gotquat.v).x)<admit_error) && \
57 (fabs(U2(expectedquat.v).y-U2(gotquat.v).y)<admit_error) && \
58 (fabs(U3(expectedquat.v).z-U3(gotquat.v).z)<admit_error) && \
59 (fabs(expectedquat.s-gotquat.s)<admit_error), \
60 "Expected Quaternion %f %f %f %f , Got Quaternion %f %f %f %f\n", \
61 expectedquat.s,U1(expectedquat.v).x,U2(expectedquat.v).y,U3(expectedquat.v).z, \
62 gotquat.s,U1(gotquat.v).x,U2(gotquat.v).y,U3(gotquat.v).z);
63
64 #define expect_vec(expectedvec,gotvec) \
65 ok( ((fabs(U1(expectedvec).x-U1(gotvec).x)<admit_error)&&(fabs(U2(expectedvec).y-U2(gotvec).y)<admit_error)&&(fabs(U3(expectedvec).z-U3(gotvec).z)<admit_error)), \
66 "Expected Vector= (%f, %f, %f)\n , Got Vector= (%f, %f, %f)\n", \
67 U1(expectedvec).x,U2(expectedvec).y,U3(expectedvec).z, U1(gotvec).x, U2(gotvec).y, U3(gotvec).z);
68
VectorTest(void)69 static void VectorTest(void)
70 {
71 D3DVALUE mod,par,theta;
72 D3DVECTOR e,r,u,v,w,axis,casnul,norm,ray,self;
73
74 U1(u).x=2.0f; U2(u).y=2.0f; U3(u).z=1.0f;
75 U1(v).x=4.0f; U2(v).y=4.0f; U3(v).z=0.0f;
76
77
78 /*______________________VectorAdd_________________________________*/
79 D3DRMVectorAdd(&r,&u,&v);
80 U1(e).x=6.0f; U2(e).y=6.0f; U3(e).z=1.0f;
81 expect_vec(e,r);
82
83 U1(self).x=9.0f; U2(self).y=18.0f; U3(self).z=27.0f;
84 D3DRMVectorAdd(&self,&self,&u);
85 U1(e).x=11.0f; U2(e).y=20.0f; U3(e).z=28.0f;
86 expect_vec(e,self);
87
88 /*_______________________VectorSubtract__________________________*/
89 D3DRMVectorSubtract(&r,&u,&v);
90 U1(e).x=-2.0f; U2(e).y=-2.0f; U3(e).z=1.0f;
91 expect_vec(e,r);
92
93 U1(self).x=9.0f; U2(self).y=18.0f; U3(self).z=27.0f;
94 D3DRMVectorSubtract(&self,&self,&u);
95 U1(e).x=7.0f; U2(e).y=16.0f; U3(e).z=26.0f;
96 expect_vec(e,self);
97
98 /*_______________________VectorCrossProduct_______________________*/
99 D3DRMVectorCrossProduct(&r,&u,&v);
100 U1(e).x=-4.0f; U2(e).y=4.0f; U3(e).z=0.0f;
101 expect_vec(e,r);
102
103 U1(self).x=9.0f; U2(self).y=18.0f; U3(self).z=27.0f;
104 D3DRMVectorCrossProduct(&self,&self,&u);
105 U1(e).x=-36.0f; U2(e).y=45.0f; U3(e).z=-18.0f;
106 expect_vec(e,self);
107
108 /*_______________________VectorDotProduct__________________________*/
109 mod=D3DRMVectorDotProduct(&u,&v);
110 ok((mod == 16.0f), "Expected 16.0f, Got %f\n", mod);
111
112 /*_______________________VectorModulus_____________________________*/
113 mod=D3DRMVectorModulus(&u);
114 ok((mod == 3.0f), "Expected 3.0f, Got %f\n", mod);
115
116 /*_______________________VectorNormalize___________________________*/
117 D3DRMVectorNormalize(&u);
118 U1(e).x=2.0f/3.0f; U2(e).y=2.0f/3.0f; U3(e).z=1.0f/3.0f;
119 expect_vec(e,u);
120
121 /* If u is the NULL vector, MSDN says that the return vector is NULL. In fact, the returned vector is (1,0,0). The following test case prove it. */
122
123 U1(casnul).x=0.0f; U2(casnul).y=0.0f; U3(casnul).z=0.0f;
124 D3DRMVectorNormalize(&casnul);
125 U1(e).x=1.0f; U2(e).y=0.0f; U3(e).z=0.0f;
126 expect_vec(e,casnul);
127
128 /*____________________VectorReflect_________________________________*/
129 U1(ray).x=3.0f; U2(ray).y=-4.0f; U3(ray).z=5.0f;
130 U1(norm).x=1.0f; U2(norm).y=-2.0f; U3(norm).z=6.0f;
131 U1(e).x=79.0f; U2(e).y=-160.0f; U3(e).z=487.0f;
132 D3DRMVectorReflect(&r,&ray,&norm);
133 expect_vec(e,r);
134
135 /*_______________________VectorRotate_______________________________*/
136 U1(w).x=3.0f; U2(w).y=4.0f; U3(w).z=0.0f;
137 U1(axis).x=0.0f; U2(axis).y=0.0f; U3(axis).z=1.0f;
138 theta=2.0f*PI/3.0f;
139 D3DRMVectorRotate(&r,&w,&axis,theta);
140 U1(e).x=-0.3f-0.4f*sqrtf(3.0f); U2(e).y=0.3f*sqrtf(3.0f)-0.4f; U3(e).z=0.0f;
141 expect_vec(e,r);
142
143 /* The same formula gives D3DRMVectorRotate, for theta in [-PI/2;+PI/2] or not. The following test proves this fact.*/
144 theta=-PI/4.0f;
145 D3DRMVectorRotate(&r,&w,&axis,theta);
146 U1(e).x=1.4f/sqrtf(2.0f); U2(e).y=0.2f/sqrtf(2.0f); U3(e).z=0.0f;
147 expect_vec(e,r);
148
149 theta=PI/8.0f;
150 D3DRMVectorRotate(&self,&self,&axis,theta);
151 U1(e).x=0.989950; U2(e).y=0.141421f; U3(e).z=0.0f;
152 expect_vec(e,r);
153
154 /*_______________________VectorScale__________________________*/
155 par=2.5f;
156 D3DRMVectorScale(&r,&v,par);
157 U1(e).x=10.0f; U2(e).y=10.0f; U3(e).z=0.0f;
158 expect_vec(e,r);
159
160 U1(self).x=9.0f; U2(self).y=18.0f; U3(self).z=27.0f;
161 D3DRMVectorScale(&self,&self,2);
162 U1(e).x=18.0f; U2(e).y=36.0f; U3(e).z=54.0f;
163 expect_vec(e,self);
164 }
165
MatrixTest(void)166 static void MatrixTest(void)
167 {
168 D3DRMQUATERNION q;
169 D3DRMMATRIX4D exp,mat;
170
171 exp[0][0]=-49.0f; exp[0][1]=4.0f; exp[0][2]=22.0f; exp[0][3]=0.0f;
172 exp[1][0]=20.0f; exp[1][1]=-39.0f; exp[1][2]=20.0f; exp[1][3]=0.0f;
173 exp[2][0]=10.0f; exp[2][1]=28.0f; exp[2][2]=-25.0f; exp[2][3]=0.0f;
174 exp[3][0]=0.0f; exp[3][1]=0.0f; exp[3][2]=0.0f; exp[3][3]=1.0f;
175 q.s=1.0f; U1(q.v).x=2.0f; U2(q.v).y=3.0f; U3(q.v).z=4.0f;
176
177 D3DRMMatrixFromQuaternion(mat,&q);
178 expect_mat(exp,mat);
179 }
180
QuaternionTest(void)181 static void QuaternionTest(void)
182 {
183 D3DVECTOR axis;
184 D3DVALUE par,theta;
185 D3DRMQUATERNION q,q1,q1final,q2,q2final,r;
186
187 /*_________________QuaternionFromRotation___________________*/
188 U1(axis).x=1.0f; U2(axis).y=1.0f; U3(axis).z=1.0f;
189 theta=2.0f*PI/3.0f;
190 D3DRMQuaternionFromRotation(&r,&axis,theta);
191 q.s=0.5f; U1(q.v).x=0.5f; U2(q.v).y=0.5f; U3(q.v).z=0.5f;
192 expect_quat(q,r);
193
194 /*_________________QuaternionSlerp_________________________*/
195 /* If the angle of the two quaternions is in ]PI/2;3PI/2[, QuaternionSlerp
196 * interpolates between the first quaternion and the opposite of the second one.
197 * The test proves this fact. */
198 par=0.31f;
199 q1.s=1.0f; U1(q1.v).x=2.0f; U2(q1.v).y=3.0f; U3(q1.v).z=50.0f;
200 q2.s=-4.0f; U1(q2.v).x=6.0f; U2(q2.v).y=7.0f; U3(q2.v).z=8.0f;
201 /* The angle between q1 and q2 is in [-PI/2,PI/2]. So, one interpolates between q1 and q2. */
202 q.s = -0.55f; U1(q.v).x=3.24f; U2(q.v).y=4.24f; U3(q.v).z=36.98f;
203 D3DRMQuaternionSlerp(&r,&q1,&q2,par);
204 expect_quat(q,r);
205
206 q1.s=1.0f; U1(q1.v).x=2.0f; U2(q1.v).y=3.0f; U3(q1.v).z=50.0f;
207 q2.s=-94.0f; U1(q2.v).x=6.0f; U2(q2.v).y=7.0f; U3(q2.v).z=-8.0f;
208 /* The angle between q1 and q2 is not in [-PI/2,PI/2]. So, one interpolates between q1 and -q2. */
209 q.s=29.83f; U1(q.v).x=-0.48f; U2(q.v).y=-0.10f; U3(q.v).z=36.98f;
210 D3DRMQuaternionSlerp(&r,&q1,&q2,par);
211 expect_quat(q,r);
212
213 /* Test the spherical interpolation part */
214 q1.s=0.1f; U1(q1.v).x=0.2f; U2(q1.v).y=0.3f; U3(q1.v).z=0.4f;
215 q2.s=0.5f; U1(q2.v).x=0.6f; U2(q2.v).y=0.7f; U3(q2.v).z=0.8f;
216 q.s = 0.243943f; U1(q.v).x = 0.351172f; U2(q.v).y = 0.458401f; U3(q.v).z = 0.565629f;
217
218 q1final=q1;
219 q2final=q2;
220 D3DRMQuaternionSlerp(&r,&q1,&q2,par);
221 expect_quat(q,r);
222
223 /* Test to show that the input quaternions are not changed */
224 expect_quat(q1,q1final);
225 expect_quat(q2,q2final);
226 }
227
ColorTest(void)228 static void ColorTest(void)
229 {
230 D3DCOLOR color, expected_color, got_color;
231 D3DVALUE expected, got, red, green, blue, alpha;
232
233 /*___________D3DRMCreateColorRGB_________________________*/
234 red=0.8f;
235 green=0.3f;
236 blue=0.55f;
237 expected_color=0xffcc4c8c;
238 got_color=D3DRMCreateColorRGB(red,green,blue);
239 ok((expected_color==got_color),"Expected color=%x, Got color=%x\n",expected_color,got_color);
240
241 /*___________D3DRMCreateColorRGBA________________________*/
242 red=0.1f;
243 green=0.4f;
244 blue=0.7f;
245 alpha=0.58f;
246 expected_color=0x931966b2;
247 got_color=D3DRMCreateColorRGBA(red,green,blue,alpha);
248 ok((expected_color==got_color),"Expected color=%x, Got color=%x\n",expected_color,got_color);
249
250 /* if a component is <0 then, then one considers this component as 0. The following test proves this fact (test only with the red component). */
251 red=-0.88f;
252 green=0.4f;
253 blue=0.6f;
254 alpha=0.41f;
255 expected_color=0x68006699;
256 got_color=D3DRMCreateColorRGBA(red,green,blue,alpha);
257 ok((expected_color==got_color),"Expected color=%x, Got color=%x\n",expected_color,got_color);
258
259 /* if a component is >1 then, then one considers this component as 1. The following test proves this fact (test only with the red component). */
260 red=2.37f;
261 green=0.4f;
262 blue=0.6f;
263 alpha=0.41f;
264 expected_color=0x68ff6699;
265 got_color=D3DRMCreateColorRGBA(red,green,blue,alpha);
266 ok((expected_color==got_color),"Expected color=%x, Got color=%x\n",expected_color,got_color);
267
268 /*___________D3DRMColorGetAlpha_________________________*/
269 color=0x0e4921bf;
270 expected=14.0f/255.0f;
271 got=D3DRMColorGetAlpha(color);
272 ok((fabs(expected-got)<admit_error),"Expected=%f, Got=%f\n",expected,got);
273
274 /*___________D3DRMColorGetBlue__________________________*/
275 color=0xc82a1455;
276 expected=1.0f/3.0f;
277 got=D3DRMColorGetBlue(color);
278 ok((fabs(expected-got)<admit_error),"Expected=%f, Got=%f\n",expected,got);
279
280 /*___________D3DRMColorGetGreen_________________________*/
281 color=0xad971203;
282 expected=6.0f/85.0f;
283 got=D3DRMColorGetGreen(color);
284 ok((fabs(expected-got)<admit_error),"Expected=%f, Got=%f\n",expected,got);
285
286 /*___________D3DRMColorGetRed__________________________*/
287 color=0xb62d7a1c;
288 expected=3.0f/17.0f;
289 got=D3DRMColorGetRed(color);
290 ok((fabs(expected-got)<admit_error),"Expected=%f, Got=%f\n",expected,got);
291 }
292
START_TEST(vector)293 START_TEST(vector)
294 {
295 VectorTest();
296 MatrixTest();
297 QuaternionTest();
298 ColorTest();
299 }
300