Export the solver functionality as an external API.

To be used for initializing the multigrid solver.

1 files changed, 3 insertions, 2 deletions

diff --git a/src/ms_solve.c b/src/ms_solve.c
index 90502be..1bf87a2 100644
--- a/src/ms_solve.c
+++ b/src/ms_solve.c
@@ -19,6 +19,7 @@
 #include "common.h"
 
 #include <errno.h>
+#include <limits.h>
 #include <math.h>
 #include <stdint.h>
 #include <stdio.h>
@@ -492,7 +493,7 @@ fail:
 int ms_solver_init(MSSolver **pctx,
                     cGH *cctkGH,
                     int basis_order_r, int basis_order_z,
-                    double sf, double filter_power, double input_filter_power)
+                    double outer_bound, double filter_power, double input_filter_power)
 {
     MSSolver *ctx;
     MSSolverPriv *s;
@@ -535,7 +536,7 @@ int ms_solver_init(MSSolver **pctx,
 
     //FIXME
     scale_factor = 1.0;
-    scale_factor = (64.0 / ctx->basis[0]->colloc_point(s->colloc_grid_order[0], ctx->nb_colloc_points[0] - 1));
+    scale_factor = (outer_bound / ctx->basis[0]->colloc_point(s->colloc_grid_order[0], ctx->nb_colloc_points[0] - 1));
     fprintf(stderr, "scale factor %16.16g\n", scale_factor);
 
     init_opencl(ctx);