lu factorization

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;
 }