/* * 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