Use the outer bound as the maximum eval radius.

1 files changed, 2 insertions, 1 deletions

diff --git a/src/qms.c b/src/qms.c
index 4be70c6..54cd657 100644
--- a/src/qms.c
+++ b/src/qms.c
@@ -243,7 +243,8 @@ static int context_init(cGH *cctkGH)
 
     qms->gh = cctkGH;
 
-    qms->max_radius = 60.0;
+    //qms->max_radius = 60.0;
+    qms->max_radius = outer_bound;
 
     ret = qms_solver_init(&qms->solver, cctkGH, basis_order_r, basis_order_z,
                           outer_bound, filter_power, 0.0, ccz4);