From c2ef7737536bf1f8c81fcfae95c0183b21c9753f Mon Sep 17 00:00:00 2001 From: rtk0c Date: Fri, 3 Jun 2022 23:25:43 -0700 Subject: Changeset: 62 Branch comment: [] [WIP] Initial migration --- .../source/test/gtx/gtx_matrix_factorisation.cpp | 105 --------------------- 1 file changed, 105 deletions(-) delete mode 100644 3rdparty/glm/source/test/gtx/gtx_matrix_factorisation.cpp (limited to '3rdparty/glm/source/test/gtx/gtx_matrix_factorisation.cpp') diff --git a/3rdparty/glm/source/test/gtx/gtx_matrix_factorisation.cpp b/3rdparty/glm/source/test/gtx/gtx_matrix_factorisation.cpp deleted file mode 100644 index 6771dba..0000000 --- a/3rdparty/glm/source/test/gtx/gtx_matrix_factorisation.cpp +++ /dev/null @@ -1,105 +0,0 @@ -#define GLM_ENABLE_EXPERIMENTAL -#include -#include -#include - -template -int test_qr(glm::mat m) -{ - int Error = 0; - - T const epsilon = static_cast(1e-10); - - glm::mat<(C < R ? C : R), R, T, Q> q(-999); - glm::mat r(-999); - - glm::qr_decompose(m, q, r); - - //Test if q*r really equals the input matrix - glm::mat tm = q*r; - glm::mat err = tm - m; - - for (glm::length_t i = 0; i < C; i++) - for (glm::length_t j = 0; j < R; j++) - Error += glm::abs(err[i][j]) > epsilon ? 1 : 0; - - //Test if the columns of q are orthonormal - for (glm::length_t i = 0; i < (C < R ? C : R); i++) - { - Error += (length(q[i]) - 1) > epsilon ? 1 : 0; - - for (glm::length_t j = 0; j epsilon ? 1 : 0; - } - - //Test if the matrix r is upper triangular - for (glm::length_t i = 0; i < C; i++) - for (glm::length_t j = i + 1; j < (C < R ? C : R); j++) - Error += glm::epsilonEqual(r[i][j], static_cast(0), glm::epsilon()) ? 0 : 1; - - return Error; -} - -template -int test_rq(glm::mat m) -{ - int Error = 0; - - T const epsilon = static_cast(1e-10); - - glm::mat q(-999); - glm::mat<(C < R ? C : R), R, T, Q> r(-999); - - glm::rq_decompose(m, r, q); - - //Test if q*r really equals the input matrix - glm::mat tm = r*q; - glm::mat err = tm - m; - - for (glm::length_t i = 0; i < C; i++) - for (glm::length_t j = 0; j < R; j++) - Error += glm::abs(err[i][j]) > epsilon ? 1 : 0; - - //Test if the rows of q are orthonormal - glm::mat<(C < R ? C : R), C, T, Q> tq = transpose(q); - - for (glm::length_t i = 0; i < (C < R ? C : R); i++) - { - Error += (length(tq[i]) - 1) > epsilon ? 1 : 0; - - for (glm::length_t j = 0; j epsilon ? 1 : 0; - } - - //Test if the matrix r is upper triangular - for (glm::length_t i = 0; i < (C < R ? C : R); i++) - for (glm::length_t j = R - (C < R ? C : R) + i + 1; j < R; j++) - Error += glm::epsilonEqual(r[i][j], static_cast(0), glm::epsilon()) ? 0 : 1; - - return Error; -} - -int main() -{ - int Error = 0; - - //Test QR square - Error += test_qr(glm::dmat3(12.0, 6.0, -4.0, -51.0, 167.0, 24.0, 4.0, -68.0, -41.0)) ? 1 : 0; - - //Test RQ square - Error += test_rq(glm::dmat3(12.0, 6.0, -4.0, -51.0, 167.0, 24.0, 4.0, -68.0, -41.0)) ? 1 : 0; - - //Test QR triangular 1 - Error += test_qr(glm::dmat3x4(12.0, 6.0, -4.0, -51.0, 167.0, 24.0, 4.0, -68.0, -41.0, 7.0, 2.0, 15.0)) ? 1 : 0; - - //Test QR triangular 2 - Error += test_qr(glm::dmat4x3(12.0, 6.0, -4.0, -51.0, 167.0, 24.0, 4.0, -68.0, -41.0, 7.0, 2.0, 15.0)) ? 1 : 0; - - //Test RQ triangular 1 : Fails at the triangular test - Error += test_rq(glm::dmat3x4(12.0, 6.0, -4.0, -51.0, 167.0, 24.0, 4.0, -68.0, -41.0, 7.0, 2.0, 15.0)) ? 1 : 0; - - //Test QR triangular 2 - Error += test_rq(glm::dmat4x3(12.0, 6.0, -4.0, -51.0, 167.0, 24.0, 4.0, -68.0, -41.0, 7.0, 2.0, 15.0)) ? 1 : 0; - - return Error; -} -- cgit v1.2.3-70-g09d2