/** * Test matrix subtraction. */ void test_m_subtract() { precision_t data_A[][2] = { { 1, 0 }, { 2, 4 } }; precision_t data_B[][2] = { { 5, 9 }, { 2, 1 } }; matrix_t *A = m_initialize(2, 2); matrix_t *B = m_initialize(2, 2); fill_matrix_data(A, data_A); fill_matrix_data(B, data_B); printf("A = \n"); m_fprint(stdout, A); printf("B = \n"); m_fprint(stdout, B); m_subtract(A, B); printf("A - B = \n"); m_fprint(stdout, A); m_free(A); m_free(B); }
void subtract_matrices(matrix_t *m1, matrix_t *m2, matrix_t *m_result) { if (m_subtract(m1, m2, m_result) != SUCCESS) printf("Could not subtract matrices. Check sizes\n"); else print_matrix("Result:", m_result); }