#include "tests.h" #include #include "matrix.h" void test_LR() { Matrix *A = matrix_from_str("[2 3 1;4 1 -5;-1 2 3]"); Matrix *b = matrix_from_str("[-1;2;3]"); Matrix **LR = matrix_LR(A, b); Matrix *L = LR[0]; Matrix *R = LR[1]; Matrix *y = matrix_forwardel(L, b); Matrix *x = matrix_backsubst(R, y); assert(matrix_eq(matrix_mult(A, x), b, 0.01)); } void run_tests() { test_LR(); }