add framwork

src/init.c

@@ -3,6 +3,13 @@
 #include <gsl/gsl_matrix.h>
 #include "common.h"
 
+typedef struct {
+  int capacity;
+  int size;
+  double *x;
+  double *y;
+} SampleData;
+
 // ODE
 // du/dphi = v
 // dv/dphi = 3Mu^2 - u
@@ -92,3 +99,9 @@ double chi(double b){
   }
   return find_root_5ord(y_previous[0], y[0], y_previous[1], y[1], phi_previous, phi);
 }
+
+static void refine_interval(double xl, double xr, double yl, double yr, SampleData *sample){}
+
+static int init(double bmax, double rela_err_limit){
+  return 0;
+}