summaryrefslogtreecommitdiff
path: root/eval_k_rtheta.c
diff options
context:
space:
mode:
Diffstat (limited to 'eval_k_rtheta.c')
-rw-r--r--eval_k_rtheta.c28
1 files changed, 6 insertions, 22 deletions
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;
}
generated by cgit v1.2.3 (git 2.39.1) at 2024-05-06 10:55:35 +0000