Incoming wave initial data.incoming_wave

5 files changed, 102 insertions, 71 deletions

diff --git a/Makefile b/Makefile
index be6287a..1933eb6 100644
--- a/Makefile
+++ b/Makefile
@@ -16,7 +16,8 @@ OBJS = basis.o               \
        threadpool.o          \
 
 TESTPROGS = nlsolve          \
-            pssolve
+            pssolve          \
+            td_solve
 
 TESTPROGS := $(addprefix tests/,$(TESTPROGS))
 
diff --git a/eq_gen/constraint_gen.py b/eq_gen/constraint_gen.py
index b7033e8..64003f0 100644
--- a/eq_gen/constraint_gen.py
+++ b/eq_gen/constraint_gen.py
@@ -67,18 +67,18 @@ for diff0 in range(psi_funcs.shape[0]):
 r, theta, phi = coords
 
 x, a, L, r_0 = sympy.symbols('x a L r_0')
-r_0 = 20
 kappa_p = sympy.sqrt(143 / sympy.pi) / 12288
 c_n = (((225 / (8 * sympy.pi)) ** 2)) ** (sympy.Rational(1, 4))
-Im2 = sympy.Piecewise((a * kappa_p * (L ** 5) * ((1 - ((x / L) ** 2)) ** 6), sympy.Abs(x) < L),
-                      (0, True))
+#Im2 = sympy.Piecewise((a * kappa_p * (L ** 5) * ((1 - ((x / L) ** 2)) ** 6), sympy.Abs(x) < L),
+#                      (0, True))
+Im2 = a * kappa_p * (L ** 5) * ((1 - ((x / L) ** 2)) ** 6)
 #Im2 = a * c_n * (L ** 5) * sympy.exp(-(x ** 2) / (L ** 2))
 Im1 = Im2.diff(x)
 I0  = Im1.diff(x)
 Ip1 = I0.diff(x)
 Ip2 = Ip1.diff(x)
 
-v = r# - r_0
+v = r - r_0
 u = -r - r_0
 
 Im2v = Im2.subs(x, v)
@@ -93,8 +93,6 @@ I0u  = I0.subs (x, u)
 Ip1u = Ip1.subs(x, u)
 Ip2u = Ip2.subs(x, u)
 
-Ip1v = sympy.Function('Ip1v')(v)
-Ip2v = Ip1v.diff(v)
 eta      = ((Ip2v / r - 2 * Ip1v / (r ** 2)) * (sympy.sin(theta) ** 2)).simplify()
 K_rtheta = ((Ip2v / (r ** 2) - 3 * Ip1v / (r ** 3) + 6 * I0v / (r ** 4) - 6 * Im1v / (r ** 5)) * r * sympy.sin(2 * theta)).simplify()
 
