From 788f68c3bbf7c74b717d8894be4610d39d0d97f2 Mon Sep 17 00:00:00 2001
From: Anton Khirnov <anton@khirnov.net>
Date: Mon, 23 Apr 2018 08:45:07 +0200
Subject: Add the simple time-antisym initial data family.

---
 teukolsky_data.h | 24 ++++++++++++++++++++++++
 1 file changed, 24 insertions(+)

(limited to 'teukolsky_data.h')

diff --git a/teukolsky_data.h b/teukolsky_data.h
index e0877cd..d95dbfe 100644
--- a/teukolsky_data.h
+++ b/teukolsky_data.h
@@ -36,6 +36,28 @@
  * Finally, free the solver context with td_context_free().
  */
 
+/**
+ * Identifies the specific initial data family to use.
+ */
+enum TDFamily {
+    /**
+     * The time-antisymmetric initial data used in Abrahams&Evans PhysRevD v49,n8 (1994).
+     * Conformally flat spatial metric.
+     *  r                / x        x  3 \        /    x  2  \
+     * K  = -60√(2/π) a | ---  -  (---)   | * exp| - (---)    | sin(2θ)
+     *  θ                \ L        L    /        \    L     /
+     */
+    TD_FAMILY_AE_TIME_ANTISYM = 0,
+    /**
+     * Simpler time-antisymmetric initial data.
+     * Conformally flat spatial metric.
+     *  r       x       /    x  2  \
+     * K  = a  ---  exp| - (---)    | sin(2θ)
+     *  θ       L       \    L     /
+     */
+    TD_FAMILY_SIMPLE_TIME_ANTISYM,
+};
+
 typedef struct TDContext {
     /**
      *  Solver internals, not to be accessed by the caller
@@ -105,6 +127,8 @@ typedef struct TDContext {
     double *coeffs[3];
 
     unsigned int solution_branch;
+
+    enum TDFamily family;
 } TDContext;
 
 /**
-- 
cgit v1.2.3