(FPCore (a b angle x-scale y-scale)
:precision binary64
(*
180.0
(/
(atan
(/
(-
(-
(/
(/
(+
(pow (* a (cos (* (/ angle 180.0) PI))) 2.0)
(pow (* b (sin (* (/ angle 180.0) PI))) 2.0))
y-scale)
y-scale)
(/
(/
(+
(pow (* a (sin (* (/ angle 180.0) PI))) 2.0)
(pow (* b (cos (* (/ angle 180.0) PI))) 2.0))
x-scale)
x-scale))
(sqrt
(+
(pow
(-
(/
(/
(+
(pow (* a (sin (* (/ angle 180.0) PI))) 2.0)
(pow (* b (cos (* (/ angle 180.0) PI))) 2.0))
x-scale)
x-scale)
(/
(/
(+
(pow (* a (cos (* (/ angle 180.0) PI))) 2.0)
(pow (* b (sin (* (/ angle 180.0) PI))) 2.0))
y-scale)
y-scale))
2.0)
(pow
(/
(/
(*
(*
(* 2.0 (- (pow b 2.0) (pow a 2.0)))
(sin (* (/ angle 180.0) PI)))
(cos (* (/ angle 180.0) PI)))
x-scale)
y-scale)
2.0))))
(/
(/
(*
(* (* 2.0 (- (pow b 2.0) (pow a 2.0))) (sin (* (/ angle 180.0) PI)))
(cos (* (/ angle 180.0) PI)))
x-scale)
y-scale)))
PI)))(FPCore (a b angle x-scale y-scale)
:precision binary64
(let* ((t_0 (* PI (* 0.005555555555555556 angle)))
(t_1 (cos t_0))
(t_2 (* (/ y-scale x-scale) (/ t_1 x-scale)))
(t_3 (* 0.005555555555555556 (* PI angle)))
(t_4 (sin t_0))
(t_5 (cos t_3)))
(if (<= a -9.6e-129)
(*
180.0
(/
(atan (/ (* x-scale (* (/ y-scale x-scale) (/ (- t_1) x-scale))) t_4))
PI))
(if (<= a -4e-225)
(*
180.0
(/
(atan (* x-scale (/ (* y-scale (sin t_3)) (* (pow x-scale 2.0) t_5))))
PI))
(if (<= a -5.1e-306)
(*
180.0
(/
(atan
(*
x-scale
(/
(* y-scale (- t_5))
(*
(pow x-scale 2.0)
(cbrt (pow (sin (* angle (* PI 0.005555555555555556))) 3.0))))))
PI))
(if (<= a 1.1e-269)
(*
180.0
(/ (atan (* -180.0 (/ y-scale (* x-scale (* PI angle))))) PI))
(if (<= a 1.5e-183)
(*
180.0
(/
(atan
(fma
-180.0
(/ y-scale (* PI (* x-scale angle)))
(*
0.001851851851851852
(/ (* y-scale PI) (/ x-scale angle)))))
PI))
(if (<= a 5.6e+143)
(* 180.0 (/ (atan (* x-scale (* t_2 (/ -1.0 t_4)))) PI))
(* 180.0 (/ (atan (* x-scale (* t_2 (/ 1.0 t_4)))) PI))))))))))double code(double a, double b, double angle, double x_45_scale, double y_45_scale) {
return 180.0 * (atan(((((((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) - (((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)) - sqrt((pow(((((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)), 2.0) + pow((((((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)))) / (((((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))) / ((double) M_PI));
}
double code(double a, double b, double angle, double x_45_scale, double y_45_scale) {
double t_0 = ((double) M_PI) * (0.005555555555555556 * angle);
double t_1 = cos(t_0);
double t_2 = (y_45_scale / x_45_scale) * (t_1 / x_45_scale);
double t_3 = 0.005555555555555556 * (((double) M_PI) * angle);
double t_4 = sin(t_0);
double t_5 = cos(t_3);
double tmp;
if (a <= -9.6e-129) {
tmp = 180.0 * (atan(((x_45_scale * ((y_45_scale / x_45_scale) * (-t_1 / x_45_scale))) / t_4)) / ((double) M_PI));
} else if (a <= -4e-225) {
tmp = 180.0 * (atan((x_45_scale * ((y_45_scale * sin(t_3)) / (pow(x_45_scale, 2.0) * t_5)))) / ((double) M_PI));
} else if (a <= -5.1e-306) {
tmp = 180.0 * (atan((x_45_scale * ((y_45_scale * -t_5) / (pow(x_45_scale, 2.0) * cbrt(pow(sin((angle * (((double) M_PI) * 0.005555555555555556))), 3.0)))))) / ((double) M_PI));
} else if (a <= 1.1e-269) {
tmp = 180.0 * (atan((-180.0 * (y_45_scale / (x_45_scale * (((double) M_PI) * angle))))) / ((double) M_PI));
} else if (a <= 1.5e-183) {
tmp = 180.0 * (atan(fma(-180.0, (y_45_scale / (((double) M_PI) * (x_45_scale * angle))), (0.001851851851851852 * ((y_45_scale * ((double) M_PI)) / (x_45_scale / angle))))) / ((double) M_PI));
} else if (a <= 5.6e+143) {
tmp = 180.0 * (atan((x_45_scale * (t_2 * (-1.0 / t_4)))) / ((double) M_PI));
} else {
tmp = 180.0 * (atan((x_45_scale * (t_2 * (1.0 / t_4)))) / ((double) M_PI));
}
return tmp;
}
function code(a, b, angle, x_45_scale, y_45_scale) return Float64(180.0 * Float64(atan(Float64(Float64(Float64(Float64(Float64(Float64((Float64(a * cos(Float64(Float64(angle / 180.0) * pi))) ^ 2.0) + (Float64(b * sin(Float64(Float64(angle / 180.0) * pi))) ^ 2.0)) / y_45_scale) / y_45_scale) - Float64(Float64(Float64((Float64(a * sin(Float64(Float64(angle / 180.0) * pi))) ^ 2.0) + (Float64(b * cos(Float64(Float64(angle / 180.0) * pi))) ^ 2.0)) / x_45_scale) / x_45_scale)) - sqrt(Float64((Float64(Float64(Float64(Float64((Float64(a * sin(Float64(Float64(angle / 180.0) * pi))) ^ 2.0) + (Float64(b * cos(Float64(Float64(angle / 180.0) * pi))) ^ 2.0)) / x_45_scale) / x_45_scale) - Float64(Float64(Float64((Float64(a * cos(Float64(Float64(angle / 180.0) * pi))) ^ 2.0) + (Float64(b * sin(Float64(Float64(angle / 180.0) * pi))) ^ 2.0)) / y_45_scale) / y_45_scale)) ^ 2.0) + (Float64(Float64(Float64(Float64(Float64(2.0 * Float64((b ^ 2.0) - (a ^ 2.0))) * sin(Float64(Float64(angle / 180.0) * pi))) * cos(Float64(Float64(angle / 180.0) * pi))) / x_45_scale) / y_45_scale) ^ 2.0)))) / Float64(Float64(Float64(Float64(Float64(2.0 * Float64((b ^ 2.0) - (a ^ 2.0))) * sin(Float64(Float64(angle / 180.0) * pi))) * cos(Float64(Float64(angle / 180.0) * pi))) / x_45_scale) / y_45_scale))) / pi)) end
function code(a, b, angle, x_45_scale, y_45_scale) t_0 = Float64(pi * Float64(0.005555555555555556 * angle)) t_1 = cos(t_0) t_2 = Float64(Float64(y_45_scale / x_45_scale) * Float64(t_1 / x_45_scale)) t_3 = Float64(0.005555555555555556 * Float64(pi * angle)) t_4 = sin(t_0) t_5 = cos(t_3) tmp = 0.0 if (a <= -9.6e-129) tmp = Float64(180.0 * Float64(atan(Float64(Float64(x_45_scale * Float64(Float64(y_45_scale / x_45_scale) * Float64(Float64(-t_1) / x_45_scale))) / t_4)) / pi)); elseif (a <= -4e-225) tmp = Float64(180.0 * Float64(atan(Float64(x_45_scale * Float64(Float64(y_45_scale * sin(t_3)) / Float64((x_45_scale ^ 2.0) * t_5)))) / pi)); elseif (a <= -5.1e-306) tmp = Float64(180.0 * Float64(atan(Float64(x_45_scale * Float64(Float64(y_45_scale * Float64(-t_5)) / Float64((x_45_scale ^ 2.0) * cbrt((sin(Float64(angle * Float64(pi * 0.005555555555555556))) ^ 3.0)))))) / pi)); elseif (a <= 1.1e-269) tmp = Float64(180.0 * Float64(atan(Float64(-180.0 * Float64(y_45_scale / Float64(x_45_scale * Float64(pi * angle))))) / pi)); elseif (a <= 1.5e-183) tmp = Float64(180.0 * Float64(atan(fma(-180.0, Float64(y_45_scale / Float64(pi * Float64(x_45_scale * angle))), Float64(0.001851851851851852 * Float64(Float64(y_45_scale * pi) / Float64(x_45_scale / angle))))) / pi)); elseif (a <= 5.6e+143) tmp = Float64(180.0 * Float64(atan(Float64(x_45_scale * Float64(t_2 * Float64(-1.0 / t_4)))) / pi)); else tmp = Float64(180.0 * Float64(atan(Float64(x_45_scale * Float64(t_2 * Float64(1.0 / t_4)))) / pi)); end return tmp end
code[a_, b_, angle_, x$45$scale_, y$45$scale_] := N[(180.0 * N[(N[ArcTan[N[(N[(N[(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] - 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[Sqrt[N[(N[Power[N[(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] - 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], 2.0], $MachinePrecision] + N[Power[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], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $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]], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision]
code[a_, b_, angle_, x$45$scale_, y$45$scale_] := Block[{t$95$0 = N[(Pi * N[(0.005555555555555556 * angle), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[Cos[t$95$0], $MachinePrecision]}, Block[{t$95$2 = N[(N[(y$45$scale / x$45$scale), $MachinePrecision] * N[(t$95$1 / x$45$scale), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(0.005555555555555556 * N[(Pi * angle), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$4 = N[Sin[t$95$0], $MachinePrecision]}, Block[{t$95$5 = N[Cos[t$95$3], $MachinePrecision]}, If[LessEqual[a, -9.6e-129], N[(180.0 * N[(N[ArcTan[N[(N[(x$45$scale * N[(N[(y$45$scale / x$45$scale), $MachinePrecision] * N[((-t$95$1) / x$45$scale), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / t$95$4), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, -4e-225], N[(180.0 * N[(N[ArcTan[N[(x$45$scale * N[(N[(y$45$scale * N[Sin[t$95$3], $MachinePrecision]), $MachinePrecision] / N[(N[Power[x$45$scale, 2.0], $MachinePrecision] * t$95$5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, -5.1e-306], N[(180.0 * N[(N[ArcTan[N[(x$45$scale * N[(N[(y$45$scale * (-t$95$5)), $MachinePrecision] / N[(N[Power[x$45$scale, 2.0], $MachinePrecision] * N[Power[N[Power[N[Sin[N[(angle * N[(Pi * 0.005555555555555556), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], 3.0], $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 1.1e-269], N[(180.0 * N[(N[ArcTan[N[(-180.0 * N[(y$45$scale / N[(x$45$scale * N[(Pi * angle), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 1.5e-183], N[(180.0 * N[(N[ArcTan[N[(-180.0 * N[(y$45$scale / N[(Pi * N[(x$45$scale * angle), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(0.001851851851851852 * N[(N[(y$45$scale * Pi), $MachinePrecision] / N[(x$45$scale / angle), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 5.6e+143], N[(180.0 * N[(N[ArcTan[N[(x$45$scale * N[(t$95$2 * N[(-1.0 / t$95$4), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision], N[(180.0 * N[(N[ArcTan[N[(x$45$scale * N[(t$95$2 * N[(1.0 / t$95$4), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision]), $MachinePrecision]]]]]]]]]]]]]
180 \cdot \frac{\tan^{-1} \left(\frac{\left(\frac{\frac{{\left(a \cdot \cos \left(\frac{angle}{180} \cdot \pi\right)\right)}^{2} + {\left(b \cdot \sin \left(\frac{angle}{180} \cdot \pi\right)\right)}^{2}}{y-scale}}{y-scale} - \frac{\frac{{\left(a \cdot \sin \left(\frac{angle}{180} \cdot \pi\right)\right)}^{2} + {\left(b \cdot \cos \left(\frac{angle}{180} \cdot \pi\right)\right)}^{2}}{x-scale}}{x-scale}\right) - \sqrt{{\left(\frac{\frac{{\left(a \cdot \sin \left(\frac{angle}{180} \cdot \pi\right)\right)}^{2} + {\left(b \cdot \cos \left(\frac{angle}{180} \cdot \pi\right)\right)}^{2}}{x-scale}}{x-scale} - \frac{\frac{{\left(a \cdot \cos \left(\frac{angle}{180} \cdot \pi\right)\right)}^{2} + {\left(b \cdot \sin \left(\frac{angle}{180} \cdot \pi\right)\right)}^{2}}{y-scale}}{y-scale}\right)}^{2} + {\left(\frac{\frac{\left(\left(2 \cdot \left({b}^{2} - {a}^{2}\right)\right) \cdot \sin \left(\frac{angle}{180} \cdot \pi\right)\right) \cdot \cos \left(\frac{angle}{180} \cdot \pi\right)}{x-scale}}{y-scale}\right)}^{2}}}{\frac{\frac{\left(\left(2 \cdot \left({b}^{2} - {a}^{2}\right)\right) \cdot \sin \left(\frac{angle}{180} \cdot \pi\right)\right) \cdot \cos \left(\frac{angle}{180} \cdot \pi\right)}{x-scale}}{y-scale}}\right)}{\pi}
\begin{array}{l}
t_0 := \pi \cdot \left(0.005555555555555556 \cdot angle\right)\\
t_1 := \cos t_0\\
t_2 := \frac{y-scale}{x-scale} \cdot \frac{t_1}{x-scale}\\
t_3 := 0.005555555555555556 \cdot \left(\pi \cdot angle\right)\\
t_4 := \sin t_0\\
t_5 := \cos t_3\\
\mathbf{if}\;a \leq -9.6 \cdot 10^{-129}:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} \left(\frac{x-scale \cdot \left(\frac{y-scale}{x-scale} \cdot \frac{-t_1}{x-scale}\right)}{t_4}\right)}{\pi}\\
\mathbf{elif}\;a \leq -4 \cdot 10^{-225}:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} \left(x-scale \cdot \frac{y-scale \cdot \sin t_3}{{x-scale}^{2} \cdot t_5}\right)}{\pi}\\
\mathbf{elif}\;a \leq -5.1 \cdot 10^{-306}:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} \left(x-scale \cdot \frac{y-scale \cdot \left(-t_5\right)}{{x-scale}^{2} \cdot \sqrt[3]{{\sin \left(angle \cdot \left(\pi \cdot 0.005555555555555556\right)\right)}^{3}}}\right)}{\pi}\\
\mathbf{elif}\;a \leq 1.1 \cdot 10^{-269}:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} \left(-180 \cdot \frac{y-scale}{x-scale \cdot \left(\pi \cdot angle\right)}\right)}{\pi}\\
\mathbf{elif}\;a \leq 1.5 \cdot 10^{-183}:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} \left(\mathsf{fma}\left(-180, \frac{y-scale}{\pi \cdot \left(x-scale \cdot angle\right)}, 0.001851851851851852 \cdot \frac{y-scale \cdot \pi}{\frac{x-scale}{angle}}\right)\right)}{\pi}\\
\mathbf{elif}\;a \leq 5.6 \cdot 10^{+143}:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} \left(x-scale \cdot \left(t_2 \cdot \frac{-1}{t_4}\right)\right)}{\pi}\\
\mathbf{else}:\\
\;\;\;\;180 \cdot \frac{\tan^{-1} \left(x-scale \cdot \left(t_2 \cdot \frac{1}{t_4}\right)\right)}{\pi}\\
\end{array}



Bits error versus a



Bits error versus b



Bits error versus angle



Bits error versus x-scale



Bits error versus y-scale
if a < -9.59999999999999954e-129Initial program 56.6
Simplified53.6
Taylor expanded in x-scale around 0 52.4
Simplified52.4
Taylor expanded in angle around inf 52.3
Taylor expanded in a around 0 40.3
Applied egg-rr37.5
if -9.59999999999999954e-129 < a < -3.9999999999999998e-225Initial program 53.4
Simplified53.2
Taylor expanded in x-scale around 0 48.3
Simplified48.3
Taylor expanded in angle around inf 48.5
Taylor expanded in a around inf 43.5
if -3.9999999999999998e-225 < a < -5.09999999999999972e-306Initial program 51.8
Simplified52.2
Taylor expanded in x-scale around 0 48.3
Simplified48.3
Taylor expanded in angle around inf 48.5
Taylor expanded in a around 0 27.7
Applied egg-rr29.7
if -5.09999999999999972e-306 < a < 1.09999999999999992e-269Initial program 52.2
Simplified51.3
Taylor expanded in x-scale around 0 48.2
Simplified48.2
Taylor expanded in angle around inf 48.5
Taylor expanded in a around 0 28.3
Taylor expanded in angle around 0 25.8
if 1.09999999999999992e-269 < a < 1.4999999999999999e-183Initial program 52.9
Simplified51.5
Taylor expanded in x-scale around 0 47.2
Simplified47.2
Taylor expanded in angle around inf 47.4
Taylor expanded in a around 0 28.8
Taylor expanded in angle around 0 31.1
Simplified31.1
if 1.4999999999999999e-183 < a < 5.59999999999999996e143Initial program 52.1
Simplified49.0
Taylor expanded in x-scale around 0 46.0
Simplified46.0
Taylor expanded in angle around inf 46.2
Taylor expanded in a around 0 36.5
Applied egg-rr34.0
if 5.59999999999999996e143 < a Initial program 63.7
Simplified63.0
Taylor expanded in x-scale around 0 63.5
Simplified63.5
Taylor expanded in angle around inf 63.5
Taylor expanded in a around 0 44.5
Applied egg-rr36.4
Final simplification35.5
herbie shell --seed 2022162
(FPCore (a b angle x-scale y-scale)
:name "raw-angle from scale-rotated-ellipse"
:precision binary64
(* 180.0 (/ (atan (/ (- (- (/ (/ (+ (pow (* a (cos (* (/ angle 180.0) PI))) 2.0) (pow (* b (sin (* (/ angle 180.0) PI))) 2.0)) y-scale) y-scale) (/ (/ (+ (pow (* a (sin (* (/ angle 180.0) PI))) 2.0) (pow (* b (cos (* (/ angle 180.0) PI))) 2.0)) x-scale) x-scale)) (sqrt (+ (pow (- (/ (/ (+ (pow (* a (sin (* (/ angle 180.0) PI))) 2.0) (pow (* b (cos (* (/ angle 180.0) PI))) 2.0)) x-scale) x-scale) (/ (/ (+ (pow (* a (cos (* (/ angle 180.0) PI))) 2.0) (pow (* b (sin (* (/ angle 180.0) PI))) 2.0)) y-scale) y-scale)) 2.0) (pow (/ (/ (* (* (* 2.0 (- (pow b 2.0) (pow a 2.0))) (sin (* (/ angle 180.0) PI))) (cos (* (/ angle 180.0) PI))) x-scale) y-scale) 2.0)))) (/ (/ (* (* (* 2.0 (- (pow b 2.0) (pow a 2.0))) (sin (* (/ angle 180.0) PI))) (cos (* (/ angle 180.0) PI))) x-scale) y-scale))) PI)))