@@ -103,7 +101,6 @@ K_rtheta = ((Ip2v / (r ** 2) - 3 * Ip1v / (r ** 3) + 6 * I0v / (r ** 4) - 6 * Im
 #K_rtheta = (((Ip2u + Ip2v) / (r ** 2) + 3 * (Ip1u - Ip1v) / (r ** 3) +
 #    6 * (I0u + I0v) / (r ** 4) + 6 * (Im1u - Im1v) / (r ** 5)) * r * sympy.sin(2 * theta)).simplify()
 
-eta = sympy.Function('eta')(r, theta)
 e23eta = sympy.exp(2 * eta / 3)
 metric = sympy.MutableDenseNDimArray([
     [e23eta,                 0,                                           0],
@@ -143,8 +140,6 @@ R_scalar = sympy.Integer(0)
 for i in range(3):
     for j in range(3):
         R_scalar += metric_u[i, j] * Ricci[i, j]
-sympy.pretty_print(R_scalar.simplify(), num_columns = 200)
-sys.exit(0)
 
 curv_m = sympy.MutableDenseNDimArray([
     [ K_rr,                                                 K_rtheta, 0],
@@ -168,29 +163,32 @@ def replace_derivatives(expr, replace):
         expr = expr.subs(replace)
     return expr
 
-R_scalar_val  = replace_derivatives(R_scalar,   replace_diff).subs(replace_val).simplify()
-Dcurv_m_r_val = replace_derivatives(Dcurv_m[0], replace_diff).subs(replace_val).simplify()
-Dcurv_m_t_val = replace_derivatives(Dcurv_m[1], replace_diff).subs(replace_val).simplify()
-
-funcs = []
-
-res = sympy.MatrixSymbol('res', 3, 3)
-funcs.append(('eval_R_scalar',  R_scalar_val))
-funcs.append(('eval_Dcurv_m_r', Dcurv_m_r_val))
-funcs.append(('eval_Dcurv_m_t', Dcurv_m_t_val))
-
-argument_sequence = (
-    a, L, r, theta,
-    *filter(lambda x: x is not None, psi_vals),
-    *filter(lambda x: x is not None, K_rr_vals),
-    *filter(lambda x: x is not None, K_phiphi_vals),
-)
-#funcs = [
-#    ('eval_K_rtheta', K_rtheta),
-#    ('eval_dK_rtheta_r', K_rtheta.diff(r)),
-#    ('eval_dK_rtheta_theta', K_rtheta.diff(theta)),
-#]
-#argument_sequence = (a, L, r, theta)
+#sys.stderr.write('before\n')
+#R_scalar_val  = replace_derivatives(R_scalar,   replace_diff).subs(replace_val).simplify()
+#sys.stderr.write('1\n')
+#Dcurv_m_r_val = replace_derivatives(Dcurv_m[0], replace_diff).subs(replace_val).simplify()
+#sys.stderr.write('2\n')
+#Dcurv_m_t_val = replace_derivatives(Dcurv_m[1], replace_diff).subs(replace_val).simplify()
+#sys.stderr.write('3\n')
+
+#funcs = []
+#
+#funcs.append(('eval_R_scalar',  R_scalar_val))
+#funcs.append(('eval_Dcurv_m_r', Dcurv_m_r_val))
+#funcs.append(('eval_Dcurv_m_t', Dcurv_m_t_val))
+#
+#argument_sequence = (
+#    a, L, r_0, r, theta,
+#    *filter(lambda x: x is not None, psi_vals),
+#    *filter(lambda x: x is not None, K_rr_vals),
+#    *filter(lambda x: x is not None, K_phiphi_vals),
+#)
+funcs = [
+    ('eval_K_rtheta', K_rtheta),
+    ('eval_dK_rtheta_r', K_rtheta.diff(r)),
+    ('eval_dK_rtheta_theta', K_rtheta.diff(theta)),
+]
+argument_sequence = (a, L, r_0, r, theta)
 gen = cg.codegen(funcs, language = 'C', argument_sequence = argument_sequence)
 sys.stdout.write(gen[1][1])
 sys.stdout.write(gen[0][1])
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;
 
 }
diff --git a/eval_k_rtheta.c b/eval_k_rtheta.c
index 0f55c9c..754a463 100644
--- a/eval_k_rtheta.c
+++ b/eval_k_rtheta.c
@@ -1,42 +1,26 @@
-/******************************************************************************
- *                      Code generated with sympy 1.1.1                       *
- *                                                                            *
- *              See http://www.sympy.org/ for more information.               *
- *                                                                            *
- *                       This file is part of 'project'                       *
- ******************************************************************************/
-
-
-#ifndef PROJECT__EVAL_K_RTHETA__H
-#define PROJECT__EVAL_K_RTHETA__H
-
-double eval_K_rtheta(double a, double L, double r_0, double x_0, double x_1);
-double eval_dK_rtheta_r(double a, double L, double r_0, double x_0, double x_1);
-double eval_dK_rtheta_theta(double a, double L, double r_0, double x_0, double x_1);
-
-#endif
-
-
+static inline
 double eval_K_rtheta(double a, double L, double r_0, double x_0, double x_1) {
 
    double eval_K_rtheta_result;
-   eval_K_rtheta_result = 15*sqrt(2)*a*(-3*pow(L, 6)*r_0*exp(4*r_0*x_0/pow(L, 2)) + 3*pow(L, 6)*r_0 + 6*pow(L, 4)*pow(r_0, 2)*x_0*exp(4*r_0*x_0/pow(L, 2)) + 6*pow(L, 4)*pow(r_0, 2)*x_0 - 3*pow(L, 4)*r_0*pow(x_0, 2)*exp(4*r_0*x_0/pow(L, 2)) + 3*pow(L, 4)*r_0*pow(x_0, 2) - 6*pow(L, 2)*pow(r_0, 3)*pow(x_0, 2)*exp(4*r_0*x_0/pow(L, 2)) + 6*pow(L, 2)*pow(r_0, 3)*pow(x_0, 2) + 6*pow(L, 2)*pow(r_0, 2)*pow(x_0, 3)*exp(4*r_0*x_0/pow(L, 2)) + 6*pow(L, 2)*pow(r_0, 2)*pow(x_0, 3) + 6*pow(L, 2)*r_0*pow(x_0, 4)*exp(4*r_0*x_0/pow(L, 2)) - 6*pow(L, 2)*r_0*pow(x_0, 4) - 6*pow(L, 2)*pow(x_0, 5)*exp(4*r_0*x_0/pow(L, 2)) - 6*pow(L, 2)*pow(x_0, 5) + 4*pow(r_0, 4)*pow(x_0, 3)*exp(4*r_0*x_0/pow(L, 2)) + 4*pow(r_0, 4)*pow(x_0, 3) - 16*pow(r_0, 3)*pow(x_0, 4)*exp(4*r_0*x_0/pow(L, 2)) + 16*pow(r_0, 3)*pow(x_0, 4) + 24*pow(r_0, 2)*pow(x_0, 5)*exp(4*r_0*x_0/pow(L, 2)) + 24*pow(r_0, 2)*pow(x_0, 5) - 16*r_0*pow(x_0, 6)*exp(4*r_0*x_0/pow(L, 2)) + 16*r_0*pow(x_0, 6) + 4*pow(x_0, 7)*exp(4*r_0*x_0/pow(L, 2)) + 4*pow(x_0, 7))*exp(-(pow(r_0, 2) + 2*r_0*x_0 + pow(x_0, 2))/pow(L, 2))*sin(2*x_1)/(sqrt(M_PI)*pow(L, 3)*pow(x_0, 4));
+   eval_K_rtheta_result = (3.0L/512.0L)*sqrt(143)*a*(5*pow(x_0, 3)*pow(pow(L, 2) - pow(r_0 - x_0, 2), 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)) + 5*pow(x_0, 2)*pow(pow(L, 2) - pow(r_0 - x_0, 2), 3)*(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), 4) - pow(pow(L, 2) - pow(r_0 - x_0, 2), 5)*(r_0 - x_0))*sin(2*x_1)/(sqrt(M_PI)*pow(L, 7)*pow(x_0, 4));
    return eval_K_rtheta_result;
 
 }
 
+static inline
 double eval_dK_rtheta_r(double a, double L, double r_0, double x_0, double x_1) {
 
    double eval_dK_rtheta_r_result;
-   eval_dK_rtheta_r_result = 15*sqrt(2)*a*(-6*pow(L, 4)*pow(r_0, 2)*exp(4*r_0*x_0/pow(L, 2)) + 6*pow(L, 4)*pow(r_0, 2) - 6*pow(L, 4)*r_0*x_0*exp(4*r_0*x_0/pow(L, 2)) + 6*pow(L, 4)*r_0*x_0 + 12*pow(L, 2)*pow(r_0, 3)*x_0*exp(4*r_0*x_0/pow(L, 2)) + 12*pow(L, 2)*pow(r_0, 3)*x_0 + 6*pow(L, 2)*pow(r_0, 2)*pow(x_0, 2)*exp(4*r_0*x_0/pow(L, 2)) + 18*pow(L, 2)*pow(r_0, 2)*pow(x_0, 2) + 24*pow(L, 2)*r_0*pow(x_0, 3)*exp(4*r_0*x_0/pow(L, 2)) - 24*pow(L, 2)*r_0*pow(x_0, 3) - 30*pow(L, 2)*pow(x_0, 4)*exp(4*r_0*x_0/pow(L, 2)) - 30*pow(L, 2)*pow(x_0, 4) - 12*pow(r_0, 4)*pow(x_0, 2)*exp(4*r_0*x_0/pow(L, 2)) + 12*pow(r_0, 4)*pow(x_0, 2) - 40*pow(r_0, 3)*pow(x_0, 3)*exp(4*r_0*x_0/pow(L, 2)) + 64*pow(r_0, 3)*pow(x_0, 3) + 144*pow(r_0, 2)*pow(x_0, 4)*exp(4*r_0*x_0/pow(L, 2)) + 120*pow(r_0, 2)*pow(x_0, 4) - 120*r_0*pow(x_0, 5)*exp(4*r_0*x_0/pow(L, 2)) + 96*r_0*pow(x_0, 5) + 28*pow(x_0, 6)*exp(4*r_0*x_0/pow(L, 2)) + 28*pow(x_0, 6) + 16*pow(r_0, 5)*pow(x_0, 3)*exp(4*r_0*x_0/pow(L, 2))/pow(L, 2) - 64*pow(r_0, 4)*pow(x_0, 4)*exp(4*r_0*x_0/pow(L, 2))/pow(L, 2) + 96*pow(r_0, 3)*pow(x_0, 5)*exp(4*r_0*x_0/pow(L, 2))/pow(L, 2) - 64*pow(r_0, 2)*pow(x_0, 6)*exp(4*r_0*x_0/pow(L, 2))/pow(L, 2) + 16*r_0*pow(x_0, 7)*exp(4*r_0*x_0/pow(L, 2))/pow(L, 2))*exp(-(pow(r_0, 2) + 2*r_0*x_0 + pow(x_0, 2))/pow(L, 2))*sin(2*x_1)/(sqrt(M_PI)*pow(L, 3)*pow(x_0, 4)) - 60*sqrt(2)*a*(-3*pow(L, 6)*r_0*exp(4*r_0*x_0/pow(L, 2)) + 3*pow(L, 6)*r_0 + 6*pow(L, 4)*pow(r_0, 2)*x_0*exp(4*r_0*x_0/pow(L, 2)) + 6*pow(L, 4)*pow(r_0, 2)*x_0 - 3*pow(L, 4)*r_0*pow(x_0, 2)*exp(4*r_0*x_0/pow(L, 2)) + 3*pow(L, 4)*r_0*pow(x_0, 2) - 6*pow(L, 2)*pow(r_0, 3)*pow(x_0, 2)*exp(4*r_0*x_0/pow(L, 2)) + 6*pow(L, 2)*pow(r_0, 3)*pow(x_0, 2) + 6*pow(L, 2)*pow(r_0, 2)*pow(x_0, 3)*exp(4*r_0*x_0/pow(L, 2)) + 6*pow(L, 2)*pow(r_0, 2)*pow(x_0, 3) + 6*pow(L, 2)*r_0*pow(x_0, 4)*exp(4*r_0*x_0/pow(L, 2)) - 6*pow(L, 2)*r_0*pow(x_0, 4) - 6*pow(L, 2)*pow(x_0, 5)*exp(4*r_0*x_0/pow(L, 2)) - 6*pow(L, 2)*pow(x_0, 5) + 4*pow(r_0, 4)*pow(x_0, 3)*exp(4*r_0*x_0/pow(L, 2)) + 4*pow(r_0, 4)*pow(x_0, 3) - 16*pow(r_0, 3)*pow(x_0, 4)*exp(4*r_0*x_0/pow(L, 2)) + 16*pow(r_0, 3)*pow(x_0, 4) + 24*pow(r_0, 2)*pow(x_0, 5)*exp(4*r_0*x_0/pow(L, 2)) + 24*pow(r_0, 2)*pow(x_0, 5) - 16*r_0*pow(x_0, 6)*exp(4*r_0*x_0/pow(L, 2)) + 16*r_0*pow(x_0, 6) + 4*pow(x_0, 7)*exp(4*r_0*x_0/pow(L, 2)) + 4*pow(x_0, 7))*exp(-(pow(r_0, 2) + 2*r_0*x_0 + pow(x_0, 2))/pow(L, 2))*sin(2*x_1)/(sqrt(M_PI)*pow(L, 3)*pow(x_0, 5)) - 15*sqrt(2)*a*(2*r_0 + 2*x_0)*(-3*pow(L, 6)*r_0*exp(4*r_0*x_0/pow(L, 2)) + 3*pow(L, 6)*r_0 + 6*pow(L, 4)*pow(r_0, 2)*x_0*exp(4*r_0*x_0/pow(L, 2)) + 6*pow(L, 4)*pow(r_0, 2)*x_0 - 3*pow(L, 4)*r_0*pow(x_0, 2)*exp(4*r_0*x_0/pow(L, 2)) + 3*pow(L, 4)*r_0*pow(x_0, 2) - 6*pow(L, 2)*pow(r_0, 3)*pow(x_0, 2)*exp(4*r_0*x_0/pow(L, 2)) + 6*pow(L, 2)*pow(r_0, 3)*pow(x_0, 2) + 6*pow(L, 2)*pow(r_0, 2)*pow(x_0, 3)*exp(4*r_0*x_0/pow(L, 2)) + 6*pow(L, 2)*pow(r_0, 2)*pow(x_0, 3) + 6*pow(L, 2)*r_0*pow(x_0, 4)*exp(4*r_0*x_0/pow(L, 2)) - 6*pow(L, 2)*r_0*pow(x_0, 4) - 6*pow(L, 2)*pow(x_0, 5)*exp(4*r_0*x_0/pow(L, 2)) - 6*pow(L, 2)*pow(x_0, 5) + 4*pow(r_0, 4)*pow(x_0, 3)*exp(4*r_0*x_0/pow(L, 2)) + 4*pow(r_0, 4)*pow(x_0, 3) - 16*pow(r_0, 3)*pow(x_0, 4)*exp(4*r_0*x_0/pow(L, 2)) + 16*pow(r_0, 3)*pow(x_0, 4) + 24*pow(r_0, 2)*pow(x_0, 5)*exp(4*r_0*x_0/pow(L, 2)) + 24*pow(r_0, 2)*pow(x_0, 5) - 16*r_0*pow(x_0, 6)*exp(4*r_0*x_0/pow(L, 2)) + 16*r_0*pow(x_0, 6) + 4*pow(x_0, 7)*exp(4*r_0*x_0/pow(L, 2)) + 4*pow(x_0, 7))*exp(-(pow(r_0, 2) + 2*r_0*x_0 + pow(x_0, 2))/pow(L, 2))*sin(2*x_1)/(sqrt(M_PI)*pow(L, 5)*pow(x_0, 4));
+   eval_dK_rtheta_r_result = (3.0L/512.0L)*sqrt(143)*a*(5*pow(x_0, 3)*pow(pow(L, 2) - pow(r_0 - x_0, 2), 2)*(-16*(pow(L, 2) - pow(r_0 - x_0, 2))*(-2*r_0 + 2*x_0) + (pow(L, 2) - pow(r_0 - x_0, 2))*(4*r_0 - 4*x_0) - 64*pow(r_0 - x_0, 3) - 16*pow(r_0 - x_0, 2)*(2*r_0 - 2*x_0)) + 5*pow(x_0, 3)*(pow(L, 2) - pow(r_0 - x_0, 2))*(4*r_0 - 4*x_0)*(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)) - 5*pow(x_0, 2)*pow(pow(L, 2) - pow(r_0 - x_0, 2), 3)*(3*pow(L, 2) - 11*pow(r_0 - x_0, 2)) + 5*pow(x_0, 2)*pow(pow(L, 2) - pow(r_0 - x_0, 2), 3)*(r_0 - x_0)*(22*r_0 - 22*x_0) + 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))*(r_0 - x_0)*(6*r_0 - 6*x_0) + 15*pow(x_0, 2)*pow(pow(L, 2) - pow(r_0 - x_0, 2), 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)) - x_0*(pow(L, 2) - 11*pow(r_0 - x_0, 2))*pow(pow(L, 2) - pow(r_0 - x_0, 2), 3)*(8*r_0 - 8*x_0) - x_0*pow(pow(L, 2) - pow(r_0 - x_0, 2), 4)*(22*r_0 - 22*x_0) + 10*x_0*pow(pow(L, 2) - pow(r_0 - x_0, 2), 3)*(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), 4) + pow(pow(L, 2) - pow(r_0 - x_0, 2), 5) - pow(pow(L, 2) - pow(r_0 - x_0, 2), 4)*(r_0 - x_0)*(10*r_0 - 10*x_0))*sin(2*x_1)/(sqrt(M_PI)*pow(L, 7)*pow(x_0, 4)) - 3.0L/128.0L*sqrt(143)*a*(5*pow(x_0, 3)*pow(pow(L, 2) - pow(r_0 - x_0, 2), 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)) + 5*pow(x_0, 2)*pow(pow(L, 2) - pow(r_0 - x_0, 2), 3)*(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), 4) - pow(pow(L, 2) - pow(r_0 - x_0, 2), 5)*(r_0 - x_0))*sin(2*x_1)/(sqrt(M_PI)*pow(L, 7)*pow(x_0, 5));
    return eval_dK_rtheta_r_result;
 
 }
 
