aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorAnton Khirnov <anton@khirnov.net>2020-01-30 13:50:18 +0100
committerAnton Khirnov <anton@khirnov.net>2020-01-30 13:50:18 +0100
commit8e19356b524f865c7d104a19f98a98adc5154592 (patch)
tree1d149c38e589fbda2a7eaa2bfea17bc149f32fbf
parente2745670894cac8ca73c6bee5c1b4cc2c8c94a3e (diff)
teukolsky_data: rename TD_FAMILY_SIMPLE_TIME_ANTISYM ->TD_FAMILY_TIME_ANTISYM_LINEAR
For consistency with previous commit.
-rw-r--r--td_constraints.c16
-rw-r--r--teukolsky_data.h5
-rw-r--r--teukolsky_data.py2
3 files changed, 12 insertions, 11 deletions
diff --git a/td_constraints.c b/td_constraints.c
index 73d48a7..25adde1 100644
--- a/td_constraints.c
+++ b/td_constraints.c
@@ -191,31 +191,31 @@ static const TDFamilyDef time_antisym_cubic = {
.constraint_eval = constraint_funcs_confflat,
};
-static double k_rtheta_eval_simple_time_antisym(TDConstraintEvalContext *ctx,
+static double k_rtheta_eval_time_antisym_linear(TDConstraintEvalContext *ctx,
double r, double theta)
{
const double r2 = SQR(r);
return ctx->amplitude * r * exp(-r2) * sin(2.0 * theta);
}
-static double k_rtheta_eval_dr_simple_time_antisym(TDConstraintEvalContext *ctx,
+static double k_rtheta_eval_dr_time_antisym_linear(TDConstraintEvalContext *ctx,
double r, double theta)
{
const double r2 = SQR(r);
return ctx->amplitude * (1.0 - 2.0 * r2) * exp(-r2) * sin(2.0 * theta);
}
-static double k_rtheta_eval_dt_simple_time_antisym(TDConstraintEvalContext *ctx,
+static double k_rtheta_eval_dt_time_antisym_linear(TDConstraintEvalContext *ctx,
double r, double theta)
{
const double r2 = SQR(r);
return 2.0 * ctx->amplitude * r * exp(-r2) * cos(2.0 * theta);
}
-static const TDFamilyDef simple_time_antisym = {
- .eval_krt = k_rtheta_eval_simple_time_antisym,
- .eval_krt_dr = k_rtheta_eval_dr_simple_time_antisym,
- .eval_krt_dt = k_rtheta_eval_dt_simple_time_antisym,
+static const TDFamilyDef time_antisym_linear = {
+ .eval_krt = k_rtheta_eval_time_antisym_linear,
+ .eval_krt_dr = k_rtheta_eval_dr_time_antisym_linear,
+ .eval_krt_dt = k_rtheta_eval_dt_time_antisym_linear,
//.a_converge = 1.018918628476058,
.a_converge = 1.0189186,
.a_diverge = 1.018918628476968,
@@ -224,7 +224,7 @@ static const TDFamilyDef simple_time_antisym = {
static const TDFamilyDef *td_families[] = {
[TD_FAMILY_TIME_ANTISYM_CUBIC] = &time_antisym_cubic,
- [TD_FAMILY_SIMPLE_TIME_ANTISYM] = &simple_time_antisym,
+ [TD_FAMILY_TIME_ANTISYM_LINEAR] = &time_antisym_linear,
};
double tdi_constraint_eval_k_rtheta(TDConstraintEvalContext *ctx, double r, double theta)
diff --git a/teukolsky_data.h b/teukolsky_data.h
index 943407b..0394260 100644
--- a/teukolsky_data.h
+++ b/teukolsky_data.h
@@ -51,13 +51,14 @@ enum TDFamily {
*/
TD_FAMILY_TIME_ANTISYM_CUBIC = 0,
/**
- * Simpler time-antisymmetric initial data.
+ * Time-antisymmetric initial data with a linear radial term multiplied by
+ * an exponential and a quadrupole angular term.
* Conformally flat spatial metric.
* r x / x 2 \
* K = a --- exp| - (---) | sin(2θ)
* θ L \ L /
*/
- TD_FAMILY_SIMPLE_TIME_ANTISYM,
+ TD_FAMILY_TIME_ANTISYM_LINEAR,
};
typedef struct TDContext {
diff --git a/teukolsky_data.py b/teukolsky_data.py
index 2b6198e..a819e5c 100644
--- a/teukolsky_data.py
+++ b/teukolsky_data.py
@@ -22,7 +22,7 @@ import numpy as np
class TeukolskyData(object):
TD_FAMILY_TIME_ANTISYM_CUBIC = 0
- TD_FAMILY_SIMPLE_TIME_ANTISYM = 1
+ TD_FAMILY_TIME_ANTISYM_LINEAR = 1
coeffs = None