Make the basis scaling factor configurable.

1 files changed, 5 insertions, 1 deletions

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;
 }