Make the basis scaling factor configurable.

2 files changed, 10 insertions, 1 deletions

diff --git a/param.ccl b/param.ccl
index b490a37..64f2063 100644
--- a/param.ccl
+++ b/param.ccl
@@ -18,3 +18,8 @@ CCTK_INT cheb_order "order"
 {
     1: :: ""
 } 20
+
+CCTK_REAL scale_factor "Scaling factor L for the SB basis"
+{
+    0: :: ""
+} 3.0
diff --git a/src/brill.c b/src/brill.c
index be9c7f3..48912c4 100644
--- a/src/brill.c
+++ b/src/brill.c
@@ -45,7 +45,9 @@ typedef struct BasisSet {
  * They are symmetric wrt beginning and decay as 1/x in infinity.
  */
 
-#define SCALE_FACTOR 1
+static CCTK_REAL scale_factor;
+
+#define SCALE_FACTOR scale_factor
 
 static long double sb_even_eval(long double coord, int idx)
 {
@@ -229,6 +231,8 @@ static BrillDataContext *init_bd(cGH *cctkGH, float amplitude, int cheb_order)
 
     bd->nb_coeffs = bd->nb_coeffs_x * bd->nb_coeffs_z;
 
+    scale_factor = sf;
+
     return bd;
 }