From 234e722c5abfb15c9850b4224cabbb3e2d0c3657 Mon Sep 17 00:00:00 2001 From: Anton Khirnov Date: Fri, 9 Feb 2018 15:34:25 +0100 Subject: Add the remaining parts for teukolsky waves initial data --- td_constraints.h | 43 +++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 43 insertions(+) create mode 100644 td_constraints.h (limited to 'td_constraints.h') diff --git a/td_constraints.h b/td_constraints.h new file mode 100644 index 0000000..44f5257 --- /dev/null +++ b/td_constraints.h @@ -0,0 +1,43 @@ +/* + * Copyright 2017 Anton Khirnov + * + * This program is free software: you can redistribute it and/or modify + * it under the terms of the GNU General Public License as published by + * the Free Software Foundation, either version 3 of the License, or + * (at your option) any later version. + * + * This program is distributed in the hope that it will be useful, + * but WITHOUT ANY WARRANTY; without even the implied warranty of + * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the + * GNU General Public License for more details. + * + * You should have received a copy of the GNU General Public License + * along with this program. If not, see . + */ + +#ifndef TEUKOLSKY_DATA_CONSTRAINTS_H +#define TEUKOLSKY_DATA_CONSTRAINTS_H + +#include "pssolve.h" + +enum TDConstraintEq { + TD_CONSTRAINT_EQ_HAM = 0, + TD_CONSTRAINT_EQ_MOM_0, + TD_CONSTRAINT_EQ_MOM_1, + TD_CONSTRAINT_EQ_NB, +}; + +enum TDConstraintVar { + TD_CONSTRAINT_VAR_CONFFAC = 0, + TD_CONSTRAINT_VAR_CURV_RR, + TD_CONSTRAINT_VAR_CURV_PHIPHI, + TD_CONSTRAINT_VAR_NB, +}; + +int tdi_constraints_eval(enum TDConstraintEq eq, + const double a, + unsigned int nb_coords[2], double *coords[2], + double *vars[TD_CONSTRAINT_VAR_NB][PSSOLVE_DIFF_ORDER_NB], + double *out); + +#endif // TEUKOLSKY_DATA_CONSTRAINTS_H -- cgit v1.2.3