From 85ee8c05d6953bb361e678aae8fa965c27b8bd0f Mon Sep 17 00:00:00 2001 From: Anton Khirnov Date: Fri, 20 Apr 2018 09:39:20 +0200 Subject: td_constraints: drop unused debugging cruft --- td_constraints.c | 132 ------------------------------------------------------- 1 file changed, 132 deletions(-) diff --git a/td_constraints.c b/td_constraints.c index 9ef6a80..6795891 100644 --- a/td_constraints.c +++ b/td_constraints.c @@ -21,88 +21,9 @@ #include "pssolve.h" #include "td_constraints.h" -#define I_C sqrt(225.0 / (8.0 * M_PI)) -static double I_m2(double x, double a) -{ - return a * I_C * exp(-SQR(x)); -} -static double I_m1(double x, double a) -{ - return a * I_C * (-2.0 * x * exp(-SQR(x))); -} -static double I_0(double x, double a) -{ - return a * I_C * (-2.0 + 4.0 * SQR(x)) * exp(-SQR(x)); -} -static double I_p1(double x, double a) -{ - return a * I_C * (12.0 * x - 8.0 * x * SQR(x)) * exp(-SQR(x)); -} -static double I_p2(double x, double a) -{ - return a * I_C * (12.0 - 48.0 * SQR(x) + 16.0 * SQR(SQR(x))) * exp(-SQR(x)); -} -static double I_p3(double x, double a) -{ - return a * I_C * (-120.0 * x + 144.0 * x * SQR(x) - 32.0 * x * SQR(SQR(x))) * exp(-SQR(x)); -} - -static double K_rtheta(double a, double r, double theta) -{ - double r2 = r * r; - double r3 = r2 * r; - - // equal to the below - return a * I_C * exp(-r2) * sin(2 * theta) * (32.0 * r3 - 48.0 * r); - - //return ((I_p2(-r, a) + I_p2(r, a)) * r3 + 3.0 * (I_p1(-r, a) - I_p1(r, a)) * r2 + - // 6.0 * (I_0(-r, a) + I_0(r, a)) * r + 6.0 * (I_m1(-r, a) - I_m1(r, a))) * sin(2.0 * theta) / r4; -} - -static double K_rtheta_dtheta(double a, double r, double theta) -{ - return 2.0 * cos(2.0 * theta) * K_rtheta(a, r, theta) / sin(2.0 * theta); -} - -static double K_rtheta_dr(double a, double r, double theta) -{ - double r2 = r * r; - return -2 * r * K_rtheta(a, r, theta) + a * I_C * exp(-r2) * sin(2 * theta) * (3.0 * 32.0 * r2 - 48.0); - //return ((-I_p3(-r, a) + I_p3(r, a)) / r + 3.0 * (-I_p2(-r, a) - I_p2(r, a)) / (SQR(r)) + - // 6.0 * (-I_p1(-r, a) + I_p1(r, a)) / (r * SQR(r)) + 6.0 * (-I_0(-r, a) - I_0(r, a)) / (SQR(SQR(r))) + - // (I_p2(-r, a) + I_p2(r, a)) / (-SQR(r)) + 3.0 * (I_p1(-r, a) - I_p1(r, a)) * (-2.0 / (r * SQR(r))) + - // 6.0 * (I_0(-r, a) + I_0(r, a)) * (-3.0 / (SQR(SQR(r)))) + 6.0 * (I_m1(-r, a) - I_m1(r, a)) * (-4.0 / (r * SQR(SQR(r))))) * sin(2.0 * theta); -} - #include "eval_constraints.c" #include "eval_k_rtheta.c" -//typedef void (*eval_metric_func)(double a, double L, double r_0, double x_0, double x_1, double psi_val_00, double psi_val_01, double psi_val_02, double psi_val_10, double psi_val_11, double psi_val_20, double *res); -// -//static const eval_metric_func eval_dmetric[3] = { -// eval_dmetric_0, -// eval_dmetric_1, -// eval_dmetric_2, -//}; -// -//static const eval_metric_func eval_d2metric[3][3] = { -// { -// eval_d2metric_00, -// eval_d2metric_01, -// eval_d2metric_02, -// }, -// { -// eval_d2metric_10, -// eval_d2metric_11, -// eval_d2metric_12, -// }, -// { -// eval_d2metric_20, -// eval_d2metric_21, -// eval_d2metric_22, -// }, -//}; - #define EQNAME constraint_eq_ham #define EQUATION TD_CONSTRAINT_EQ_HAM #include "td_constraints_template.c" @@ -131,65 +52,12 @@ void [TD_CONSTRAINT_EQ_MOM_1] = constraint_eq_mom_1, }; -static double q_gundlach(const double a, double rho, double z) -{ - return a * SQR(rho) * exp(-(SQR(rho) + SQR(z))); -} - -static double dq_rho_gundlach(const double a, double rho, double z) -{ - return a * 2 * rho * exp(-SQR(rho) - SQR(z)) * (1.0 - rho * (rho)); -} - -static double d2q_rho_gundlach(const double a, double rho, double z) -{ - double rho2 = SQR(rho); - return a * 2 * exp(-SQR(rho) - SQR(z)) * (1.0 - 4.0 * rho * (rho) - rho2 + 2.0 * rho2 * SQR(rho)); -} - -static double d2q_rho_z_gundlach(const double a, double rho, double z) -{ - return a * 4 * z * rho * exp(-SQR(rho) - SQR(z)) * (rho * (rho) - 1.0); -} - -static double dq_z_gundlach(const double a, double rho, double z) -{ - return -a * 2 * z * SQR(rho) * exp(-SQR(rho) - SQR(z)); -} - -static double d2q_z_gundlach(const double a, double rho, double z) -{ - return a * 2 * SQR(rho) * exp(-SQR(rho) - SQR(z)) * (2 * SQR(z) - 1.0); -} - - 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) { - //for (int idx1 = 0; idx1 < nb_coords[1]; idx1++) { - // const double theta = coords[1][idx1]; - - // for (int idx0 = 0; idx0 < nb_coords[0]; idx0++) { - // const double r = coords[0][idx0]; - // const double rho = r * sin(theta); - // const double z = r * cos(theta); - - // const double ct = cos(theta); - // const double st = sin(theta); - - // const int idx = idx1 * nb_coords[0] + idx0; - - // const double laplace = vars[eq][PSSOLVE_DIFF_ORDER_20][idx] + - // vars[eq][PSSOLVE_DIFF_ORDER_10][idx] * (2.0 / r) + - // vars[eq][PSSOLVE_DIFF_ORDER_02][idx] * (1.0 / SQR(r)) + - // vars[eq][PSSOLVE_DIFF_ORDER_01][idx] * (ct / (SQR(r) * st)); - - // out[idx] = laplace + 0.25 * (d2q_z_gundlach(a, rho, z) + d2q_rho_gundlach(a, rho, z)) * (1.0 + vars[eq][PSSOLVE_DIFF_ORDER_00][idx]); - // } - //} constraint_funcs[eq](a, nb_coords, coords, vars, out); return 0; } -- cgit v1.2.3