MatrixTests.c 6.2 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223
  1. #include "../Tests.h"
  2. #include "core/Matrix.h"
  3. typedef CoreVector3 V3;
  4. #define CV3(a, b, c) (&(V3){{a, b, c}})
  5. #define CV30 CV3(0.0f, 0.0f, 0.0f)
  6. static void testInit() {
  7. Matrix m = UNIT_MATRIX;
  8. const float* data = (float*)&m;
  9. for(int i = 0; i < 16; i++) {
  10. int x = i % 4;
  11. int y = i / 4;
  12. TEST_FLOAT(x == y, data[i], 0.0f);
  13. }
  14. }
  15. static void testTranspose() {
  16. Matrix m;
  17. float* data = (float*)&m;
  18. for(int i = 0; i < 16; i++) {
  19. data[i] = (float)(i + 1);
  20. }
  21. Matrix t = m;
  22. transposeMatrix(&t);
  23. Matrix m2 = t;
  24. transposeMatrix(&m2);
  25. const float* mp = (float*)&m;
  26. const float* tp = (float*)&t;
  27. for(int x = 0; x < 4; x++) {
  28. for(int y = 0; y < 4; y++) {
  29. TEST_FLOAT(mp[y * 4 + x], tp[x * 4 + y], 0.0f);
  30. }
  31. }
  32. const float* mp2 = (float*)&m2;
  33. for(int i = 0; i < 16; i++) {
  34. TEST_FLOAT(mp[i], mp2[i], 0.0f);
  35. }
  36. }
  37. static void testScale() {
  38. Matrix m = UNIT_MATRIX;
  39. scaleMatrix(&m, CV3(2.0f, 3.0f, 4.0f));
  40. TEST_V3(CV3(-8.0f, 18.0f, 28.0f),
  41. mulMatrixV3(CV30, &m, CV3(-4.0f, 6.0f, 7.0f)));
  42. }
  43. static void testUniformScale() {
  44. Matrix m = UNIT_MATRIX;
  45. scaleMatrixF(&m, 2.0f);
  46. TEST_V3(CV3(-8.0f, 12.0f, 14.0f),
  47. mulMatrixV3(CV30, &m, CV3(-4.0f, 6.0f, 7.0f)));
  48. }
  49. static void testTranslateX() {
  50. Matrix m = UNIT_MATRIX;
  51. translateMatrixX(&m, 5.0f);
  52. TEST_V3(CV3(1.0f, 6.0f, 7.0f),
  53. mulMatrixV3(CV30, &m, CV3(-4.0f, 6.0f, 7.0f)));
  54. }
  55. static void testTranslateY() {
  56. Matrix m = UNIT_MATRIX;
  57. translateMatrixY(&m, 6.0f);
  58. TEST_V3(CV3(-4.0f, 12.0f, 7.0f),
  59. mulMatrixV3(CV30, &m, CV3(-4.0f, 6.0f, 7.0f)));
  60. }
  61. static void testTranslateZ() {
  62. Matrix m = UNIT_MATRIX;
  63. translateMatrixZ(&m, 7.0f);
  64. TEST_V3(CV3(-4.0f, 6.0f, 14.0f),
  65. mulMatrixV3(CV30, &m, CV3(-4.0f, 6.0f, 7.0f)));
  66. }
  67. static void testTranslate() {
  68. Matrix m = UNIT_MATRIX;
  69. translateMatrix(&m, CV3(1.0f, 2.0f, 3.0f));
  70. TEST_V3(CV3(-3.0f, 8.0f, 10.0f),
  71. mulMatrixV3(CV30, &m, CV3(-4.0f, 6.0f, 7.0f)));
  72. }
  73. static void testTranslateTo() {
  74. char buffer[1024];
  75. Matrix m;
  76. for(int i = 0; i < 16; i++) {
  77. ((float*)&m)[i] = (float)i + 1.0f;
  78. }
  79. translateMatrixTo(&m, CV3(6.0f, 8.0f, 9.0f));
  80. toStringMatrix(&m, buffer, sizeof(buffer));
  81. TEST_STRING("[[1.000, 0.000, 0.000, 6.000], [0.000, 1.000, 0.000, 8.000], "
  82. "[0.000, 0.000, 1.000, 9.000], [0.000, 0.000, 0.000, 1.000]]",
  83. buffer);
  84. }
  85. static void testCombination() {
  86. Matrix m = UNIT_MATRIX;
  87. scaleMatrixF(&m, 2.0f);
  88. translateMatrixX(&m, 1.0f);
  89. translateMatrixY(&m, 2.0f);
  90. translateMatrixZ(&m, 3.0f);
  91. translateMatrix(&m, CV3(-4.0f, 2.0f, 3.0f));
  92. scaleMatrix(&m, CV3(2.0f, 3.0f, 4.0f));
  93. scaleMatrixF(&m, 0.5f);
  94. TEST_V3(CV3(-1.0f, 9.0f, 16.0f),
  95. mulMatrixV3(CV30, &m, CV3(1.0f, 1.0f, 1.0f)));
  96. }
  97. static void testMatrixCombination() {
  98. Matrix a = UNIT_MATRIX;
  99. scaleMatrixF(&a, 2.0f);
  100. translateMatrix(&a, CV3(1.0f, 2.0f, 3.0f));
  101. Matrix b = UNIT_MATRIX;
  102. scaleMatrixF(&b, 3.0f);
  103. translateMatrix(&b, CV3(1.0f, 1.0f, 1.0f));
  104. Matrix c = UNIT_MATRIX;
  105. translateMatrix(&c, CV3(-1.0f, -2.0f, -3.0f));
  106. mulSetMatrix(&c, mulMatrix(&ZERO_MATRIX, &b, &a));
  107. TEST_V3(CV3(9.0f, 11.0f, 13.0f),
  108. mulMatrixV3(CV30, &c, CV3(1.0f, 1.0f, 1.0f)));
  109. }
  110. static void testRotateX() {
  111. Matrix m = UNIT_MATRIX;
  112. rotateMatrixX(&m, 90.0f);
  113. TEST_V3(CV3(1.0f, 0.0f, 0.0f),
  114. mulMatrixV3(CV30, &m, CV3(1.0f, 0.0f, 0.0f)));
  115. TEST_V3(CV3(0.0f, 0.0f, 1.0f),
  116. mulMatrixV3(CV30, &m, CV3(0.0f, 1.0f, 0.0f)));
  117. TEST_V3(CV3(0.0f, -1.0f, 0.0f),
  118. mulMatrixV3(CV30, &m, CV3(0.0f, 0.0f, 1.0f)));
  119. }
  120. static void testRotateY() {
  121. Matrix m = UNIT_MATRIX;
  122. rotateMatrixY(&m, 90.0f);
  123. TEST_V3(CV3(0.0f, 0.0f, -1.0f),
  124. mulMatrixV3(CV30, &m, CV3(1.0f, 0.0f, 0.0f)));
  125. TEST_V3(CV3(0.0f, 1.0f, 0.0f),
  126. mulMatrixV3(CV30, &m, CV3(0.0f, 1.0f, 0.0f)));
  127. TEST_V3(CV3(1.0f, 0.0f, 0.0f),
  128. mulMatrixV3(CV30, &m, CV3(0.0f, 0.0f, 1.0f)));
  129. }
  130. static void testRotateZ() {
  131. Matrix m = UNIT_MATRIX;
  132. rotateMatrixZ(&m, 90.0f);
  133. TEST_V3(CV3(0.0f, 1.0f, 0.0f),
  134. mulMatrixV3(CV30, &m, CV3(1.0f, 0.0f, 0.0f)));
  135. TEST_V3(CV3(-1.0f, 0.0f, 0.0f),
  136. mulMatrixV3(CV30, &m, CV3(0.0f, 1.0f, 0.0f)));
  137. TEST_V3(CV3(0.0f, 0.0f, 1.0f),
  138. mulMatrixV3(CV30, &m, CV3(0.0f, 0.0f, 1.0f)));
  139. }
  140. static void testQuaternionMatrix() {
  141. Quaternion q1 = UNIT_QUATERNION;
  142. axisAngleQ(&q1, CV3(1.0f, 0.0f, 0.0f), 48.0f);
  143. Quaternion q2 = UNIT_QUATERNION;
  144. axisAngleQ(&q2, CV3(0.0f, 1.0f, 0.0f), 52.0f);
  145. Quaternion q3 = UNIT_QUATERNION;
  146. axisAngleQ(&q3, CV3(0.0f, 0.0f, 1.0f), 60.0f);
  147. Matrix m = UNIT_MATRIX;
  148. translateMatrix(&m, CV3(1.0f, 2.0f, 3.0f));
  149. rotateMatrix(&m, &q1);
  150. rotateMatrix(&m, &q2);
  151. rotateMatrix(&m, &q3);
  152. translateMatrix(&m, CV3(1.0f, 2.0f, 3.0f));
  153. Matrix check = UNIT_MATRIX;
  154. translateMatrix(&check, CV3(1.0f, 2.0f, 3.0f));
  155. rotateMatrixX(&check, 48.0f);
  156. rotateMatrixY(&check, 52.0f);
  157. rotateMatrixZ(&check, 60.0f);
  158. translateMatrix(&check, CV3(1.0f, 2.0f, 3.0f));
  159. for(int i = 0; i < 16; i++) {
  160. TEST_FLOAT(((float*)&check)[i], ((float*)&m)[i], 0.0001f);
  161. }
  162. }
  163. static void testToString() {
  164. Matrix m;
  165. for(int i = 0; i < 16; i++) {
  166. ((float*)&m)[i] = (float)i + 1.0f;
  167. }
  168. char buffer[1024];
  169. size_t n = toStringMatrix(&m, buffer, sizeof(buffer));
  170. TEST_SIZE(127, n);
  171. TEST_STRING(
  172. "[[1.000, 2.000, 3.000, 4.000], [5.000, 6.000, 7.000, 8.000], "
  173. "[9.000, 10.000, 11.000, 12.000], [13.000, 14.000, 15.000, 16.000]]",
  174. buffer);
  175. n = toStringMatrix(&m, buffer, 20);
  176. TEST_SIZE(127, n);
  177. TEST_STRING("[[1.000, 2.000, 3.0", buffer);
  178. }
  179. void testMatrix() {
  180. testInit();
  181. testTranspose();
  182. testScale();
  183. testUniformScale();
  184. testTranslateX();
  185. testTranslateY();
  186. testTranslateZ();
  187. testTranslate();
  188. testTranslateTo();
  189. testCombination();
  190. testMatrixCombination();
  191. testRotateX();
  192. testRotateY();
  193. testRotateZ();
  194. testQuaternionMatrix();
  195. testToString();
  196. }