diff options
Diffstat (limited to 'test/main.c')
| -rw-r--r-- | test/main.c | 106 |
1 files changed, 61 insertions, 45 deletions
diff --git a/test/main.c b/test/main.c index 798bfbf..3551af2 100644 --- a/test/main.c +++ b/test/main.c @@ -5,56 +5,72 @@ double f(double x) { return (x - 1) / (x + 1); } int main() { - printf("smaceps(): %.10e\n", smaceps()); - printf("dmaceps(): %.10e\n", dmaceps()); - - Array_double *v = InitArray(double, {3, 1, -4, 1, 5, -9, 3}); - Array_double *w = InitArray(double, {-2, 7, 1, -8, -2, 8, 5}); - - char v_s[256]; - char w_s[256]; - format_vector_into(v, v_s); - format_vector_into(w, w_s); - - printf("v: %s\n", v_s); - printf("w: %s\n", w_s); - printf("l1_norm(v): %f\n", l1_norm(v)); - printf("l2_norm(v): %f\n", l2_norm(v)); - printf("linf_norm(v): %f\n", linf_norm(v)); - - printf("l1_dist(v, w): %f\n", l1_distance(v, w)); - printf("l2_dist(v, w): %f\n", l2_distance(v, w)); - printf("linf_dist(v, w): %f\n", linf_distance(v, w)); - - double h = 0.001; - printf("approx f'(1) w/ c.d.: %f\n", central_derivative_at(&f, 1, h)); - printf("approx f'(1) w/ fw.d.: %f\n", forward_derivative_at(&f, 1, h)); - printf("approx f'(1) w/ bw.d.: %f\n", backward_derivative_at(&f, 1, h)); - - v = InitArray(double, {1, 2, 3, 4, 5}); - w = InitArray(double, {2, 3, 4, 5, 6}); - format_vector_into(v, v_s); - format_vector_into(w, w_s); - printf("v: %s\n", v_s); - printf("w: %s\n", w_s); - - Line *line = least_squares_lin_reg(v, w); - printf("least_squares_lin_reg(v, w): (%f)x + %f\n", line->m, line->a); - - v = InitArray(double, {1, 2, 3, 4, 5, 6, 7}); - w = InitArray(double, {0.5, 3, 2, 3.5, 5, 6, 7.5}); - format_vector_into(v, v_s); - format_vector_into(w, w_s); - printf("v: %s\n", v_s); - printf("w: %s\n", w_s); - line = least_squares_lin_reg(v, w); - printf("least_squares_lin_reg(v, w): (%f)x + %f\n", line->m, line->a); + // printf("smaceps(): %.10e\n", smaceps()); + // printf("dmaceps(): %.10e\n", dmaceps()); + // + // Array_double *v = InitArray(double, {3, 1, -4, 1, 5, -9, 3}); + // Array_double *w = InitArray(double, {-2, 7, 1, -8, -2, 8, 5}); + // + // char v_s[256]; + // char w_s[256]; + // format_vector_into(v, v_s); + // format_vector_into(w, w_s); + // + // printf("v: %s\n", v_s); + // printf("w: %s\n", w_s); + // printf("l1_norm(v): %f\n", l1_norm(v)); + // printf("l2_norm(v): %f\n", l2_norm(v)); + // printf("linf_norm(v): %f\n", linf_norm(v)); + // + // printf("l1_dist(v, w): %f\n", l1_distance(v, w)); + // printf("l2_dist(v, w): %f\n", l2_distance(v, w)); + // printf("linf_dist(v, w): %f\n", linf_distance(v, w)); + // + // double h = 0.001; + // printf("approx f'(1) w/ c.d.: %f\n", central_derivative_at(&f, 1, h)); + // printf("approx f'(1) w/ fw.d.: %f\n", forward_derivative_at(&f, 1, h)); + // printf("approx f'(1) w/ bw.d.: %f\n", backward_derivative_at(&f, 1, h)); + // + // v = InitArray(double, {1, 2, 3, 4, 5}); + // w = InitArray(double, {2, 3, 4, 5, 6}); + // format_vector_into(v, v_s); + // format_vector_into(w, w_s); + // printf("v: %s\n", v_s); + // printf("w: %s\n", w_s); + // + // Line *line = least_squares_lin_reg(v, w); + // printf("least_squares_lin_reg(v, w): (%f)x + %f\n", line->m, line->a); + // + // v = InitArray(double, {1, 2, 3, 4, 5, 6, 7}); + // w = InitArray(double, {0.5, 3, 2, 3.5, 5, 6, 7.5}); + // + // format_vector_into(v, v_s); + // format_vector_into(w, w_s); + // printf("v: %s\n", v_s); + // printf("w: %s\n", w_s); + // line = least_squares_lin_reg(v, w); + // printf("least_squares_lin_reg(v, w): (%f)x + %f\n", line->m, line->a); char m_s[256]; Matrix_double *m = InitMatrixWithSize(double, 8, 8, 0.0); put_identity_diagonal(m); + for (int i = 0; i < 8; i++) { + for (int j = 0; j < 8; j++) { + m->data[i]->data[j] = (i + 1.0) + j * 3 + (rand() % 12); + } + } + format_matrix_into(m, m_s); - printf("%s", m_s); + printf("%s\n", m_s); + + Matrix_double **u_l = put_lu_decomp(m); + Matrix_double *u = u_l[0]; + Matrix_double *l = u_l[1]; + + format_matrix_into(u, m_s); + printf("%s\n", m_s); + format_matrix_into(l, m_s); + printf("%s\n", m_s); return 0; } |