From 62020d51e5b966f31b6bc8aa0f6ee073b8f0984b Mon Sep 17 00:00:00 2001
From: TheOnePath <ethan.sc@closedless.xyz>
Date: Wed, 11 Sep 2024 14:23:57 +0100
Subject: [PATCH] Removed functions

---
 src/demo.c | 51 ++++-----------------------------------------------
 1 file changed, 4 insertions(+), 47 deletions(-)

diff --git a/src/demo.c b/src/demo.c
index 54fc9e1..64f5184 100644
--- a/src/demo.c
+++ b/src/demo.c
@@ -6,56 +6,12 @@
 #include "kernels.h"
 #include "solve.h"
 #include "gp.h"
+#include "linalg.h"
+#include "matio.h"
 
 #define RNG_MAX 10
 #define X_MAX 10
 
-int linspace ( double start, double end, int length, gsl_matrix *out ) {
-  if ( out->size2 != length )
-    return 1;
-
-  double by = end / length;
-  for ( int i = 0 ; i < length ; i++ )
-    (void)gsl_matrix_set(out, 0, i, i*by);
-
-  return 0;
-}
-
-void print_matrix ( gsl_matrix *matrix ) {
-  printf("\n");
-  for ( int i = 0 ; i < matrix->size1 ; i++ ) {
-    for ( int j = 0 ; j < matrix->size2 ; j++ ) {
-      printf(" %.6e ", gsl_matrix_get(matrix, i, j));
-    }
-    printf("\n");
-  }
-  printf("\n");
-}
-
-void print_flatten ( gsl_matrix *matrix ) {
-  printf("\n");
-  if ( matrix->size1 == 1 ) {
-    for ( int i = 0 ; i < matrix->size2 ; i++ ) {
-      if ( i == matrix->size2-1 ) {
-        printf("%g\n", gsl_matrix_get(matrix, 0, i));
-        return;
-      }
-      printf("%g, ", gsl_matrix_get(matrix, 0, i));
-    }
-    printf("\n");
-    return;
-  }
-
-  for ( int i = 0 ; i < matrix->size1 ; i++ ) {
-    if ( i == matrix->size1-1 ) {
-      printf("%g\n", gsl_matrix_get(matrix, i, 0));
-      return;
-    }
-    printf("%g, ", gsl_matrix_get(matrix, i, 0));
-  }
-  printf("\n");
-}
-
 
 int main(void)
 {
@@ -83,10 +39,11 @@ int main(void)
   };
   gp_model model = gpfit(X, Y, Xs, info, noise, &inv_chol);
 
-  print_flatten(X);
+/*  print_flatten(X);
   print_flatten(Y);
   print_flatten(Xs);
   print_flatten(model.meanf);
+*/
 
   /* free all allocations made */
   gsl_rng_free(rng);