+static inline
 double eval_dK_rtheta_theta(double a, double L, double r_0, double x_0, double x_1) {
 
    double eval_dK_rtheta_theta_result;
-   eval_dK_rtheta_theta_result = 30*sqrt(2)*a*(-3*pow(L, 6)*r_0*exp(4*r_0*x_0/pow(L, 2)) + 3*pow(L, 6)*r_0 + 6*pow(L, 4)*pow(r_0, 2)*x_0*exp(4*r_0*x_0/pow(L, 2)) + 6*pow(L, 4)*pow(r_0, 2)*x_0 - 3*pow(L, 4)*r_0*pow(x_0, 2)*exp(4*r_0*x_0/pow(L, 2)) + 3*pow(L, 4)*r_0*pow(x_0, 2) - 6*pow(L, 2)*pow(r_0, 3)*pow(x_0, 2)*exp(4*r_0*x_0/pow(L, 2)) + 6*pow(L, 2)*pow(r_0, 3)*pow(x_0, 2) + 6*pow(L, 2)*pow(r_0, 2)*pow(x_0, 3)*exp(4*r_0*x_0/pow(L, 2)) + 6*pow(L, 2)*pow(r_0, 2)*pow(x_0, 3) + 6*pow(L, 2)*r_0*pow(x_0, 4)*exp(4*r_0*x_0/pow(L, 2)) - 6*pow(L, 2)*r_0*pow(x_0, 4) - 6*pow(L, 2)*pow(x_0, 5)*exp(4*r_0*x_0/pow(L, 2)) - 6*pow(L, 2)*pow(x_0, 5) + 4*pow(r_0, 4)*pow(x_0, 3)*exp(4*r_0*x_0/pow(L, 2)) + 4*pow(r_0, 4)*pow(x_0, 3) - 16*pow(r_0, 3)*pow(x_0, 4)*exp(4*r_0*x_0/pow(L, 2)) + 16*pow(r_0, 3)*pow(x_0, 4) + 24*pow(r_0, 2)*pow(x_0, 5)*exp(4*r_0*x_0/pow(L, 2)) + 24*pow(r_0, 2)*pow(x_0, 5) - 16*r_0*pow(x_0, 6)*exp(4*r_0*x_0/pow(L, 2)) + 16*r_0*pow(x_0, 6) + 4*pow(x_0, 7)*exp(4*r_0*x_0/pow(L, 2)) + 4*pow(x_0, 7))*exp(-(pow(r_0, 2) + 2*r_0*x_0 + pow(x_0, 2))/pow(L, 2))*cos(2*x_1)/(sqrt(M_PI)*pow(L, 3)*pow(x_0, 4));
+   eval_dK_rtheta_theta_result = (3.0L/256.0L)*sqrt(143)*a*(5*pow(x_0, 3)*pow(pow(L, 2) - pow(r_0 - x_0, 2), 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)) + 5*pow(x_0, 2)*pow(pow(L, 2) - pow(r_0 - x_0, 2), 3)*(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), 4) - pow(pow(L, 2) - pow(r_0 - x_0, 2), 5)*(r_0 - x_0))*cos(2*x_1)/(sqrt(M_PI)*pow(L, 7)*pow(x_0, 4));
    return eval_dK_rtheta_theta_result;
 
 }
