double code(double a, double b, double angle, double x_45_scale, double y_45_scale) {
return ((((((2.0 * (pow(b, 2.0) - pow(a, 2.0))) * sin(((angle / 180.0) * ((double) M_PI)))) * cos(((angle / 180.0) * ((double) M_PI)))) / x_45_scale) / y_45_scale) * (((((2.0 * (pow(b, 2.0) - pow(a, 2.0))) * sin(((angle / 180.0) * ((double) M_PI)))) * cos(((angle / 180.0) * ((double) M_PI)))) / x_45_scale) / y_45_scale)) - ((4.0 * (((pow((a * sin(((angle / 180.0) * ((double) M_PI)))), 2.0) + pow((b * cos(((angle / 180.0) * ((double) M_PI)))), 2.0)) / x_45_scale) / x_45_scale)) * (((pow((a * cos(((angle / 180.0) * ((double) M_PI)))), 2.0) + pow((b * sin(((angle / 180.0) * ((double) M_PI)))), 2.0)) / y_45_scale) / y_45_scale));
}
↓
double code(double a, double b, double angle, double x_45_scale, double y_45_scale) {
double t_0 = (angle / 180.0) * ((double) M_PI);
double t_1 = sin(t_0);
double t_2 = cos(t_0);
double t_3 = ((((2.0 * (pow(b, 2.0) - pow(a, 2.0))) * t_1) * t_2) / x_45_scale) / y_45_scale;
double t_4 = (t_3 * t_3) - ((4.0 * (((pow((a * t_1), 2.0) + pow((b * t_2), 2.0)) / x_45_scale) / x_45_scale)) * (((pow((a * t_2), 2.0) + pow((b * t_1), 2.0)) / y_45_scale) / y_45_scale));
double t_5 = angle * (((double) M_PI) * 0.005555555555555556);
double t_6 = sin(t_5);
double t_7 = log(exp(t_6));
double t_8 = cos(t_5);
double t_9 = pow((a * t_6), 2.0) + pow((b * t_8), 2.0);
double t_10 = pow((a * t_8), 2.0);
double tmp;
if (t_4 <= 0.0) {
tmp = t_4;
} else if (t_4 <= 5e-32) {
tmp = log(exp(fma(((t_10 + pow((b * t_6), 2.0)) / (y_45_scale * y_45_scale)), (((-4.0 / x_45_scale) * t_9) / x_45_scale), pow(((((b + a) / y_45_scale) * ((b - a) / x_45_scale)) * (t_6 * (2.0 * t_8))), 2.0))));
} else {
tmp = fma(((t_9 / (x_45_scale * y_45_scale)) * ((t_10 + pow((b * t_7), 2.0)) / (x_45_scale * y_45_scale))), -4.0, (((4.0 / x_45_scale) * (pow((t_7 * fma(b, b, -(a * a))), 2.0) / (x_45_scale * (y_45_scale / t_8)))) * (t_8 / y_45_scale)));
}
return tmp;
}
code[a_, b_, angle_, x$45$scale_, y$45$scale_] := N[(N[(N[(N[(N[(N[(N[(2.0 * N[(N[Power[b, 2.0], $MachinePrecision] - N[Power[a, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Sin[N[(N[(angle / 180.0), $MachinePrecision] * Pi), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[Cos[N[(N[(angle / 180.0), $MachinePrecision] * Pi), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / x$45$scale), $MachinePrecision] / y$45$scale), $MachinePrecision] * N[(N[(N[(N[(N[(2.0 * N[(N[Power[b, 2.0], $MachinePrecision] - N[Power[a, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Sin[N[(N[(angle / 180.0), $MachinePrecision] * Pi), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[Cos[N[(N[(angle / 180.0), $MachinePrecision] * Pi), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / x$45$scale), $MachinePrecision] / y$45$scale), $MachinePrecision]), $MachinePrecision] - N[(N[(4.0 * N[(N[(N[(N[Power[N[(a * N[Sin[N[(N[(angle / 180.0), $MachinePrecision] * Pi), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision] + N[Power[N[(b * N[Cos[N[(N[(angle / 180.0), $MachinePrecision] * Pi), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] / x$45$scale), $MachinePrecision] / x$45$scale), $MachinePrecision]), $MachinePrecision] * N[(N[(N[(N[Power[N[(a * N[Cos[N[(N[(angle / 180.0), $MachinePrecision] * Pi), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision] + N[Power[N[(b * N[Sin[N[(N[(angle / 180.0), $MachinePrecision] * Pi), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] / y$45$scale), $MachinePrecision] / y$45$scale), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
↓
code[a_, b_, angle_, x$45$scale_, y$45$scale_] := Block[{t$95$0 = N[(N[(angle / 180.0), $MachinePrecision] * Pi), $MachinePrecision]}, Block[{t$95$1 = N[Sin[t$95$0], $MachinePrecision]}, Block[{t$95$2 = N[Cos[t$95$0], $MachinePrecision]}, Block[{t$95$3 = N[(N[(N[(N[(N[(2.0 * N[(N[Power[b, 2.0], $MachinePrecision] - N[Power[a, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * t$95$1), $MachinePrecision] * t$95$2), $MachinePrecision] / x$45$scale), $MachinePrecision] / y$45$scale), $MachinePrecision]}, Block[{t$95$4 = N[(N[(t$95$3 * t$95$3), $MachinePrecision] - N[(N[(4.0 * N[(N[(N[(N[Power[N[(a * t$95$1), $MachinePrecision], 2.0], $MachinePrecision] + N[Power[N[(b * t$95$2), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] / x$45$scale), $MachinePrecision] / x$45$scale), $MachinePrecision]), $MachinePrecision] * N[(N[(N[(N[Power[N[(a * t$95$2), $MachinePrecision], 2.0], $MachinePrecision] + N[Power[N[(b * t$95$1), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] / y$45$scale), $MachinePrecision] / y$45$scale), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$5 = N[(angle * N[(Pi * 0.005555555555555556), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$6 = N[Sin[t$95$5], $MachinePrecision]}, Block[{t$95$7 = N[Log[N[Exp[t$95$6], $MachinePrecision]], $MachinePrecision]}, Block[{t$95$8 = N[Cos[t$95$5], $MachinePrecision]}, Block[{t$95$9 = N[(N[Power[N[(a * t$95$6), $MachinePrecision], 2.0], $MachinePrecision] + N[Power[N[(b * t$95$8), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$10 = N[Power[N[(a * t$95$8), $MachinePrecision], 2.0], $MachinePrecision]}, If[LessEqual[t$95$4, 0.0], t$95$4, If[LessEqual[t$95$4, 5e-32], N[Log[N[Exp[N[(N[(N[(t$95$10 + N[Power[N[(b * t$95$6), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] / N[(y$45$scale * y$45$scale), $MachinePrecision]), $MachinePrecision] * N[(N[(N[(-4.0 / x$45$scale), $MachinePrecision] * t$95$9), $MachinePrecision] / x$45$scale), $MachinePrecision] + N[Power[N[(N[(N[(N[(b + a), $MachinePrecision] / y$45$scale), $MachinePrecision] * N[(N[(b - a), $MachinePrecision] / x$45$scale), $MachinePrecision]), $MachinePrecision] * N[(t$95$6 * N[(2.0 * t$95$8), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision], N[(N[(N[(t$95$9 / N[(x$45$scale * y$45$scale), $MachinePrecision]), $MachinePrecision] * N[(N[(t$95$10 + N[Power[N[(b * t$95$7), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] / N[(x$45$scale * y$45$scale), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * -4.0 + N[(N[(N[(4.0 / x$45$scale), $MachinePrecision] * N[(N[Power[N[(t$95$7 * N[(b * b + (-N[(a * a), $MachinePrecision])), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision] / N[(x$45$scale * N[(y$45$scale / t$95$8), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(t$95$8 / y$45$scale), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]]]]]]]]]