d3dx9_36: Implement D3DXCreateMesh and initial ID3DXMesh methods.
[wine] / dlls / d3dx9_36 / tests / math.c
1 /*
2  * Copyright 2008 David Adam
3  * Copyright 2008 Luis Busquets
4  * Copyright 2008 Philip Nilsson
5  * Copyright 2008 Henri Verbeet
6  *
7  * This library is free software; you can redistribute it and/or
8  * modify it under the terms of the GNU Lesser General Public
9  * License as published by the Free Software Foundation; either
10  * version 2.1 of the License, or (at your option) any later version.
11  *
12  * This library is distributed in the hope that it will be useful,
13  * but WITHOUT ANY WARRANTY; without even the implied warranty of
14  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
15  * Lesser General Public License for more details.
16  *
17  * You should have received a copy of the GNU Lesser General Public
18  * License along with this library; if not, write to the Free Software
19  * Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301, USA
20  */
21
22 #include "wine/test.h"
23 #include "d3dx9.h"
24
25 #define ARRAY_SIZE 5
26
27 #define admitted_error 0.0001f
28
29 #define relative_error(exp, out) ((exp == 0.0f) ? fabs(exp - out) : (fabs(1.0f - out/ exp) ))
30
31 #define expect_color(expectedcolor,gotcolor) ok((relative_error(expectedcolor.r, gotcolor.r)<admitted_error)&&(relative_error(expectedcolor.g, gotcolor.g)<admitted_error)&&(relative_error(expectedcolor.b, gotcolor.b)<admitted_error)&&(relative_error(expectedcolor.a, gotcolor.a)<admitted_error),"Expected Color= (%f, %f, %f, %f)\n , Got Color= (%f, %f, %f, %f)\n", expectedcolor.r, expectedcolor.g, expectedcolor.b, expectedcolor.a, gotcolor.r, gotcolor.g, gotcolor.b, gotcolor.a);
32
33 static inline BOOL compare_matrix(const D3DXMATRIX *m1, const D3DXMATRIX *m2)
34 {
35     int i, j;
36
37     for (i = 0; i < 4; ++i)
38     {
39         for (j = 0; j < 4; ++j)
40         {
41             if (relative_error(U(*m1).m[i][j], U(*m2).m[i][j]) > admitted_error)
42                 return FALSE;
43         }
44     }
45
46     return TRUE;
47 }
48
49 #define expect_mat(expectedmat, gotmat) \
50 do { \
51     const D3DXMATRIX *__m1 = (expectedmat); \
52     const D3DXMATRIX *__m2 = (gotmat); \
53     ok(compare_matrix(__m1, __m2), "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" \
54             "Got matrix=\n(%f,%f,%f,%f\n %f,%f,%f,%f\n %f,%f,%f,%f\n %f,%f,%f,%f)\n", \
55             U(*__m1).m[0][0], U(*__m1).m[0][1], U(*__m1).m[0][2], U(*__m1).m[0][3], \
56             U(*__m1).m[1][0], U(*__m1).m[1][1], U(*__m1).m[1][2], U(*__m1).m[1][3], \
57             U(*__m1).m[2][0], U(*__m1).m[2][1], U(*__m1).m[2][2], U(*__m1).m[2][3], \
58             U(*__m1).m[3][0], U(*__m1).m[3][1], U(*__m1).m[3][2], U(*__m1).m[3][3], \
59             U(*__m2).m[0][0], U(*__m2).m[0][1], U(*__m2).m[0][2], U(*__m2).m[0][3], \
60             U(*__m2).m[1][0], U(*__m2).m[1][1], U(*__m2).m[1][2], U(*__m2).m[1][3], \
61             U(*__m2).m[2][0], U(*__m2).m[2][1], U(*__m2).m[2][2], U(*__m2).m[2][3], \
62             U(*__m2).m[3][0], U(*__m2).m[3][1], U(*__m2).m[3][2], U(*__m2).m[3][3]); \
63 } while(0)
64
65 #define compare_rotation(exp, got) \
66     ok(relative_error(exp.w, got.w) < admitted_error && \
67        relative_error(exp.x, got.x) < admitted_error && \
68        relative_error(exp.y, got.y) < admitted_error && \
69        relative_error(exp.z, got.z) < admitted_error, \
70        "Expected rotation = (%f, %f, %f, %f), \
71         got rotation = (%f, %f, %f, %f)\n", \
72         exp.w, exp.x, exp.y, exp.z, got.w, got.x, got.y, got.z)
73
74 #define compare_scale(exp, got) \
75     ok(relative_error(exp.x, got.x) < admitted_error && \
76        relative_error(exp.y, got.y) < admitted_error && \
77        relative_error(exp.z, got.z) < admitted_error, \
78        "Expected scale = (%f, %f, %f), \
79         got scale = (%f, %f, %f)\n", \
80         exp.x, exp.y, exp.z, got.x, got.y, got.z)
81
82 #define compare_translation(exp, got) \
83     ok(relative_error(exp.x, got.x) < admitted_error && \
84        relative_error(exp.y, got.y) < admitted_error && \
85        relative_error(exp.z, got.z) < admitted_error, \
86        "Expected translation = (%f, %f, %f), \
87         got translation = (%f, %f, %f)\n", \
88         exp.x, exp.y, exp.z, got.x, got.y, got.z)
89
90 #define compare_vectors(exp, out) \
91     for (i = 0; i < ARRAY_SIZE + 2; ++i) { \
92         ok(relative_error(exp[i].x, out[i].x) < admitted_error && \
93            relative_error(exp[i].y, out[i].y) < admitted_error && \
94            relative_error(exp[i].z, out[i].z) < admitted_error && \
95            relative_error(exp[i].w, out[i].w) < admitted_error, \
96             "Got (%f, %f, %f, %f), expected (%f, %f, %f, %f) for index %d.\n", \
97             out[i].x, out[i].y, out[i].z, out[i].w, \
98             exp[i].x, exp[i].y, exp[i].z, exp[i].w, \
99             i); \
100     }
101
102 #define compare_planes(exp, out) \
103     for (i = 0; i < ARRAY_SIZE + 2; ++i) { \
104         ok(relative_error(exp[i].a, out[i].a) < admitted_error && \
105            relative_error(exp[i].b, out[i].b) < admitted_error && \
106            relative_error(exp[i].c, out[i].c) < admitted_error && \
107            relative_error(exp[i].d, out[i].d) < admitted_error, \
108             "Got (%f, %f, %f, %f), expected (%f, %f, %f, %f) for index %d.\n", \
109             out[i].a, out[i].b, out[i].c, out[i].d, \
110             exp[i].a, exp[i].b, exp[i].c, exp[i].d, \
111             i); \
112     }
113
114 #define expect_plane(expectedplane,gotplane) ok((relative_error(expectedplane.a, gotplane.a)<admitted_error)&&(relative_error(expectedplane.b, gotplane.b)<admitted_error)&&(relative_error(expectedplane.c, gotplane.c)<admitted_error)&&(relative_error(expectedplane.d, gotplane.d)<admitted_error),"Expected Plane= (%f, %f, %f, %f)\n , Got Plane= (%f, %f, %f, %f)\n", expectedplane.a, expectedplane.b, expectedplane.c, expectedplane.d, gotplane.a, gotplane.b, gotplane.c, gotplane.d);
115
116 #define expect_vec(expectedvec,gotvec) ok((relative_error(expectedvec.x, gotvec.x)<admitted_error)&&(relative_error(expectedvec.y, gotvec.y)<admitted_error),"Expected Vector= (%f, %f)\n , Got Vector= (%f, %f)\n", expectedvec.x, expectedvec.y, gotvec.x, gotvec.y);
117
118 #define expect_vec3(expectedvec,gotvec) ok((relative_error(expectedvec.x, gotvec.x)<admitted_error)&&(relative_error(expectedvec.y, gotvec.y)<admitted_error)&&(relative_error(expectedvec.z, gotvec.z)<admitted_error),"Expected Vector= (%f, %f, %f)\n , Got Vector= (%f, %f, %f)\n", expectedvec.x, expectedvec.y, expectedvec.z, gotvec.x, gotvec.y, gotvec.z);
119
120 #define expect_vec4(expectedvec,gotvec) ok((relative_error(expectedvec.x, gotvec.x)<admitted_error)&&(relative_error(expectedvec.y, gotvec.y)<admitted_error)&&(relative_error(expectedvec.z, gotvec.z)<admitted_error)&&(relative_error(expectedvec.w, gotvec.w)<admitted_error),"Expected Vector= (%f, %f, %f, %f)\n , Got Vector= (%f, %f, %f, %f)\n", expectedvec.x, expectedvec.y, expectedvec.z, expectedvec.w, gotvec.x, gotvec.y, gotvec.z, gotvec.w);
121
122
123 static void D3DXColorTest(void)
124 {
125     D3DXCOLOR color, color1, color2, expected, got;
126     LPD3DXCOLOR funcpointer;
127     FLOAT scale;
128
129     color.r = 0.2f; color.g = 0.75f; color.b = 0.41f; color.a = 0.93f;
130     color1.r = 0.6f; color1.g = 0.55f; color1.b = 0.23f; color1.a = 0.82f;
131     color2.r = 0.3f; color2.g = 0.5f; color2.b = 0.76f; color2.a = 0.11f;
132
133     scale = 0.3f;
134
135 /*_______________D3DXColorAdd________________*/
136     expected.r = 0.9f; expected.g = 1.05f; expected.b = 0.99f, expected.a = 0.93f;
137     D3DXColorAdd(&got,&color1,&color2);
138     expect_color(expected,got);
139     /* Test the NULL case */
140     funcpointer = D3DXColorAdd(&got,NULL,&color2);
141     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
142     funcpointer = D3DXColorAdd(NULL,NULL,&color2);
143     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
144     funcpointer = D3DXColorAdd(NULL,NULL,NULL);
145     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
146
147 /*_______________D3DXColorAdjustContrast______*/
148     expected.r = 0.41f; expected.g = 0.575f; expected.b = 0.473f, expected.a = 0.93f;
149     D3DXColorAdjustContrast(&got,&color,scale);
150     expect_color(expected,got);
151
152 /*_______________D3DXColorAdjustSaturation______*/
153     expected.r = 0.486028f; expected.g = 0.651028f; expected.b = 0.549028f, expected.a = 0.93f;
154     D3DXColorAdjustSaturation(&got,&color,scale);
155     expect_color(expected,got);
156
157 /*_______________D3DXColorLerp________________*/
158     expected.r = 0.32f; expected.g = 0.69f; expected.b = 0.356f; expected.a = 0.897f;
159     D3DXColorLerp(&got,&color,&color1,scale);
160     expect_color(expected,got);
161     /* Test the NULL case */
162     funcpointer = D3DXColorLerp(&got,NULL,&color1,scale);
163     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
164     funcpointer = D3DXColorLerp(NULL,NULL,&color1,scale);
165     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
166     funcpointer = D3DXColorLerp(NULL,NULL,NULL,scale);
167     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
168
169 /*_______________D3DXColorModulate________________*/
170     expected.r = 0.18f; expected.g = 0.275f; expected.b = 0.1748f; expected.a = 0.0902f;
171     D3DXColorModulate(&got,&color1,&color2);
172     expect_color(expected,got);
173     /* Test the NULL case */
174     funcpointer = D3DXColorModulate(&got,NULL,&color2);
175     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
176     funcpointer = D3DXColorModulate(NULL,NULL,&color2);
177     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
178     funcpointer = D3DXColorModulate(NULL,NULL,NULL);
179     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
180
181 /*_______________D3DXColorNegative________________*/
182     expected.r = 0.8f; expected.g = 0.25f; expected.b = 0.59f; expected.a = 0.93f;
183     D3DXColorNegative(&got,&color);
184     expect_color(got,expected);
185     /* Test the greater than 1 case */
186     color1.r = 0.2f; color1.g = 1.75f; color1.b = 0.41f; color1.a = 0.93f;
187     expected.r = 0.8f; expected.g = -0.75f; expected.b = 0.59f; expected.a = 0.93f;
188     D3DXColorNegative(&got,&color1);
189     expect_color(got,expected);
190     /* Test the negative case */
191     color1.r = 0.2f; color1.g = -0.75f; color1.b = 0.41f; color1.a = 0.93f;
192     expected.r = 0.8f; expected.g = 1.75f; expected.b = 0.59f; expected.a = 0.93f;
193     D3DXColorNegative(&got,&color1);
194     expect_color(got,expected);
195     /* Test the NULL case */
196     funcpointer = D3DXColorNegative(&got,NULL);
197     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
198     funcpointer = D3DXColorNegative(NULL,NULL);
199     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
200
201 /*_______________D3DXColorScale________________*/
202     expected.r = 0.06f; expected.g = 0.225f; expected.b = 0.123f; expected.a = 0.279f;
203     D3DXColorScale(&got,&color,scale);
204     expect_color(expected,got);
205     /* Test the NULL case */
206     funcpointer = D3DXColorScale(&got,NULL,scale);
207     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
208     funcpointer = D3DXColorScale(NULL,NULL,scale);
209     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
210
211 /*_______________D3DXColorSubtract_______________*/
212     expected.r = -0.1f; expected.g = 0.25f; expected.b = -0.35f, expected.a = 0.82f;
213     D3DXColorSubtract(&got,&color,&color2);
214     expect_color(expected,got);
215     /* Test the NULL case */
216     funcpointer = D3DXColorSubtract(&got,NULL,&color2);
217     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
218     funcpointer = D3DXColorSubtract(NULL,NULL,&color2);
219     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
220     funcpointer = D3DXColorSubtract(NULL,NULL,NULL);
221     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
222 }
223
224 static void D3DXFresnelTest(void)
225 {
226     FLOAT expected, got;
227
228     expected = 0.089187;
229     got = D3DXFresnelTerm(0.5f,1.5);
230     ok(relative_error(got, expected) < admitted_error, "Expected: %f, Got: %f\n", expected, got);
231 }
232
233 static void D3DXMatrixTest(void)
234 {
235     D3DXMATRIX expectedmat, gotmat, mat, mat2, mat3;
236     LPD3DXMATRIX funcpointer;
237     D3DXPLANE plane;
238     D3DXQUATERNION q, r;
239     D3DXVECTOR3 at, axis, eye, last;
240     D3DXVECTOR4 light;
241     BOOL expected, got;
242     FLOAT angle, determinant, expectedfloat, gotfloat;
243
244     U(mat).m[0][1] = 5.0f; U(mat).m[0][2] = 7.0f; U(mat).m[0][3] = 8.0f;
245     U(mat).m[1][0] = 11.0f; U(mat).m[1][2] = 16.0f; U(mat).m[1][3] = 33.0f;
246     U(mat).m[2][0] = 19.0f; U(mat).m[2][1] = -21.0f; U(mat).m[2][3] = 43.0f;
247     U(mat).m[3][0] = 2.0f; U(mat).m[3][1] = 3.0f; U(mat).m[3][2] = -4.0f;
248     U(mat).m[0][0] = 10.0f; U(mat).m[1][1] = 20.0f; U(mat).m[2][2] = 30.0f;
249     U(mat).m[3][3] = -40.0f;
250
251     U(mat2).m[0][0] = 1.0f; U(mat2).m[1][0] = 2.0f; U(mat2).m[2][0] = 3.0f;
252     U(mat2).m[3][0] = 4.0f; U(mat2).m[0][1] = 5.0f; U(mat2).m[1][1] = 6.0f;
253     U(mat2).m[2][1] = 7.0f; U(mat2).m[3][1] = 8.0f; U(mat2).m[0][2] = -8.0f;
254     U(mat2).m[1][2] = -7.0f; U(mat2).m[2][2] = -6.0f; U(mat2).m[3][2] = -5.0f;
255     U(mat2).m[0][3] = -4.0f; U(mat2).m[1][3] = -3.0f; U(mat2).m[2][3] = -2.0f;
256     U(mat2).m[3][3] = -1.0f;
257
258     plane.a = -3.0f; plane.b = -1.0f; plane.c = 4.0f; plane.d = 7.0f;
259
260     q.x = 1.0f; q.y = -4.0f; q.z =7.0f; q.w = -11.0f;
261     r.x = 0.87f; r.y = 0.65f; r.z =0.43f; r.w= 0.21f;
262
263     at.x = -2.0f; at.y = 13.0f; at.z = -9.0f;
264     axis.x = 1.0f; axis.y = -3.0f; axis.z = 7.0f;
265     eye.x = 8.0f; eye.y = -5.0f; eye.z = 5.75f;
266     last.x = 9.7f; last.y = -8.6; last.z = 1.3f;
267
268     light.x = 9.6f; light.y = 8.5f; light.z = 7.4; light.w = 6.3;
269
270     angle = D3DX_PI/3.0f;
271
272 /*____________D3DXMatrixAffineTransformation______*/
273     U(expectedmat).m[0][0] = -459.239990f; U(expectedmat).m[0][1] = -576.719971f; U(expectedmat).m[0][2] = -263.440002f; U(expectedmat).m[0][3] = 0.0f;
274     U(expectedmat).m[1][0] = 519.760010f; U(expectedmat).m[1][1] = -352.440002f; U(expectedmat).m[1][2] = -277.679993f; U(expectedmat).m[1][3] = 0.0f;
275     U(expectedmat).m[2][0] = 363.119995f; U(expectedmat).m[2][1] = -121.040001f; U(expectedmat).m[2][2] = -117.479996f; U(expectedmat).m[2][3] = 0.0f;
276     U(expectedmat).m[3][0] = -1239.0f; U(expectedmat).m[3][1] = 667.0f; U(expectedmat).m[3][2] = 567.0f; U(expectedmat).m[3][3] = 1.0f;
277     D3DXMatrixAffineTransformation(&gotmat,3.56f,&at,&q,&axis);
278     expect_mat(&expectedmat, &gotmat);
279 /* Test the NULL case */
280     U(expectedmat).m[0][0] = -459.239990f; U(expectedmat).m[0][1] = -576.719971f; U(expectedmat).m[0][2] = -263.440002f; U(expectedmat).m[0][3] = 0.0f;
281     U(expectedmat).m[1][0] = 519.760010f; U(expectedmat).m[1][1] = -352.440002f; U(expectedmat).m[1][2] = -277.679993f; U(expectedmat).m[1][3] = 0.0f;
282     U(expectedmat).m[2][0] = 363.119995f; U(expectedmat).m[2][1] = -121.040001f; U(expectedmat).m[2][2] = -117.479996f; U(expectedmat).m[2][3] = 0.0f;
283     U(expectedmat).m[3][0] = 1.0f; U(expectedmat).m[3][1] = -3.0f; U(expectedmat).m[3][2] = 7.0f; U(expectedmat).m[3][3] = 1.0f;
284     D3DXMatrixAffineTransformation(&gotmat,3.56f,NULL,&q,&axis);
285     expect_mat(&expectedmat, &gotmat);
286
287 /*____________D3DXMatrixfDeterminant_____________*/
288     expectedfloat = -147888.0f;
289     gotfloat = D3DXMatrixDeterminant(&mat);
290     ok(relative_error(gotfloat, expectedfloat ) < admitted_error, "Expected: %f, Got: %f\n", expectedfloat, gotfloat);
291
292 /*____________D3DXMatrixInverse______________*/
293     U(expectedmat).m[0][0] = 16067.0f/73944.0f; U(expectedmat).m[0][1] = -10165.0f/147888.0f; U(expectedmat).m[0][2] = -2729.0f/147888.0f; U(expectedmat).m[0][3] = -1631.0f/49296.0f;
294     U(expectedmat).m[1][0] = -565.0f/36972.0f; U(expectedmat).m[1][1] = 2723.0f/73944.0f; U(expectedmat).m[1][2] = -1073.0f/73944.0f; U(expectedmat).m[1][3] = 289.0f/24648.0f;
295     U(expectedmat).m[2][0] = -389.0f/2054.0f; U(expectedmat).m[2][1] = 337.0f/4108.0f; U(expectedmat).m[2][2] = 181.0f/4108.0f; U(expectedmat).m[2][3] = 317.0f/4108.0f;
296     U(expectedmat).m[3][0] = 163.0f/5688.0f; U(expectedmat).m[3][1] = -101.0f/11376.0f; U(expectedmat).m[3][2] = -73.0f/11376.0f; U(expectedmat).m[3][3] = -127.0f/3792.0f;
297     expectedfloat = -147888.0f;
298     D3DXMatrixInverse(&gotmat,&determinant,&mat);
299     expect_mat(&expectedmat, &gotmat);
300     ok(relative_error( determinant, expectedfloat ) < admitted_error, "Expected: %f, Got: %f\n", expectedfloat, determinant);
301     funcpointer = D3DXMatrixInverse(&gotmat,NULL,&mat2);
302     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
303
304 /*____________D3DXMatrixIsIdentity______________*/
305     expected = FALSE;
306     memset(&mat3, 0, sizeof(mat3));
307     got = D3DXMatrixIsIdentity(&mat3);
308     ok(expected == got, "Expected : %d, Got : %d\n", expected, got);
309     D3DXMatrixIdentity(&mat3);
310     expected = TRUE;
311     got = D3DXMatrixIsIdentity(&mat3);
312     ok(expected == got, "Expected : %d, Got : %d\n", expected, got);
313     U(mat3).m[0][0] = 0.000009f;
314     expected = FALSE;
315     got = D3DXMatrixIsIdentity(&mat3);
316     ok(expected == got, "Expected : %d, Got : %d\n", expected, got);
317     /* Test the NULL case */
318     expected = FALSE;
319     got = D3DXMatrixIsIdentity(NULL);
320     ok(expected == got, "Expected : %d, Got : %d\n", expected, got);
321
322 /*____________D3DXMatrixLookatLH_______________*/
323     U(expectedmat).m[0][0] = -0.822465f; U(expectedmat).m[0][1] = -0.409489f; U(expectedmat).m[0][2] = -0.394803f; U(expectedmat).m[0][3] = 0.0f;
324     U(expectedmat).m[1][0] = -0.555856f; U(expectedmat).m[1][1] = 0.431286f; U(expectedmat).m[1][2] = 0.710645f; U(expectedmat).m[1][3] = 0.0f;
325     U(expectedmat).m[2][0] = -0.120729f; U(expectedmat).m[2][1] = 0.803935f; U(expectedmat).m[2][2] = -0.582335f; U(expectedmat).m[2][3] = 0.0f;
326     U(expectedmat).m[3][0] = 4.494634f; U(expectedmat).m[3][1] = 0.809719f; U(expectedmat).m[3][2] = 10.060076f; U(expectedmat).m[3][3] = 1.0f;
327     D3DXMatrixLookAtLH(&gotmat,&eye,&at,&axis);
328     expect_mat(&expectedmat, &gotmat);
329
330 /*____________D3DXMatrixLookatRH_______________*/
331     U(expectedmat).m[0][0] = 0.822465f; U(expectedmat).m[0][1] = -0.409489f; U(expectedmat).m[0][2] = 0.394803f; U(expectedmat).m[0][3] = 0.0f;
332     U(expectedmat).m[1][0] = 0.555856f; U(expectedmat).m[1][1] = 0.431286f; U(expectedmat).m[1][2] = -0.710645f; U(expectedmat).m[1][3] = 0.0f;
333     U(expectedmat).m[2][0] = 0.120729f; U(expectedmat).m[2][1] = 0.803935f; U(expectedmat).m[2][2] = 0.582335f; U(expectedmat).m[2][3] = 0.0f;
334     U(expectedmat).m[3][0] = -4.494634f; U(expectedmat).m[3][1] = 0.809719f; U(expectedmat).m[3][2] = -10.060076f; U(expectedmat).m[3][3] = 1.0f;
335     D3DXMatrixLookAtRH(&gotmat,&eye,&at,&axis);
336     expect_mat(&expectedmat, &gotmat);
337
338 /*____________D3DXMatrixMultiply______________*/
339     U(expectedmat).m[0][0] = 73.0f; U(expectedmat).m[0][1] = 193.0f; U(expectedmat).m[0][2] = -197.0f; U(expectedmat).m[0][3] = -77.0f;
340     U(expectedmat).m[1][0] = 231.0f; U(expectedmat).m[1][1] = 551.0f; U(expectedmat).m[1][2] = -489.0f; U(expectedmat).m[1][3] = -169.0;
341     U(expectedmat).m[2][0] = 239.0f; U(expectedmat).m[2][1] = 523.0f; U(expectedmat).m[2][2] = -400.0f; U(expectedmat).m[2][3] = -116.0f;
342     U(expectedmat).m[3][0] = -164.0f; U(expectedmat).m[3][1] = -320.0f; U(expectedmat).m[3][2] = 187.0f; U(expectedmat).m[3][3] = 31.0f;
343     D3DXMatrixMultiply(&gotmat,&mat,&mat2);
344     expect_mat(&expectedmat, &gotmat);
345
346 /*____________D3DXMatrixMultiplyTranspose____*/
347     U(expectedmat).m[0][0] = 73.0f; U(expectedmat).m[0][1] = 231.0f; U(expectedmat).m[0][2] = 239.0f; U(expectedmat).m[0][3] = -164.0f;
348     U(expectedmat).m[1][0] = 193.0f; U(expectedmat).m[1][1] = 551.0f; U(expectedmat).m[1][2] = 523.0f; U(expectedmat).m[1][3] = -320.0;
349     U(expectedmat).m[2][0] = -197.0f; U(expectedmat).m[2][1] = -489.0f; U(expectedmat).m[2][2] = -400.0f; U(expectedmat).m[2][3] = 187.0f;
350     U(expectedmat).m[3][0] = -77.0f; U(expectedmat).m[3][1] = -169.0f; U(expectedmat).m[3][2] = -116.0f; U(expectedmat).m[3][3] = 31.0f;
351     D3DXMatrixMultiplyTranspose(&gotmat,&mat,&mat2);
352     expect_mat(&expectedmat, &gotmat);
353
354 /*____________D3DXMatrixOrthoLH_______________*/
355     U(expectedmat).m[0][0] = 0.8f; U(expectedmat).m[0][1] = 0.0f; U(expectedmat).m[0][2] = 0.0f; U(expectedmat).m[0][3] = 0.0f;
356     U(expectedmat).m[1][0] = 0.0f; U(expectedmat).m[1][1] = 0.270270f; U(expectedmat).m[1][2] = 0.0f; U(expectedmat).m[1][3] = 0.0f;
357     U(expectedmat).m[2][0] = 0.0f; U(expectedmat).m[2][1] = 0.0f; U(expectedmat).m[2][2] = -0.151515f; U(expectedmat).m[2][3] = 0.0f;
358     U(expectedmat).m[3][0] = 0.0f; U(expectedmat).m[3][1] = 0.0f; U(expectedmat).m[3][2] = -0.484848f; U(expectedmat).m[3][3] = 1.0f;
359     D3DXMatrixOrthoLH(&gotmat, 2.5f, 7.4f, -3.2f, -9.8f);
360     expect_mat(&expectedmat, &gotmat);
361
362 /*____________D3DXMatrixOrthoOffCenterLH_______________*/
363     U(expectedmat).m[0][0] = 3.636364f; U(expectedmat).m[0][1] = 0.0f; U(expectedmat).m[0][2] = 0.0f; U(expectedmat).m[0][3] = 0.0f;
364     U(expectedmat).m[1][0] = 0.0f; U(expectedmat).m[1][1] = 0.180180f; U(expectedmat).m[1][2] = 0.0; U(expectedmat).m[1][3] = 0.0f;
365     U(expectedmat).m[2][0] = 0.0f; U(expectedmat).m[2][1] = 0.0f; U(expectedmat).m[2][2] = -0.045662f; U(expectedmat).m[2][3] = 0.0f;
366     U(expectedmat).m[3][0] = -1.727272f; U(expectedmat).m[3][1] = -0.567568f; U(expectedmat).m[3][2] = 0.424658f; U(expectedmat).m[3][3] = 1.0f;
367     D3DXMatrixOrthoOffCenterLH(&gotmat, 0.2f, 0.75f, -2.4f, 8.7f, 9.3, -12.6);
368     expect_mat(&expectedmat, &gotmat);
369
370 /*____________D3DXMatrixOrthoOffCenterRH_______________*/
371     U(expectedmat).m[0][0] = 3.636364f; U(expectedmat).m[0][1] = 0.0f; U(expectedmat).m[0][2] = 0.0f; U(expectedmat).m[0][3] = 0.0f;
372     U(expectedmat).m[1][0] = 0.0f; U(expectedmat).m[1][1] = 0.180180f; U(expectedmat).m[1][2] = 0.0; U(expectedmat).m[1][3] = 0.0f;
373     U(expectedmat).m[2][0] = 0.0f; U(expectedmat).m[2][1] = 0.0f; U(expectedmat).m[2][2] = 0.045662f; U(expectedmat).m[2][3] = 0.0f;
374     U(expectedmat).m[3][0] = -1.727272f; U(expectedmat).m[3][1] = -0.567568f; U(expectedmat).m[3][2] = 0.424658f; U(expectedmat).m[3][3] = 1.0f;
375     D3DXMatrixOrthoOffCenterRH(&gotmat, 0.2f, 0.75f, -2.4f, 8.7f, 9.3, -12.6);
376     expect_mat(&expectedmat, &gotmat);
377
378 /*____________D3DXMatrixOrthoRH_______________*/
379     U(expectedmat).m[0][0] = 0.8f; U(expectedmat).m[0][1] = 0.0f; U(expectedmat).m[0][2] = 0.0f; U(expectedmat).m[0][3] = 0.0f;
380     U(expectedmat).m[1][0] = 0.0f; U(expectedmat).m[1][1] = 0.270270f; U(expectedmat).m[1][2] = 0.0f; U(expectedmat).m[1][3] = 0.0f;
381     U(expectedmat).m[2][0] = 0.0f; U(expectedmat).m[2][1] = 0.0f; U(expectedmat).m[2][2] = 0.151515f; U(expectedmat).m[2][3] = 0.0f;
382     U(expectedmat).m[3][0] = 0.0f; U(expectedmat).m[3][1] = 0.0f; U(expectedmat).m[3][2] = -0.484848f; U(expectedmat).m[3][3] = 1.0f;
383     D3DXMatrixOrthoRH(&gotmat, 2.5f, 7.4f, -3.2f, -9.8f);
384     expect_mat(&expectedmat, &gotmat);
385
386 /*____________D3DXMatrixPerspectiveFovLH_______________*/
387     U(expectedmat).m[0][0] = 13.288858f; U(expectedmat).m[0][1] = 0.0f; U(expectedmat).m[0][2] = 0.0f; U(expectedmat).m[0][3] = 0.0f;
388     U(expectedmat).m[1][0] = 0.0f; U(expectedmat).m[1][1] = 9.966644f; U(expectedmat).m[1][2] = 0.0; U(expectedmat).m[1][3] = 0.0f;
389     U(expectedmat).m[2][0] = 0.0f; U(expectedmat).m[2][1] = 0.0f; U(expectedmat).m[2][2] = 0.783784f; U(expectedmat).m[2][3] = 1.0f;
390     U(expectedmat).m[3][0] = 0.0f; U(expectedmat).m[3][1] = 0.0f; U(expectedmat).m[3][2] = 1.881081f; U(expectedmat).m[3][3] = 0.0f;
391     D3DXMatrixPerspectiveFovLH(&gotmat, 0.2f, 0.75f, -2.4f, 8.7f);
392     expect_mat(&expectedmat, &gotmat);
393
394 /*____________D3DXMatrixPerspectiveFovRH_______________*/
395     U(expectedmat).m[0][0] = 13.288858f; U(expectedmat).m[0][1] = 0.0f; U(expectedmat).m[0][2] = 0.0f; U(expectedmat).m[0][3] = 0.0f;
396     U(expectedmat).m[1][0] = 0.0f; U(expectedmat).m[1][1] = 9.966644f; U(expectedmat).m[1][2] = 0.0; U(expectedmat).m[1][3] = 0.0f;
397     U(expectedmat).m[2][0] = 0.0f; U(expectedmat).m[2][1] = 0.0f; U(expectedmat).m[2][2] = -0.783784f; U(expectedmat).m[2][3] = -1.0f;
398     U(expectedmat).m[3][0] = 0.0f; U(expectedmat).m[3][1] = 0.0f; U(expectedmat).m[3][2] = 1.881081f; U(expectedmat).m[3][3] = 0.0f;
399     D3DXMatrixPerspectiveFovRH(&gotmat, 0.2f, 0.75f, -2.4f, 8.7f);
400     expect_mat(&expectedmat, &gotmat);
401
402 /*____________D3DXMatrixPerspectiveLH_______________*/
403     U(expectedmat).m[0][0] = -24.0f; U(expectedmat).m[0][1] = 0.0f; U(expectedmat).m[0][2] = 0.0f; U(expectedmat).m[0][3] = 0.0f;
404     U(expectedmat).m[1][0] = 0.0f; U(expectedmat).m[1][1] = -6.4f; U(expectedmat).m[1][2] = 0.0; U(expectedmat).m[1][3] = 0.0f;
405     U(expectedmat).m[2][0] = 0.0f; U(expectedmat).m[2][1] = 0.0f; U(expectedmat).m[2][2] = 0.783784f; U(expectedmat).m[2][3] = 1.0f;
406     U(expectedmat).m[3][0] = 0.0f; U(expectedmat).m[3][1] = 0.0f; U(expectedmat).m[3][2] = 1.881081f; U(expectedmat).m[3][3] = 0.0f;
407     D3DXMatrixPerspectiveLH(&gotmat, 0.2f, 0.75f, -2.4f, 8.7f);
408     expect_mat(&expectedmat, &gotmat);
409
410 /*____________D3DXMatrixPerspectiveOffCenterLH_______________*/
411     U(expectedmat).m[0][0] = 11.636364f; U(expectedmat).m[0][1] = 0.0f; U(expectedmat).m[0][2] = 0.0f; U(expectedmat).m[0][3] = 0.0f;
412     U(expectedmat).m[1][0] = 0.0f; U(expectedmat).m[1][1] = 0.576577f; U(expectedmat).m[1][2] = 0.0; U(expectedmat).m[1][3] = 0.0f;
413     U(expectedmat).m[2][0] = -1.727273f; U(expectedmat).m[2][1] = -0.567568f; U(expectedmat).m[2][2] = 0.840796f; U(expectedmat).m[2][3] = 1.0f;
414     U(expectedmat).m[3][0] = 0.0f; U(expectedmat).m[3][1] = 0.0f; U(expectedmat).m[3][2] = -2.690547f; U(expectedmat).m[3][3] = 0.0f;
415     D3DXMatrixPerspectiveOffCenterLH(&gotmat, 0.2f, 0.75f, -2.4f, 8.7f, 3.2f, -16.9f);
416     expect_mat(&expectedmat, &gotmat);
417
418 /*____________D3DXMatrixPerspectiveOffCenterRH_______________*/
419     U(expectedmat).m[0][0] = 11.636364f; U(expectedmat).m[0][1] = 0.0f; U(expectedmat).m[0][2] = 0.0f; U(expectedmat).m[0][3] = 0.0f;
420     U(expectedmat).m[1][0] = 0.0f; U(expectedmat).m[1][1] = 0.576577f; U(expectedmat).m[1][2] = 0.0; U(expectedmat).m[1][3] = 0.0f;
421     U(expectedmat).m[2][0] = 1.727273f; U(expectedmat).m[2][1] = 0.567568f; U(expectedmat).m[2][2] = -0.840796f; U(expectedmat).m[2][3] = -1.0f;
422     U(expectedmat).m[3][0] = 0.0f; U(expectedmat).m[3][1] = 0.0f; U(expectedmat).m[3][2] = -2.690547f; U(expectedmat).m[3][3] = 0.0f;
423     D3DXMatrixPerspectiveOffCenterRH(&gotmat, 0.2f, 0.75f, -2.4f, 8.7f, 3.2f, -16.9f);
424     expect_mat(&expectedmat, &gotmat);
425
426 /*____________D3DXMatrixPerspectiveRH_______________*/
427     U(expectedmat).m[0][0] = -24.0f; U(expectedmat).m[0][1] = -0.0f; U(expectedmat).m[0][2] = 0.0f; U(expectedmat).m[0][3] = 0.0f;
428     U(expectedmat).m[1][0] = 0.0f; U(expectedmat).m[1][1] = -6.4f; U(expectedmat).m[1][2] = 0.0; U(expectedmat).m[1][3] = 0.0f;
429     U(expectedmat).m[2][0] = 0.0f; U(expectedmat).m[2][1] = 0.0f; U(expectedmat).m[2][2] = -0.783784f; U(expectedmat).m[2][3] = -1.0f;
430     U(expectedmat).m[3][0] = 0.0f; U(expectedmat).m[3][1] = 0.0f; U(expectedmat).m[3][2] = 1.881081f; U(expectedmat).m[3][3] = 0.0f;
431     D3DXMatrixPerspectiveRH(&gotmat, 0.2f, 0.75f, -2.4f, 8.7f);
432     expect_mat(&expectedmat, &gotmat);
433
434 /*____________D3DXMatrixReflect______________*/
435     U(expectedmat).m[0][0] = 0.307692f; U(expectedmat).m[0][1] = -0.230769f; U(expectedmat).m[0][2] = 0.923077f; U(expectedmat).m[0][3] = 0.0f;
436     U(expectedmat).m[1][0] = -0.230769; U(expectedmat).m[1][1] = 0.923077f; U(expectedmat).m[1][2] = 0.307693f; U(expectedmat).m[1][3] = 0.0f;
437     U(expectedmat).m[2][0] = 0.923077f; U(expectedmat).m[2][1] = 0.307693f; U(expectedmat).m[2][2] = -0.230769f; U(expectedmat).m[2][3] = 0.0f;
438     U(expectedmat).m[3][0] = 1.615385f; U(expectedmat).m[3][1] = 0.538462f; U(expectedmat).m[3][2] = -2.153846f; U(expectedmat).m[3][3] = 1.0f;
439     D3DXMatrixReflect(&gotmat,&plane);
440     expect_mat(&expectedmat, &gotmat);
441
442 /*____________D3DXMatrixRotationAxis_____*/
443     U(expectedmat).m[0][0] = 0.508475f; U(expectedmat).m[0][1] = 0.763805f; U(expectedmat).m[0][2] = 0.397563f; U(expectedmat).m[0][3] = 0.0f;
444     U(expectedmat).m[1][0] = -0.814652f; U(expectedmat).m[1][1] = 0.576271f; U(expectedmat).m[1][2] = -0.065219f; U(expectedmat).m[1][3] = 0.0f;
445     U(expectedmat).m[2][0] = -0.278919f; U(expectedmat).m[2][1] = -0.290713f; U(expectedmat).m[2][2] = 0.915254f; U(expectedmat).m[2][3] = 0.0f;
446     U(expectedmat).m[3][0] = 0.0f; U(expectedmat).m[3][1] = 0.0f; U(expectedmat).m[3][2] = 0.0f; U(expectedmat).m[3][3] = 1.0f;
447     D3DXMatrixRotationAxis(&gotmat,&axis,angle);
448     expect_mat(&expectedmat, &gotmat);
449
450 /*____________D3DXMatrixRotationQuaternion______________*/
451     U(expectedmat).m[0][0] = -129.0f; U(expectedmat).m[0][1] = -162.0f; U(expectedmat).m[0][2] = -74.0f; U(expectedmat).m[0][3] = 0.0f;
452     U(expectedmat).m[1][0] = 146.0f; U(expectedmat).m[1][1] = -99.0f; U(expectedmat).m[1][2] = -78.0f; U(expectedmat).m[1][3] = 0.0f;
453     U(expectedmat).m[2][0] = 102.0f; U(expectedmat).m[2][1] = -34.0f; U(expectedmat).m[2][2] = -33.0f; U(expectedmat).m[2][3] = 0.0f;
454     U(expectedmat).m[3][0] = 0.0f; U(expectedmat).m[3][1] = 0.0f; U(expectedmat).m[3][2] = 0.0f; U(expectedmat).m[3][3] = 1.0f;
455     D3DXMatrixRotationQuaternion(&gotmat,&q);
456     expect_mat(&expectedmat, &gotmat);
457
458 /*____________D3DXMatrixRotationX______________*/
459     U(expectedmat).m[0][0] = 1.0f; U(expectedmat).m[0][1] = 0.0f; U(expectedmat).m[0][2] = 0.0f; U(expectedmat).m[0][3] = 0.0f;
460     U(expectedmat).m[1][0] = 0.0; U(expectedmat).m[1][1] = 0.5f; U(expectedmat).m[1][2] = sqrt(3.0f)/2.0f; U(expectedmat).m[1][3] = 0.0f;
461     U(expectedmat).m[2][0] = 0.0f; U(expectedmat).m[2][1] = -sqrt(3.0f)/2.0f; U(expectedmat).m[2][2] = 0.5f; U(expectedmat).m[2][3] = 0.0f;
462     U(expectedmat).m[3][0] = 0.0f; U(expectedmat).m[3][1] = 0.0f; U(expectedmat).m[3][2] = 0.0f; U(expectedmat).m[3][3] = 1.0f;
463     D3DXMatrixRotationX(&gotmat,angle);
464     expect_mat(&expectedmat, &gotmat);
465
466 /*____________D3DXMatrixRotationY______________*/
467     U(expectedmat).m[0][0] = 0.5f; U(expectedmat).m[0][1] = 0.0f; U(expectedmat).m[0][2] = -sqrt(3.0f)/2.0f; U(expectedmat).m[0][3] = 0.0f;
468     U(expectedmat).m[1][0] = 0.0; U(expectedmat).m[1][1] = 1.0f; U(expectedmat).m[1][2] = 0.0f; U(expectedmat).m[1][3] = 0.0f;
469     U(expectedmat).m[2][0] = sqrt(3.0f)/2.0f; U(expectedmat).m[2][1] = 0.0f; U(expectedmat).m[2][2] = 0.5f; U(expectedmat).m[2][3] = 0.0f;
470     U(expectedmat).m[3][0] = 0.0f; U(expectedmat).m[3][1] = 0.0f; U(expectedmat).m[3][2] = 0.0f; U(expectedmat).m[3][3] = 1.0f;
471     D3DXMatrixRotationY(&gotmat,angle);
472     expect_mat(&expectedmat, &gotmat);
473
474 /*____________D3DXMatrixRotationYawPitchRoll____*/
475     U(expectedmat).m[0][0] = 0.888777f; U(expectedmat).m[0][1] = 0.091875f; U(expectedmat).m[0][2] = -0.449037f; U(expectedmat).m[0][3] = 0.0f;
476     U(expectedmat).m[1][0] = 0.351713f; U(expectedmat).m[1][1] = 0.491487f; U(expectedmat).m[1][2] = 0.796705f; U(expectedmat).m[1][3] = 0.0f;
477     U(expectedmat).m[2][0] = 0.293893f; U(expectedmat).m[2][1] = -0.866025f; U(expectedmat).m[2][2] = 0.404509f; U(expectedmat).m[2][3] = 0.0f;
478     U(expectedmat).m[3][0] = 0.0f; U(expectedmat).m[3][1] = 0.0f; U(expectedmat).m[3][2] = 0.0f; U(expectedmat).m[3][3] = 1.0f;
479     D3DXMatrixRotationYawPitchRoll(&gotmat, 3.0f*angle/5.0f, angle, 3.0f*angle/17.0f);
480     expect_mat(&expectedmat, &gotmat);
481
482 /*____________D3DXMatrixRotationZ______________*/
483     U(expectedmat).m[0][0] = 0.5f; U(expectedmat).m[0][1] = sqrt(3.0f)/2.0f; U(expectedmat).m[0][2] = 0.0f; U(expectedmat).m[0][3] = 0.0f;
484     U(expectedmat).m[1][0] = -sqrt(3.0f)/2.0f; U(expectedmat).m[1][1] = 0.5f; U(expectedmat).m[1][2] = 0.0f; U(expectedmat).m[1][3] = 0.0f;
485     U(expectedmat).m[2][0] = 0.0f; U(expectedmat).m[2][1] = 0.0f; U(expectedmat).m[2][2] = 1.0f; U(expectedmat).m[2][3] = 0.0f;
486     U(expectedmat).m[3][0] = 0.0f; U(expectedmat).m[3][1] = 0.0f; U(expectedmat).m[3][2] = 0.0f; U(expectedmat).m[3][3] = 1.0f;
487     D3DXMatrixRotationZ(&gotmat,angle);
488     expect_mat(&expectedmat, &gotmat);
489
490 /*____________D3DXMatrixScaling______________*/
491     U(expectedmat).m[0][0] = 0.69f; U(expectedmat).m[0][1] = 0.0f; U(expectedmat).m[0][2] = 0.0f; U(expectedmat).m[0][3] = 0.0f;
492     U(expectedmat).m[1][0] = 0.0; U(expectedmat).m[1][1] = 0.53f; U(expectedmat).m[1][2] = 0.0f; U(expectedmat).m[1][3] = 0.0f;
493     U(expectedmat).m[2][0] = 0.0f; U(expectedmat).m[2][1] = 0.0f; U(expectedmat).m[2][2] = 4.11f; U(expectedmat).m[2][3] = 0.0f;
494     U(expectedmat).m[3][0] = 0.0f; U(expectedmat).m[3][1] = 0.0f; U(expectedmat).m[3][2] = 0.0f; U(expectedmat).m[3][3] = 1.0f;
495     D3DXMatrixScaling(&gotmat,0.69f,0.53f,4.11f);
496     expect_mat(&expectedmat, &gotmat);
497
498 /*____________D3DXMatrixShadow______________*/
499     U(expectedmat).m[0][0] = 12.786773f; U(expectedmat).m[0][1] = 5.000961f; U(expectedmat).m[0][2] = 4.353778f; U(expectedmat).m[0][3] = 3.706595f;
500     U(expectedmat).m[1][0] = 1.882715; U(expectedmat).m[1][1] = 8.805615f; U(expectedmat).m[1][2] = 1.451259f; U(expectedmat).m[1][3] = 1.235532f;
501     U(expectedmat).m[2][0] = -7.530860f; U(expectedmat).m[2][1] = -6.667949f; U(expectedmat).m[2][2] = 1.333590f; U(expectedmat).m[2][3] = -4.942127f;
502     U(expectedmat).m[3][0] = -13.179006f; U(expectedmat).m[3][1] = -11.668910f; U(expectedmat).m[3][2] = -10.158816f; U(expectedmat).m[3][3] = -1.510094f;
503     D3DXMatrixShadow(&gotmat,&light,&plane);
504     expect_mat(&expectedmat, &gotmat);
505
506 /*____________D3DXMatrixTransformation______________*/
507     U(expectedmat).m[0][0] = -0.2148f; U(expectedmat).m[0][1] = 1.3116f; U(expectedmat).m[0][2] = 0.4752f; U(expectedmat).m[0][3] = 0.0f;
508     U(expectedmat).m[1][0] = 0.9504f; U(expectedmat).m[1][1] = -0.8836f; U(expectedmat).m[1][2] = 0.9244f; U(expectedmat).m[1][3] = 0.0f;
509     U(expectedmat).m[2][0] = 1.0212f; U(expectedmat).m[2][1] = 0.1936f; U(expectedmat).m[2][2] = -1.3588f; U(expectedmat).m[2][3] = 0.0f;
510     U(expectedmat).m[3][0] = 18.2985f; U(expectedmat).m[3][1] = -29.624001f; U(expectedmat).m[3][2] = 15.683499f; U(expectedmat).m[3][3] = 1.0f;
511     D3DXMatrixTransformation(&gotmat,&at,&q,NULL,&eye,&r,&last);
512     expect_mat(&expectedmat, &gotmat);
513
514 /*____________D3DXMatrixTranslation______________*/
515     U(expectedmat).m[0][0] = 1.0f; U(expectedmat).m[0][1] = 0.0f; U(expectedmat).m[0][2] = 0.0f; U(expectedmat).m[0][3] = 0.0f;
516     U(expectedmat).m[1][0] = 0.0; U(expectedmat).m[1][1] = 1.0f; U(expectedmat).m[1][2] = 0.0f; U(expectedmat).m[1][3] = 0.0f;
517     U(expectedmat).m[2][0] = 0.0f; U(expectedmat).m[2][1] = 0.0f; U(expectedmat).m[2][2] = 1.0f; U(expectedmat).m[2][3] = 0.0f;
518     U(expectedmat).m[3][0] = 0.69f; U(expectedmat).m[3][1] = 0.53f; U(expectedmat).m[3][2] = 4.11f; U(expectedmat).m[3][3] = 1.0f;
519     D3DXMatrixTranslation(&gotmat,0.69f,0.53f,4.11f);
520     expect_mat(&expectedmat, &gotmat);
521
522 /*____________D3DXMatrixTranspose______________*/
523     U(expectedmat).m[0][0] = 10.0f; U(expectedmat).m[0][1] = 11.0f; U(expectedmat).m[0][2] = 19.0f; U(expectedmat).m[0][3] = 2.0f;
524     U(expectedmat).m[1][0] = 5.0; U(expectedmat).m[1][1] = 20.0f; U(expectedmat).m[1][2] = -21.0f; U(expectedmat).m[1][3] = 3.0f;
525     U(expectedmat).m[2][0] = 7.0f; U(expectedmat).m[2][1] = 16.0f; U(expectedmat).m[2][2] = 30.f; U(expectedmat).m[2][3] = -4.0f;
526     U(expectedmat).m[3][0] = 8.0f; U(expectedmat).m[3][1] = 33.0f; U(expectedmat).m[3][2] = 43.0f; U(expectedmat).m[3][3] = -40.0f;
527     D3DXMatrixTranspose(&gotmat,&mat);
528     expect_mat(&expectedmat, &gotmat);
529 }
530
531 static void D3DXPlaneTest(void)
532 {
533     D3DXMATRIX mat;
534     D3DXPLANE expectedplane, gotplane, nulplane, plane;
535     D3DXVECTOR3 expectedvec, gotvec, vec1, vec2, vec3;
536     LPD3DXVECTOR3 funcpointer;
537     D3DXVECTOR4 vec;
538     FLOAT expected, got;
539
540     U(mat).m[0][1] = 5.0f; U(mat).m[0][2] = 7.0f; U(mat).m[0][3] = 8.0f;
541     U(mat).m[1][0] = 11.0f; U(mat).m[1][2] = 16.0f; U(mat).m[1][3] = 33.0f;
542     U(mat).m[2][0] = 19.0f; U(mat).m[2][1] = -21.0f; U(mat).m[2][3] = 43.0f;
543     U(mat).m[3][0] = 2.0f; U(mat).m[3][1] = 3.0f; U(mat).m[3][2] = -4.0f;
544     U(mat).m[0][0] = 10.0f; U(mat).m[1][1] = 20.0f; U(mat).m[2][2] = 30.0f;
545     U(mat).m[3][3] = -40.0f;
546
547     plane.a = -3.0f; plane.b = -1.0f; plane.c = 4.0f; plane.d = 7.0f;
548
549     vec.x = 2.0f; vec.y = 5.0f; vec.z = -6.0f; vec.w = 11.0f;
550
551 /*_______________D3DXPlaneDot________________*/
552     expected = 42.0f;
553     got = D3DXPlaneDot(&plane,&vec),
554     ok( expected == got, "Expected : %f, Got : %f\n",expected, got);
555     expected = 0.0f;
556     got = D3DXPlaneDot(NULL,&vec),
557     ok( expected == got, "Expected : %f, Got : %f\n",expected, got);
558     expected = 0.0f;
559     got = D3DXPlaneDot(NULL,NULL),
560     ok( expected == got, "Expected : %f, Got : %f\n",expected, got);
561
562 /*_______________D3DXPlaneDotCoord________________*/
563     expected = -28.0f;
564     got = D3DXPlaneDotCoord(&plane,&vec),
565     ok( expected == got, "Expected : %f, Got : %f\n",expected, got);
566     expected = 0.0f;
567     got = D3DXPlaneDotCoord(NULL,&vec),
568     ok( expected == got, "Expected : %f, Got : %f\n",expected, got);
569     expected = 0.0f;
570     got = D3DXPlaneDotCoord(NULL,NULL),
571     ok( expected == got, "Expected : %f, Got : %f\n",expected, got);
572
573 /*_______________D3DXPlaneDotNormal______________*/
574     expected = -35.0f;
575     got = D3DXPlaneDotNormal(&plane,&vec),
576     ok( expected == got, "Expected : %f, Got : %f\n",expected, got);
577     expected = 0.0f;
578     got = D3DXPlaneDotNormal(NULL,&vec),
579     ok( expected == got, "Expected : %f, Got : %f\n",expected, got);
580     expected = 0.0f;
581     got = D3DXPlaneDotNormal(NULL,NULL),
582     ok( expected == got, "Expected : %f, Got : %f\n",expected, got);
583
584 /*_______________D3DXPlaneFromPointNormal_______*/
585     vec1.x = 11.0f; vec1.y = 13.0f; vec1.z = 15.0f;
586     vec2.x = 17.0f; vec2.y = 31.0f; vec2.z = 24.0f;
587     expectedplane.a = 17.0f; expectedplane.b = 31.0f; expectedplane.c = 24.0f; expectedplane.d = -950.0f;
588     D3DXPlaneFromPointNormal(&gotplane,&vec1,&vec2);
589     expect_plane(expectedplane, gotplane);
590
591 /*_______________D3DXPlaneFromPoints_______*/
592     vec1.x = 1.0f; vec1.y = 2.0f; vec1.z = 3.0f;
593     vec2.x = 1.0f; vec2.y = -6.0f; vec2.z = -5.0f;
594     vec3.x = 83.0f; vec3.y = 74.0f; vec3.z = 65.0f;
595     expectedplane.a = 0.085914f; expectedplane.b = -0.704492f; expectedplane.c = 0.704492f; expectedplane.d = -0.790406f;
596     D3DXPlaneFromPoints(&gotplane,&vec1,&vec2,&vec3);
597     expect_plane(expectedplane, gotplane);
598
599 /*_______________D3DXPlaneIntersectLine___________*/
600     vec1.x = 9.0f; vec1.y = 6.0f; vec1.z = 3.0f;
601     vec2.x = 2.0f; vec2.y = 5.0f; vec2.z = 8.0f;
602     expectedvec.x = 20.0f/3.0f; expectedvec.y = 17.0f/3.0f; expectedvec.z = 14.0f/3.0f;
603     D3DXPlaneIntersectLine(&gotvec,&plane,&vec1,&vec2);
604     expect_vec3(expectedvec, gotvec);
605     /* Test a parallel line */
606     vec1.x = 11.0f; vec1.y = 13.0f; vec1.z = 15.0f;
607     vec2.x = 17.0f; vec2.y = 31.0f; vec2.z = 24.0f;
608     expectedvec.x = 20.0f/3.0f; expectedvec.y = 17.0f/3.0f; expectedvec.z = 14.0f/3.0f;
609     funcpointer = D3DXPlaneIntersectLine(&gotvec,&plane,&vec1,&vec2);
610     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
611
612 /*_______________D3DXPlaneNormalize______________*/
613     expectedplane.a = -3.0f/sqrt(26.0f); expectedplane.b = -1.0f/sqrt(26.0f); expectedplane.c = 4.0f/sqrt(26.0f); expectedplane.d = 7.0/sqrt(26.0f);
614     D3DXPlaneNormalize(&gotplane, &plane);
615     expect_plane(expectedplane, gotplane);
616     nulplane.a = 0.0; nulplane.b = 0.0f, nulplane.c = 0.0f; nulplane.d = 0.0f;
617     expectedplane.a = 0.0f; expectedplane.b = 0.0f; expectedplane.c = 0.0f; expectedplane.d = 0.0f;
618     D3DXPlaneNormalize(&gotplane, &nulplane);
619     expect_plane(expectedplane, gotplane);
620     nulplane.a = 0.0; nulplane.b = 0.0f, nulplane.c = 0.0f; nulplane.d = 4.3f;
621     expectedplane.a = 0.0f; expectedplane.b = 0.0f; expectedplane.c = 0.0f; expectedplane.d = 0.0f;
622     D3DXPlaneNormalize(&gotplane, &nulplane);
623     expect_plane(expectedplane, gotplane);
624
625 /*_______________D3DXPlaneTransform____________*/
626     expectedplane.a = 49.0f; expectedplane.b = -98.0f; expectedplane.c = 55.0f; expectedplane.d = -165.0f;
627     D3DXPlaneTransform(&gotplane,&plane,&mat);
628     expect_plane(expectedplane, gotplane);
629 }
630
631 static void D3DXQuaternionTest(void)
632 {
633     D3DXMATRIX mat;
634     D3DXQUATERNION expectedquat, gotquat, Nq, Nq1, nul, smallq, smallr, q, r, s, t, u;
635     LPD3DXQUATERNION funcpointer;
636     D3DXVECTOR3 axis, expectedvec;
637     FLOAT angle, expected, got, scale, scale2;
638     BOOL expectedbool, gotbool;
639
640     nul.x = 0.0f; nul.y = 0.0f; nul.z = 0.0f; nul.w = 0.0f;
641     q.x = 1.0f, q.y = 2.0f; q.z = 4.0f; q.w = 10.0f;
642     r.x = -3.0f; r.y = 4.0f; r.z = -5.0f; r.w = 7.0;
643     t.x = -1111.0f, t.y = 111.0f; t.z = -11.0f; t.w = 1.0f;
644     u.x = 91.0f; u.y = - 82.0f; u.z = 7.3f; u.w = -6.4f;
645     smallq.x = 0.1f; smallq.y = 0.2f; smallq.z= 0.3f; smallq.w = 0.4f;
646     smallr.x = 0.5f; smallr.y = 0.6f; smallr.z= 0.7f; smallr.w = 0.8f;
647
648     scale = 0.3f;
649     scale2 = 0.78f;
650
651 /*_______________D3DXQuaternionBaryCentric________________________*/
652     expectedquat.x = -867.444458; expectedquat.y = 87.851111f; expectedquat.z = -9.937778f; expectedquat.w = 3.235555f;
653     D3DXQuaternionBaryCentric(&gotquat,&q,&r,&t,scale,scale2);
654     expect_vec4(expectedquat,gotquat);
655
656 /*_______________D3DXQuaternionConjugate________________*/
657     expectedquat.x = -1.0f; expectedquat.y = -2.0f; expectedquat.z = -4.0f; expectedquat.w = 10.0f;
658     D3DXQuaternionConjugate(&gotquat,&q);
659     expect_vec4(expectedquat,gotquat);
660     /* Test the NULL case */
661     funcpointer = D3DXQuaternionConjugate(&gotquat,NULL);
662     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
663     funcpointer = D3DXQuaternionConjugate(NULL,NULL);
664     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
665
666 /*_______________D3DXQuaternionDot______________________*/
667     expected = 55.0f;
668     got = D3DXQuaternionDot(&q,&r);
669     ok(relative_error(got, expected ) < admitted_error, "Expected: %f, Got: %f\n", expected, got);
670     /* Tests the case NULL */
671     expected=0.0f;
672     got = D3DXQuaternionDot(NULL,&r);
673     ok(relative_error(got, expected ) < admitted_error, "Expected: %f, Got: %f\n", expected, got);
674     expected=0.0f;
675     got = D3DXQuaternionDot(NULL,NULL);
676     ok(relative_error(got, expected ) < admitted_error, "Expected: %f, Got: %f\n", expected, got);
677
678 /*_______________D3DXQuaternionExp______________________________*/
679     expectedquat.x = -0.216382f; expectedquat.y = -0.432764f; expectedquat.z = -0.8655270f; expectedquat.w = -0.129449f;
680     D3DXQuaternionExp(&gotquat,&q);
681     expect_vec4(expectedquat,gotquat);
682     /* Test the null quaternion */
683     expectedquat.x = 0.0f; expectedquat.y = 0.0f; expectedquat.z = 0.0f; expectedquat.w = 1.0f;
684     D3DXQuaternionExp(&gotquat,&nul);
685     expect_vec4(expectedquat,gotquat);
686     /* Test the case where the norm of the quaternion is <1 */
687     Nq1.x = 0.2f; Nq1.y = 0.1f; Nq1.z = 0.3; Nq1.w= 0.9f;
688     expectedquat.x = 0.195366; expectedquat.y = 0.097683f; expectedquat.z = 0.293049f; expectedquat.w = 0.930813f;
689     D3DXQuaternionExp(&gotquat,&Nq1);
690     expect_vec4(expectedquat,gotquat);
691
692 /*_______________D3DXQuaternionIdentity________________*/
693     expectedquat.x = 0.0f; expectedquat.y = 0.0f; expectedquat.z = 0.0f; expectedquat.w = 1.0f;
694     D3DXQuaternionIdentity(&gotquat);
695     expect_vec4(expectedquat,gotquat);
696     /* Test the NULL case */
697     funcpointer = D3DXQuaternionIdentity(NULL);
698     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
699
700 /*_______________D3DXQuaternionInverse________________________*/
701     expectedquat.x = -1.0f/121.0f; expectedquat.y = -2.0f/121.0f; expectedquat.z = -4.0f/121.0f; expectedquat.w = 10.0f/121.0f;
702     D3DXQuaternionInverse(&gotquat,&q);
703     expect_vec4(expectedquat,gotquat);
704
705     expectedquat.x = 1.0f; expectedquat.y = 2.0f; expectedquat.z = 4.0f; expectedquat.w = 10.0f;
706     D3DXQuaternionInverse(&gotquat,&gotquat);
707     expect_vec4(expectedquat,gotquat);
708
709
710 /*_______________D3DXQuaternionIsIdentity________________*/
711     s.x = 0.0f; s.y = 0.0f; s.z = 0.0f; s.w = 1.0f;
712     expectedbool = TRUE;
713     gotbool = D3DXQuaternionIsIdentity(&s);
714     ok( expectedbool == gotbool, "Expected boolean : %d, Got bool : %d\n", expectedbool, gotbool);
715     s.x = 2.3f; s.y = -4.2f; s.z = 1.2f; s.w=0.2f;
716     expectedbool = FALSE;
717     gotbool = D3DXQuaternionIsIdentity(&q);
718     ok( expectedbool == gotbool, "Expected boolean : %d, Got bool : %d\n", expectedbool, gotbool);
719     /* Test the NULL case */
720     gotbool = D3DXQuaternionIsIdentity(NULL);
721     ok(gotbool == FALSE, "Expected boolean: %d, Got boolean: %d\n", FALSE, gotbool);
722
723 /*_______________D3DXQuaternionLength__________________________*/
724    expected = 11.0f;
725    got = D3DXQuaternionLength(&q);
726    ok(relative_error(got, expected ) < admitted_error, "Expected: %f, Got: %f\n", expected, got);
727    /* Tests the case NULL */
728     expected=0.0f;
729     got = D3DXQuaternionLength(NULL);
730     ok(relative_error(got, expected ) < admitted_error, "Expected: %f, Got: %f\n", expected, got);
731
732 /*_______________D3DXQuaternionLengthSq________________________*/
733     expected = 121.0f;
734     got = D3DXQuaternionLengthSq(&q);
735     ok(relative_error(got, expected ) < admitted_error, "Expected: %f, Got: %f\n", expected, got);
736     /* Tests the case NULL */
737     expected=0.0f;
738     got = D3DXQuaternionLengthSq(NULL);
739     ok(relative_error(got, expected ) < admitted_error, "Expected: %f, Got: %f\n", expected, got);
740
741 /*_______________D3DXQuaternionLn______________________________*/
742     expectedquat.x = 1.0f; expectedquat.y = 2.0f; expectedquat.z = 4.0f; expectedquat.w = 0.0f;
743     D3DXQuaternionLn(&gotquat,&q);
744     expect_vec4(expectedquat,gotquat);
745     expectedquat.x = -3.0f; expectedquat.y = 4.0f; expectedquat.z = -5.0f; expectedquat.w = 0.0f;
746     D3DXQuaternionLn(&gotquat,&r);
747     expect_vec4(expectedquat,gotquat);
748     Nq.x = 1.0f/11.0f; Nq.y = 2.0f/11.0f; Nq.z = 4.0f/11.0f; Nq.w=10.0f/11.0f;
749     expectedquat.x = 0.093768f; expectedquat.y = 0.187536f; expectedquat.z = 0.375073f; expectedquat.w = 0.0f;
750     D3DXQuaternionLn(&gotquat,&Nq);
751     expect_vec4(expectedquat,gotquat);
752     /* Test the cas where the norm of the quaternion is <1 */
753     Nq1.x = 0.2f; Nq1.y = 0.1f; Nq1.z = 0.3; Nq1.w= 0.9f;
754     expectedquat.x = 0.206945f; expectedquat.y = 0.103473f; expectedquat.z = 0.310418f; expectedquat.w = 0.0f;
755     D3DXQuaternionLn(&gotquat,&Nq1);
756     todo_wine{ expect_vec4(expectedquat,gotquat) };
757
758 /*_______________D3DXQuaternionMultiply________________________*/
759     expectedquat.x = 3.0f; expectedquat.y = 61.0f; expectedquat.z = -32.0f; expectedquat.w = 85.0f;
760     D3DXQuaternionMultiply(&gotquat,&q,&r);
761     expect_vec4(expectedquat,gotquat);
762
763 /*_______________D3DXQuaternionNormalize________________________*/
764     expectedquat.x = 1.0f/11.0f; expectedquat.y = 2.0f/11.0f; expectedquat.z = 4.0f/11.0f; expectedquat.w = 10.0f/11.0f;
765     D3DXQuaternionNormalize(&gotquat,&q);
766     expect_vec4(expectedquat,gotquat);
767
768 /*_______________D3DXQuaternionRotationAxis___________________*/
769     axis.x = 2.0f; axis.y = 7.0; axis.z = 13.0f;
770     angle = D3DX_PI/3.0f;
771     expectedquat.x = 0.067116; expectedquat.y = 0.234905f; expectedquat.z = 0.436251f; expectedquat.w = 0.866025f;
772     D3DXQuaternionRotationAxis(&gotquat,&axis,angle);
773     expect_vec4(expectedquat,gotquat);
774  /* Test the nul quaternion */
775     axis.x = 0.0f; axis.y = 0.0; axis.z = 0.0f;
776     expectedquat.x = 0.0f; expectedquat.y = 0.0f; expectedquat.z = 0.0f; expectedquat.w = 0.866025f;
777     D3DXQuaternionRotationAxis(&gotquat,&axis,angle);
778     expect_vec4(expectedquat,gotquat);
779
780 /*_______________D3DXQuaternionRotationMatrix___________________*/
781     /* test when the trace is >0 */
782     U(mat).m[0][1] = 5.0f; U(mat).m[0][2] = 7.0f; U(mat).m[0][3] = 8.0f;
783     U(mat).m[1][0] = 11.0f; U(mat).m[1][2] = 16.0f; U(mat).m[1][3] = 33.0f;
784     U(mat).m[2][0] = 19.0f; U(mat).m[2][1] = -21.0f; U(mat).m[2][3] = 43.0f;
785     U(mat).m[3][0] = 2.0f; U(mat).m[3][1] = 3.0f; U(mat).m[3][2] = -4.0f;
786     U(mat).m[0][0] = 10.0f; U(mat).m[1][1] = 20.0f; U(mat).m[2][2] = 30.0f;
787     U(mat).m[3][3] = 48.0f;
788     expectedquat.x = 2.368682f; expectedquat.y = 0.768221f; expectedquat.z = -0.384111f; expectedquat.w = 3.905125f;
789     D3DXQuaternionRotationMatrix(&gotquat,&mat);
790     expect_vec4(expectedquat,gotquat);
791     /* test the case when the greater element is (2,2) */
792     U(mat).m[0][1] = 5.0f; U(mat).m[0][2] = 7.0f; U(mat).m[0][3] = 8.0f;
793     U(mat).m[1][0] = 11.0f; U(mat).m[1][2] = 16.0f; U(mat).m[1][3] = 33.0f;
794     U(mat).m[2][0] = 19.0f; U(mat).m[2][1] = -21.0f; U(mat).m[2][3] = 43.0f;
795     U(mat).m[3][0] = 2.0f; U(mat).m[3][1] = 3.0f; U(mat).m[3][2] = -4.0f;
796     U(mat).m[0][0] = -10.0f; U(mat).m[1][1] = -60.0f; U(mat).m[2][2] = 40.0f;
797     U(mat).m[3][3] = 48.0f;
798     expectedquat.x = 1.233905f; expectedquat.y = -0.237290f; expectedquat.z = 5.267827f; expectedquat.w = -0.284747f;
799     D3DXQuaternionRotationMatrix(&gotquat,&mat);
800     expect_vec4(expectedquat,gotquat);
801     /* test the case when the greater element is (1,1) */
802     U(mat).m[0][1] = 5.0f; U(mat).m[0][2] = 7.0f; U(mat).m[0][3] = 8.0f;
803     U(mat).m[1][0] = 11.0f; U(mat).m[1][2] = 16.0f; U(mat).m[1][3] = 33.0f;
804     U(mat).m[2][0] = 19.0f; U(mat).m[2][1] = -21.0f; U(mat).m[2][3] = 43.0f;
805     U(mat).m[3][0] = 2.0f; U(mat).m[3][1] = 3.0f; U(mat).m[3][2] = -4.0f;
806     U(mat).m[0][0] = -10.0f; U(mat).m[1][1] = 60.0f; U(mat).m[2][2] = -80.0f;
807     U(mat).m[3][3] = 48.0f;
808     expectedquat.x = 0.651031f; expectedquat.y = 6.144103f; expectedquat.z = -0.203447f; expectedquat.w = 0.488273f;
809     D3DXQuaternionRotationMatrix(&gotquat,&mat);
810     expect_vec4(expectedquat,gotquat);
811     /* test the case when the trace is near 0 in a matrix which is not a rotation */
812     U(mat).m[0][1] = 5.0f; U(mat).m[0][2] = 7.0f; U(mat).m[0][3] = 8.0f;
813     U(mat).m[1][0] = 11.0f; U(mat).m[1][2] = 16.0f; U(mat).m[1][3] = 33.0f;
814     U(mat).m[2][0] = 19.0f; U(mat).m[2][1] = -21.0f; U(mat).m[2][3] = 43.0f;
815     U(mat).m[3][0] = 2.0f; U(mat).m[3][1] = 3.0f; U(mat).m[3][2] = -4.0f;
816     U(mat).m[0][0] = -10.0f; U(mat).m[1][1] = 10.0f; U(mat).m[2][2] = -0.9f;
817     U(mat).m[3][3] = 48.0f;
818     expectedquat.x = 1.709495f; expectedquat.y = 2.339872f; expectedquat.z = -0.534217f; expectedquat.w = 1.282122f;
819     D3DXQuaternionRotationMatrix(&gotquat,&mat);
820     expect_vec4(expectedquat,gotquat);
821     /* test the case when the trace is 0.49 in a matrix which is not a rotation */
822     U(mat).m[0][1] = 5.0f; U(mat).m[0][2] = 7.0f; U(mat).m[0][3] = 8.0f;
823     U(mat).m[1][0] = 11.0f; U(mat).m[1][2] = 16.0f; U(mat).m[1][3] = 33.0f;
824     U(mat).m[2][0] = 19.0f; U(mat).m[2][1] = -21.0f; U(mat).m[2][3] = 43.0f;
825     U(mat).m[3][0] = 2.0f; U(mat).m[3][1] = 3.0f; U(mat).m[3][2] = -4.0f;
826     U(mat).m[0][0] = -10.0f; U(mat).m[1][1] = 10.0f; U(mat).m[2][2] = -0.51f;
827     U(mat).m[3][3] = 48.0f;
828     expectedquat.x = 1.724923f; expectedquat.y = 2.318944f; expectedquat.z = -0.539039f; expectedquat.w = 1.293692f;
829     D3DXQuaternionRotationMatrix(&gotquat,&mat);
830     expect_vec4(expectedquat,gotquat);
831     /* test the case when the trace is 0.51 in a matrix which is not a rotation */
832     U(mat).m[0][1] = 5.0f; U(mat).m[0][2] = 7.0f; U(mat).m[0][3] = 8.0f;
833     U(mat).m[1][0] = 11.0f; U(mat).m[1][2] = 16.0f; U(mat).m[1][3] = 33.0f;
834     U(mat).m[2][0] = 19.0f; U(mat).m[2][1] = -21.0f; U(mat).m[2][3] = 43.0f;
835     U(mat).m[3][0] = 2.0f; U(mat).m[3][1] = 3.0f; U(mat).m[3][2] = -4.0f;
836     U(mat).m[0][0] = -10.0f; U(mat).m[1][1] = 10.0f; U(mat).m[2][2] = -0.49f;
837     U(mat).m[3][3] = 48.0f;
838     expectedquat.x = 1.725726f; expectedquat.y = 2.317865f; expectedquat.z = -0.539289f; expectedquat.w = 1.294294f;
839     D3DXQuaternionRotationMatrix(&gotquat,&mat);
840     expect_vec4(expectedquat,gotquat);
841     /* test the case when the trace is 0.99 in a matrix which is not a rotation */
842     U(mat).m[0][1] = 5.0f; U(mat).m[0][2] = 7.0f; U(mat).m[0][3] = 8.0f;
843     U(mat).m[1][0] = 11.0f; U(mat).m[1][2] = 16.0f; U(mat).m[1][3] = 33.0f;
844     U(mat).m[2][0] = 19.0f; U(mat).m[2][1] = -21.0f; U(mat).m[2][3] = 43.0f;
845     U(mat).m[3][0] = 2.0f; U(mat).m[3][1] = 3.0f; U(mat).m[3][2] = -4.0f;
846     U(mat).m[0][0] = -10.0f; U(mat).m[1][1] = 10.0f; U(mat).m[2][2] = -0.01f;
847     U(mat).m[3][3] = 48.0f;
848     expectedquat.x = 1.745328f; expectedquat.y = 2.291833f; expectedquat.z = -0.545415f; expectedquat.w = 1.308996f;
849     D3DXQuaternionRotationMatrix(&gotquat,&mat);
850     expect_vec4(expectedquat,gotquat);
851     /* test the case when the trace is 1.0 in a matrix which is not a rotation */
852     U(mat).m[0][1] = 5.0f; U(mat).m[0][2] = 7.0f; U(mat).m[0][3] = 8.0f;
853     U(mat).m[1][0] = 11.0f; U(mat).m[1][2] = 16.0f; U(mat).m[1][3] = 33.0f;
854     U(mat).m[2][0] = 19.0f; U(mat).m[2][1] = -21.0f; U(mat).m[2][3] = 43.0f;
855     U(mat).m[3][0] = 2.0f; U(mat).m[3][1] = 3.0f; U(mat).m[3][2] = -4.0f;
856     U(mat).m[0][0] = -10.0f; U(mat).m[1][1] = 10.0f; U(mat).m[2][2] = 0.0f;
857     U(mat).m[3][3] = 48.0f;
858     expectedquat.x = 1.745743f; expectedquat.y = 2.291288f; expectedquat.z = -0.545545f; expectedquat.w = 1.309307f;
859     D3DXQuaternionRotationMatrix(&gotquat,&mat);
860     expect_vec4(expectedquat,gotquat);
861     /* test the case when the trace is 1.01 in a matrix which is not a rotation */
862     U(mat).m[0][1] = 5.0f; U(mat).m[0][2] = 7.0f; U(mat).m[0][3] = 8.0f;
863     U(mat).m[1][0] = 11.0f; U(mat).m[1][2] = 16.0f; U(mat).m[1][3] = 33.0f;
864     U(mat).m[2][0] = 19.0f; U(mat).m[2][1] = -21.0f; U(mat).m[2][3] = 43.0f;
865     U(mat).m[3][0] = 2.0f; U(mat).m[3][1] = 3.0f; U(mat).m[3][2] = -4.0f;
866     U(mat).m[0][0] = -10.0f; U(mat).m[1][1] = 10.0f; U(mat).m[2][2] = 0.01f;
867     U(mat).m[3][3] = 48.0f;
868     expectedquat.x = 18.408188f; expectedquat.y = 5.970223f; expectedquat.z = -2.985111f; expectedquat.w = 0.502494f;
869     D3DXQuaternionRotationMatrix(&gotquat,&mat);
870     expect_vec4(expectedquat,gotquat);
871     /* test the case when the trace is 1.5 in a matrix which is not a rotation */
872     U(mat).m[0][1] = 5.0f; U(mat).m[0][2] = 7.0f; U(mat).m[0][3] = 8.0f;
873     U(mat).m[1][0] = 11.0f; U(mat).m[1][2] = 16.0f; U(mat).m[1][3] = 33.0f;
874     U(mat).m[2][0] = 19.0f; U(mat).m[2][1] = -21.0f; U(mat).m[2][3] = 43.0f;
875     U(mat).m[3][0] = 2.0f; U(mat).m[3][1] = 3.0f; U(mat).m[3][2] = -4.0f;
876     U(mat).m[0][0] = -10.0f; U(mat).m[1][1] = 10.0f; U(mat).m[2][2] = 0.5f;
877     U(mat).m[3][3] = 48.0f;
878     expectedquat.x = 15.105186f; expectedquat.y = 4.898980f; expectedquat.z = -2.449490f; expectedquat.w = 0.612372f;
879     D3DXQuaternionRotationMatrix(&gotquat,&mat);
880     expect_vec4(expectedquat,gotquat);
881     /* test the case when the trace is 1.7 in a matrix which is not a rotation */
882     U(mat).m[0][1] = 5.0f; U(mat).m[0][2] = 7.0f; U(mat).m[0][3] = 8.0f;
883     U(mat).m[1][0] = 11.0f; U(mat).m[1][2] = 16.0f; U(mat).m[1][3] = 33.0f;
884     U(mat).m[2][0] = 19.0f; U(mat).m[2][1] = -21.0f; U(mat).m[2][3] = 43.0f;
885     U(mat).m[3][0] = 2.0f; U(mat).m[3][1] = 3.0f; U(mat).m[3][2] = -4.0f;
886     U(mat).m[0][0] = -10.0f; U(mat).m[1][1] = 10.0f; U(mat).m[2][2] = 0.70f;
887     U(mat).m[3][3] = 48.0f;
888     expectedquat.x = 14.188852f; expectedquat.y = 4.601790f; expectedquat.z = -2.300895f; expectedquat.w = 0.651920f;
889     D3DXQuaternionRotationMatrix(&gotquat,&mat);
890     expect_vec4(expectedquat,gotquat);
891     /* test the case when the trace is 1.99 in a matrix which is not a rotation */
892     U(mat).m[0][1] = 5.0f; U(mat).m[0][2] = 7.0f; U(mat).m[0][3] = 8.0f;
893     U(mat).m[1][0] = 11.0f; U(mat).m[1][2] = 16.0f; U(mat).m[1][3] = 33.0f;
894     U(mat).m[2][0] = 19.0f; U(mat).m[2][1] = -21.0f; U(mat).m[2][3] = 43.0f;
895     U(mat).m[3][0] = 2.0f; U(mat).m[3][1] = 3.0f; U(mat).m[3][2] = -4.0f;
896     U(mat).m[0][0] = -10.0f; U(mat).m[1][1] = 10.0f; U(mat).m[2][2] = 0.99f;
897     U(mat).m[3][3] = 48.0f;
898     expectedquat.x = 13.114303f; expectedquat.y = 4.253287f; expectedquat.z = -2.126644f; expectedquat.w = 0.705337f;
899     D3DXQuaternionRotationMatrix(&gotquat,&mat);
900     expect_vec4(expectedquat,gotquat);
901     /* test the case when the trace is 2.0 in a matrix which is not a rotation */
902     U(mat).m[0][1] = 5.0f; U(mat).m[0][2] = 7.0f; U(mat).m[0][3] = 8.0f;
903     U(mat).m[1][0] = 11.0f; U(mat).m[1][2] = 16.0f; U(mat).m[1][3] = 33.0f;
904     U(mat).m[2][0] = 19.0f; U(mat).m[2][1] = -21.0f; U(mat).m[2][3] = 43.0f;
905     U(mat).m[3][0] = 2.0f; U(mat).m[3][1] = 3.0f; U(mat).m[3][2] = -4.0f;
906     U(mat).m[0][0] = -10.0f; U(mat).m[1][1] = 10.0f; U(mat).m[2][2] = 2.0f;
907     U(mat).m[3][3] = 48.0f;
908     expectedquat.x = 10.680980f; expectedquat.y = 3.464102f; expectedquat.z = -1.732051f; expectedquat.w = 0.866025f;
909     D3DXQuaternionRotationMatrix(&gotquat,&mat);
910     expect_vec4(expectedquat,gotquat);
911
912 /*_______________D3DXQuaternionRotationYawPitchRoll__________*/
913     expectedquat.x = 0.303261f; expectedquat.y = 0.262299f; expectedquat.z = 0.410073f; expectedquat.w = 0.819190f;
914     D3DXQuaternionRotationYawPitchRoll(&gotquat,D3DX_PI/4.0f,D3DX_PI/11.0f,D3DX_PI/3.0f);
915     expect_vec4(expectedquat,gotquat);
916
917 /*_______________D3DXQuaternionSlerp________________________*/
918     expectedquat.x = -0.2f; expectedquat.y = 2.6f; expectedquat.z = 1.3f; expectedquat.w = 9.1f;
919     D3DXQuaternionSlerp(&gotquat,&q,&r,scale);
920     expect_vec4(expectedquat,gotquat);
921     expectedquat.x = 334.0f; expectedquat.y = -31.9f; expectedquat.z = 6.1f; expectedquat.w = 6.7f;
922     D3DXQuaternionSlerp(&gotquat,&q,&t,scale);
923     expect_vec4(expectedquat,gotquat);
924     expectedquat.x = 0.239485f; expectedquat.y = 0.346580f; expectedquat.z = 0.453676f; expectedquat.w = 0.560772f;
925     D3DXQuaternionSlerp(&gotquat,&smallq,&smallr,scale);
926     expect_vec4(expectedquat,gotquat);
927
928 /*_______________D3DXQuaternionSquad________________________*/
929     expectedquat.x = -156.296f; expectedquat.y = 30.242f; expectedquat.z = -2.5022f; expectedquat.w = 7.3576f;
930     D3DXQuaternionSquad(&gotquat,&q,&r,&t,&u,scale);
931     expect_vec4(expectedquat,gotquat);
932
933 /*_______________D3DXQuaternionToAxisAngle__________________*/
934     Nq.x = 1.0f/22.0f; Nq.y = 2.0f/22.0f; Nq.z = 4.0f/22.0f; Nq.w = 10.0f/22.0f;
935     expectedvec.x = 1.0f/22.0f; expectedvec.y = 2.0f/22.0f; expectedvec.z = 4.0f/22.0f;
936     expected = 2.197869f;
937     D3DXQuaternionToAxisAngle(&Nq,&axis,&angle);
938     expect_vec3(expectedvec,axis);
939     ok(relative_error(angle,  expected ) < admitted_error, "Expected: %f, Got: %f\n", expected, angle);
940     /* Test if |w|>1.0f */
941     expectedvec.x = 1.0f; expectedvec.y = 2.0f; expectedvec.z = 4.0f;
942     expected = 0.0f;
943     D3DXQuaternionToAxisAngle(&q,&axis,&angle);
944     expect_vec3(expectedvec,axis);
945     /* Test the null quaternion */
946     expectedvec.x = 0.0f; expectedvec.y = 0.0f; expectedvec.z = 0.0f;
947     expected = 3.141593f;
948     D3DXQuaternionToAxisAngle(&nul,&axis,&angle);
949     expect_vec3(expectedvec,axis);
950     ok(relative_error(angle, expected ) < admitted_error, "Expected: %f, Got: %f\n", expected, angle);
951 }
952
953 static void D3DXVector2Test(void)
954 {
955     D3DXVECTOR2 expectedvec, gotvec, nul, u, v, w, x;
956     LPD3DXVECTOR2 funcpointer;
957     D3DXVECTOR4 expectedtrans, gottrans;
958     D3DXMATRIX mat;
959     FLOAT coeff1, coeff2, expected, got, scale;
960
961     nul.x = 0.0f; nul.y = 0.0f;
962     u.x = 3.0f; u.y = 4.0f;
963     v.x = -7.0f; v.y = 9.0f;
964     w.x = 4.0f; w.y = -3.0f;
965     x.x = 2.0f; x.y = -11.0f;
966
967     U(mat).m[0][0] = 1.0f; U(mat).m[0][1] = 2.0f; U(mat).m[0][2] = 3.0f; U(mat).m[0][3] = 4.0f;
968     U(mat).m[1][0] = 5.0f; U(mat).m[1][1] = 6.0f; U(mat).m[1][2] = 7.0f; U(mat).m[1][3] = 8.0f;
969     U(mat).m[2][0] = 9.0f; U(mat).m[2][1] = 10.0f; U(mat).m[2][2] = 11.0f; U(mat).m[2][3] = 12.0f;
970     U(mat).m[3][0] = 13.0f; U(mat).m[3][1] = 14.0f; U(mat).m[3][2] = 15.0f; U(mat).m[3][3] = 16.0f;
971
972     coeff1 = 2.0f; coeff2 = 5.0f;
973     scale = -6.5f;
974
975 /*_______________D3DXVec2Add__________________________*/
976     expectedvec.x = -4.0f; expectedvec.y = 13.0f;
977     D3DXVec2Add(&gotvec,&u,&v);
978     expect_vec(expectedvec,gotvec);
979     /* Tests the case NULL */
980     funcpointer = D3DXVec2Add(&gotvec,NULL,&v);
981     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
982     funcpointer = D3DXVec2Add(NULL,NULL,NULL);
983     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
984
985 /*_______________D3DXVec2BaryCentric___________________*/
986     expectedvec.x = -12.0f; expectedvec.y = -21.0f;
987     D3DXVec2BaryCentric(&gotvec,&u,&v,&w,coeff1,coeff2);
988     expect_vec(expectedvec,gotvec);
989
990 /*_______________D3DXVec2CatmullRom____________________*/
991     expectedvec.x = 5820.25f; expectedvec.y = -3654.5625f;
992     D3DXVec2CatmullRom(&gotvec,&u,&v,&w,&x,scale);
993     expect_vec(expectedvec,gotvec);
994
995 /*_______________D3DXVec2CCW__________________________*/
996    expected = 55.0f;
997    got = D3DXVec2CCW(&u,&v);
998    ok(relative_error(got, expected ) < admitted_error, "Expected: %f, Got: %f\n", expected, got);
999    /* Tests the case NULL */
1000     expected=0.0f;
1001     got = D3DXVec2CCW(NULL,&v);
1002     ok(relative_error(got, expected ) < admitted_error, "Expected: %f, Got: %f\n", expected, got);
1003     expected=0.0f;
1004     got = D3DXVec2CCW(NULL,NULL);
1005     ok(relative_error(got, expected ) < admitted_error, "Expected: %f, Got: %f\n", expected, got);
1006
1007 /*_______________D3DXVec2Dot__________________________*/
1008     expected = 15.0f;
1009     got = D3DXVec2Dot(&u,&v);
1010     ok(relative_error(got,  expected ) < admitted_error, "Expected: %f, Got: %f\n", expected, got);
1011     /* Tests the case NULL */
1012     expected=0.0f;
1013     got = D3DXVec2Dot(NULL,&v);
1014     ok(relative_error(got, expected ) < admitted_error, "Expected: %f, Got: %f\n", expected, got);
1015     expected=0.0f;
1016     got = D3DXVec2Dot(NULL,NULL);
1017     ok(relative_error(got, expected ) < admitted_error, "Expected: %f, Got: %f\n", expected, got);
1018
1019 /*_______________D3DXVec2Hermite__________________________*/
1020     expectedvec.x = 2604.625f; expectedvec.y = -4533.0f;
1021     D3DXVec2Hermite(&gotvec,&u,&v,&w,&x,scale);
1022     expect_vec(expectedvec,gotvec);
1023
1024 /*_______________D3DXVec2Length__________________________*/
1025    expected = 5.0f;
1026    got = D3DXVec2Length(&u);
1027    ok(relative_error(got, expected ) < admitted_error, "Expected: %f, Got: %f\n", expected, got);
1028    /* Tests the case NULL */
1029     expected=0.0f;
1030     got = D3DXVec2Length(NULL);
1031     ok(relative_error(got, expected ) < admitted_error, "Expected: %f, Got: %f\n", expected, got);
1032
1033 /*_______________D3DXVec2LengthSq________________________*/
1034    expected = 25.0f;
1035    got = D3DXVec2LengthSq(&u);
1036    ok(relative_error(got, expected ) < admitted_error, "Expected: %f, Got: %f\n", expected, got);
1037    /* Tests the case NULL */
1038     expected=0.0f;
1039     got = D3DXVec2LengthSq(NULL);
1040     ok(relative_error(got, expected ) < admitted_error, "Expected: %f, Got: %f\n", expected, got);
1041
1042 /*_______________D3DXVec2Lerp__________________________*/
1043    expectedvec.x = 68.0f; expectedvec.y = -28.5f;
1044    D3DXVec2Lerp(&gotvec,&u,&v,scale);
1045    expect_vec(expectedvec,gotvec);
1046    /* Tests the case NULL */
1047     funcpointer = D3DXVec2Lerp(&gotvec,NULL,&v,scale);
1048     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
1049     funcpointer = D3DXVec2Lerp(NULL,NULL,NULL,scale);
1050     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
1051
1052 /*_______________D3DXVec2Maximize__________________________*/
1053    expectedvec.x = 3.0f; expectedvec.y = 9.0f;
1054    D3DXVec2Maximize(&gotvec,&u,&v);
1055    expect_vec(expectedvec,gotvec);
1056    /* Tests the case NULL */
1057     funcpointer = D3DXVec2Maximize(&gotvec,NULL,&v);
1058     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
1059     funcpointer = D3DXVec2Maximize(NULL,NULL,NULL);
1060     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
1061
1062 /*_______________D3DXVec2Minimize__________________________*/
1063     expectedvec.x = -7.0f; expectedvec.y = 4.0f;
1064     D3DXVec2Minimize(&gotvec,&u,&v);
1065     expect_vec(expectedvec,gotvec);
1066     /* Tests the case NULL */
1067     funcpointer = D3DXVec2Minimize(&gotvec,NULL,&v);
1068     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
1069     funcpointer = D3DXVec2Minimize(NULL,NULL,NULL);
1070     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
1071
1072 /*_______________D3DXVec2Normalize_________________________*/
1073     expectedvec.x = 0.6f; expectedvec.y = 0.8f;
1074     D3DXVec2Normalize(&gotvec,&u);
1075     expect_vec(expectedvec,gotvec);
1076     /* Test the nul vector */
1077     expectedvec.x = 0.0f; expectedvec.y = 0.0f;
1078     D3DXVec2Normalize(&gotvec,&nul);
1079     expect_vec(expectedvec,gotvec);
1080
1081 /*_______________D3DXVec2Scale____________________________*/
1082     expectedvec.x = -19.5f; expectedvec.y = -26.0f;
1083     D3DXVec2Scale(&gotvec,&u,scale);
1084     expect_vec(expectedvec,gotvec);
1085     /* Tests the case NULL */
1086     funcpointer = D3DXVec2Scale(&gotvec,NULL,scale);
1087     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
1088     funcpointer = D3DXVec2Scale(NULL,NULL,scale);
1089     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
1090
1091 /*_______________D3DXVec2Subtract__________________________*/
1092    expectedvec.x = 10.0f; expectedvec.y = -5.0f;
1093    D3DXVec2Subtract(&gotvec,&u,&v);
1094    expect_vec(expectedvec,gotvec);
1095    /* Tests the case NULL */
1096     funcpointer = D3DXVec2Subtract(&gotvec,NULL,&v);
1097     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
1098     funcpointer = D3DXVec2Subtract(NULL,NULL,NULL);
1099     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
1100
1101 /*_______________D3DXVec2Transform_______________________*/
1102     expectedtrans.x = 36.0f; expectedtrans.y = 44.0f; expectedtrans.z = 52.0f; expectedtrans.w = 60.0f;
1103     D3DXVec2Transform(&gottrans,&u,&mat);
1104     expect_vec4(expectedtrans,gottrans);
1105
1106 /*_______________D3DXVec2TransformCoord_______________________*/
1107     expectedvec.x = 0.6f; expectedvec.y = 11.0f/15.0f;
1108     D3DXVec2TransformCoord(&gotvec,&u,&mat);
1109     expect_vec(expectedvec,gotvec);
1110
1111  /*_______________D3DXVec2TransformNormal______________________*/
1112     expectedvec.x = 23.0f; expectedvec.y = 30.0f;
1113     D3DXVec2TransformNormal(&gotvec,&u,&mat);
1114     expect_vec(expectedvec,gotvec);
1115 }
1116
1117 static void D3DXVector3Test(void)
1118 {
1119     D3DVIEWPORT9 viewport;
1120     D3DXVECTOR3 expectedvec, gotvec, nul, u, v, w, x;
1121     LPD3DXVECTOR3 funcpointer;
1122     D3DXVECTOR4 expectedtrans, gottrans;
1123     D3DXMATRIX mat, projection, view, world;
1124     FLOAT coeff1, coeff2, expected, got, scale;
1125
1126     nul.x = 0.0f; nul.y = 0.0f; nul.z = 0.0f;
1127     u.x = 9.0f; u.y = 6.0f; u.z = 2.0f;
1128     v.x = 2.0f; v.y = -3.0f; v.z = -4.0;
1129     w.x = 3.0f; w.y = -5.0f; w.z = 7.0f;
1130     x.x = 4.0f; x.y = 1.0f; x.z = 11.0f;
1131
1132     viewport.Width = 800; viewport.MinZ = 0.2f; viewport.X = 10;
1133     viewport.Height = 680; viewport.MaxZ = 0.9f; viewport.Y = 5;
1134
1135     U(mat).m[0][0] = 1.0f; U(mat).m[0][1] = 2.0f; U(mat).m[0][2] = 3.0f; U(mat).m[0][3] = 4.0f;
1136     U(mat).m[1][0] = 5.0f; U(mat).m[1][1] = 6.0f; U(mat).m[1][2] = 7.0f; U(mat).m[1][3] = 8.0f;
1137     U(mat).m[2][0] = 9.0f; U(mat).m[2][1] = 10.0f; U(mat).m[2][2] = 11.0f; U(mat).m[2][3] = 12.0f;
1138     U(mat).m[3][0] = 13.0f; U(mat).m[3][1] = 14.0f; U(mat).m[3][2] = 15.0f; U(mat).m[3][3] = 16.0f;
1139
1140     U(view).m[0][1] = 5.0f; U(view).m[0][2] = 7.0f; U(view).m[0][3] = 8.0f;
1141     U(view).m[1][0] = 11.0f; U(view).m[1][2] = 16.0f; U(view).m[1][3] = 33.0f;
1142     U(view).m[2][0] = 19.0f; U(view).m[2][1] = -21.0f; U(view).m[2][3] = 43.0f;
1143     U(view).m[3][0] = 2.0f; U(view).m[3][1] = 3.0f; U(view).m[3][2] = -4.0f;
1144     U(view).m[0][0] = 10.0f; U(view).m[1][1] = 20.0f; U(view).m[2][2] = 30.0f;
1145     U(view).m[3][3] = -40.0f;
1146
1147     U(world).m[0][0] = 21.0f; U(world).m[0][1] = 2.0f; U(world).m[0][2] = 3.0f; U(world).m[0][3] = 4.0;
1148     U(world).m[1][0] = 5.0f; U(world).m[1][1] = 23.0f; U(world).m[1][2] = 7.0f; U(world).m[1][3] = 8.0f;
1149     U(world).m[2][0] = -8.0f; U(world).m[2][1] = -7.0f; U(world).m[2][2] = 25.0f; U(world).m[2][3] = -5.0f;
1150     U(world).m[3][0] = -4.0f; U(world).m[3][1] = -3.0f; U(world).m[3][2] = -2.0f; U(world).m[3][3] = 27.0f;
1151
1152     coeff1 = 2.0f; coeff2 = 5.0f;
1153     scale = -6.5f;
1154
1155 /*_______________D3DXVec3Add__________________________*/
1156     expectedvec.x = 11.0f; expectedvec.y = 3.0f; expectedvec.z = -2.0f;
1157     D3DXVec3Add(&gotvec,&u,&v);
1158     expect_vec3(expectedvec,gotvec);
1159     /* Tests the case NULL */
1160     funcpointer = D3DXVec3Add(&gotvec,NULL,&v);
1161     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
1162     funcpointer = D3DXVec3Add(NULL,NULL,NULL);
1163     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
1164
1165 /*_______________D3DXVec3BaryCentric___________________*/
1166     expectedvec.x = -35.0f; expectedvec.y = -67.0; expectedvec.z = 15.0f;
1167     D3DXVec3BaryCentric(&gotvec,&u,&v,&w,coeff1,coeff2);
1168
1169     expect_vec3(expectedvec,gotvec);
1170
1171 /*_______________D3DXVec3CatmullRom____________________*/
1172     expectedvec.x = 1458.0f; expectedvec.y = 22.1875f; expectedvec.z = 4141.375f;
1173     D3DXVec3CatmullRom(&gotvec,&u,&v,&w,&x,scale);
1174     expect_vec3(expectedvec,gotvec);
1175
1176 /*_______________D3DXVec3Cross________________________*/
1177     expectedvec.x = -18.0f; expectedvec.y = 40.0f; expectedvec.z = -39.0f;
1178     D3DXVec3Cross(&gotvec,&u,&v);
1179     expect_vec3(expectedvec,gotvec);
1180     /* Tests the case NULL */
1181     funcpointer = D3DXVec3Cross(&gotvec,NULL,&v);
1182     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
1183     funcpointer = D3DXVec3Cross(NULL,NULL,NULL);
1184     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
1185
1186 /*_______________D3DXVec3Dot__________________________*/
1187     expected = -8.0f;
1188     got = D3DXVec3Dot(&u,&v);
1189     ok(relative_error(got, expected ) < admitted_error, "Expected: %f, Got: %f\n", expected, got);
1190     /* Tests the case NULL */
1191     expected=0.0f;
1192     got = D3DXVec3Dot(NULL,&v);
1193     ok(relative_error(got, expected ) < admitted_error, "Expected: %f, Got: %f\n", expected, got);
1194     expected=0.0f;
1195     got = D3DXVec3Dot(NULL,NULL);
1196     ok(relative_error(got, expected ) < admitted_error, "Expected: %f, Got: %f\n", expected, got);
1197
1198 /*_______________D3DXVec3Hermite__________________________*/
1199     expectedvec.x = -6045.75f; expectedvec.y = -6650.0f; expectedvec.z = 1358.875f;
1200     D3DXVec3Hermite(&gotvec,&u,&v,&w,&x,scale);
1201     expect_vec3(expectedvec,gotvec);
1202
1203 /*_______________D3DXVec3Length__________________________*/
1204    expected = 11.0f;
1205    got = D3DXVec3Length(&u);
1206    ok(relative_error(got, expected ) < admitted_error, "Expected: %f, Got: %f\n", expected, got);
1207    /* Tests the case NULL */
1208     expected=0.0f;
1209     got = D3DXVec3Length(NULL);
1210     ok(relative_error(got, expected ) < admitted_error, "Expected: %f, Got: %f\n", expected, got);
1211
1212 /*_______________D3DXVec3LengthSq________________________*/
1213     expected = 121.0f;
1214     got = D3DXVec3LengthSq(&u);
1215     ok(relative_error(got, expected ) < admitted_error, "Expected: %f, Got: %f\n", expected, got);
1216    /* Tests the case NULL */
1217     expected=0.0f;
1218     got = D3DXVec3LengthSq(NULL);
1219     ok(relative_error(got, expected ) < admitted_error, "Expected: %f, Got: %f\n", expected, got);
1220
1221 /*_______________D3DXVec3Lerp__________________________*/
1222     expectedvec.x = 54.5f; expectedvec.y = 64.5f, expectedvec.z = 41.0f ;
1223     D3DXVec3Lerp(&gotvec,&u,&v,scale);
1224     expect_vec3(expectedvec,gotvec);
1225     /* Tests the case NULL */
1226     funcpointer = D3DXVec3Lerp(&gotvec,NULL,&v,scale);
1227     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
1228     funcpointer = D3DXVec3Lerp(NULL,NULL,NULL,scale);
1229     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
1230
1231 /*_______________D3DXVec3Maximize__________________________*/
1232     expectedvec.x = 9.0f; expectedvec.y = 6.0f; expectedvec.z = 2.0f;
1233     D3DXVec3Maximize(&gotvec,&u,&v);
1234     expect_vec3(expectedvec,gotvec);
1235     /* Tests the case NULL */
1236     funcpointer = D3DXVec3Maximize(&gotvec,NULL,&v);
1237     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
1238     funcpointer = D3DXVec3Maximize(NULL,NULL,NULL);
1239     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
1240
1241 /*_______________D3DXVec3Minimize__________________________*/
1242     expectedvec.x = 2.0f; expectedvec.y = -3.0f; expectedvec.z = -4.0f;
1243     D3DXVec3Minimize(&gotvec,&u,&v);
1244     expect_vec3(expectedvec,gotvec);
1245     /* Tests the case NULL */
1246     funcpointer = D3DXVec3Minimize(&gotvec,NULL,&v);
1247     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
1248     funcpointer = D3DXVec3Minimize(NULL,NULL,NULL);
1249     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
1250
1251 /*_______________D3DXVec3Normalize_________________________*/
1252     expectedvec.x = 9.0f/11.0f; expectedvec.y = 6.0f/11.0f; expectedvec.z = 2.0f/11.0f;
1253     D3DXVec3Normalize(&gotvec,&u);
1254     expect_vec3(expectedvec,gotvec);
1255     /* Test the nul vector */
1256     expectedvec.x = 0.0f; expectedvec.y = 0.0f; expectedvec.z = 0.0f;
1257     D3DXVec3Normalize(&gotvec,&nul);
1258     expect_vec3(expectedvec,gotvec);
1259
1260 /*_______________D3DXVec3Project_________________________*/
1261     expectedvec.x = 1135.721924f; expectedvec.y = 147.086914f; expectedvec.z = 0.153412f;
1262     D3DXMatrixPerspectiveFovLH(&projection,D3DX_PI/4.0f,20.0f/17.0f,1.0f,1000.0f);
1263     D3DXVec3Project(&gotvec,&u,&viewport,&projection,&view,&world);
1264     expect_vec3(expectedvec,gotvec);
1265
1266 /*_______________D3DXVec3Scale____________________________*/
1267     expectedvec.x = -58.5f; expectedvec.y = -39.0f; expectedvec.z = -13.0f;
1268     D3DXVec3Scale(&gotvec,&u,scale);
1269     expect_vec3(expectedvec,gotvec);
1270     /* Tests the case NULL */
1271     funcpointer = D3DXVec3Scale(&gotvec,NULL,scale);
1272     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
1273     funcpointer = D3DXVec3Scale(NULL,NULL,scale);
1274     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
1275
1276 /*_______________D3DXVec3Subtract_______________________*/
1277     expectedvec.x = 7.0f; expectedvec.y = 9.0f; expectedvec.z = 6.0f;
1278     D3DXVec3Subtract(&gotvec,&u,&v);
1279     expect_vec3(expectedvec,gotvec);
1280     /* Tests the case NULL */
1281     funcpointer = D3DXVec3Subtract(&gotvec,NULL,&v);
1282     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
1283     funcpointer = D3DXVec3Subtract(NULL,NULL,NULL);
1284     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
1285
1286 /*_______________D3DXVec3Transform_______________________*/
1287     expectedtrans.x = 70.0f; expectedtrans.y = 88.0f; expectedtrans.z = 106.0f; expectedtrans.w = 124.0f;
1288     D3DXVec3Transform(&gottrans,&u,&mat);
1289     expect_vec4(expectedtrans,gottrans);
1290
1291 /*_______________D3DXVec3TransformCoord_______________________*/
1292     expectedvec.x = 70.0f/124.0f; expectedvec.y = 88.0f/124.0f; expectedvec.z = 106.0f/124.0f;
1293     D3DXVec3TransformCoord(&gotvec,&u,&mat);
1294     expect_vec3(expectedvec,gotvec);
1295
1296 /*_______________D3DXVec3TransformNormal______________________*/
1297     expectedvec.x = 57.0f; expectedvec.y = 74.0f; expectedvec.z = 91.0f;
1298     D3DXVec3TransformNormal(&gotvec,&u,&mat);
1299     expect_vec3(expectedvec,gotvec);
1300
1301 /*_______________D3DXVec3Unproject_________________________*/
1302     expectedvec.x = -2.913411f; expectedvec.y = 1.593215f; expectedvec.z = 0.380724f;
1303     D3DXMatrixPerspectiveFovLH(&projection,D3DX_PI/4.0f,20.0f/17.0f,1.0f,1000.0f);
1304     D3DXVec3Unproject(&gotvec,&u,&viewport,&projection,&view,&world);
1305     expect_vec3(expectedvec,gotvec);
1306 }
1307
1308 static void D3DXVector4Test(void)
1309 {
1310     D3DXVECTOR4 expectedvec, gotvec, u, v, w, x;
1311     LPD3DXVECTOR4 funcpointer;
1312     D3DXVECTOR4 expectedtrans, gottrans;
1313     D3DXMATRIX mat;
1314     FLOAT coeff1, coeff2, expected, got, scale;
1315
1316     u.x = 1.0f; u.y = 2.0f; u.z = 4.0f; u.w = 10.0;
1317     v.x = -3.0f; v.y = 4.0f; v.z = -5.0f; v.w = 7.0;
1318     w.x = 4.0f; w.y =6.0f; w.z = -2.0f; w.w = 1.0f;
1319     x.x = 6.0f; x.y = -7.0f; x.z =8.0f; x.w = -9.0f;
1320
1321     U(mat).m[0][0] = 1.0f; U(mat).m[0][1] = 2.0f; U(mat).m[0][2] = 3.0f; U(mat).m[0][3] = 4.0f;
1322     U(mat).m[1][0] = 5.0f; U(mat).m[1][1] = 6.0f; U(mat).m[1][2] = 7.0f; U(mat).m[1][3] = 8.0f;
1323     U(mat).m[2][0] = 9.0f; U(mat).m[2][1] = 10.0f; U(mat).m[2][2] = 11.0f; U(mat).m[2][3] = 12.0f;
1324     U(mat).m[3][0] = 13.0f; U(mat).m[3][1] = 14.0f; U(mat).m[3][2] = 15.0f; U(mat).m[3][3] = 16.0f;
1325
1326     coeff1 = 2.0f; coeff2 = 5.0;
1327     scale = -6.5f;
1328
1329 /*_______________D3DXVec4Add__________________________*/
1330     expectedvec.x = -2.0f; expectedvec.y = 6.0f; expectedvec.z = -1.0f; expectedvec.w = 17.0f;
1331     D3DXVec4Add(&gotvec,&u,&v);
1332     expect_vec4(expectedvec,gotvec);
1333     /* Tests the case NULL */
1334     funcpointer = D3DXVec4Add(&gotvec,NULL,&v);
1335     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
1336     funcpointer = D3DXVec4Add(NULL,NULL,NULL);
1337     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
1338
1339 /*_______________D3DXVec4BaryCentric____________________*/
1340     expectedvec.x = 8.0f; expectedvec.y = 26.0; expectedvec.z =  -44.0f; expectedvec.w = -41.0f;
1341     D3DXVec4BaryCentric(&gotvec,&u,&v,&w,coeff1,coeff2);
1342     expect_vec4(expectedvec,gotvec);
1343
1344 /*_______________D3DXVec4CatmullRom____________________*/
1345     expectedvec.x = 2754.625f; expectedvec.y = 2367.5625f; expectedvec.z = 1060.1875f; expectedvec.w = 131.3125f;
1346     D3DXVec4CatmullRom(&gotvec,&u,&v,&w,&x,scale);
1347     expect_vec4(expectedvec,gotvec);
1348
1349 /*_______________D3DXVec4Cross_________________________*/
1350     expectedvec.x = 390.0f; expectedvec.y = -393.0f; expectedvec.z = -316.0f; expectedvec.w = 166.0f;
1351     D3DXVec4Cross(&gotvec,&u,&v,&w);
1352     expect_vec4(expectedvec,gotvec);
1353
1354 /*_______________D3DXVec4Dot__________________________*/
1355     expected = 55.0f;
1356     got = D3DXVec4Dot(&u,&v);
1357     ok(relative_error(got, expected ) < admitted_error, "Expected: %f, Got: %f\n", expected, got);
1358     /* Tests the case NULL */
1359     expected=0.0f;
1360     got = D3DXVec4Dot(NULL,&v);
1361     ok(relative_error(got, expected ) < admitted_error, "Expected: %f, Got: %f\n", expected, got);
1362     expected=0.0f;
1363     got = D3DXVec4Dot(NULL,NULL);
1364     ok(relative_error(got, expected ) < admitted_error, "Expected: %f, Got: %f\n", expected, got);
1365
1366 /*_______________D3DXVec4Hermite_________________________*/
1367     expectedvec.x = 1224.625f; expectedvec.y = 3461.625f; expectedvec.z = -4758.875f; expectedvec.w = -5781.5f;
1368     D3DXVec4Hermite(&gotvec,&u,&v,&w,&x,scale);
1369     expect_vec4(expectedvec,gotvec);
1370
1371 /*_______________D3DXVec4Length__________________________*/
1372    expected = 11.0f;
1373    got = D3DXVec4Length(&u);
1374    ok(relative_error(got, expected ) < admitted_error, "Expected: %f, Got: %f\n", expected, got);
1375    /* Tests the case NULL */
1376     expected=0.0f;
1377     got = D3DXVec4Length(NULL);
1378     ok(relative_error(got, expected ) < admitted_error, "Expected: %f, Got: %f\n", expected, got);
1379
1380 /*_______________D3DXVec4LengthSq________________________*/
1381     expected = 121.0f;
1382     got = D3DXVec4LengthSq(&u);
1383     ok(relative_error(got, expected ) < admitted_error, "Expected: %f, Got: %f\n", expected, got);
1384     /* Tests the case NULL */
1385     expected=0.0f;
1386     got = D3DXVec4LengthSq(NULL);
1387     ok(relative_error(got, expected ) < admitted_error, "Expected: %f, Got: %f\n", expected, got);
1388
1389 /*_______________D3DXVec4Lerp__________________________*/
1390     expectedvec.x = 27.0f; expectedvec.y = -11.0f; expectedvec.z = 62.5;  expectedvec.w = 29.5;
1391     D3DXVec4Lerp(&gotvec,&u,&v,scale);
1392     expect_vec4(expectedvec,gotvec);
1393     /* Tests the case NULL */
1394     funcpointer = D3DXVec4Lerp(&gotvec,NULL,&v,scale);
1395     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
1396     funcpointer = D3DXVec4Lerp(NULL,NULL,NULL,scale);
1397     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
1398
1399 /*_______________D3DXVec4Maximize__________________________*/
1400     expectedvec.x = 1.0f; expectedvec.y = 4.0f; expectedvec.z = 4.0f; expectedvec.w = 10.0;
1401     D3DXVec4Maximize(&gotvec,&u,&v);
1402     expect_vec4(expectedvec,gotvec);
1403     /* Tests the case NULL */
1404     funcpointer = D3DXVec4Maximize(&gotvec,NULL,&v);
1405     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
1406     funcpointer = D3DXVec4Maximize(NULL,NULL,NULL);
1407     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
1408
1409 /*_______________D3DXVec4Minimize__________________________*/
1410     expectedvec.x = -3.0f; expectedvec.y = 2.0f; expectedvec.z = -5.0f; expectedvec.w = 7.0;
1411     D3DXVec4Minimize(&gotvec,&u,&v);
1412     expect_vec4(expectedvec,gotvec);
1413     /* Tests the case NULL */
1414     funcpointer = D3DXVec4Minimize(&gotvec,NULL,&v);
1415     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
1416     funcpointer = D3DXVec4Minimize(NULL,NULL,NULL);
1417     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
1418
1419 /*_______________D3DXVec4Normalize_________________________*/
1420     expectedvec.x = 1.0f/11.0f; expectedvec.y = 2.0f/11.0f; expectedvec.z = 4.0f/11.0f; expectedvec.w = 10.0f/11.0f;
1421     D3DXVec4Normalize(&gotvec,&u);
1422     expect_vec4(expectedvec,gotvec);
1423
1424 /*_______________D3DXVec4Scale____________________________*/
1425     expectedvec.x = -6.5f; expectedvec.y = -13.0f; expectedvec.z = -26.0f; expectedvec.w = -65.0f;
1426     D3DXVec4Scale(&gotvec,&u,scale);
1427     expect_vec4(expectedvec,gotvec);
1428     /* Tests the case NULL */
1429     funcpointer = D3DXVec4Scale(&gotvec,NULL,scale);
1430     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
1431     funcpointer = D3DXVec4Scale(NULL,NULL,scale);
1432     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
1433
1434 /*_______________D3DXVec4Subtract__________________________*/
1435     expectedvec.x = 4.0f; expectedvec.y = -2.0f; expectedvec.z = 9.0f; expectedvec.w = 3.0f;
1436     D3DXVec4Subtract(&gotvec,&u,&v);
1437     expect_vec4(expectedvec,gotvec);
1438     /* Tests the case NULL */
1439     funcpointer = D3DXVec4Subtract(&gotvec,NULL,&v);
1440     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
1441     funcpointer = D3DXVec4Subtract(NULL,NULL,NULL);
1442     ok(funcpointer == NULL, "Expected: %p, Got: %p\n", NULL, funcpointer);
1443
1444 /*_______________D3DXVec4Transform_______________________*/
1445     expectedtrans.x = 177.0f; expectedtrans.y = 194.0f; expectedtrans.z = 211.0f; expectedtrans.w = 228.0f;
1446     D3DXVec4Transform(&gottrans,&u,&mat);
1447     expect_vec4(expectedtrans,gottrans);
1448 }
1449
1450 static void test_matrix_stack(void)
1451 {
1452     ID3DXMatrixStack *stack;
1453     ULONG refcount;
1454     HRESULT hr;
1455
1456     const D3DXMATRIX mat1 = {{{
1457          1.0f,  2.0f,  3.0f,  4.0f,
1458          5.0f,  6.0f,  7.0f,  8.0f,
1459          9.0f, 10.0f, 11.0f, 12.0f,
1460         13.0f, 14.0f, 15.0f, 16.0f
1461     }}};
1462
1463     const D3DXMATRIX mat2 = {{{
1464         17.0f, 18.0f, 19.0f, 20.0f,
1465         21.0f, 22.0f, 23.0f, 24.0f,
1466         25.0f, 26.0f, 27.0f, 28.0f,
1467         29.0f, 30.0f, 31.0f, 32.0f
1468     }}};
1469
1470     hr = D3DXCreateMatrixStack(0, &stack);
1471     ok(SUCCEEDED(hr), "Failed to create a matrix stack, hr %#x\n", hr);
1472     if (FAILED(hr)) return;
1473
1474     ok(D3DXMatrixIsIdentity(ID3DXMatrixStack_GetTop(stack)),
1475             "The top of an empty matrix stack should be an identity matrix\n");
1476
1477     hr = ID3DXMatrixStack_Pop(stack);
1478     ok(SUCCEEDED(hr), "Pop failed, hr %#x\n", hr);
1479
1480     hr = ID3DXMatrixStack_Push(stack);
1481     ok(SUCCEEDED(hr), "Push failed, hr %#x\n", hr);
1482     ok(D3DXMatrixIsIdentity(ID3DXMatrixStack_GetTop(stack)), "The top should be an identity matrix\n");
1483
1484     hr = ID3DXMatrixStack_Push(stack);
1485     ok(SUCCEEDED(hr), "Push failed, hr %#x\n", hr);
1486
1487     hr = ID3DXMatrixStack_LoadMatrix(stack, &mat1);
1488     ok(SUCCEEDED(hr), "LoadMatrix failed, hr %#x\n", hr);
1489     expect_mat(&mat1, ID3DXMatrixStack_GetTop(stack));
1490
1491     hr = ID3DXMatrixStack_Push(stack);
1492     ok(SUCCEEDED(hr), "Push failed, hr %#x\n", hr);
1493     expect_mat(&mat1, ID3DXMatrixStack_GetTop(stack));
1494
1495     hr = ID3DXMatrixStack_LoadMatrix(stack, &mat2);
1496     ok(SUCCEEDED(hr), "LoadMatrix failed, hr %#x\n", hr);
1497     expect_mat(&mat2, ID3DXMatrixStack_GetTop(stack));
1498
1499     hr = ID3DXMatrixStack_Push(stack);
1500     ok(SUCCEEDED(hr), "Push failed, hr %#x\n", hr);
1501     expect_mat(&mat2, ID3DXMatrixStack_GetTop(stack));
1502
1503     hr = ID3DXMatrixStack_LoadIdentity(stack);
1504     ok(SUCCEEDED(hr), "LoadIdentity failed, hr %#x\n", hr);
1505     ok(D3DXMatrixIsIdentity(ID3DXMatrixStack_GetTop(stack)), "The top should be an identity matrix\n");
1506
1507     hr = ID3DXMatrixStack_Pop(stack);
1508     ok(SUCCEEDED(hr), "Pop failed, hr %#x\n", hr);
1509     expect_mat(&mat2, ID3DXMatrixStack_GetTop(stack));
1510
1511     hr = ID3DXMatrixStack_Pop(stack);
1512     ok(SUCCEEDED(hr), "Pop failed, hr %#x\n", hr);
1513     expect_mat(&mat1, ID3DXMatrixStack_GetTop(stack));
1514
1515     hr = ID3DXMatrixStack_Pop(stack);
1516     ok(SUCCEEDED(hr), "Pop failed, hr %#x\n", hr);
1517     ok(D3DXMatrixIsIdentity(ID3DXMatrixStack_GetTop(stack)), "The top should be an identity matrix\n");
1518
1519     hr = ID3DXMatrixStack_Pop(stack);
1520     ok(SUCCEEDED(hr), "Pop failed, hr %#x\n", hr);
1521     ok(D3DXMatrixIsIdentity(ID3DXMatrixStack_GetTop(stack)), "The top should be an identity matrix\n");
1522
1523     refcount = ID3DXMatrixStack_Release(stack);
1524     ok(!refcount, "Matrix stack has %u references left.\n", refcount);
1525 }
1526
1527 static void test_Matrix_AffineTransformation2D(void)
1528 {
1529     D3DXMATRIX exp_mat, got_mat;
1530     D3DXVECTOR2 center, position;
1531     FLOAT angle, scale;
1532
1533     center.x = 3.0f;
1534     center.y = 4.0f;
1535
1536     position.x = -6.0f;
1537     position.y = 7.0f;
1538
1539     angle = D3DX_PI/3.0f;
1540
1541     scale = 20.0f;
1542
1543     U(exp_mat).m[0][0] = 10.0f;
1544     U(exp_mat).m[1][0] = -17.320507f;
1545     U(exp_mat).m[2][0] = 0.0f;
1546     U(exp_mat).m[3][0] = -1.035898f;
1547     U(exp_mat).m[0][1] = 17.320507f;
1548     U(exp_mat).m[1][1] = 10.0f;
1549     U(exp_mat).m[2][1] = 0.0f;
1550     U(exp_mat).m[3][1] = 6.401924f;
1551     U(exp_mat).m[0][2] = 0.0f;
1552     U(exp_mat).m[1][2] = 0.0f;
1553     U(exp_mat).m[2][2] = 1.0f;
1554     U(exp_mat).m[3][2] = 0.0f;
1555     U(exp_mat).m[0][3] = 0.0f;
1556     U(exp_mat).m[1][3] = 0.0f;
1557     U(exp_mat).m[2][3] = 0.0f;
1558     U(exp_mat).m[3][3] = 1.0f;
1559
1560     D3DXMatrixAffineTransformation2D(&got_mat, scale, &center, angle, &position);
1561
1562     expect_mat(&exp_mat, &got_mat);
1563
1564 /*______________*/
1565
1566     center.x = 3.0f;
1567     center.y = 4.0f;
1568
1569     angle = D3DX_PI/3.0f;
1570
1571     scale = 20.0f;
1572
1573     U(exp_mat).m[0][0] = 10.0f;
1574     U(exp_mat).m[1][0] = -17.320507f;
1575     U(exp_mat).m[2][0] = 0.0f;
1576     U(exp_mat).m[3][0] = 4.964102f;
1577     U(exp_mat).m[0][1] = 17.320507f;
1578     U(exp_mat).m[1][1] = 10.0f;
1579     U(exp_mat).m[2][1] = 0.0f;
1580     U(exp_mat).m[3][1] = -0.598076f;
1581     U(exp_mat).m[0][2] = 0.0f;
1582     U(exp_mat).m[1][2] = 0.0f;
1583     U(exp_mat).m[2][2] = 1.0f;
1584     U(exp_mat).m[3][2] = 0.0f;
1585     U(exp_mat).m[0][3] = 0.0f;
1586     U(exp_mat).m[1][3] = 0.0f;
1587     U(exp_mat).m[2][3] = 0.0f;
1588     U(exp_mat).m[3][3] = 1.0f;
1589
1590     D3DXMatrixAffineTransformation2D(&got_mat, scale, &center, angle, NULL);
1591
1592     expect_mat(&exp_mat, &got_mat);
1593
1594 /*______________*/
1595
1596     position.x = -6.0f;
1597     position.y = 7.0f;
1598
1599     angle = D3DX_PI/3.0f;
1600
1601     scale = 20.0f;
1602
1603     U(exp_mat).m[0][0] = 10.0f;
1604     U(exp_mat).m[1][0] = -17.320507f;
1605     U(exp_mat).m[2][0] = 0.0f;
1606     U(exp_mat).m[3][0] = -6.0f;
1607     U(exp_mat).m[0][1] = 17.320507f;
1608     U(exp_mat).m[1][1] = 10.0f;
1609     U(exp_mat).m[2][1] = 0.0f;
1610     U(exp_mat).m[3][1] = 7.0f;
1611     U(exp_mat).m[0][2] = 0.0f;
1612     U(exp_mat).m[1][2] = 0.0f;
1613     U(exp_mat).m[2][2] = 1.0f;
1614     U(exp_mat).m[3][2] = 0.0f;
1615     U(exp_mat).m[0][3] = 0.0f;
1616     U(exp_mat).m[1][3] = 0.0f;
1617     U(exp_mat).m[2][3] = 0.0f;
1618     U(exp_mat).m[3][3] = 1.0f;
1619
1620     D3DXMatrixAffineTransformation2D(&got_mat, scale, NULL, angle, &position);
1621
1622     expect_mat(&exp_mat, &got_mat);
1623
1624 /*______________*/
1625
1626     angle = 5.0f * D3DX_PI/4.0f;
1627
1628     scale = -20.0f;
1629
1630     U(exp_mat).m[0][0] = 14.142133f;
1631     U(exp_mat).m[1][0] = -14.142133f;
1632     U(exp_mat).m[2][0] = 0.0f;
1633     U(exp_mat).m[3][0] = 0.0f;
1634     U(exp_mat).m[0][1] = 14.142133;
1635     U(exp_mat).m[1][1] = 14.142133f;
1636     U(exp_mat).m[2][1] = 0.0f;
1637     U(exp_mat).m[3][1] = 0.0f;
1638     U(exp_mat).m[0][2] = 0.0f;
1639     U(exp_mat).m[1][2] = 0.0f;
1640     U(exp_mat).m[2][2] = 1.0f;
1641     U(exp_mat).m[3][2] = 0.0f;
1642     U(exp_mat).m[0][3] = 0.0f;
1643     U(exp_mat).m[1][3] = 0.0f;
1644     U(exp_mat).m[2][3] = 0.0f;
1645     U(exp_mat).m[3][3] = 1.0f;
1646
1647     D3DXMatrixAffineTransformation2D(&got_mat, scale, NULL, angle, NULL);
1648
1649     expect_mat(&exp_mat, &got_mat);
1650 }
1651
1652 static void test_Matrix_Decompose(void)
1653 {
1654     D3DXMATRIX pm;
1655     D3DXQUATERNION exp_rotation, got_rotation;
1656     D3DXVECTOR3 exp_scale, got_scale, exp_translation, got_translation;
1657     HRESULT hr;
1658
1659 /*___________*/
1660
1661     U(pm).m[0][0] = -0.9238790f;
1662     U(pm).m[1][0] = -0.2705984f;
1663     U(pm).m[2][0] = 0.2705984f;
1664     U(pm).m[3][0] = -5.0f;
1665     U(pm).m[0][1] = 0.2705984f;
1666     U(pm).m[1][1] = 0.03806049f;
1667     U(pm).m[2][1] = 0.9619395f;
1668     U(pm).m[3][1] = 0.0f;
1669     U(pm).m[0][2] = -0.2705984f;
1670     U(pm).m[1][2] = 0.9619395f;
1671     U(pm).m[2][2] = 0.03806049f;
1672     U(pm).m[3][2] = 10.0f;
1673     U(pm).m[0][3] = 0.0f;
1674     U(pm).m[1][3] = 0.0f;
1675     U(pm).m[2][3] = 0.0f;
1676     U(pm).m[3][3] = 1.0f;
1677
1678     exp_scale.x = 1.0f;
1679     exp_scale.y = 1.0f;
1680     exp_scale.z = 1.0f;
1681
1682     exp_rotation.w = 0.195091f;
1683     exp_rotation.x = 0.0f;
1684     exp_rotation.y = 0.693520f;
1685     exp_rotation.z = 0.693520f;
1686
1687     exp_translation.x = -5.0f;
1688     exp_translation.y = 0.0f;
1689     exp_translation.z = 10.0f;
1690
1691     D3DXMatrixDecompose(&got_scale, &got_rotation, &got_translation, &pm);
1692
1693     compare_scale(exp_scale, got_scale);
1694     compare_rotation(exp_rotation, got_rotation);
1695     compare_translation(exp_translation, got_translation);
1696
1697 /*_________*/
1698
1699     U(pm).m[0][0] = -2.255813f;
1700     U(pm).m[1][0] = 1.302324f;
1701     U(pm).m[2][0] = 1.488373f;
1702     U(pm).m[3][0] = 1.0f;
1703     U(pm).m[0][1] = 1.302327f;
1704     U(pm).m[1][1] = -0.7209296f;
1705     U(pm).m[2][1] = 2.60465f;
1706     U(pm).m[3][1] = 2.0f;
1707     U(pm).m[0][2] = 1.488371f;
1708     U(pm).m[1][2] = 2.604651f;
1709     U(pm).m[2][2] = -0.02325551f;
1710     U(pm).m[3][2] = 3.0f;
1711     U(pm).m[0][3] = 0.0f;
1712     U(pm).m[1][3] = 0.0f;
1713     U(pm).m[2][3] = 0.0f;
1714     U(pm).m[3][3] = 1.0f;
1715
1716     exp_scale.x = 3.0f;
1717     exp_scale.y = 3.0f;
1718     exp_scale.z = 3.0f;
1719
1720     exp_rotation.w = 0.0;
1721     exp_rotation.x = 0.352180f;
1722     exp_rotation.y = 0.616316f;
1723     exp_rotation.z = 0.704361f;
1724
1725     exp_translation.x = 1.0f;
1726     exp_translation.y = 2.0f;
1727     exp_translation.z = 3.0f;
1728
1729     D3DXMatrixDecompose(&got_scale, &got_rotation, &got_translation, &pm);
1730
1731     compare_scale(exp_scale, got_scale);
1732     compare_rotation(exp_rotation, got_rotation);
1733     compare_translation(exp_translation, got_translation);
1734
1735 /*_____________*/
1736
1737     U(pm).m[0][0] = 2.427051f;
1738     U(pm).m[1][0] = 0.0f;
1739     U(pm).m[2][0] = 1.763355f;
1740     U(pm).m[3][0] = 5.0f;
1741     U(pm).m[0][1] = 0.0f;
1742     U(pm).m[1][1] = 3.0f;
1743     U(pm).m[2][1] = 0.0f;
1744     U(pm).m[3][1] = 5.0f;
1745     U(pm).m[0][2] = -1.763355f;
1746     U(pm).m[1][2] = 0.0f;
1747     U(pm).m[2][2] = 2.427051f;
1748     U(pm).m[3][2] = 5.0f;
1749     U(pm).m[0][3] = 0.0f;
1750     U(pm).m[1][3] = 0.0f;
1751     U(pm).m[2][3] = 0.0f;
1752     U(pm).m[3][3] = 1.0f;
1753
1754     exp_scale.x = 3.0f;
1755     exp_scale.y = 3.0f;
1756     exp_scale.z = 3.0f;
1757
1758     exp_rotation.w = 0.951057f;
1759     exp_rotation.x = 0.0f;
1760     exp_rotation.y = 0.309017f;
1761     exp_rotation.z = 0.0f;
1762
1763     exp_translation.x = 5.0f;
1764     exp_translation.y = 5.0f;
1765     exp_translation.z = 5.0f;
1766
1767     D3DXMatrixDecompose(&got_scale, &got_rotation, &got_translation, &pm);
1768
1769     compare_scale(exp_scale, got_scale);
1770     compare_rotation(exp_rotation, got_rotation);
1771     compare_translation(exp_translation, got_translation);
1772
1773 /*_____________*/
1774
1775     U(pm).m[0][0] = -0.9238790f;
1776     U(pm).m[1][0] = -0.2705984f;
1777     U(pm).m[2][0] = 0.2705984f;
1778     U(pm).m[3][0] = -5.0f;
1779     U(pm).m[0][1] = 0.2705984f;
1780     U(pm).m[1][1] = 0.03806049f;
1781     U(pm).m[2][1] = 0.9619395f;
1782     U(pm).m[3][1] = 0.0f;
1783     U(pm).m[0][2] = -0.2705984f;
1784     U(pm).m[1][2] = 0.9619395f;
1785     U(pm).m[2][2] = 0.03806049f;
1786     U(pm).m[3][2] = 10.0f;
1787     U(pm).m[0][3] = 0.0f;
1788     U(pm).m[1][3] = 0.0f;
1789     U(pm).m[2][3] = 0.0f;
1790     U(pm).m[3][3] = 1.0f;
1791
1792     exp_scale.x = 1.0f;
1793     exp_scale.y = 1.0f;
1794     exp_scale.z = 1.0f;
1795
1796     exp_rotation.w = 0.195091f;
1797     exp_rotation.x = 0.0f;
1798     exp_rotation.y = 0.693520f;
1799     exp_rotation.z = 0.693520f;
1800
1801     exp_translation.x = -5.0f;
1802     exp_translation.y = 0.0f;
1803     exp_translation.z = 10.0f;
1804
1805     D3DXMatrixDecompose(&got_scale, &got_rotation, &got_translation, &pm);
1806
1807     compare_scale(exp_scale, got_scale);
1808     compare_rotation(exp_rotation, got_rotation);
1809     compare_translation(exp_translation, got_translation);
1810
1811 /*__________*/
1812
1813     U(pm).m[0][0] = -0.9238790f;
1814     U(pm).m[1][0] = -0.5411968f;
1815     U(pm).m[2][0] = 0.8117952f;
1816     U(pm).m[3][0] = -5.0f;
1817     U(pm).m[0][1] = 0.2705984f;
1818     U(pm).m[1][1] = 0.07612098f;
1819     U(pm).m[2][1] = 2.8858185f;
1820     U(pm).m[3][1] = 0.0f;
1821     U(pm).m[0][2] = -0.2705984f;
1822     U(pm).m[1][2] = 1.9238790f;
1823     U(pm).m[2][2] = 0.11418147f;
1824     U(pm).m[3][2] = 10.0f;
1825     U(pm).m[0][3] = 0.0f;
1826     U(pm).m[1][3] = 0.0f;
1827     U(pm).m[2][3] = 0.0f;
1828     U(pm).m[3][3] = 1.0f;
1829
1830     exp_scale.x = 1.0f;
1831     exp_scale.y = 2.0f;
1832     exp_scale.z = 3.0f;
1833
1834     exp_rotation.w = 0.195091f;
1835     exp_rotation.x = 0.0f;
1836     exp_rotation.y = 0.693520f;
1837     exp_rotation.z = 0.693520f;
1838
1839     exp_translation.x = -5.0f;
1840     exp_translation.y = 0.0f;
1841     exp_translation.z = 10.0f;
1842
1843     D3DXMatrixDecompose(&got_scale, &got_rotation, &got_translation, &pm);
1844
1845     compare_scale(exp_scale, got_scale);
1846     compare_rotation(exp_rotation, got_rotation);
1847     compare_translation(exp_translation, got_translation);
1848
1849 /*__________*/
1850
1851     U(pm).m[0][0] = 0.7156004f;
1852     U(pm).m[1][0] = -0.5098283f;
1853     U(pm).m[2][0] = -0.4774843f;
1854     U(pm).m[3][0] = -5.0f;
1855     U(pm).m[0][1] = -0.6612288f;
1856     U(pm).m[1][1] = -0.7147621f;
1857     U(pm).m[2][1] = -0.2277977f;
1858     U(pm).m[3][1] = 0.0f;
1859     U(pm).m[0][2] = -0.2251499f;
1860     U(pm).m[1][2] = 0.4787385f;
1861     U(pm).m[2][2] = -0.8485972f;
1862     U(pm).m[3][2] = 10.0f;
1863     U(pm).m[0][3] = 0.0f;
1864     U(pm).m[1][3] = 0.0f;
1865     U(pm).m[2][3] = 0.0f;
1866     U(pm).m[3][3] = 1.0f;
1867
1868     exp_scale.x = 1.0f;
1869     exp_scale.y = 1.0f;
1870     exp_scale.z = 1.0f;
1871
1872     exp_rotation.w = 0.195091f;
1873     exp_rotation.x = 0.905395f;
1874     exp_rotation.y = -0.323355f;
1875     exp_rotation.z = -0.194013f;
1876
1877     exp_translation.x = -5.0f;
1878     exp_translation.y = 0.0f;
1879     exp_translation.z = 10.0f;
1880
1881     D3DXMatrixDecompose(&got_scale, &got_rotation, &got_translation, &pm);
1882
1883     compare_scale(exp_scale, got_scale);
1884     compare_rotation(exp_rotation, got_rotation);
1885     compare_translation(exp_translation, got_translation);
1886
1887 /*_____________*/
1888
1889     U(pm).m[0][0] = 0.06554436f;
1890     U(pm).m[1][0] = -0.6873012f;
1891     U(pm).m[2][0] = 0.7234092f;
1892     U(pm).m[3][0] = -5.0f;
1893     U(pm).m[0][1] = -0.9617381f;
1894     U(pm).m[1][1] = -0.2367795f;
1895     U(pm).m[2][1] = -0.1378230f;
1896     U(pm).m[3][1] = 0.0f;
1897     U(pm).m[0][2] = 0.2660144f;
1898     U(pm).m[1][2] = -0.6866967f;
1899     U(pm).m[2][2] = -0.6765233f;
1900     U(pm).m[3][2] = 10.0f;
1901     U(pm).m[0][3] = 0.0f;
1902     U(pm).m[1][3] = 0.0f;
1903     U(pm).m[2][3] = 0.0f;
1904     U(pm).m[3][3] = 1.0f;
1905
1906     exp_scale.x = 1.0f;
1907     exp_scale.y = 1.0f;
1908     exp_scale.z = 1.0f;
1909
1910     exp_rotation.w = -0.195091f;
1911     exp_rotation.x = 0.703358f;
1912     exp_rotation.y = -0.586131f;
1913     exp_rotation.z = 0.351679f;
1914
1915     exp_translation.x = -5.0f;
1916     exp_translation.y = 0.0f;
1917     exp_translation.z = 10.0f;
1918
1919     D3DXMatrixDecompose(&got_scale, &got_rotation, &got_translation, &pm);
1920
1921     compare_scale(exp_scale, got_scale);
1922     compare_rotation(exp_rotation, got_rotation);
1923     compare_translation(exp_translation, got_translation);
1924
1925 /*_________*/
1926
1927     U(pm).m[0][0] = 7.121047f;
1928     U(pm).m[1][0] = -5.883487f;
1929     U(pm).m[2][0] = 11.81843f;
1930     U(pm).m[3][0] = -5.0f;
1931     U(pm).m[0][1] = 5.883487f;
1932     U(pm).m[1][1] = -10.60660f;
1933     U(pm).m[2][1] = -8.825232f;
1934     U(pm).m[3][1] = 0.0f;
1935     U(pm).m[0][2] = 11.81843f;
1936     U(pm).m[1][2] = 8.8252320f;
1937     U(pm).m[2][2] = -2.727645f;
1938     U(pm).m[3][2] = 2.0f;
1939     U(pm).m[0][3] = 0.0f;
1940     U(pm).m[1][3] = 0.0f;
1941     U(pm).m[2][3] = 0.0f;
1942     U(pm).m[3][3] = 1.0f;
1943
1944     exp_scale.x = 15.0f;
1945     exp_scale.y = 15.0f;
1946     exp_scale.z = 15.0f;
1947
1948     exp_rotation.w = 0.382684f;
1949     exp_rotation.x = 0.768714f;
1950     exp_rotation.y = 0.0f;
1951     exp_rotation.z = 0.512476f;
1952
1953     exp_translation.x = -5.0f;
1954     exp_translation.y = 0.0f;
1955     exp_translation.z = 2.0f;
1956
1957     D3DXMatrixDecompose(&got_scale, &got_rotation, &got_translation, &pm);
1958
1959     compare_scale(exp_scale, got_scale);
1960     compare_rotation(exp_rotation, got_rotation);
1961     compare_translation(exp_translation, got_translation);
1962
1963 /*__________*/
1964
1965     U(pm).m[0][0] = 0.0f;
1966     U(pm).m[1][0] = 4.0f;
1967     U(pm).m[2][0] = 5.0f;
1968     U(pm).m[3][0] = -5.0f;
1969     U(pm).m[0][1] = 0.0f;
1970     U(pm).m[1][1] = -10.60660f;
1971     U(pm).m[2][1] = -8.825232f;
1972     U(pm).m[3][1] = 6.0f;
1973     U(pm).m[0][2] = 0.0f;
1974     U(pm).m[1][2] = 8.8252320f;
1975     U(pm).m[2][2] = 2.727645;
1976     U(pm).m[3][2] = 3.0f;
1977     U(pm).m[0][3] = 0.0f;
1978     U(pm).m[1][3] = 0.0f;
1979     U(pm).m[2][3] = 0.0f;
1980     U(pm).m[3][3] = 1.0f;
1981
1982     hr = D3DXMatrixDecompose(&got_scale, &got_rotation, &got_translation, &pm);
1983     ok(hr == D3DERR_INVALIDCALL, "Expected D3DERR_INVALIDCALL, got %x\n", hr);
1984 }
1985
1986 static void test_Matrix_Transformation2D(void)
1987 {
1988     D3DXMATRIX exp_mat, got_mat;
1989     D3DXVECTOR2 rot_center, sca, sca_center, trans;
1990     FLOAT rot, sca_rot;
1991
1992     rot_center.x = 3.0f;
1993     rot_center.y = 4.0f;
1994
1995     sca.x = 12.0f;
1996     sca.y = -3.0f;
1997
1998     sca_center.x = 9.0f;
1999     sca_center.y = -5.0f;
2000
2001     trans.x = -6.0f;
2002     trans.y = 7.0f;
2003
2004     rot = D3DX_PI/3.0f;
2005
2006     sca_rot = 5.0f*D3DX_PI/4.0f;
2007
2008     U(exp_mat).m[0][0] = -4.245192f;
2009     U(exp_mat).m[1][0] = -0.147116f;
2010     U(exp_mat).m[2][0] = 0.0f;
2011     U(exp_mat).m[3][0] = 45.265373f;
2012     U(exp_mat).m[0][1] = 7.647113f;
2013     U(exp_mat).m[1][1] = 8.745192f;
2014     U(exp_mat).m[2][1] = 0.0f;
2015     U(exp_mat).m[3][1] = -13.401899f;
2016     U(exp_mat).m[0][2] = 0.0f;
2017     U(exp_mat).m[1][2] = 0.0f;
2018     U(exp_mat).m[2][2] = 1.0f;
2019     U(exp_mat).m[3][2] = 0.0f;
2020     U(exp_mat).m[0][3] = 0.0f;
2021     U(exp_mat).m[1][3] = 0.0f;
2022     U(exp_mat).m[2][3] = 0.0f;
2023     U(exp_mat).m[3][3] = 1.0f;
2024
2025     D3DXMatrixTransformation2D(&got_mat, &sca_center, sca_rot, &sca, &rot_center, rot, &trans);
2026
2027     expect_mat(&exp_mat, &got_mat);
2028
2029 /*_________*/
2030
2031     sca_center.x = 9.0f;
2032     sca_center.y = -5.0f;
2033
2034     trans.x = -6.0f;
2035     trans.y = 7.0f;
2036
2037     rot = D3DX_PI/3.0f;
2038
2039     sca_rot = 5.0f*D3DX_PI/4.0f;
2040
2041     U(exp_mat).m[0][0] = 0.50f;
2042     U(exp_mat).m[1][0] = -0.866025f;
2043     U(exp_mat).m[2][0] = 0.0f;
2044     U(exp_mat).m[3][0] = -6.0f;
2045     U(exp_mat).m[0][1] = 0.866025f;
2046     U(exp_mat).m[1][1] = 0.50f;
2047     U(exp_mat).m[2][1] = 0.0f;
2048     U(exp_mat).m[3][1] = 7.0f;
2049     U(exp_mat).m[0][2] = 0.0f;
2050     U(exp_mat).m[1][2] = 0.0f;
2051     U(exp_mat).m[2][2] = 1.0f;
2052     U(exp_mat).m[3][2] = 0.0f;
2053     U(exp_mat).m[0][3] = 0.0f;
2054     U(exp_mat).m[1][3] = 0.0f;
2055     U(exp_mat).m[2][3] = 0.0f;
2056     U(exp_mat).m[3][3] = 1.0f;
2057
2058     D3DXMatrixTransformation2D(&got_mat, &sca_center, sca_rot, NULL, NULL, rot, &trans);
2059
2060     expect_mat(&exp_mat, &got_mat);
2061
2062 /*_________*/
2063
2064     U(exp_mat).m[0][0] = 0.50f;
2065     U(exp_mat).m[1][0] = -0.866025f;
2066     U(exp_mat).m[2][0] = 0.0f;
2067     U(exp_mat).m[3][0] = 0.0f;
2068     U(exp_mat).m[0][1] = 0.866025f;
2069     U(exp_mat).m[1][1] = 0.50f;
2070     U(exp_mat).m[2][1] = 0.0f;
2071     U(exp_mat).m[3][1] = 0.0f;
2072     U(exp_mat).m[0][2] = 0.0f;
2073     U(exp_mat).m[1][2] = 0.0f;
2074     U(exp_mat).m[2][2] = 1.0f;
2075     U(exp_mat).m[3][2] = 0.0f;
2076     U(exp_mat).m[0][3] = 0.0f;
2077     U(exp_mat).m[1][3] = 0.0f;
2078     U(exp_mat).m[2][3] = 0.0f;
2079     U(exp_mat).m[3][3] = 1.0f;
2080
2081     D3DXMatrixTransformation2D(&got_mat, NULL, sca_rot, NULL, NULL, rot, NULL);
2082
2083     expect_mat(&exp_mat, &got_mat);
2084 }
2085
2086 static void test_D3DXVec_Array(void)
2087 {
2088     unsigned int i;
2089     D3DVIEWPORT9 viewport;
2090     D3DXMATRIX mat, projection, view, world;
2091     D3DXVECTOR4 inp_vec[ARRAY_SIZE];
2092     D3DXVECTOR4 out_vec[ARRAY_SIZE + 2];
2093     D3DXVECTOR4 exp_vec[ARRAY_SIZE + 2];
2094     D3DXPLANE inp_plane[ARRAY_SIZE];
2095     D3DXPLANE out_plane[ARRAY_SIZE + 2];
2096     D3DXPLANE exp_plane[ARRAY_SIZE + 2];
2097
2098     viewport.Width = 800; viewport.MinZ = 0.2f; viewport.X = 10;
2099     viewport.Height = 680; viewport.MaxZ = 0.9f; viewport.Y = 5;
2100
2101     for (i = 0; i < ARRAY_SIZE + 2; ++i) {
2102         out_vec[i].x = out_vec[i].y = out_vec[i].z = out_vec[i].w = 0.0f;
2103         exp_vec[i].x = exp_vec[i].y = exp_vec[i].z = exp_vec[i].w = 0.0f;
2104         out_plane[i].a = out_plane[i].b = out_plane[i].c = out_plane[i].d = 0.0f;
2105         exp_plane[i].a = exp_plane[i].b = exp_plane[i].c = exp_plane[i].d = 0.0f;
2106     }
2107
2108     for (i = 0; i < ARRAY_SIZE; ++i) {
2109         inp_plane[i].a = inp_plane[i].c = inp_vec[i].x = inp_vec[i].z = i;
2110         inp_plane[i].b = inp_plane[i].d = inp_vec[i].y = inp_vec[i].w = ARRAY_SIZE - i;
2111     }
2112
2113     U(mat).m[0][0] = 1.0f; U(mat).m[0][1] = 2.0f; U(mat).m[0][2] = 3.0f; U(mat).m[0][3] = 4.0f;
2114     U(mat).m[1][0] = 5.0f; U(mat).m[1][1] = 6.0f; U(mat).m[1][2] = 7.0f; U(mat).m[1][3] = 8.0f;
2115     U(mat).m[2][0] = 9.0f; U(mat).m[2][1] = 10.0f; U(mat).m[2][2] = 11.0f; U(mat).m[2][3] = 12.0f;
2116     U(mat).m[3][0] = 13.0f; U(mat).m[3][1] = 14.0f; U(mat).m[3][2] = 15.0f; U(mat).m[3][3] = 16.0f;
2117
2118     D3DXMatrixPerspectiveFovLH(&projection,D3DX_PI/4.0f,20.0f/17.0f,1.0f,1000.0f);
2119
2120     U(view).m[0][1] = 5.0f; U(view).m[0][2] = 7.0f; U(view).m[0][3] = 8.0f;
2121     U(view).m[1][0] = 11.0f; U(view).m[1][2] = 16.0f; U(view).m[1][3] = 33.0f;
2122     U(view).m[2][0] = 19.0f; U(view).m[2][1] = -21.0f; U(view).m[2][3] = 43.0f;
2123     U(view).m[3][0] = 2.0f; U(view).m[3][1] = 3.0f; U(view).m[3][2] = -4.0f;
2124     U(view).m[0][0] = 10.0f; U(view).m[1][1] = 20.0f; U(view).m[2][2] = 30.0f;
2125     U(view).m[3][3] = -40.0f;
2126
2127     U(world).m[0][0] = 21.0f; U(world).m[0][1] = 2.0f; U(world).m[0][2] = 3.0f; U(world).m[0][3] = 4.0;
2128     U(world).m[1][0] = 5.0f; U(world).m[1][1] = 23.0f; U(world).m[1][2] = 7.0f; U(world).m[1][3] = 8.0f;
2129     U(world).m[2][0] = -8.0f; U(world).m[2][1] = -7.0f; U(world).m[2][2] = 25.0f; U(world).m[2][3] = -5.0f;
2130     U(world).m[3][0] = -4.0f; U(world).m[3][1] = -3.0f; U(world).m[3][2] = -2.0f; U(world).m[3][3] = 27.0f;
2131
2132     /* D3DXVec2TransformCoordArray */
2133     exp_vec[1].x = 0.678571f; exp_vec[1].y = 0.785714f;
2134     exp_vec[2].x = 0.653846f; exp_vec[2].y = 0.769231f;
2135     exp_vec[3].x = 0.625f;    exp_vec[3].y = 0.75f;
2136     exp_vec[4].x = 0.590909f; exp_vec[4].y = 8.0f/11.0f;
2137     exp_vec[5].x = 0.55f;     exp_vec[5].y = 0.7f;
2138     D3DXVec2TransformCoordArray((D3DXVECTOR2*)(out_vec + 1), sizeof(D3DXVECTOR4), (D3DXVECTOR2*)inp_vec, sizeof(D3DXVECTOR4), &mat, ARRAY_SIZE);
2139     compare_vectors(exp_vec, out_vec);
2140
2141     /* D3DXVec2TransformNormalArray */
2142     exp_vec[1].x = 25.0f; exp_vec[1].y = 30.0f;
2143     exp_vec[2].x = 21.0f; exp_vec[2].y = 26.0f;
2144     exp_vec[3].x = 17.0f; exp_vec[3].y = 22.0f;
2145     exp_vec[4].x = 13.0f; exp_vec[4].y = 18.0f;
2146     exp_vec[5].x =  9.0f; exp_vec[5].y = 14.0f;
2147     D3DXVec2TransformNormalArray((D3DXVECTOR2*)(out_vec + 1), sizeof(D3DXVECTOR4), (D3DXVECTOR2*)inp_vec, sizeof(D3DXVECTOR4), &mat, ARRAY_SIZE);
2148     compare_vectors(exp_vec, out_vec);
2149
2150     /* D3DXVec3TransformCoordArray */
2151     exp_vec[1].x = 0.678571f; exp_vec[1].y = 0.785714f;  exp_vec[1].z = 0.892857f;
2152     exp_vec[2].x = 0.671875f; exp_vec[2].y = 0.78125f;   exp_vec[2].z = 0.890625f;
2153     exp_vec[3].x = 6.0f/9.0f; exp_vec[3].y = 7.0f/9.0f;  exp_vec[3].z = 8.0f/9.0f;
2154     exp_vec[4].x = 0.6625f;   exp_vec[4].y = 0.775f;     exp_vec[4].z = 0.8875f;
2155     exp_vec[5].x = 0.659091f; exp_vec[5].y = 0.772727f;  exp_vec[5].z = 0.886364f;
2156     D3DXVec3TransformCoordArray((D3DXVECTOR3*)(out_vec + 1), sizeof(D3DXVECTOR4), (D3DXVECTOR3*)inp_vec, sizeof(D3DXVECTOR4), &mat, ARRAY_SIZE);
2157     compare_vectors(exp_vec, out_vec);
2158
2159     /* D3DXVec3TransformNormalArray */
2160     exp_vec[1].x = 25.0f; exp_vec[1].y = 30.0f; exp_vec[1].z = 35.0f;
2161     exp_vec[2].x = 30.0f; exp_vec[2].y = 36.0f; exp_vec[2].z = 42.0f;
2162     exp_vec[3].x = 35.0f; exp_vec[3].y = 42.0f; exp_vec[3].z = 49.0f;
2163     exp_vec[4].x = 40.0f; exp_vec[4].y = 48.0f; exp_vec[4].z = 56.0f;
2164     exp_vec[5].x = 45.0f; exp_vec[5].y = 54.0f; exp_vec[5].z = 63.0f;
2165     D3DXVec3TransformNormalArray((D3DXVECTOR3*)(out_vec + 1), sizeof(D3DXVECTOR4), (D3DXVECTOR3*)inp_vec, sizeof(D3DXVECTOR4), &mat, ARRAY_SIZE);
2166     compare_vectors(exp_vec, out_vec);
2167
2168     /* D3DXVec3ProjectArray */
2169     exp_vec[1].x = 1089.554199f; exp_vec[1].y = -226.590622f; exp_vec[1].z = 0.215273f;
2170     exp_vec[2].x = 1068.903320f; exp_vec[2].y =  103.085129f; exp_vec[2].z = 0.183050f;
2171     exp_vec[3].x = 1051.778931f; exp_vec[3].y =  376.462250f; exp_vec[3].z = 0.156329f;
2172     exp_vec[4].x = 1037.348877f; exp_vec[4].y =  606.827393f; exp_vec[4].z = 0.133813f;
2173     exp_vec[5].x = 1025.023560f; exp_vec[5].y =  803.591248f; exp_vec[5].z = 0.114581f;
2174     D3DXVec3ProjectArray((D3DXVECTOR3*)(out_vec + 1), sizeof(D3DXVECTOR4), (CONST D3DXVECTOR3*)inp_vec, sizeof(D3DXVECTOR4), &viewport, &projection, &view, &world, ARRAY_SIZE);
2175     compare_vectors(exp_vec, out_vec);
2176
2177     /* D3DXVec3UnprojectArray */
2178     exp_vec[1].x = -6.124031f; exp_vec[1].y = 3.225360f; exp_vec[1].z = 0.620571f;
2179     exp_vec[2].x = -3.807109f; exp_vec[2].y = 2.046579f; exp_vec[2].z = 0.446894f;
2180     exp_vec[3].x = -2.922839f; exp_vec[3].y = 1.596689f; exp_vec[3].z = 0.380609f;
2181     exp_vec[4].x = -2.456225f; exp_vec[4].y = 1.359290f; exp_vec[4].z = 0.345632f;
2182     exp_vec[5].x = -2.167897f; exp_vec[5].y = 1.212597f; exp_vec[5].z = 0.324019f;
2183     D3DXVec3UnprojectArray((D3DXVECTOR3*)(out_vec + 1), sizeof(D3DXVECTOR4), (CONST D3DXVECTOR3*)inp_vec, sizeof(D3DXVECTOR4), &viewport, &projection, &view, &world, ARRAY_SIZE);
2184     compare_vectors(exp_vec, out_vec);
2185
2186     /* D3DXVec2TransformArray */
2187     exp_vec[1].x = 38.0f; exp_vec[1].y = 44.0f; exp_vec[1].z = 50.0f; exp_vec[1].w = 56.0f;
2188     exp_vec[2].x = 34.0f; exp_vec[2].y = 40.0f; exp_vec[2].z = 46.0f; exp_vec[2].w = 52.0f;
2189     exp_vec[3].x = 30.0f; exp_vec[3].y = 36.0f; exp_vec[3].z = 42.0f; exp_vec[3].w = 48.0f;
2190     exp_vec[4].x = 26.0f; exp_vec[4].y = 32.0f; exp_vec[4].z = 38.0f; exp_vec[4].w = 44.0f;
2191     exp_vec[5].x = 22.0f; exp_vec[5].y = 28.0f; exp_vec[5].z = 34.0f; exp_vec[5].w = 40.0f;
2192     D3DXVec2TransformArray(out_vec + 1, sizeof(D3DXVECTOR4), (D3DXVECTOR2*)inp_vec, sizeof(D3DXVECTOR4), &mat, ARRAY_SIZE);
2193     compare_vectors(exp_vec, out_vec);
2194
2195     /* D3DXVec3TransformArray */
2196     exp_vec[1].x = 38.0f; exp_vec[1].y = 44.0f; exp_vec[1].z = 50.0f; exp_vec[1].w = 56.0f;
2197     exp_vec[2].x = 43.0f; exp_vec[2].y = 50.0f; exp_vec[2].z = 57.0f; exp_vec[2].w = 64.0f;
2198     exp_vec[3].x = 48.0f; exp_vec[3].y = 56.0f; exp_vec[3].z = 64.0f; exp_vec[3].w = 72.0f;
2199     exp_vec[4].x = 53.0f; exp_vec[4].y = 62.0f; exp_vec[4].z = 71.0f; exp_vec[4].w = 80.0f;
2200     exp_vec[5].x = 58.0f; exp_vec[5].y = 68.0f; exp_vec[5].z = 78.0f; exp_vec[5].w = 88.0f;
2201     D3DXVec3TransformArray(out_vec + 1, sizeof(D3DXVECTOR4), (D3DXVECTOR3*)inp_vec, sizeof(D3DXVECTOR4), &mat, ARRAY_SIZE);
2202     compare_vectors(exp_vec, out_vec);
2203
2204     /* D3DXVec4TransformArray */
2205     exp_vec[1].x = 90.0f; exp_vec[1].y = 100.0f; exp_vec[1].z = 110.0f; exp_vec[1].w = 120.0f;
2206     exp_vec[2].x = 82.0f; exp_vec[2].y = 92.0f;  exp_vec[2].z = 102.0f; exp_vec[2].w = 112.0f;
2207     exp_vec[3].x = 74.0f; exp_vec[3].y = 84.0f;  exp_vec[3].z = 94.0f;  exp_vec[3].w = 104.0f;
2208     exp_vec[4].x = 66.0f; exp_vec[4].y = 76.0f;  exp_vec[4].z = 86.0f;  exp_vec[4].w = 96.0f;
2209     exp_vec[5].x = 58.0f; exp_vec[5].y = 68.0f;  exp_vec[5].z = 78.0f;  exp_vec[5].w = 88.0f;
2210     D3DXVec4TransformArray(out_vec + 1, sizeof(D3DXVECTOR4), inp_vec, sizeof(D3DXVECTOR4), &mat, ARRAY_SIZE);
2211     compare_vectors(exp_vec, out_vec);
2212
2213     /* D3DXPlaneTransformArray */
2214     exp_plane[1].a = 90.0f; exp_plane[1].b = 100.0f; exp_plane[1].c = 110.0f; exp_plane[1].d = 120.0f;
2215     exp_plane[2].a = 82.0f; exp_plane[2].b = 92.0f;  exp_plane[2].c = 102.0f; exp_plane[2].d = 112.0f;
2216     exp_plane[3].a = 74.0f; exp_plane[3].b = 84.0f;  exp_plane[3].c = 94.0f;  exp_plane[3].d = 104.0f;
2217     exp_plane[4].a = 66.0f; exp_plane[4].b = 76.0f;  exp_plane[4].c = 86.0f;  exp_plane[4].d = 96.0f;
2218     exp_plane[5].a = 58.0f; exp_plane[5].b = 68.0f;  exp_plane[5].c = 78.0f;  exp_plane[5].d = 88.0f;
2219     D3DXPlaneTransformArray(out_plane + 1, sizeof(D3DXPLANE), inp_plane, sizeof(D3DXPLANE), &mat, ARRAY_SIZE);
2220     compare_planes(exp_plane, out_plane);
2221 }
2222
2223 START_TEST(math)
2224 {
2225     D3DXColorTest();
2226     D3DXFresnelTest();
2227     D3DXMatrixTest();
2228     D3DXPlaneTest();
2229     D3DXQuaternionTest();
2230     D3DXVector2Test();
2231     D3DXVector3Test();
2232     D3DXVector4Test();
2233     test_matrix_stack();
2234     test_Matrix_AffineTransformation2D();
2235     test_Matrix_Decompose();
2236     test_Matrix_Transformation2D();
2237     test_D3DXVec_Array();
2238 }