Add the simple time-antisym initial data family.

1 files changed, 7 insertions, 5 deletions

diff --git a/init.c b/init.c
index a7d54a8..ff995b6 100644
--- a/init.c
+++ b/init.c
@@ -249,14 +249,14 @@ static int teukolsky_init_check_options(TDContext *td)
     return 0;
 }
 
-static int constraint_eval_alloc(TDContext *td, double amplitude,
+static int constraint_eval_alloc(const TDContext *td, double amplitude,
                                  TDConstraintEvalContext **pce)
 {
     TDPriv *priv = td->priv;
     TDConstraintEvalContext *ce;
     int ret;
 
-    ret = tdi_constraint_eval_alloc(&ce, TD_FAMILY_AE_TIME_ANTISYM);
+    ret = tdi_constraint_eval_alloc(&ce, td->family);
     if (ret < 0) {
         tdi_log(&priv->logger, 0, "Error allocating the constraints evaluator\n");
         return ret;
@@ -405,8 +405,8 @@ int td_solve(TDContext *td, double *coeffs_init[3])
             tdi_constraint_eval_free(&ce);
             if (ret == -EDOM) {
                 inverse_step = 0.5 * inverse_step;
-                if (fabs(inverse_step) < 1e-2)
-                    return ret;
+                //if (fabs(inverse_step) < 1e-2)
+                //    return ret;
                 continue;
             } else if (ret < 0)
                 return ret;
@@ -446,7 +446,9 @@ TDContext *td_context_alloc(void)
 
     td->nb_threads = 1;
 
-    td->amplitude   = 1.0;
+    td->family          = TD_FAMILY_SIMPLE_TIME_ANTISYM;
+    td->solution_branch = 0;
+    td->amplitude       = 1.0;
 
     td->nb_coeffs[0] = 32;
     td->nb_coeffs[1] = 16;