diff --git a/td_constraints_template.c b/td_constraints_template.c
index df09290..5770617 100644
--- a/td_constraints_template.c
+++ b/td_constraints_template.c
@@ -16,8 +16,8 @@
  * along with this program.  If not, see <http://www.gnu.org/licenses/>.
  */
 
-#define L 1.0
-#define R_0 0.0
+#define L 1
+#define R_0 5.0
 
 static void EQNAME(double a,
                    unsigned int nb_coords[2], double *coords[2],
@@ -53,9 +53,14 @@ static void EQNAME(double a,
             //const double K_rtheta_val = K_rtheta(a, r, theta);
             //const double dK_rtheta_r_val = K_rtheta_dr(a, r, theta);
             //const double dK_rtheta_theta_val = K_rtheta_dtheta(a, r, theta);
-            const double K_rtheta_val        = eval_K_rtheta(a, L, R_0, coord0, coord1);
-            const double dK_rtheta_r_val     = eval_dK_rtheta_r(a, L, R_0, coord0, coord1);
-            const double dK_rtheta_theta_val = eval_dK_rtheta_theta(a, L, R_0, coord0, coord1);
+            double K_rtheta_val        = eval_K_rtheta(a, L, R_0, coord0, coord1);
+            double dK_rtheta_r_val     = eval_dK_rtheta_r(a, L, R_0, coord0, coord1);
+            double dK_rtheta_theta_val = eval_dK_rtheta_theta(a, L, R_0, coord0, coord1);
+            if (fabs(r - R_0) > L) {
+                K_rtheta_val = 0.0;
+                dK_rtheta_r_val     = 0.0;
+                dK_rtheta_theta_val = 0.0;
+            }
 
             const double Km[3][3] = {{ K_rr,                               K_rtheta_val,       0 },
                                      { K_rtheta_val / SQR(r), -(K_rr + K_phiphi), 0 },
@@ -157,8 +162,14 @@ static void EQNAME(double a,
 #endif
 
             if (EQUATION == TD_CONSTRAINT_EQ_HAM) {
-                double Rscal = eval_R_scalar(a, L, coord0, coord1, psi, dpsi_01, dpsi_02, dpsi_10, dpsi_11, dpsi_20,
-                                             K_rr, dK_rr_01, dK_rr_10, K_phiphi, dK_phiphi_01, dK_phiphi_10);
+                double Rscal;
+                if (fabs(r - R_0) < L) {
+                    Rscal = eval_R_scalar(a, L, R_0, coord0, coord1, psi, dpsi_01, dpsi_02, dpsi_10, dpsi_11, dpsi_20,
+                                          K_rr, dK_rr_01, dK_rr_10, K_phiphi, dK_phiphi_01, dK_phiphi_10);
+                } else {
+                    Rscal = eval_R_scalar1(a, L, R_0, coord0, coord1, psi, dpsi_01, dpsi_02, dpsi_10, dpsi_11, dpsi_20,
+                                          K_rr, dK_rr_01, dK_rr_10, K_phiphi, dK_phiphi_01, dK_phiphi_10);
+                }
                 double K2 = 0.0;
                 for (int i = 0; i < 3; i++)
                     for (int j = 0; j < 3; j++)
@@ -166,11 +177,21 @@ static void EQNAME(double a,
 
                 out[idx] = Rscal - K2;
             } else if (EQUATION == TD_CONSTRAINT_EQ_MOM_0) {
-                out[idx] = eval_Dcurv_m_r(a, L, coord0, coord1, psi, dpsi_01, dpsi_02, dpsi_10, dpsi_11, dpsi_20,
-                                          K_rr, dK_rr_01, dK_rr_10, K_phiphi, dK_phiphi_01, dK_phiphi_10);
+                if (fabs(r - R_0) < L) {
+                    out[idx] = eval_Dcurv_m_r(a, L, R_0, coord0, coord1, psi, dpsi_01, dpsi_02, dpsi_10, dpsi_11, dpsi_20,
+                                              K_rr, dK_rr_01, dK_rr_10, K_phiphi, dK_phiphi_01, dK_phiphi_10);
+                } else {
+                    out[idx] = eval_Dcurv_m_r1(a, L, R_0, coord0, coord1, psi, dpsi_01, dpsi_02, dpsi_10, dpsi_11, dpsi_20,
+                                              K_rr, dK_rr_01, dK_rr_10, K_phiphi, dK_phiphi_01, dK_phiphi_10);
+                }
             } else {
-                out[idx] = eval_Dcurv_m_t(a, L, coord0, coord1, psi, dpsi_01, dpsi_02, dpsi_10, dpsi_11, dpsi_20,
-                                          K_rr, dK_rr_01, dK_rr_10, K_phiphi, dK_phiphi_01, dK_phiphi_10);
+                if (fabs(r - R_0) < L) {
+                    out[idx] = eval_Dcurv_m_t(a, L, R_0, coord0, coord1, psi, dpsi_01, dpsi_02, dpsi_10, dpsi_11, dpsi_20,
+                                              K_rr, dK_rr_01, dK_rr_10, K_phiphi, dK_phiphi_01, dK_phiphi_10);
+                } else {
+                    out[idx] = eval_Dcurv_m_t1(a, L, R_0, coord0, coord1, psi, dpsi_01, dpsi_02, dpsi_10, dpsi_11, dpsi_20,
+                                              K_rr, dK_rr_01, dK_rr_10, K_phiphi, dK_phiphi_01, dK_phiphi_10);
+                }
             }
 #if 0
             } else {