/*
 * Copyright (c), Recep Aslantas.
 *
 * MIT License (MIT), http://opensource.org/licenses/MIT
 * Full license can be found in the LICENSE file
 */

#include "test_common.h"

#define m 3
#define n 3

TEST_IMPL(mat3_identity) {
  mat3 m1 = GLM_MAT3_IDENTITY_INIT;
  mat3 m2 = GLM_MAT3_IDENTITY_INIT;
  mat3 m3;
  int  i, j;
  
  /* test identity matrix multiplication */
  glmc_mat3_mul(m1, m2, m3);
  for (i = 0; i < m; i++) {
    for (j = 0; j < n; j++) {
      if (i == j) {
        ASSERT(glm_eq(m3[i][j], 1.0f))
      } else {
        ASSERT(glm_eq(m3[i][j], 0.0f))
      }
    }
  }
  
  TEST_SUCCESS
}

TEST_IMPL(mat3_mul) {
  mat3 m1 = GLM_MAT3_IDENTITY_INIT;
  mat3 m2 = GLM_MAT3_IDENTITY_INIT;
  mat3 m3;
  mat3 m4 = GLM_MAT3_ZERO_INIT;
  int  i, j, k;

  /* test random matrices */
  /* random matrices */
  test_rand_mat3(m1);
  test_rand_mat3(m2);

  glmc_mat3_mul(m1, m2, m3);
  for (i = 0; i < m; i++) {
    for (j = 0; j < n; j++) {
      for (k = 0; k < m; k++)
        /* column-major */
        m4[i][j] += m1[k][j] * m2[i][k];
    }
  }
  
  ASSERTIFY(test_assert_mat3_eq(m3, m4))
  
  TEST_SUCCESS
}

TEST_IMPL(mat3_inv) {
  mat3 m1 = GLM_MAT3_IDENTITY_INIT;
  mat3 m2 = GLM_MAT3_IDENTITY_INIT;
  mat3 m3;
  int  i;
  
  for (i = 0; i < 100000; i++) {
    test_rand_mat3(m1);
    test_rand_mat3(m2);

    /* test inverse precise */
    glmc_mat3_inv(m1, m2);
    glmc_mat3_inv(m2, m3);
    ASSERTIFY(test_assert_mat3_eq(m1, m3))
  }
  
  TEST_SUCCESS
}