diff options
Diffstat (limited to 'eval_constraints.c')
| -rw-r--r-- | eval_constraints.c | 37 |
1 files changed, 32 insertions, 5 deletions
diff --git a/eval_constraints.c b/eval_constraints.c index 71d78bf..42246dd 100644 --- a/eval_constraints.c +++ b/eval_constraints.c @@ -1,5 +1,32 @@ static inline -double eval_R_scalar(double a, double L, 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 K_rr_val_00, double K_rr_val_01, double K_rr_val_10, double K_phiphi_val_00, double K_phiphi_val_01, double K_phiphi_val_10) { +double eval_R_scalar(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 K_rr_val_00, double K_rr_val_01, double K_rr_val_10, double K_phiphi_val_00, double K_phiphi_val_01, double K_phiphi_val_10) { + + double eval_R_scalar_result; + eval_R_scalar_result = -1.0L/294912.0L*(514800*pow(M_PI, 131.0L/2.0L)*pow(L, 16)*pow(a, 2)*psi_val_00*pow(r_0, 2)*pow(sin(x_1), 3) + 257400*pow(M_PI, 131.0L/2.0L)*pow(L, 16)*pow(a, 2)*psi_val_00*r_0*x_0*pow(sin(x_1), 5) - 514800*pow(M_PI, 131.0L/2.0L)*pow(L, 16)*pow(a, 2)*psi_val_00*r_0*x_0*pow(sin(x_1), 3) - 96525*pow(M_PI, 131.0L/2.0L)*pow(L, 16)*pow(a, 2)*psi_val_00*pow(x_0, 2)*pow(sin(x_1), 5) + 128700*pow(M_PI, 131.0L/2.0L)*pow(L, 16)*pow(a, 2)*psi_val_00*pow(x_0, 2)*pow(sin(x_1), 3) + 184320*sqrt(143)*pow(M_PI, 66)*pow(L, 15)*a*psi_val_00*r_0*pow(x_0, 2)*pow(sin(x_1), 3) - 115200*sqrt(143)*pow(M_PI, 66)*pow(L, 15)*a*psi_val_00*r_0*pow(x_0, 2)*sin(x_1) - 86400*sqrt(143)*pow(M_PI, 66)*pow(L, 15)*a*psi_val_00*pow(x_0, 3)*pow(sin(x_1), 3) + 57600*sqrt(143)*pow(M_PI, 66)*pow(L, 15)*a*psi_val_00*pow(x_0, 3)*sin(x_1) - 184320*sqrt(143)*pow(M_PI, 66)*pow(L, 15)*a*psi_val_01*r_0*pow(x_0, 2)*pow(sin(x_1), 2)*cos(x_1) + 92160*sqrt(143)*pow(M_PI, 66)*pow(L, 15)*a*psi_val_01*pow(x_0, 3)*pow(sin(x_1), 2)*cos(x_1) + 184320*sqrt(143)*pow(M_PI, 66)*pow(L, 15)*a*psi_val_10*r_0*pow(x_0, 3)*pow(sin(x_1), 3) - 46080*sqrt(143)*pow(M_PI, 66)*pow(L, 15)*a*psi_val_10*pow(x_0, 4)*pow(sin(x_1), 3) - 6864000*pow(M_PI, 131.0L/2.0L)*pow(L, 14)*pow(a, 2)*psi_val_00*pow(r_0, 4)*pow(sin(x_1), 3) - 6864000*pow(M_PI, 131.0L/2.0L)*pow(L, 14)*pow(a, 2)*psi_val_00*pow(r_0, 3)*x_0*pow(sin(x_1), 5) + 13728000*pow(M_PI, 131.0L/2.0L)*pow(L, 14)*pow(a, 2)*psi_val_00*pow(r_0, 3)*x_0*pow(sin(x_1), 3) + 3861000*pow(M_PI, 131.0L/2.0L)*pow(L, 14)*pow(a, 2)*psi_val_00*pow(r_0, 2)*pow(x_0, 2)*pow(sin(x_1), 5) - 5148000*pow(M_PI, 131.0L/2.0L)*pow(L, 14)*pow(a, 2)*psi_val_00*pow(r_0, 2)*pow(x_0, 2)*pow(sin(x_1), 3) + 5148000*pow(M_PI, 131.0L/2.0L)*pow(L, 14)*pow(a, 2)*psi_val_00*r_0*pow(x_0, 3)*pow(sin(x_1), 5) - 3432000*pow(M_PI, 131.0L/2.0L)*pow(L, 14)*pow(a, 2)*psi_val_00*r_0*pow(x_0, 3)*pow(sin(x_1), 3) - 2145000*pow(M_PI, 131.0L/2.0L)*pow(L, 14)*pow(a, 2)*psi_val_00*pow(x_0, 4)*pow(sin(x_1), 5) + 1716000*pow(M_PI, 131.0L/2.0L)*pow(L, 14)*pow(a, 2)*psi_val_00*pow(x_0, 4)*pow(sin(x_1), 3) + 2359296*pow(M_PI, 133.0L/2.0L)*pow(L, 14)*psi_val_01*pow(x_0, 4)*cos(x_1) + 2359296*pow(M_PI, 133.0L/2.0L)*pow(L, 14)*psi_val_02*pow(x_0, 4)*sin(x_1) + 4718592*pow(M_PI, 133.0L/2.0L)*pow(L, 14)*psi_val_10*pow(x_0, 5)*sin(x_1) + 2359296*pow(M_PI, 133.0L/2.0L)*pow(L, 14)*psi_val_20*pow(x_0, 6)*sin(x_1) - 1228800*sqrt(143)*pow(M_PI, 66)*pow(L, 13)*a*psi_val_00*pow(r_0, 3)*pow(x_0, 2)*pow(sin(x_1), 3) + 768000*sqrt(143)*pow(M_PI, 66)*pow(L, 13)*a*psi_val_00*pow(r_0, 3)*pow(x_0, 2)*sin(x_1) + 1728000*sqrt(143)*pow(M_PI, 66)*pow(L, 13)*a*psi_val_00*pow(r_0, 2)*pow(x_0, 3)*pow(sin(x_1), 3) - 1152000*sqrt(143)*pow(M_PI, 66)*pow(L, 13)*a*psi_val_00*pow(r_0, 2)*pow(x_0, 3)*sin(x_1) - 268800*sqrt(143)*pow(M_PI, 66)*pow(L, 13)*a*psi_val_00*pow(x_0, 5)*pow(sin(x_1), 3) + 384000*sqrt(143)*pow(M_PI, 66)*pow(L, 13)*a*psi_val_00*pow(x_0, 5)*sin(x_1) + 1228800*sqrt(143)*pow(M_PI, 66)*pow(L, 13)*a*psi_val_01*pow(r_0, 3)*pow(x_0, 2)*pow(sin(x_1), 2)*cos(x_1) - 1843200*sqrt(143)*pow(M_PI, 66)*pow(L, 13)*a*psi_val_01*pow(r_0, 2)*pow(x_0, 3)*pow(sin(x_1), 2)*cos(x_1) + 614400*sqrt(143)*pow(M_PI, 66)*pow(L, 13)*a*psi_val_01*pow(x_0, 5)*pow(sin(x_1), 2)*cos(x_1) - 1228800*sqrt(143)*pow(M_PI, 66)*pow(L, 13)*a*psi_val_10*pow(r_0, 3)*pow(x_0, 3)*pow(sin(x_1), 3) + 921600*sqrt(143)*pow(M_PI, 66)*pow(L, 13)*a*psi_val_10*pow(r_0, 2)*pow(x_0, 4)*pow(sin(x_1), 3) + 307200*sqrt(143)*pow(M_PI, 66)*pow(L, 13)*a*psi_val_10*pow(x_0, 6)*pow(sin(x_1), 3) + 37294400*pow(M_PI, 131.0L/2.0L)*pow(L, 12)*pow(a, 2)*psi_val_00*pow(r_0, 6)*pow(sin(x_1), 3) + 55941600*pow(M_PI, 131.0L/2.0L)*pow(L, 12)*pow(a, 2)*psi_val_00*pow(r_0, 5)*x_0*pow(sin(x_1), 5) - 111883200*pow(M_PI, 131.0L/2.0L)*pow(L, 12)*pow(a, 2)*psi_val_00*pow(r_0, 5)*x_0*pow(sin(x_1), 3) - 52123500*pow(M_PI, 131.0L/2.0L)*pow(L, 12)*pow(a, 2)*psi_val_00*pow(r_0, 4)*pow(x_0, 2)*pow(sin(x_1), 5) + 69498000*pow(M_PI, 131.0L/2.0L)*pow(L, 12)*pow(a, 2)*psi_val_00*pow(r_0, 4)*pow(x_0, 2)*pow(sin(x_1), 3) - 142428000*pow(M_PI, 131.0L/2.0L)*pow(L, 12)*pow(a, 2)*psi_val_00*pow(r_0, 3)*pow(x_0, 3)*pow(sin(x_1), 5) + 94952000*pow(M_PI, 131.0L/2.0L)*pow(L, 12)*pow(a, 2)*psi_val_00*pow(r_0, 3)*pow(x_0, 3)*pow(sin(x_1), 3) + 232089000*pow(M_PI, 131.0L/2.0L)*pow(L, 12)*pow(a, 2)*psi_val_00*pow(r_0, 2)*pow(x_0, 4)*pow(sin(x_1), 5) - 142428000*pow(M_PI, 131.0L/2.0L)*pow(L, 12)*pow(a, 2)*psi_val_00*pow(r_0, 2)*pow(x_0, 4)*pow(sin(x_1), 3) - 108108000*pow(M_PI, 131.0L/2.0L)*pow(L, 12)*pow(a, 2)*psi_val_00*r_0*pow(x_0, 5)*pow(sin(x_1), 5) + 57657600*pow(M_PI, 131.0L/2.0L)*pow(L, 12)*pow(a, 2)*psi_val_00*r_0*pow(x_0, 5)*pow(sin(x_1), 3) + 14628900*pow(M_PI, 131.0L/2.0L)*pow(L, 12)*pow(a, 2)*psi_val_00*pow(x_0, 6)*pow(sin(x_1), 5) - 5090800*pow(M_PI, 131.0L/2.0L)*pow(L, 12)*pow(a, 2)*psi_val_00*pow(x_0, 6)*pow(sin(x_1), 3) + 2580480*sqrt(143)*pow(M_PI, 66)*pow(L, 11)*a*psi_val_00*pow(r_0, 5)*pow(x_0, 2)*pow(sin(x_1), 3) - 1612800*sqrt(143)*pow(M_PI, 66)*pow(L, 11)*a*psi_val_00*pow(r_0, 5)*pow(x_0, 2)*sin(x_1) - 6048000*sqrt(143)*pow(M_PI, 66)*pow(L, 11)*a*psi_val_00*pow(r_0, 4)*pow(x_0, 3)*pow(sin(x_1), 3) + 4032000*sqrt(143)*pow(M_PI, 66)*pow(L, 11)*a*psi_val_00*pow(r_0, 4)*pow(x_0, 3)*sin(x_1) + 5644800*sqrt(143)*pow(M_PI, 66)*pow(L, 11)*a*psi_val_00*pow(r_0, 2)*pow(x_0, 5)*pow(sin(x_1), 3) - 8064000*sqrt(143)*pow(M_PI, 66)*pow(L, 11)*a*psi_val_00*pow(r_0, 2)*pow(x_0, 5)*sin(x_1) + 8064000*sqrt(143)*pow(M_PI, 66)*pow(L, 11)*a*psi_val_00*r_0*pow(x_0, 6)*sin(x_1) - 2177280*sqrt(143)*pow(M_PI, 66)*pow(L, 11)*a*psi_val_00*pow(x_0, 7)*pow(sin(x_1), 3) - 2419200*sqrt(143)*pow(M_PI, 66)*pow(L, 11)*a*psi_val_00*pow(x_0, 7)*sin(x_1) - 2580480*sqrt(143)*pow(M_PI, 66)*pow(L, 11)*a*psi_val_01*pow(r_0, 5)*pow(x_0, 2)*pow(sin(x_1), 2)*cos(x_1) + 6451200*sqrt(143)*pow(M_PI, 66)*pow(L, 11)*a*psi_val_01*pow(r_0, 4)*pow(x_0, 3)*pow(sin(x_1), 2)*cos(x_1) - 12902400*sqrt(143)*pow(M_PI, 66)*pow(L, 11)*a*psi_val_01*pow(r_0, 2)*pow(x_0, 5)*pow(sin(x_1), 2)*cos(x_1) + 12902400*sqrt(143)*pow(M_PI, 66)*pow(L, 11)*a*psi_val_01*r_0*pow(x_0, 6)*pow(sin(x_1), 2)*cos(x_1) - 3870720*sqrt(143)*pow(M_PI, 66)*pow(L, 11)*a*psi_val_01*pow(x_0, 7)*pow(sin(x_1), 2)*cos(x_1) + 2580480*sqrt(143)*pow(M_PI, 66)*pow(L, 11)*a*psi_val_10*pow(r_0, 5)*pow(x_0, 3)*pow(sin(x_1), 3) - 3225600*sqrt(143)*pow(M_PI, 66)*pow(L, 11)*a*psi_val_10*pow(r_0, 4)*pow(x_0, 4)*pow(sin(x_1), 3) - 6451200*sqrt(143)*pow(M_PI, 66)*pow(L, 11)*a*psi_val_10*pow(r_0, 2)*pow(x_0, 6)*pow(sin(x_1), 3) + 12902400*sqrt(143)*pow(M_PI, 66)*pow(L, 11)*a*psi_val_10*r_0*pow(x_0, 7)*pow(sin(x_1), 3) - 5806080*sqrt(143)*pow(M_PI, 66)*pow(L, 11)*a*psi_val_10*pow(x_0, 8)*pow(sin(x_1), 3) - 108451200*pow(M_PI, 131.0L/2.0L)*pow(L, 10)*pow(a, 2)*psi_val_00*pow(r_0, 8)*pow(sin(x_1), 3) - 216902400*pow(M_PI, 131.0L/2.0L)*pow(L, 10)*pow(a, 2)*psi_val_00*pow(r_0, 7)*x_0*pow(sin(x_1), 5) + 433804800*pow(M_PI, 131.0L/2.0L)*pow(L, 10)*pow(a, 2)*psi_val_00*pow(r_0, 7)*x_0*pow(sin(x_1), 3) + 286486200*pow(M_PI, 131.0L/2.0L)*pow(L, 10)*pow(a, 2)*psi_val_00*pow(r_0, 6)*pow(x_0, 2)*pow(sin(x_1), 5) - 381981600*pow(M_PI, 131.0L/2.0L)*pow(L, 10)*pow(a, 2)*psi_val_00*pow(r_0, 6)*pow(x_0, 2)*pow(sin(x_1), 3) + 1117116000*pow(M_PI, 131.0L/2.0L)*pow(L, 10)*pow(a, 2)*psi_val_00*pow(r_0, 5)*pow(x_0, 3)*pow(sin(x_1), 5) - 744744000*pow(M_PI, 131.0L/2.0L)*pow(L, 10)*pow(a, 2)*psi_val_00*pow(r_0, 5)*pow(x_0, 3)*pow(sin(x_1), 3) - 3012009000*pow(M_PI, 131.0L/2.0L)*pow(L, 10)*pow(a, 2)*psi_val_00*pow(r_0, 4)*pow(x_0, 4)*pow(sin(x_1), 5) + 1861860000*pow(M_PI, 131.0L/2.0L)*pow(L, 10)*pow(a, 2)*psi_val_00*pow(r_0, 4)*pow(x_0, 4)*pow(sin(x_1), 3) + 2738736000*pow(M_PI, 131.0L/2.0L)*pow(L, 10)*pow(a, 2)*psi_val_00*pow(r_0, 3)*pow(x_0, 5)*pow(sin(x_1), 5) - 1470268800*pow(M_PI, 131.0L/2.0L)*pow(L, 10)*pow(a, 2)*psi_val_00*pow(r_0, 3)*pow(x_0, 5)*pow(sin(x_1), 3) - 1057656600*pow(M_PI, 131.0L/2.0L)*pow(L, 10)*pow(a, 2)*psi_val_00*pow(r_0, 2)*pow(x_0, 6)*pow(sin(x_1), 5) + 343543200*pow(M_PI, 131.0L/2.0L)*pow(L, 10)*pow(a, 2)*psi_val_00*pow(r_0, 2)*pow(x_0, 6)*pow(sin(x_1), 3) + 160960800*pow(M_PI, 131.0L/2.0L)*pow(L, 10)*pow(a, 2)*psi_val_00*r_0*pow(x_0, 7)*pow(sin(x_1), 5) + 122865600*pow(M_PI, 131.0L/2.0L)*pow(L, 10)*pow(a, 2)*psi_val_00*r_0*pow(x_0, 7)*pow(sin(x_1), 3) - 16731000*pow(M_PI, 131.0L/2.0L)*pow(L, 10)*pow(a, 2)*psi_val_00*pow(x_0, 8)*pow(sin(x_1), 5) - 56628000*pow(M_PI, 131.0L/2.0L)*pow(L, 10)*pow(a, 2)*psi_val_00*pow(x_0, 8)*pow(sin(x_1), 3) - 2211840*sqrt(143)*pow(M_PI, 66)*pow(L, 9)*a*psi_val_00*pow(r_0, 7)*pow(x_0, 2)*pow(sin(x_1), 3) + 1382400*sqrt(143)*pow(M_PI, 66)*pow(L, 9)*a*psi_val_00*pow(r_0, 7)*pow(x_0, 2)*sin(x_1) + 7257600*sqrt(143)*pow(M_PI, 66)*pow(L, 9)*a*psi_val_00*pow(r_0, 6)*pow(x_0, 3)*pow(sin(x_1), 3) - 4838400*sqrt(143)*pow(M_PI, 66)*pow(L, 9)*a*psi_val_00*pow(r_0, 6)*pow(x_0, 3)*sin(x_1) - 16934400*sqrt(143)*pow(M_PI, 66)*pow(L, 9)*a*psi_val_00*pow(r_0, 4)*pow(x_0, 5)*pow(sin(x_1), 3) + 24192000*sqrt(143)*pow(M_PI, 66)*pow(L, 9)*a*psi_val_00*pow(r_0, 4)*pow(x_0, 5)*sin(x_1) - 48384000*sqrt(143)*pow(M_PI, 66)*pow(L, 9)*a*psi_val_00*pow(r_0, 3)*pow(x_0, 6)*sin(x_1) + 39191040*sqrt(143)*pow(M_PI, 66)*pow(L, 9)*a*psi_val_00*pow(r_0, 2)*pow(x_0, 7)*pow(sin(x_1), 3) + 43545600*sqrt(143)*pow(M_PI, 66)*pow(L, 9)*a*psi_val_00*pow(r_0, 2)*pow(x_0, 7)*sin(x_1) - 38707200*sqrt(143)*pow(M_PI, 66)*pow(L, 9)*a*psi_val_00*r_0*pow(x_0, 8)*pow(sin(x_1), 3) - 19353600*sqrt(143)*pow(M_PI, 66)*pow(L, 9)*a*psi_val_00*r_0*pow(x_0, 8)*sin(x_1) + 11404800*sqrt(143)*pow(M_PI, 66)*pow(L, 9)*a*psi_val_00*pow(x_0, 9)*pow(sin(x_1), 3) + 3456000*sqrt(143)*pow(M_PI, 66)*pow(L, 9)*a*psi_val_00*pow(x_0, 9)*sin(x_1) + 2211840*sqrt(143)*pow(M_PI, 66)*pow(L, 9)*a*psi_val_01*pow(r_0, 7)*pow(x_0, 2)*pow(sin(x_1), 2)*cos(x_1) - 7741440*sqrt(143)*pow(M_PI, 66)*pow(L, 9)*a*psi_val_01*pow(r_0, 6)*pow(x_0, 3)*pow(sin(x_1), 2)*cos(x_1) + 38707200*sqrt(143)*pow(M_PI, 66)*pow(L, 9)*a*psi_val_01*pow(r_0, 4)*pow(x_0, 5)*pow(sin(x_1), 2)*cos(x_1) - 77414400*sqrt(143)*pow(M_PI, 66)*pow(L, 9)*a*psi_val_01*pow(r_0, 3)*pow(x_0, 6)*pow(sin(x_1), 2)*cos(x_1) + 69672960*sqrt(143)*pow(M_PI, 66)*pow(L, 9)*a*psi_val_01*pow(r_0, 2)*pow(x_0, 7)*pow(sin(x_1), 2)*cos(x_1) - 30965760*sqrt(143)*pow(M_PI, 66)*pow(L, 9)*a*psi_val_01*r_0*pow(x_0, 8)*pow(sin(x_1), 2)*cos(x_1) + 5529600*sqrt(143)*pow(M_PI, 66)*pow(L, 9)*a*psi_val_01*pow(x_0, 9)*pow(sin(x_1), 2)*cos(x_1) - 2211840*sqrt(143)*pow(M_PI, 66)*pow(L, 9)*a*psi_val_10*pow(r_0, 7)*pow(x_0, 3)*pow(sin(x_1), 3) + 3870720*sqrt(143)*pow(M_PI, 66)*pow(L, 9)*a*psi_val_10*pow(r_0, 6)*pow(x_0, 4)*pow(sin(x_1), 3) + 19353600*sqrt(143)*pow(M_PI, 66)*pow(L, 9)*a*psi_val_10*pow(r_0, 4)*pow(x_0, 6)*pow(sin(x_1), 3) - 77414400*sqrt(143)*pow(M_PI, 66)*pow(L, 9)*a*psi_val_10*pow(r_0, 3)*pow(x_0, 7)*pow(sin(x_1), 3) + 104509440*sqrt(143)*pow(M_PI, 66)*pow(L, 9)*a*psi_val_10*pow(r_0, 2)*pow(x_0, 8)*pow(sin(x_1), 3) - 61931520*sqrt(143)*pow(M_PI, 66)*pow(L, 9)*a*psi_val_10*r_0*pow(x_0, 9)*pow(sin(x_1), 3) + 13824000*sqrt(143)*pow(M_PI, 66)*pow(L, 9)*a*psi_val_10*pow(x_0, 10)*pow(sin(x_1), 3) + 187044000*pow(M_PI, 131.0L/2.0L)*pow(L, 8)*pow(a, 2)*psi_val_00*pow(r_0, 10)*pow(sin(x_1), 3) + 467610000*pow(M_PI, 131.0L/2.0L)*pow(L, 8)*pow(a, 2)*psi_val_00*pow(r_0, 9)*x_0*pow(sin(x_1), 5) - 935220000*pow(M_PI, 131.0L/2.0L)*pow(L, 8)*pow(a, 2)*psi_val_00*pow(r_0, 9)*x_0*pow(sin(x_1), 3) - 803667150*pow(M_PI, 131.0L/2.0L)*pow(L, 8)*pow(a, 2)*psi_val_00*pow(r_0, 8)*pow(x_0, 2)*pow(sin(x_1), 5) + 1071556200*pow(M_PI, 131.0L/2.0L)*pow(L, 8)*pow(a, 2)*psi_val_00*pow(r_0, 8)*pow(x_0, 2)*pow(sin(x_1), 3) - 3975285600*pow(M_PI, 131.0L/2.0L)*pow(L, 8)*pow(a, 2)*psi_val_00*pow(r_0, 7)*pow(x_0, 3)*pow(sin(x_1), 5) + 2650190400*pow(M_PI, 131.0L/2.0L)*pow(L, 8)*pow(a, 2)*psi_val_00*pow(r_0, 7)*pow(x_0, 3)*pow(sin(x_1), 3) + 14870255400*pow(M_PI, 131.0L/2.0L)*pow(L, 8)*pow(a, 2)*psi_val_00*pow(r_0, 6)*pow(x_0, 4)*pow(sin(x_1), 5) - 9275666400*pow(M_PI, 131.0L/2.0L)*pow(L, 8)*pow(a, 2)*psi_val_00*pow(r_0, 6)*pow(x_0, 4)*pow(sin(x_1), 3) - 19654034400*pow(M_PI, 131.0L/2.0L)*pow(L, 8)*pow(a, 2)*psi_val_00*pow(r_0, 5)*pow(x_0, 5)*pow(sin(x_1), 5) + 10695484800*pow(M_PI, 131.0L/2.0L)*pow(L, 8)*pow(a, 2)*psi_val_00*pow(r_0, 5)*pow(x_0, 5)*pow(sin(x_1), 3) + 11689177500*pow(M_PI, 131.0L/2.0L)*pow(L, 8)*pow(a, 2)*psi_val_00*pow(r_0, 4)*pow(x_0, 6)*pow(sin(x_1), 5) - 3549546000*pow(M_PI, 131.0L/2.0L)*pow(L, 8)*pow(a, 2)*psi_val_00*pow(r_0, 4)*pow(x_0, 6)*pow(sin(x_1), 3) - 2760357600*pow(M_PI, 131.0L/2.0L)*pow(L, 8)*pow(a, 2)*psi_val_00*pow(r_0, 3)*pow(x_0, 7)*pow(sin(x_1), 5) - 3893947200*pow(M_PI, 131.0L/2.0L)*pow(L, 8)*pow(a, 2)*psi_val_00*pow(r_0, 3)*pow(x_0, 7)*pow(sin(x_1), 3) + 782238600*pow(M_PI, 131.0L/2.0L)*pow(L, 8)*pow(a, 2)*psi_val_00*pow(r_0, 2)*pow(x_0, 8)*pow(sin(x_1), 5) + 4752633600*pow(M_PI, 131.0L/2.0L)*pow(L, 8)*pow(a, 2)*psi_val_00*pow(r_0, 2)*pow(x_0, 8)*pow(sin(x_1), 3) - 943628400*pow(M_PI, 131.0L/2.0L)*pow(L, 8)*pow(a, 2)*psi_val_00*r_0*pow(x_0, 9)*pow(sin(x_1), 5) - 2025909600*pow(M_PI, 131.0L/2.0L)*pow(L, 8)*pow(a, 2)*psi_val_00*r_0*pow(x_0, 9)*pow(sin(x_1), 3) + 327691650*pow(M_PI, 131.0L/2.0L)*pow(L, 8)*pow(a, 2)*psi_val_00*pow(x_0, 10)*pow(sin(x_1), 5) + 323380200*pow(M_PI, 131.0L/2.0L)*pow(L, 8)*pow(a, 2)*psi_val_00*pow(x_0, 10)*pow(sin(x_1), 3) + 675840*sqrt(143)*pow(M_PI, 66)*pow(L, 7)*a*psi_val_00*pow(r_0, 9)*pow(x_0, 2)*pow(sin(x_1), 3) - 422400*sqrt(143)*pow(M_PI, 66)*pow(L, 7)*a*psi_val_00*pow(r_0, 9)*pow(x_0, 2)*sin(x_1) - 2851200*sqrt(143)*pow(M_PI, 66)*pow(L, 7)*a*psi_val_00*pow(r_0, 8)*pow(x_0, 3)*pow(sin(x_1), 3) + 1900800*sqrt(143)*pow(M_PI, 66)*pow(L, 7)*a*psi_val_00*pow(r_0, 8)*pow(x_0, 3)*sin(x_1) + 12418560*sqrt(143)*pow(M_PI, 66)*pow(L, 7)*a*psi_val_00*pow(r_0, 6)*pow(x_0, 5)*pow(sin(x_1), 3) - 17740800*sqrt(143)*pow(M_PI, 66)*pow(L, 7)*a*psi_val_00*pow(r_0, 6)*pow(x_0, 5)*sin(x_1) + 53222400*sqrt(143)*pow(M_PI, 66)*pow(L, 7)*a*psi_val_00*pow(r_0, 5)*pow(x_0, 6)*sin(x_1) - 71850240*sqrt(143)*pow(M_PI, 66)*pow(L, 7)*a*psi_val_00*pow(r_0, 4)*pow(x_0, 7)*pow(sin(x_1), 3) - 79833600*sqrt(143)*pow(M_PI, 66)*pow(L, 7)*a*psi_val_00*pow(r_0, 4)*pow(x_0, 7)*sin(x_1) + 141926400*sqrt(143)*pow(M_PI, 66)*pow(L, 7)*a*psi_val_00*pow(r_0, 3)*pow(x_0, 8)*pow(sin(x_1), 3) + 70963200*sqrt(143)*pow(M_PI, 66)*pow(L, 7)*a*psi_val_00*pow(r_0, 3)*pow(x_0, 8)*sin(x_1) - 125452800*sqrt(143)*pow(M_PI, 66)*pow(L, 7)*a*psi_val_00*pow(r_0, 2)*pow(x_0, 9)*pow(sin(x_1), 3) - 38016000*sqrt(143)*pow(M_PI, 66)*pow(L, 7)*a*psi_val_00*pow(r_0, 2)*pow(x_0, 9)*sin(x_1) + 54743040*sqrt(143)*pow(M_PI, 66)*pow(L, 7)*a*psi_val_00*r_0*pow(x_0, 10)*pow(sin(x_1), 3) + 11404800*sqrt(143)*pow(M_PI, 66)*pow(L, 7)*a*psi_val_00*r_0*pow(x_0, 10)*sin(x_1) - 9609600*sqrt(143)*pow(M_PI, 66)*pow(L, 7)*a*psi_val_00*pow(x_0, 11)*pow(sin(x_1), 3) - 1478400*sqrt(143)*pow(M_PI, 66)*pow(L, 7)*a*psi_val_00*pow(x_0, 11)*sin(x_1) - 675840*sqrt(143)*pow(M_PI, 66)*pow(L, 7)*a*psi_val_01*pow(r_0, 9)*pow(x_0, 2)*pow(sin(x_1), 2)*cos(x_1) + 3041280*sqrt(143)*pow(M_PI, 66)*pow(L, 7)*a*psi_val_01*pow(r_0, 8)*pow(x_0, 3)*pow(sin(x_1), 2)*cos(x_1) - 28385280*sqrt(143)*pow(M_PI, 66)*pow(L, 7)*a*psi_val_01*pow(r_0, 6)*pow(x_0, 5)*pow(sin(x_1), 2)*cos(x_1) + 85155840*sqrt(143)*pow(M_PI, 66)*pow(L, 7)*a*psi_val_01*pow(r_0, 5)*pow(x_0, 6)*pow(sin(x_1), 2)*cos(x_1) - 127733760*sqrt(143)*pow(M_PI, 66)*pow(L, 7)*a*psi_val_01*pow(r_0, 4)*pow(x_0, 7)*pow(sin(x_1), 2)*cos(x_1) + 113541120*sqrt(143)*pow(M_PI, 66)*pow(L, 7)*a*psi_val_01*pow(r_0, 3)*pow(x_0, 8)*pow(sin(x_1), 2)*cos(x_1) - 60825600*sqrt(143)*pow(M_PI, 66)*pow(L, 7)*a*psi_val_01*pow(r_0, 2)*pow(x_0, 9)*pow(sin(x_1), 2)*cos(x_1) + 18247680*sqrt(143)*pow(M_PI, 66)*pow(L, 7)*a*psi_val_01*r_0*pow(x_0, 10)*pow(sin(x_1), 2)*cos(x_1) - 2365440*sqrt(143)*pow(M_PI, 66)*pow(L, 7)*a*psi_val_01*pow(x_0, 11)*pow(sin(x_1), 2)*cos(x_1) + 675840*sqrt(143)*pow(M_PI, 66)*pow(L, 7)*a*psi_val_10*pow(r_0, 9)*pow(x_0, 3)*pow(sin(x_1), 3) - 1520640*sqrt(143)*pow(M_PI, 66)*pow(L, 7)*a*psi_val_10*pow(r_0, 8)*pow(x_0, 4)*pow(sin(x_1), 3) - 14192640*sqrt(143)*pow(M_PI, 66)*pow(L, 7)*a*psi_val_10*pow(r_0, 6)*pow(x_0, 6)*pow(sin(x_1), 3) + 85155840*sqrt(143)*pow(M_PI, 66)*pow(L, 7)*a*psi_val_10*pow(r_0, 5)*pow(x_0, 7)*pow(sin(x_1), 3) - 191600640*sqrt(143)*pow(M_PI, 66)*pow(L, 7)*a*psi_val_10*pow(r_0, 4)*pow(x_0, 8)*pow(sin(x_1), 3) + 227082240*sqrt(143)*pow(M_PI, 66)*pow(L, 7)*a*psi_val_10*pow(r_0, 3)*pow(x_0, 9)*pow(sin(x_1), 3) - 152064000*sqrt(143)*pow(M_PI, 66)*pow(L, 7)*a*psi_val_10*pow(r_0, 2)*pow(x_0, 10)*pow(sin(x_1), 3) + 54743040*sqrt(143)*pow(M_PI, 66)*pow(L, 7)*a*psi_val_10*r_0*pow(x_0, 11)*pow(sin(x_1), 3) - 8279040*sqrt(143)*pow(M_PI, 66)*pow(L, 7)*a*psi_val_10*pow(x_0, 12)*pow(sin(x_1), 3) - 198140800*pow(M_PI, 131.0L/2.0L)*pow(L, 6)*pow(a, 2)*psi_val_00*pow(r_0, 12)*pow(sin(x_1), 3) - 594422400*pow(M_PI, 131.0L/2.0L)*pow(L, 6)*pow(a, 2)*psi_val_00*pow(r_0, 11)*x_0*pow(sin(x_1), 5) + 1188844800*pow(M_PI, 131.0L/2.0L)*pow(L, 6)*pow(a, 2)*psi_val_00*pow(r_0, 11)*x_0*pow(sin(x_1), 3) + 1262547000*pow(M_PI, 131.0L/2.0L)*pow(L, 6)*pow(a, 2)*psi_val_00*pow(r_0, 10)*pow(x_0, 2)*pow(sin(x_1), 5) - 1683396000*pow(M_PI, 131.0L/2.0L)*pow(L, 6)*pow(a, 2)*psi_val_00*pow(r_0, 10)*pow(x_0, 2)*pow(sin(x_1), 3) + 7442292000*pow(M_PI, 131.0L/2.0L)*pow(L, 6)*pow(a, 2)*psi_val_00*pow(r_0, 9)*pow(x_0, 3)*pow(sin(x_1), 5) - 4961528000*pow(M_PI, 131.0L/2.0L)*pow(L, 6)*pow(a, 2)*psi_val_00*pow(r_0, 9)*pow(x_0, 3)*pow(sin(x_1), 3) - 35476155000*pow(M_PI, 131.0L/2.0L)*pow(L, 6)*pow(a, 2)*psi_val_00*pow(r_0, 8)*pow(x_0, 4)*pow(sin(x_1), 5) + 22326876000*pow(M_PI, 131.0L/2.0L)*pow(L, 6)*pow(a, 2)*psi_val_00*pow(r_0, 8)*pow(x_0, 4)*pow(sin(x_1), 3) + 60540480000*pow(M_PI, 131.0L/2.0L)*pow(L, 6)*pow(a, 2)*psi_val_00*pow(r_0, 7)*pow(x_0, 5)*pow(sin(x_1), 5) - 33383750400*pow(M_PI, 131.0L/2.0L)*pow(L, 6)*pow(a, 2)*psi_val_00*pow(r_0, 7)*pow(x_0, 5)*pow(sin(x_1), 3) - 46115269200*pow(M_PI, 131.0L/2.0L)*pow(L, 6)*pow(a, 2)*psi_val_00*pow(r_0, 6)*pow(x_0, 6)*pow(sin(x_1), 5) + 12651038400*pow(M_PI, 131.0L/2.0L)*pow(L, 6)*pow(a, 2)*psi_val_00*pow(r_0, 6)*pow(x_0, 6)*pow(sin(x_1), 3) + 10738728000*pow(M_PI, 131.0L/2.0L)*pow(L, 6)*pow(a, 2)*psi_val_00*pow(r_0, 5)*pow(x_0, 7)*pow(sin(x_1), 5) + 31897008000*pow(M_PI, 131.0L/2.0L)*pow(L, 6)*pow(a, 2)*psi_val_00*pow(r_0, 5)*pow(x_0, 7)*pow(sin(x_1), 3) - 4254822000*pow(M_PI, 131.0L/2.0L)*pow(L, 6)*pow(a, 2)*psi_val_00*pow(r_0, 4)*pow(x_0, 8)*pow(sin(x_1), 5) - 59274072000*pow(M_PI, 131.0L/2.0L)*pow(L, 6)*pow(a, 2)*psi_val_00*pow(r_0, 4)*pow(x_0, 8)*pow(sin(x_1), 3) + 21106800000*pow(M_PI, 131.0L/2.0L)*pow(L, 6)*pow(a, 2)*psi_val_00*pow(r_0, 3)*pow(x_0, 9)*pow(sin(x_1), 5) + 49439104000*pow(M_PI, 131.0L/2.0L)*pow(L, 6)*pow(a, 2)*psi_val_00*pow(r_0, 3)*pow(x_0, 9)*pow(sin(x_1), 3) - 23511345000*pow(M_PI, 131.0L/2.0L)*pow(L, 6)*pow(a, 2)*psi_val_00*pow(r_0, 2)*pow(x_0, 10)*pow(sin(x_1), 5) - 23443648800*pow(M_PI, 131.0L/2.0L)*pow(L, 6)*pow(a, 2)*psi_val_00*pow(r_0, 2)*pow(x_0, 10)*pow(sin(x_1), 3) + 10677981600*pow(M_PI, 131.0L/2.0L)*pow(L, 6)*pow(a, 2)*psi_val_00*r_0*pow(x_0, 11)*pow(sin(x_1), 5) + 6134356800*pow(M_PI, 131.0L/2.0L)*pow(L, 6)*pow(a, 2)*psi_val_00*r_0*pow(x_0, 11)*pow(sin(x_1), 3) - 1816815000*pow(M_PI, 131.0L/2.0L)*pow(L, 6)*pow(a, 2)*psi_val_00*pow(x_0, 12)*pow(sin(x_1), 5) - 692692000*pow(M_PI, 131.0L/2.0L)*pow(L, 6)*pow(a, 2)*psi_val_00*pow(x_0, 12)*pow(sin(x_1), 3) + 126984000*pow(M_PI, 131.0L/2.0L)*pow(L, 4)*pow(a, 2)*psi_val_00*pow(r_0, 14)*pow(sin(x_1), 3) + 444444000*pow(M_PI, 131.0L/2.0L)*pow(L, 4)*pow(a, 2)*psi_val_00*pow(r_0, 13)*x_0*pow(sin(x_1), 5) - 888888000*pow(M_PI, 131.0L/2.0L)*pow(L, 4)*pow(a, 2)*psi_val_00*pow(r_0, 13)*x_0*pow(sin(x_1), 3) - 1127025900*pow(M_PI, 131.0L/2.0L)*pow(L, 4)*pow(a, 2)*psi_val_00*pow(r_0, 12)*pow(x_0, 2)*pow(sin(x_1), 5) + 1502701200*pow(M_PI, 131.0L/2.0L)*pow(L, 4)*pow(a, 2)*psi_val_00*pow(r_0, 12)*pow(x_0, 2)*pow(sin(x_1), 3) - 7618010400*pow(M_PI, 131.0L/2.0L)*pow(L, 4)*pow(a, 2)*psi_val_00*pow(r_0, 11)*pow(x_0, 3)*pow(sin(x_1), 5) + 5078673600*pow(M_PI, 131.0L/2.0L)*pow(L, 4)*pow(a, 2)*psi_val_00*pow(r_0, 11)*pow(x_0, 3)*pow(sin(x_1), 3) + 44000556600*pow(M_PI, 131.0L/2.0L)*pow(L, 4)*pow(a, 2)*psi_val_00*pow(r_0, 10)*pow(x_0, 4)*pow(sin(x_1), 5) - 27932704800*pow(M_PI, 131.0L/2.0L)*pow(L, 4)*pow(a, 2)*psi_val_00*pow(r_0, 10)*pow(x_0, 4)*pow(sin(x_1), 3) - 90846756000*pow(M_PI, 131.0L/2.0L)*pow(L, 4)*pow(a, 2)*psi_val_00*pow(r_0, 9)*pow(x_0, 5)*pow(sin(x_1), 5) + 50738688000*pow(M_PI, 131.0L/2.0L)*pow(L, 4)*pow(a, 2)*psi_val_00*pow(r_0, 9)*pow(x_0, 5)*pow(sin(x_1), 3) + 80445865500*pow(M_PI, 131.0L/2.0L)*pow(L, 4)*pow(a, 2)*psi_val_00*pow(r_0, 8)*pow(x_0, 6)*pow(sin(x_1), 5) - 18828810000*pow(M_PI, 131.0L/2.0L)*pow(L, 4)*pow(a, 2)*psi_val_00*pow(r_0, 8)*pow(x_0, 6)*pow(sin(x_1), 3) - 9268459200*pow(M_PI, 131.0L/2.0L)*pow(L, 4)*pow(a, 2)*psi_val_00*pow(r_0, 7)*pow(x_0, 7)*pow(sin(x_1), 5) - 100616630400*pow(M_PI, 131.0L/2.0L)*pow(L, 4)*pow(a, 2)*psi_val_00*pow(r_0, 7)*pow(x_0, 7)*pow(sin(x_1), 3) + 1632430800*pow(M_PI, 131.0L/2.0L)*pow(L, 4)*pow(a, 2)*psi_val_00*pow(r_0, 6)*pow(x_0, 8)*pow(sin(x_1), 5) + 244497052800*pow(M_PI, 131.0L/2.0L)*pow(L, 4)*pow(a, 2)*psi_val_00*pow(r_0, 6)*pow(x_0, 8)*pow(sin(x_1), 3) - 121448527200*pow(M_PI, 131.0L/2.0L)*pow(L, 4)*pow(a, 2)*psi_val_00*pow(r_0, 5)*pow(x_0, 9)*pow(sin(x_1), 5) - 300362462400*pow(M_PI, 131.0L/2.0L)*pow(L, 4)*pow(a, 2)*psi_val_00*pow(r_0, 5)*pow(x_0, 9)*pow(sin(x_1), 3) + 236559823500*pow(M_PI, 131.0L/2.0L)*pow(L, 4)*pow(a, 2)*psi_val_00*pow(r_0, 4)*pow(x_0, 10)*pow(sin(x_1), 5) + 235393158000*pow(M_PI, 131.0L/2.0L)*pow(L, 4)*pow(a, 2)*psi_val_00*pow(r_0, 4)*pow(x_0, 10)*pow(sin(x_1), 3) - 216107892000*pow(M_PI, 131.0L/2.0L)*pow(L, 4)*pow(a, 2)*psi_val_00*pow(r_0, 3)*pow(x_0, 11)*pow(sin(x_1), 5) - 122594472000*pow(M_PI, 131.0L/2.0L)*pow(L, 4)*pow(a, 2)*psi_val_00*pow(r_0, 3)*pow(x_0, 11)*pow(sin(x_1), 3) + 110282772600*pow(M_PI, 131.0L/2.0L)*pow(L, 4)*pow(a, 2)*psi_val_00*pow(r_0, 2)*pow(x_0, 12)*pow(sin(x_1), 5) + 41400559200*pow(M_PI, 131.0L/2.0L)*pow(L, 4)*pow(a, 2)*psi_val_00*pow(r_0, 2)*pow(x_0, 12)*pow(sin(x_1), 3) - 30560930400*pow(M_PI, 131.0L/2.0L)*pow(L, 4)*pow(a, 2)*psi_val_00*r_0*pow(x_0, 13)*pow(sin(x_1), 5) - 8254646400*pow(M_PI, 131.0L/2.0L)*pow(L, 4)*pow(a, 2)*psi_val_00*r_0*pow(x_0, 13)*pow(sin(x_1), 3) + 3611708100*pow(M_PI, 131.0L/2.0L)*pow(L, 4)*pow(a, 2)*psi_val_00*pow(x_0, 14)*pow(sin(x_1), 5) + 740797200*pow(M_PI, 131.0L/2.0L)*pow(L, 4)*pow(a, 2)*psi_val_00*pow(x_0, 14)*pow(sin(x_1), 3) - 45302400*pow(M_PI, 131.0L/2.0L)*pow(L, 2)*pow(a, 2)*psi_val_00*pow(r_0, 16)*pow(sin(x_1), 3) - 181209600*pow(M_PI, 131.0L/2.0L)*pow(L, 2)*pow(a, 2)*psi_val_00*pow(r_0, 15)*x_0*pow(sin(x_1), 5) + 362419200*pow(M_PI, 131.0L/2.0L)*pow(L, 2)*pow(a, 2)*psi_val_00*pow(r_0, 15)*x_0*pow(sin(x_1), 3) + 535134600*pow(M_PI, 131.0L/2.0L)*pow(L, 2)*pow(a, 2)*psi_val_00*pow(r_0, 14)*pow(x_0, 2)*pow(sin(x_1), 5) - 713512800*pow(M_PI, 131.0L/2.0L)*pow(L, 2)*pow(a, 2)*psi_val_00*pow(r_0, 14)*pow(x_0, 2)*pow(sin(x_1), 3) + 4043239200*pow(M_PI, 131.0L/2.0L)*pow(L, 2)*pow(a, 2)*psi_val_00*pow(r_0, 13)*pow(x_0, 3)*pow(sin(x_1), 5) - 2695492800*pow(M_PI, 131.0L/2.0L)*pow(L, 2)*pow(a, 2)*psi_val_00*pow(r_0, 13)*pow(x_0, 3)*pow(sin(x_1), 3) - 27371143800*pow(M_PI, 131.0L/2.0L)*pow(L, 2)*pow(a, 2)*psi_val_00*pow(r_0, 12)*pow(x_0, 4)*pow(sin(x_1), 5) + 17520703200*pow(M_PI, 131.0L/2.0L)*pow(L, 2)*pow(a, 2)*psi_val_00*pow(r_0, 12)*pow(x_0, 4)*pow(sin(x_1), 3) + 65643177600*pow(M_PI, 131.0L/2.0L)*pow(L, 2)*pow(a, 2)*psi_val_00*pow(r_0, 11)*pow(x_0, 5)*pow(sin(x_1), 5) - 37102665600*pow(M_PI, 131.0L/2.0L)*pow(L, 2)*pow(a, 2)*psi_val_00*pow(r_0, 11)*pow(x_0, 5)*pow(sin(x_1), 3) - 63443179800*pow(M_PI, 131.0L/2.0L)*pow(L, 2)*pow(a, 2)*psi_val_00*pow(r_0, 10)*pow(x_0, 6)*pow(sin(x_1), 5) + 11336925600*pow(M_PI, 131.0L/2.0L)*pow(L, 2)*pow(a, 2)*psi_val_00*pow(r_0, 10)*pow(x_0, 6)*pow(sin(x_1), 3) - 9592783200*pow(M_PI, 131.0L/2.0L)*pow(L, 2)*pow(a, 2)*psi_val_00*pow(r_0, 9)*pow(x_0, 7)*pow(sin(x_1), 5) + 132803985600*pow(M_PI, 131.0L/2.0L)*pow(L, 2)*pow(a, 2)*psi_val_00*pow(r_0, 9)*pow(x_0, 7)*pow(sin(x_1), 3) + 14295738600*pow(M_PI, 131.0L/2.0L)*pow(L, 2)*pow(a, 2)*psi_val_00*pow(r_0, 8)*pow(x_0, 8)*pow(sin(x_1), 5) - 393553274400*pow(M_PI, 131.0L/2.0L)*pow(L, 2)*pow(a, 2)*psi_val_00*pow(r_0, 8)*pow(x_0, 8)*pow(sin(x_1), 3) + 251156505600*pow(M_PI, 131.0L/2.0L)*pow(L, 2)*pow(a, 2)*psi_val_00*pow(r_0, 7)*pow(x_0, 9)*pow(sin(x_1), 5) + 634867833600*pow(M_PI, 131.0L/2.0L)*pow(L, 2)*pow(a, 2)*psi_val_00*pow(r_0, 7)*pow(x_0, 9)*pow(sin(x_1), 3) - 703979476200*pow(M_PI, 131.0L/2.0L)*pow(L, 2)*pow(a, 2)*psi_val_00*pow(r_0, 6)*pow(x_0, 10)*pow(sin(x_1), 5) - 691552461600*pow(M_PI, 131.0L/2.0L)*pow(L, 2)*pow(a, 2)*psi_val_00*pow(r_0, 6)*pow(x_0, 10)*pow(sin(x_1), 3) + 965382818400*pow(M_PI, 131.0L/2.0L)*pow(L, 2)*pow(a, 2)*psi_val_00*pow(r_0, 5)*pow(x_0, 11)*pow(sin(x_1), 5) + 537988651200*pow(M_PI, 131.0L/2.0L)*pow(L, 2)*pow(a, 2)*psi_val_00*pow(r_0, 5)*pow(x_0, 11)*pow(sin(x_1), 3) - 819409991400*pow(M_PI, 131.0L/2.0L)*pow(L, 2)*pow(a, 2)*psi_val_00*pow(r_0, 4)*pow(x_0, 12)*pow(sin(x_1), 5) - 301974472800*pow(M_PI, 131.0L/2.0L)*pow(L, 2)*pow(a, 2)*psi_val_00*pow(r_0, 4)*pow(x_0, 12)*pow(sin(x_1), 3) + 453159907200*pow(M_PI, 131.0L/2.0L)*pow(L, 2)*pow(a, 2)*psi_val_00*pow(r_0, 3)*pow(x_0, 13)*pow(sin(x_1), 5) + 120187267200*pow(M_PI, 131.0L/2.0L)*pow(L, 2)*pow(a, 2)*psi_val_00*pow(r_0, 3)*pow(x_0, 13)*pow(sin(x_1), 3) - 160362001800*pow(M_PI, 131.0L/2.0L)*pow(L, 2)*pow(a, 2)*psi_val_00*pow(r_0, 2)*pow(x_0, 14)*pow(sin(x_1), 5) - 32311936800*pow(M_PI, 131.0L/2.0L)*pow(L, 2)*pow(a, 2)*psi_val_00*pow(r_0, 2)*pow(x_0, 14)*pow(sin(x_1), 3) + 33195333600*pow(M_PI, 131.0L/2.0L)*pow(L, 2)*pow(a, 2)*psi_val_00*r_0*pow(x_0, 15)*pow(sin(x_1), 5) + 5277729600*pow(M_PI, 131.0L/2.0L)*pow(L, 2)*pow(a, 2)*psi_val_00*r_0*pow(x_0, 15)*pow(sin(x_1), 3) - 3072069000*pow(M_PI, 131.0L/2.0L)*pow(L, 2)*pow(a, 2)*psi_val_00*pow(x_0, 16)*pow(sin(x_1), 5) - 396396000*pow(M_PI, 131.0L/2.0L)*pow(L, 2)*pow(a, 2)*psi_val_00*pow(x_0, 16)*pow(sin(x_1), 3) + 6921200*pow(M_PI, 131.0L/2.0L)*pow(a, 2)*psi_val_00*pow(r_0, 18)*pow(sin(x_1), 3) + 31145400*pow(M_PI, 131.0L/2.0L)*pow(a, 2)*psi_val_00*pow(r_0, 17)*x_0*pow(sin(x_1), 5) - 62290800*pow(M_PI, 131.0L/2.0L)*pow(a, 2)*psi_val_00*pow(r_0, 17)*x_0*pow(sin(x_1), 3) - 105115725*pow(M_PI, 131.0L/2.0L)*pow(a, 2)*psi_val_00*pow(r_0, 16)*pow(x_0, 2)*pow(sin(x_1), 5) + 140154300*pow(M_PI, 131.0L/2.0L)*pow(a, 2)*psi_val_00*pow(r_0, 16)*pow(x_0, 2)*pow(sin(x_1), 3) - 872071200*pow(M_PI, 131.0L/2.0L)*pow(a, 2)*psi_val_00*pow(r_0, 15)*pow(x_0, 3)*pow(sin(x_1), 5) + 581380800*pow(M_PI, 131.0L/2.0L)*pow(a, 2)*psi_val_00*pow(r_0, 15)*pow(x_0, 3)*pow(sin(x_1), 3) + 6758551800*pow(M_PI, 131.0L/2.0L)*pow(a, 2)*psi_val_00*pow(r_0, 14)*pow(x_0, 4)*pow(sin(x_1), 5) - 4360356000*pow(M_PI, 131.0L/2.0L)*pow(a, 2)*psi_val_00*pow(r_0, 14)*pow(x_0, 4)*pow(sin(x_1), 3) - 18313495200*pow(M_PI, 131.0L/2.0L)*pow(a, 2)*psi_val_00*pow(r_0, 13)*pow(x_0, 5)*pow(sin(x_1), 5) + 10464854400*pow(M_PI, 131.0L/2.0L)*pow(a, 2)*psi_val_00*pow(r_0, 13)*pow(x_0, 5)*pow(sin(x_1), 3) + 18422504100*pow(M_PI, 131.0L/2.0L)*pow(a, 2)*psi_val_00*pow(r_0, 12)*pow(x_0, 6)*pow(sin(x_1), 5) - 1889487600*pow(M_PI, 131.0L/2.0L)*pow(a, 2)*psi_val_00*pow(r_0, 12)*pow(x_0, 6)*pow(sin(x_1), 3) + 11336925600*pow(M_PI, 131.0L/2.0L)*pow(a, 2)*psi_val_00*pow(r_0, 11)*pow(x_0, 7)*pow(sin(x_1), 5) - 61543310400*pow(M_PI, 131.0L/2.0L)*pow(a, 2)*psi_val_00*pow(r_0, 11)*pow(x_0, 7)*pow(sin(x_1), 3) - 13361376600*pow(M_PI, 131.0L/2.0L)*pow(a, 2)*psi_val_00*pow(r_0, 10)*pow(x_0, 8)*pow(sin(x_1), 5) + 213782025600*pow(M_PI, 131.0L/2.0L)*pow(a, 2)*psi_val_00*pow(r_0, 10)*pow(x_0, 8)*pow(sin(x_1), 3) - 169244103600*pow(M_PI, 131.0L/2.0L)*pow(a, 2)*psi_val_00*pow(r_0, 9)*pow(x_0, 9)*pow(sin(x_1), 5) - 425584588000*pow(M_PI, 131.0L/2.0L)*pow(a, 2)*psi_val_00*pow(r_0, 9)*pow(x_0, 9)*pow(sin(x_1), 3) + 615736771650*pow(M_PI, 131.0L/2.0L)*pow(a, 2)*psi_val_00*pow(r_0, 8)*pow(x_0, 10)*pow(sin(x_1), 5) + 592354362600*pow(M_PI, 131.0L/2.0L)*pow(a, 2)*psi_val_00*pow(r_0, 8)*pow(x_0, 10)*pow(sin(x_1), 3) - 1122355634400*pow(M_PI, 131.0L/2.0L)*pow(a, 2)*psi_val_00*pow(r_0, 7)*pow(x_0, 11)*pow(sin(x_1), 5) - 612193982400*pow(M_PI, 131.0L/2.0L)*pow(a, 2)*psi_val_00*pow(r_0, 7)*pow(x_0, 11)*pow(sin(x_1), 3) + 1329254526600*pow(M_PI, 131.0L/2.0L)*pow(a, 2)*psi_val_00*pow(r_0, 6)*pow(x_0, 12)*pow(sin(x_1), 5) + 479929850400*pow(M_PI, 131.0L/2.0L)*pow(a, 2)*psi_val_00*pow(r_0, 6)*pow(x_0, 12)*pow(sin(x_1), 3) - 1099681783200*pow(M_PI, 131.0L/2.0L)*pow(a, 2)*psi_val_00*pow(r_0, 5)*pow(x_0, 13)*pow(sin(x_1), 5) - 286039353600*pow(M_PI, 131.0L/2.0L)*pow(a, 2)*psi_val_00*pow(r_0, 5)*pow(x_0, 13)*pow(sin(x_1), 3) + 647185839300*pow(M_PI, 131.0L/2.0L)*pow(a, 2)*psi_val_00*pow(r_0, 4)*pow(x_0, 14)*pow(sin(x_1), 5) + 128007594000*pow(M_PI, 131.0L/2.0L)*pow(a, 2)*psi_val_00*pow(r_0, 4)*pow(x_0, 14)*pow(sin(x_1), 3) - 267476695200*pow(M_PI, 131.0L/2.0L)*pow(a, 2)*psi_val_00*pow(r_0, 3)*pow(x_0, 15)*pow(sin(x_1), 5) - 41776363200*pow(M_PI, 131.0L/2.0L)*pow(a, 2)*psi_val_00*pow(r_0, 3)*pow(x_0, 15)*pow(sin(x_1), 3) + 74157197400*pow(M_PI, 131.0L/2.0L)*pow(a, 2)*psi_val_00*pow(r_0, 2)*pow(x_0, 16)*pow(sin(x_1), 5) + 9405910800*pow(M_PI, 131.0L/2.0L)*pow(a, 2)*psi_val_00*pow(r_0, 2)*pow(x_0, 16)*pow(sin(x_1), 3) - 12427014600*pow(M_PI, 131.0L/2.0L)*pow(a, 2)*psi_val_00*r_0*pow(x_0, 17)*pow(sin(x_1), 5) - 1308106800*pow(M_PI, 131.0L/2.0L)*pow(a, 2)*psi_val_00*r_0*pow(x_0, 17)*pow(sin(x_1), 3) + 953827875*pow(M_PI, 131.0L/2.0L)*pow(a, 2)*psi_val_00*pow(x_0, 18)*pow(sin(x_1), 5) + 84784700*pow(M_PI, 131.0L/2.0L)*pow(a, 2)*psi_val_00*pow(x_0, 18)*pow(sin(x_1), 3))*exp(-5.0L/768.0L*sqrt(143)*a*(6*L*r_0/pow(x_0, 2) - 3*L/x_0 - 40*pow(r_0, 3)/(L*pow(x_0, 2)) + 60*pow(r_0, 2)/(L*x_0) - 20*x_0/L + 84*pow(r_0, 5)/(pow(L, 3)*pow(x_0, 2)) - 210*pow(r_0, 4)/(pow(L, 3)*x_0) + 420*pow(r_0, 2)*x_0/pow(L, 3) - 420*r_0*pow(x_0, 2)/pow(L, 3) + 126*pow(x_0, 3)/pow(L, 3) - 72*pow(r_0, 7)/(pow(L, 5)*pow(x_0, 2)) + 252*pow(r_0, 6)/(pow(L, 5)*x_0) - 1260*pow(r_0, 4)*x_0/pow(L, 5) + 2520*pow(r_0, 3)*pow(x_0, 2)/pow(L, 5) - 2268*pow(r_0, 2)*pow(x_0, 3)/pow(L, 5) + 1008*r_0*pow(x_0, 4)/pow(L, 5) - 180*pow(x_0, 5)/pow(L, 5) + 22*pow(r_0, 9)/(pow(L, 7)*pow(x_0, 2)) - 99*pow(r_0, 8)/(pow(L, 7)*x_0) + 924*pow(r_0, 6)*x_0/pow(L, 7) - 2772*pow(r_0, 5)*pow(x_0, 2)/pow(L, 7) + 4158*pow(r_0, 4)*pow(x_0, 3)/pow(L, 7) - 3696*pow(r_0, 3)*pow(x_0, 4)/pow(L, 7) + 1980*pow(r_0, 2)*pow(x_0, 5)/pow(L, 7) - 594*r_0*pow(x_0, 6)/pow(L, 7) + 77*pow(x_0, 7)/pow(L, 7))*pow(sin(x_1), 2)/sqrt(M_PI))/(pow(M_PI, 133.0L/2.0L)*pow(L, 14)*pow(psi_val_00, 5)*pow(x_0, 6)*sin(x_1)); + return eval_R_scalar_result; + +} + +static inline +double eval_Dcurv_m_r(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 K_rr_val_00, double K_rr_val_01, double K_rr_val_10, double K_phiphi_val_00, double K_phiphi_val_01, double K_phiphi_val_10) { + + double eval_Dcurv_m_r_result; + eval_Dcurv_m_r_result = (1.0L/196608.0L)*(196608*pow(M_PI, 9.0L/2.0L)*K_rr_val_10*pow(L, 14)*psi_val_00*pow(x_0, 8)*sin(x_1) + 2304*sqrt(143)*pow(M_PI, 4)*pow(L, 7)*a*psi_val_00*pow(x_0, 2)*pow(pow(L, 2) - pow(r_0 - x_0, 2), 2)*(5*pow(x_0, 3)*(16*(-pow(L, 2) + pow(r_0 - x_0, 2))*pow(r_0 - x_0, 2) + pow(pow(L, 2) - pow(r_0 - x_0, 2), 2) + 16*pow(r_0 - x_0, 4)) + 5*pow(x_0, 2)*(pow(L, 2) - pow(r_0 - x_0, 2))*(3*pow(L, 2) - 11*pow(r_0 - x_0, 2))*(r_0 - x_0) - x_0*(pow(L, 2) - 11*pow(r_0 - x_0, 2))*pow(pow(L, 2) - pow(r_0 - x_0, 2), 2) + pow(pow(L, 2) - pow(r_0 - x_0, 2), 3)*(-r_0 + x_0))*sin(x_1)*cos(2*x_1) + 256*pow(M_PI, 7.0L/2.0L)*pow(L, 7)*pow(x_0, 5)*(-K_phiphi_val_00 + K_rr_val_00)*(768*M_PI*pow(L, 7)*pow(x_0, 2)*(psi_val_00 + 2*psi_val_10*x_0) + 5*sqrt(143)*sqrt(M_PI)*a*psi_val_00*(pow(L, 2) - pow(r_0 - x_0, 2))*(-x_0*((pow(L, 2) - pow(r_0 - x_0, 2))*(12*x_0*(r_0 - x_0)*(9*pow(L, 2) - 8*pow(-r_0 + x_0, 2) - 25*pow(r_0 - x_0, 2)) + 48*(-pow(L, 2) + pow(r_0 - x_0, 2))*pow(r_0 - x_0, 2) + 3*pow(pow(L, 2) - pow(r_0 - x_0, 2), 2) - 2*(pow(L, 2) - pow(r_0 - x_0, 2))*(3*pow(L, 2) - 11*pow(r_0 - x_0, 2)) + 44*(pow(L, 2) - pow(r_0 - x_0, 2))*pow(r_0 - x_0, 2) + 4*(3*pow(L, 2) - 11*pow(r_0 - x_0, 2))*pow(r_0 - x_0, 2) + 48*pow(r_0 - x_0, 4)) + 4*(r_0 - x_0)*(3*x_0*(16*(-pow(L, 2) + pow(r_0 - x_0, 2))*pow(r_0 - x_0, 2) + pow(pow(L, 2) - pow(r_0 - x_0, 2), 2) + 16*pow(r_0 - x_0, 4)) + 2*(pow(L, 2) - pow(r_0 - x_0, 2))*(3*pow(L, 2) - 11*pow(r_0 - x_0, 2))*(r_0 - x_0))) + 2*(pow(L, 2) - pow(r_0 - x_0, 2))*(3*x_0*(16*(-pow(L, 2) + pow(r_0 - x_0, 2))*pow(r_0 - x_0, 2) + pow(pow(L, 2) - pow(r_0 - x_0, 2), 2) + 16*pow(r_0 - x_0, 4)) + 2*(pow(L, 2) - pow(r_0 - x_0, 2))*(3*pow(L, 2) - 11*pow(r_0 - x_0, 2))*(r_0 - x_0)))*pow(sin(x_1), 2))*sin(x_1) + 128*pow(M_PI, 7.0L/2.0L)*pow(L, 7)*pow(x_0, 5)*(K_phiphi_val_00 + 2*K_rr_val_00)*(1536*M_PI*pow(L, 7)*pow(x_0, 2)*(psi_val_00 + 2*psi_val_10*x_0) + 5*sqrt(143)*sqrt(M_PI)*a*psi_val_00*(pow(L, 2) - pow(r_0 - x_0, 2))*(x_0*((pow(L, 2) - pow(r_0 - x_0, 2))*(12*x_0*(r_0 - x_0)*(9*pow(L, 2) - 8*pow(-r_0 + x_0, 2) - 25*pow(r_0 - x_0, 2)) + 48*(-pow(L, 2) + pow(r_0 - x_0, 2))*pow(r_0 - x_0, 2) + 3*pow(pow(L, 2) - pow(r_0 - x_0, 2), 2) - 2*(pow(L, 2) - pow(r_0 - x_0, 2))*(3*pow(L, 2) - 11*pow(r_0 - x_0, 2)) + 44*(pow(L, 2) - pow(r_0 - x_0, 2))*pow(r_0 - x_0, 2) + 4*(3*pow(L, 2) - 11*pow(r_0 - x_0, 2))*pow(r_0 - x_0, 2) + 48*pow(r_0 - x_0, 4)) + 4*(r_0 - x_0)*(3*x_0*(16*(-pow(L, 2) + pow(r_0 - x_0, 2))*pow(r_0 - x_0, 2) + pow(pow(L, 2) - pow(r_0 - x_0, 2), 2) + 16*pow(r_0 - x_0, 4)) + 2*(pow(L, 2) - pow(r_0 - x_0, 2))*(3*pow(L, 2) - 11*pow(r_0 - x_0, 2))*(r_0 - x_0))) - 2*(pow(L, 2) - pow(r_0 - x_0, 2))*(3*x_0*(16*(-pow(L, 2) + pow(r_0 - x_0, 2))*pow(r_0 - x_0, 2) + pow(pow(L, 2) - pow(r_0 - x_0, 2), 2) + 16*pow(r_0 - x_0, 4)) + 2*(pow(L, 2) - pow(r_0 - x_0, 2))*(3*pow(L, 2) - 11*pow(r_0 - x_0, 2))*(r_0 - x_0)))*pow(sin(x_1), 2))*sin(x_1) + (3.0L/2.0L)*sqrt(143)*pow(M_PI, 7.0L/2.0L)*a*pow(pow(L, 2) - pow(r_0 - x_0, 2), 2)*(3072*sqrt(M_PI)*pow(L, 7)*psi_val_01*pow(x_0, 2) + 5*sqrt(143)*a*psi_val_00*pow(pow(L, 2) - pow(r_0 - x_0, 2), 2)*(3*x_0*(16*(-pow(L, 2) + pow(r_0 - x_0, 2))*pow(r_0 - x_0, 2) + pow(pow(L, 2) - pow(r_0 - x_0, 2), 2) + 16*pow(r_0 - x_0, 4)) + 2*(pow(L, 2) - pow(r_0 - x_0, 2))*(3*pow(L, 2) - 11*pow(r_0 - x_0, 2))*(r_0 - x_0))*sin(2*x_1))*(5*pow(x_0, 3)*(16*(-pow(L, 2) + pow(r_0 - x_0, 2))*pow(r_0 - x_0, 2) + pow(pow(L, 2) - pow(r_0 - x_0, 2), 2) + 16*pow(r_0 - x_0, 4)) + 5*pow(x_0, 2)*(pow(L, 2) - pow(r_0 - x_0, 2))*(3*pow(L, 2) - 11*pow(r_0 - x_0, 2))*(r_0 - x_0) - x_0*(pow(L, 2) - 11*pow(r_0 - x_0, 2))*pow(pow(L, 2) - pow(r_0 - x_0, 2), 2) + pow(pow(L, 2) - pow(r_0 - x_0, 2), 3)*(-r_0 + x_0))*sin(x_1)*sin(2*x_1) + 3*sqrt(143)*pow(M_PI, 7.0L/2.0L)*a*pow(pow(L, 2) - pow(r_0 - x_0, 2), 2)*(384*sqrt(M_PI)*pow(L, 7)*pow(x_0, 2)*(psi_val_00*cos(x_1) + 2*psi_val_01*sin(x_1)) - 5*sqrt(143)*a*psi_val_00*pow(pow(L, 2) - pow(r_0 - x_0, 2), 2)*(3*x_0*(16*(-pow(L, 2) + pow(r_0 - x_0, 2))*pow(r_0 - x_0, 2) + pow(pow(L, 2) - pow(r_0 - x_0, 2), 2) + 16*pow(r_0 - x_0, 4)) + 2*(pow(L, 2) - pow(r_0 - x_0, 2))*(3*pow(L, 2) - 11*pow(r_0 - x_0, 2))*(r_0 - x_0))*pow(sin(x_1), 2)*cos(x_1))*(5*pow(x_0, 3)*(16*(-pow(L, 2) + pow(r_0 - x_0, 2))*pow(r_0 - x_0, 2) + pow(pow(L, 2) - pow(r_0 - x_0, 2), 2) + 16*pow(r_0 - x_0, 4)) + 5*pow(x_0, 2)*(pow(L, 2) - pow(r_0 - x_0, 2))*(3*pow(L, 2) - 11*pow(r_0 - x_0, 2))*(r_0 - x_0) - x_0*(pow(L, 2) - 11*pow(r_0 - x_0, 2))*pow(pow(L, 2) - pow(r_0 - x_0, 2), 2) + pow(pow(L, 2) - pow(r_0 - x_0, 2), 3)*(-r_0 + x_0))*sin(2*x_1))/(pow(M_PI, 9.0L/2.0L)*pow(L, 14)*psi_val_00*pow(x_0, 8)*sin(x_1)); + return eval_Dcurv_m_r_result; + +} + +static inline +double eval_Dcurv_m_t(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 K_rr_val_00, double K_rr_val_01, double K_rr_val_10, double K_phiphi_val_00, double K_phiphi_val_01, double K_phiphi_val_10) { + + double eval_Dcurv_m_t_result; + eval_Dcurv_m_t_result = (1.0L/786432.0L)*(-786432*pow(M_PI, 5)*pow(L, 14)*psi_val_00*pow(x_0, 7)*(K_phiphi_val_01 + K_rr_val_01)*sin(x_1) - 4608*sqrt(143)*pow(M_PI, 9.0L/2.0L)*pow(L, 7)*a*psi_val_00*pow(x_0, 3)*(pow(L, 2) - pow(r_0 - x_0, 2))*(60*pow(x_0, 3)*(-3*pow(L, 2) + 11*pow(r_0 - x_0, 2))*(pow(L, 2) - pow(r_0 - x_0, 2))*(r_0 - x_0) - 20*pow(x_0, 3)*(r_0 - x_0)*(16*(-pow(L, 2) + pow(r_0 - x_0, 2))*pow(r_0 - x_0, 2) + pow(pow(L, 2) - pow(r_0 - x_0, 2), 2) + 16*pow(r_0 - x_0, 4)) + 5*pow(x_0, 2)*pow(pow(L, 2) - pow(r_0 - x_0, 2), 2)*(3*pow(L, 2) - 11*pow(r_0 - x_0, 2)) - 110*pow(x_0, 2)*pow(pow(L, 2) - pow(r_0 - x_0, 2), 2)*pow(r_0 - x_0, 2) - 30*pow(x_0, 2)*(pow(L, 2) - pow(r_0 - x_0, 2))*(3*pow(L, 2) - 11*pow(r_0 - x_0, 2))*pow(r_0 - x_0, 2) + 15*pow(x_0, 2)*(pow(L, 2) - pow(r_0 - x_0, 2))*(-pow(pow(L, 2) - pow(r_0 - x_0, 2), 2) + 16*(pow(L, 2) - pow(r_0 - x_0, 2))*pow(r_0 - x_0, 2) - 16*pow(r_0 - x_0, 4)) + 8*x_0*(pow(L, 2) - 11*pow(r_0 - x_0, 2))*pow(pow(L, 2) - pow(r_0 - x_0, 2), 2)*(r_0 - x_0) + 22*x_0*pow(pow(L, 2) - pow(r_0 - x_0, 2), 3)*(r_0 - x_0) - 10*x_0*pow(pow(L, 2) - pow(r_0 - x_0, 2), 2)*(3*pow(L, 2) - 11*pow(r_0 - x_0, 2))*(r_0 - x_0) + (pow(L, 2) - 11*pow(r_0 - x_0, 2))*pow(pow(L, 2) - pow(r_0 - x_0, 2), 3) - pow(pow(L, 2) - pow(r_0 - x_0, 2), 4) + 10*pow(pow(L, 2) - pow(r_0 - x_0, 2), 3)*pow(r_0 - x_0, 2))*sin(x_1)*sin(2*x_1) - 18432*sqrt(143)*pow(M_PI, 9.0L/2.0L)*pow(L, 7)*a*psi_val_00*pow(x_0, 2)*pow(pow(L, 2) - pow(r_0 - x_0, 2), 2)*(5*pow(x_0, 3)*(16*(-pow(L, 2) + pow(r_0 - x_0, 2))*pow(r_0 - x_0, 2) + pow(pow(L, 2) - pow(r_0 - x_0, 2), 2) + 16*pow(r_0 - x_0, 4)) + 5*pow(x_0, 2)*(pow(L, 2) - pow(r_0 - x_0, 2))*(3*pow(L, 2) - 11*pow(r_0 - x_0, 2))*(r_0 - x_0) - x_0*(pow(L, 2) - 11*pow(r_0 - x_0, 2))*pow(pow(L, 2) - pow(r_0 - x_0, 2), 2) + pow(pow(L, 2) - pow(r_0 - x_0, 2), 3)*(-r_0 + x_0))*sin(x_1)*sin(2*x_1) - 512*pow(M_PI, 9.0L/2.0L)*pow(L, 7)*pow(x_0, 5)*(K_phiphi_val_00 + 2*K_rr_val_00)*(3072*sqrt(M_PI)*pow(L, 7)*psi_val_01*pow(x_0, 2) + 5*sqrt(143)*a*psi_val_00*pow(pow(L, 2) - pow(r_0 - x_0, 2), 2)*(3*x_0*(16*(-pow(L, 2) + pow(r_0 - x_0, 2))*pow(r_0 - x_0, 2) + pow(pow(L, 2) - pow(r_0 - x_0, 2), 2) + 16*pow(r_0 - x_0, 4)) + 2*(pow(L, 2) - pow(r_0 - x_0, 2))*(3*pow(L, 2) - 11*pow(r_0 - x_0, 2))*(r_0 - x_0))*sin(2*x_1))*sin(x_1) - 2048*pow(M_PI, 9.0L/2.0L)*pow(L, 7)*pow(x_0, 5)*(2*K_phiphi_val_00 + K_rr_val_00)*(384*sqrt(M_PI)*pow(L, 7)*pow(x_0, 2)*(psi_val_00*cos(x_1) + 2*psi_val_01*sin(x_1)) - 5*sqrt(143)*a*psi_val_00*pow(pow(L, 2) - pow(r_0 - x_0, 2), 2)*(3*x_0*(16*(-pow(L, 2) + pow(r_0 - x_0, 2))*pow(r_0 - x_0, 2) + pow(pow(L, 2) - pow(r_0 - x_0, 2), 2) + 16*pow(r_0 - x_0, 4)) + 2*(pow(L, 2) - pow(r_0 - x_0, 2))*(3*pow(L, 2) - 11*pow(r_0 - x_0, 2))*(r_0 - x_0))*pow(sin(x_1), 2)*cos(x_1)) + 6*sqrt(143)*pow(M_PI, 7.0L/2.0L)*a*pow(pow(L, 2) - pow(r_0 - x_0, 2), 2)*(768*M_PI*pow(L, 7)*pow(x_0, 2)*(psi_val_00 + 2*psi_val_10*x_0) + 5*sqrt(143)*sqrt(M_PI)*a*psi_val_00*(pow(L, 2) - pow(r_0 - x_0, 2))*(-x_0*((pow(L, 2) - pow(r_0 - x_0, 2))*(-36*x_0*(-3*pow(L, 2) + 11*pow(r_0 - x_0, 2))*(r_0 - x_0) + 48*(-pow(L, 2) + pow(r_0 - x_0, 2))*pow(r_0 - x_0, 2) + 3*pow(pow(L, 2) - pow(r_0 - x_0, 2), 2) - 2*(pow(L, 2) - pow(r_0 - x_0, 2))*(3*pow(L, 2) - 11*pow(r_0 - x_0, 2)) + 44*(pow(L, 2) - pow(r_0 - x_0, 2))*pow(r_0 - x_0, 2) + 4*(3*pow(L, 2) - 11*pow(r_0 - x_0, 2))*pow(r_0 - x_0, 2) + 48*pow(r_0 - x_0, 4)) + 4*(r_0 - x_0)*(3*x_0*(16*(-pow(L, 2) + pow(r_0 - x_0, 2))*pow(r_0 - x_0, 2) + pow(pow(L, 2) - pow(r_0 - x_0, 2), 2) + 16*pow(r_0 - x_0, 4)) + 2*(pow(L, 2) - pow(r_0 - x_0, 2))*(3*pow(L, 2) - 11*pow(r_0 - x_0, 2))*(r_0 - x_0))) + 2*(pow(L, 2) - pow(r_0 - x_0, 2))*(3*x_0*(16*(-pow(L, 2) + pow(r_0 - x_0, 2))*pow(r_0 - x_0, 2) + pow(pow(L, 2) - pow(r_0 - x_0, 2), 2) + 16*pow(r_0 - x_0, 4)) + 2*(pow(L, 2) - pow(r_0 - x_0, 2))*(3*pow(L, 2) - 11*pow(r_0 - x_0, 2))*(r_0 - x_0)))*pow(sin(x_1), 2))*(5*pow(x_0, 3)*(16*(-pow(L, 2) + pow(r_0 - x_0, 2))*pow(r_0 - x_0, 2) + pow(pow(L, 2) - pow(r_0 - x_0, 2), 2) + 16*pow(r_0 - x_0, 4)) + 5*pow(x_0, 2)*(pow(L, 2) - pow(r_0 - x_0, 2))*(3*pow(L, 2) - 11*pow(r_0 - x_0, 2))*(r_0 - x_0) - x_0*(pow(L, 2) - 11*pow(r_0 - x_0, 2))*pow(pow(L, 2) - pow(r_0 - x_0, 2), 2) + pow(pow(L, 2) - pow(r_0 - x_0, 2), 3)*(-r_0 + x_0))*sin(x_1)*sin(2*x_1) + 3*sqrt(143)*pow(M_PI, 7.0L/2.0L)*a*pow(pow(L, 2) - pow(r_0 - x_0, 2), 2)*(5*pow(x_0, 3)*(16*(-pow(L, 2) + pow(r_0 - x_0, 2))*pow(r_0 - x_0, 2) + pow(pow(L, 2) - pow(r_0 - x_0, 2), 2) + 16*pow(r_0 - x_0, 4)) + 5*pow(x_0, 2)*(pow(L, 2) - pow(r_0 - x_0, 2))*(3*pow(L, 2) - 11*pow(r_0 - x_0, 2))*(r_0 - x_0) - x_0*(pow(L, 2) - 11*pow(r_0 - x_0, 2))*pow(pow(L, 2) - pow(r_0 - x_0, 2), 2) + pow(pow(L, 2) - pow(r_0 - x_0, 2), 3)*(-r_0 + x_0))*(3072*M_PI*pow(L, 7)*psi_val_10*pow(x_0, 3) + 1536*M_PI*pow(L, 7)*pow(x_0, 2)*(psi_val_00 + 2*psi_val_10*x_0) + 5*sqrt(143)*sqrt(M_PI)*a*psi_val_00*(pow(L, 2) - pow(r_0 - x_0, 2))*(x_0*((pow(L, 2) - pow(r_0 - x_0, 2))*(-36*x_0*(-3*pow(L, 2) + 11*pow(r_0 - x_0, 2))*(r_0 - x_0) + 48*(-pow(L, 2) + pow(r_0 - x_0, 2))*pow(r_0 - x_0, 2) + 3*pow(pow(L, 2) - pow(r_0 - x_0, 2), 2) - 2*(pow(L, 2) - pow(r_0 - x_0, 2))*(3*pow(L, 2) - 11*pow(r_0 - x_0, 2)) + 44*(pow(L, 2) - pow(r_0 - x_0, 2))*pow(r_0 - x_0, 2) + 4*(3*pow(L, 2) - 11*pow(r_0 - x_0, 2))*pow(r_0 - x_0, 2) + 48*pow(r_0 - x_0, 4)) + 4*(r_0 - x_0)*(3*x_0*(16*(-pow(L, 2) + pow(r_0 - x_0, 2))*pow(r_0 - x_0, 2) + pow(pow(L, 2) - pow(r_0 - x_0, 2), 2) + 16*pow(r_0 - x_0, 4)) + 2*(pow(L, 2) - pow(r_0 - x_0, 2))*(3*pow(L, 2) - 11*pow(r_0 - x_0, 2))*(r_0 - x_0))) - 2*(pow(L, 2) - pow(r_0 - x_0, 2))*(3*x_0*(16*(-pow(L, 2) + pow(r_0 - x_0, 2))*pow(r_0 - x_0, 2) + pow(pow(L, 2) - pow(r_0 - x_0, 2), 2) + 16*pow(r_0 - x_0, 4)) + 2*(pow(L, 2) - pow(r_0 - x_0, 2))*(3*pow(L, 2) - 11*pow(r_0 - x_0, 2))*(r_0 - x_0)))*pow(sin(x_1), 2) + 5*sqrt(143)*sqrt(M_PI)*a*psi_val_00*(pow(L, 2) - pow(r_0 - x_0, 2))*(x_0*((pow(L, 2) - pow(r_0 - x_0, 2))*(12*x_0*(r_0 - x_0)*(9*pow(L, 2) - 8*pow(-r_0 + x_0, 2) - 25*pow(r_0 - x_0, 2)) + 48*(-pow(L, 2) + pow(r_0 - x_0, 2))*pow(r_0 - x_0, 2) + 3*pow(pow(L, 2) - pow(r_0 - x_0, 2), 2) - 2*(pow(L, 2) - pow(r_0 - x_0, 2))*(3*pow(L, 2) - 11*pow(r_0 - x_0, 2)) + 44*(pow(L, 2) - pow(r_0 - x_0, 2))*pow(r_0 - x_0, 2) + 4*(3*pow(L, 2) - 11*pow(r_0 - x_0, 2))*pow(r_0 - x_0, 2) + 48*pow(r_0 - x_0, 4)) + 4*(r_0 - x_0)*(3*x_0*(16*(-pow(L, 2) + pow(r_0 - x_0, 2))*pow(r_0 - x_0, 2) + pow(pow(L, 2) - pow(r_0 - x_0, 2), 2) + 16*pow(r_0 - x_0, 4)) + 2*(pow(L, 2) - pow(r_0 - x_0, 2))*(3*pow(L, 2) - 11*pow(r_0 - x_0, 2))*(r_0 - x_0))) - 2*(pow(L, 2) - pow(r_0 - x_0, 2))*(3*x_0*(16*(-pow(L, 2) + pow(r_0 - x_0, 2))*pow(r_0 - x_0, 2) + pow(pow(L, 2) - pow(r_0 - x_0, 2), 2) + 16*pow(r_0 - x_0, 4)) + 2*(pow(L, 2) - pow(r_0 - x_0, 2))*(3*pow(L, 2) - 11*pow(r_0 - x_0, 2))*(r_0 - x_0)))*pow(sin(x_1), 2))*sin(x_1)*sin(2*x_1))/(pow(M_PI, 5)*pow(L, 14)*psi_val_00*pow(x_0, 7)*sin(x_1)); + return eval_Dcurv_m_t_result; + +} + +static inline +double eval_R_scalar1(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 K_rr_val_00, double K_rr_val_01, double K_rr_val_10, double K_phiphi_val_00, double K_phiphi_val_01, double K_phiphi_val_10) { double eval_R_scalar_result; eval_R_scalar_result = -(8.0*psi_val_01/tan(x_1) + 8.0*psi_val_02 + 16.0*psi_val_10*x_0 + 8.0*psi_val_20*pow(x_0, 2))/(pow(psi_val_00, 5)*pow(x_0, 2)); @@ -8,19 +35,19 @@ double eval_R_scalar(double a, double L, double x_0, double x_1, double psi_val_ } static inline -double eval_Dcurv_m_r(double a, double L, 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 K_rr_val_00, double K_rr_val_01, double K_rr_val_10, double K_phiphi_val_00, double K_phiphi_val_01, double K_phiphi_val_10) { +double eval_Dcurv_m_r1(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 K_rr_val_00, double K_rr_val_01, double K_rr_val_10, double K_phiphi_val_00, double K_phiphi_val_01, double K_phiphi_val_10) { double eval_Dcurv_m_r_result; - eval_Dcurv_m_r_result = 3.0*K_rr_val_00/x_0 + 6.0*K_rr_val_00*psi_val_10/psi_val_00 + 1.0*K_rr_val_10 - 540.0*sqrt(2)*a*exp(-pow(x_0, 2)/pow(L, 2))*cos(2*x_1)/(sqrt(M_PI)*L*x_0) - 180.0*sqrt(2)*a*exp(-pow(x_0, 2)/pow(L, 2))/(sqrt(M_PI)*L*x_0) - 1080.0*sqrt(2)*a*psi_val_01*exp(-pow(x_0, 2)/pow(L, 2))*sin(2*x_1)/(sqrt(M_PI)*L*psi_val_00*x_0) + 360.0*sqrt(2)*a*x_0*exp(-pow(x_0, 2)/pow(L, 2))*cos(2*x_1)/(sqrt(M_PI)*pow(L, 3)) + 120.0*sqrt(2)*a*x_0*exp(-pow(x_0, 2)/pow(L, 2))/(sqrt(M_PI)*pow(L, 3)) + 720.0*sqrt(2)*a*psi_val_01*x_0*exp(-pow(x_0, 2)/pow(L, 2))*sin(2*x_1)/(sqrt(M_PI)*pow(L, 3)*psi_val_00); + eval_Dcurv_m_r_result = 3.0*K_rr_val_00/x_0 + 6.0*K_rr_val_00*psi_val_10/psi_val_00 + 1.0*K_rr_val_10; return eval_Dcurv_m_r_result; } static inline -double eval_Dcurv_m_t(double a, double L, 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 K_rr_val_00, double K_rr_val_01, double K_rr_val_10, double K_phiphi_val_00, double K_phiphi_val_01, double K_phiphi_val_10) { +double eval_Dcurv_m_t1(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 K_rr_val_00, double K_rr_val_01, double K_rr_val_10, double K_phiphi_val_00, double K_phiphi_val_01, double K_phiphi_val_10) { double eval_Dcurv_m_t_result; - eval_Dcurv_m_t_result = -2.0*K_phiphi_val_00/tan(x_1) - 6.0*K_phiphi_val_00*psi_val_01/psi_val_00 - 1.0*K_phiphi_val_01 - 1.0*K_rr_val_00/tan(x_1) - 6.0*K_rr_val_00*psi_val_01/psi_val_00 - 1.0*K_rr_val_01 - 540.0*sqrt(2)*a*exp(-pow(x_0, 2)/pow(L, 2))*sin(2*x_1)/(sqrt(M_PI)*L) - 1080.0*sqrt(2)*a*psi_val_10*x_0*exp(-pow(x_0, 2)/pow(L, 2))*sin(2*x_1)/(sqrt(M_PI)*L*psi_val_00) + 960.0*sqrt(2)*a*pow(x_0, 2)*exp(-pow(x_0, 2)/pow(L, 2))*sin(2*x_1)/(sqrt(M_PI)*pow(L, 3)) + 720.0*sqrt(2)*a*psi_val_10*pow(x_0, 3)*exp(-pow(x_0, 2)/pow(L, 2))*sin(2*x_1)/(sqrt(M_PI)*pow(L, 3)*psi_val_00) - 240.0*sqrt(2)*a*pow(x_0, 4)*exp(-pow(x_0, 2)/pow(L, 2))*sin(2*x_1)/(sqrt(M_PI)*pow(L, 5)); + eval_Dcurv_m_t_result = -2.0*K_phiphi_val_00/tan(x_1) - 6.0*K_phiphi_val_00*psi_val_01/psi_val_00 - 1.0*K_phiphi_val_01 - 1.0*K_rr_val_00/tan(x_1) - 6.0*K_rr_val_00*psi_val_01/psi_val_00 - 1.0*K_rr_val_01; return eval_Dcurv_m_t_result; } |