void matrix_to_entity(MATRIX *m) {
int depth, cur, maxdepth = MIN(m->cols, m->rows) ;
for(depth = 0; depth < maxdepth; depth++) {
for (cur = depth+1; cur < maxdepth; cur++) {
FRAC *piv = matrix_get_val(m, depth, depth) ;
FRAC *tz = matrix_get_val(m, cur, depth) ;
FRAC fac ;
if(tz->num != 0) {
frac_div(&fac, piv, tz) ;
frac_negate(&fac, &fac) ;
matrix_multiply_row(m, cur, &fac) ;
matrix_add_row(m, cur, depth) ;
}
}
}
for(depth = maxdepth - 1; depth >= 0; depth--) {
for (cur = depth-1; cur >= 0; cur--) {
FRAC *piv = matrix_get_val(m, depth, depth) ;
FRAC *tz = matrix_get_val(m, cur, depth) ;
FRAC fac ;
if(tz->num != 0) {
frac_div(&fac, piv, tz) ;
frac_negate(&fac, &fac) ;
matrix_multiply_row(m, cur, &fac) ;
matrix_add_row(m, cur, depth) ;
}
}
}
for(depth = 0; depth < maxdepth; depth++) {
FRAC t, *e = matrix_get_val(m, depth, depth) ;
if (e->num != 0) {
frac_inverse(&t,e) ;
matrix_multiply_row(m, depth, &t) ;
}
}
}
static void test_fractions(CuTest *tc) {
variant a, b;
a = frac_make(120, 12000);
CuAssertIntEquals(tc, 1, a.sa[0]);
CuAssertIntEquals(tc, 100, a.sa[1]);
b = frac_make(23, 2300);
a = frac_add(a, b);
CuAssertIntEquals(tc, 1, a.sa[0]);
CuAssertIntEquals(tc, 50, a.sa[1]);
a = frac_mul(a, b);
CuAssertIntEquals(tc, 1, a.sa[0]);
CuAssertIntEquals(tc, 5000, a.sa[1]);
a = frac_div(b, b);
CuAssertIntEquals(tc, 1, a.sa[0]);
CuAssertIntEquals(tc, 1, a.sa[1]);
a = frac_sub(a, a);
CuAssertIntEquals(tc, 0, a.sa[0]);
a = frac_sub(frac_one, a);
CuAssertIntEquals(tc, 1, a.sa[0]);
CuAssertIntEquals(tc, 1, a.sa[1]);
a = frac_mul(a, frac_zero);
CuAssertIntEquals(tc, 0, a.sa[0]);
CuAssertIntEquals(tc, 1, frac_sign(frac_make(-1, -1)));
CuAssertIntEquals(tc, 1, frac_sign(frac_make(1, 1)));
CuAssertIntEquals(tc, -1, frac_sign(frac_make(-1, 1)));
CuAssertIntEquals(tc, -1, frac_sign(frac_make(1, -1)));
CuAssertIntEquals(tc, 0, frac_sign(frac_make(0, 1)));
/* we reduce large integers by calculating the gcd */
a = frac_make(480000, 3000);
CuAssertIntEquals(tc, 160, a.sa[0]);
CuAssertIntEquals(tc, 1, a.sa[1]);
/* if num is too big for a short, and the gcd is 1, we cheat: */
a = frac_make(480001, 3000);
CuAssertIntEquals(tc, 32000, a.sa[0]);
CuAssertIntEquals(tc, 200, a.sa[1]);
}
