From 0ce79392c39fb71f17ed5f59ba399ca5677aa7cf Mon Sep 17 00:00:00 2001 From: Elizabeth Hunt Date: Fri, 13 Oct 2023 09:58:25 -0600 Subject: [PATCH] format --- lib/lizfcm.a | Bin 14128 -> 14128 bytes test/main.c | 114 +++++++++++++++++++++++++-------------------------- 2 files changed, 56 insertions(+), 58 deletions(-) diff --git a/lib/lizfcm.a b/lib/lizfcm.a index 47bbc82c6ce73ecb42787a52b330983cc2e889a2..78ca79eefc769195e8b9beb2aafcacaeaa03fd2a 100644 GIT binary patch delta 77 zcmdmxw;^wW44aXmk(s&iMx`nyFw@Lp^CG4$b{N};%Ze4oF5-`3hOsXRT@rw?S1Gp% F0{}Kr6z2c{ delta 77 zcmdmxw;^wW44Z+OsgaS%Mx`nyFw@X-^CG4$b{N};%Ze4oF5-`3hOsXRT@rw?S1Gp% F0{}F~6x#p* diff --git a/test/main.c b/test/main.c index ac7fba3..75e60ee 100644 --- a/test/main.c +++ b/test/main.c @@ -5,71 +5,69 @@ double f(double x) { return (x - 1) / (x + 1); } int main() { - printf("Basic Routines\n"); - printf("smaceps(): %.10e\n", smaceps()); - printf("dmaceps(): %.10e\n", dmaceps()); - printf("========\n"); - - Array_double *v = InitArray(double, {3, 1, -4, 1, 5, -9, 3}); - char v_s[256]; - strcpy(v_s, ""); - format_vector_into(v, v_s); + printf("Basic Routines\n"); + printf("smaceps(): %.10e\n", smaceps()); + printf("dmaceps(): %.10e\n", dmaceps()); + printf("========\n"); - Array_double *w = InitArray(double, {-2, 7, 1, -8, -2, 8, 5}); - char w_s[256]; - strcpy(w_s, ""); - format_vector_into(w, w_s); - - printf("Norm, Distance\n"); - printf("v: %s", v_s); - printf("w: %s", 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)); - printf("========\n"); - - double h = 0.001; - printf("Derivative Approxs\n"); - printf("f(x) = (x-1)/(x+1)\n"); - 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)); - printf("========\n"); - - v = InitArray(double, {1, 2, 3, 4, 5}); - strcpy(v_s, ""); - format_vector_into(v, v_s); - w = InitArray(double, {2, 3, 4, 5, 6}); - strcpy(w_s, ""); - format_vector_into(w, w_s); - Line *line = least_squares_lin_reg(v, w); - printf("Least Squares\n"); - printf("v: %s", v_s); - printf("w: %s", w_s); - 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}); - strcpy(v_s, ""); - format_vector_into(v, v_s); - w = InitArray(double, {0.5, 3, 2, 3.5, 5, 6, 7.5}); - strcpy(w_s, ""); - format_vector_into(w, w_s); - printf("v: %s", v_s); - printf("w: %s", 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("========\n"); + Array_double *v = InitArray(double, {3, 1, -4, 1, 5, -9, 3}); + char v_s[256]; + strcpy(v_s, ""); + format_vector_into(v, v_s); + + Array_double *w = InitArray(double, {-2, 7, 1, -8, -2, 8, 5}); + char w_s[256]; + strcpy(w_s, ""); + format_vector_into(w, w_s); + + printf("Norm, Distance\n"); + printf("v: %s", v_s); + printf("w: %s", 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)); + printf("========\n"); + + double h = 0.001; + printf("Derivative Approxs\n"); + printf("f(x) = (x-1)/(x+1)\n"); + 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)); + printf("========\n"); + + v = InitArray(double, {1, 2, 3, 4, 5}); + strcpy(v_s, ""); + format_vector_into(v, v_s); + w = InitArray(double, {2, 3, 4, 5, 6}); + strcpy(w_s, ""); + format_vector_into(w, w_s); + Line *line = least_squares_lin_reg(v, w); + printf("Least Squares\n"); + printf("v: %s", v_s); + printf("w: %s", w_s); + 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}); + strcpy(v_s, ""); + format_vector_into(v, v_s); + w = InitArray(double, {0.5, 3, 2, 3.5, 5, 6, 7.5}); + strcpy(w_s, ""); + format_vector_into(w, w_s); + printf("v: %s", v_s); + printf("w: %s", 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("========\n"); - printf("LU Decomp\n"); char m_s[2048]; Matrix_double *m = InitMatrixWithSize(double, 8, 8, 0.0); for (int i = 0; i < 8; i++) { - for (int j = 0; j < 8; j++) { + 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("m = %s", m_s);