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