\frac{-\sqrt{\left(\left(2 \cdot \frac{4 \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)}{{\left(x-scale \cdot y-scale\right)}^{2}}\right) \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)\right) \cdot \left(\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) - \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}}\right)}}{\frac{4 \cdot \left(\left(b \cdot a\right) \cdot \left(b \cdot \left(-a\right)\right)\right)}{{\left(x-scale \cdot y-scale\right)}^{2}}}\begin{array}{l}
\mathbf{if}\;y-scale \leq -1.755251793459713 \cdot 10^{-50}:\\
\;\;\;\;\left(\frac{a \cdot \left(\sqrt{8} \cdot \cos \left(0.005555555555555556 \cdot \left(angle \cdot \pi\right)\right)\right)}{{y-scale}^{2} \cdot x-scale} \cdot {\left(y-scale \cdot x-scale\right)}^{2}\right) \cdot 0.25\\
\mathbf{elif}\;y-scale \leq 6.879951835079887 \cdot 10^{-110}:\\
\;\;\;\;-0.25 \cdot \log \left({\left(e^{\frac{\sqrt{\frac{\left(a \cdot a\right) \cdot {\cos \left(0.005555555555555556 \cdot \left(angle \cdot \pi\right)\right)}^{2}}{y-scale \cdot y-scale} \cdot \left(\frac{8}{{\left(y-scale \cdot x-scale\right)}^{2}} \cdot {\left(a \cdot b\right)}^{4}\right)}}{\left(a \cdot b\right) \cdot \left(b \cdot \left(-a\right)\right)}}\right)}^{\left({\left(y-scale \cdot x-scale\right)}^{2}\right)}\right)\\
\mathbf{else}:\\
\;\;\;\;\left(a \cdot \left(\sqrt{8} \cdot \left(\cos \left(0.005555555555555556 \cdot \left(angle \cdot \pi\right)\right) \cdot x-scale\right)\right)\right) \cdot 0.25\\
\end{array}(FPCore (a b angle x-scale y-scale)
:precision binary64
(/
(-
(sqrt
(*
(*
(* 2.0 (/ (* 4.0 (* (* b a) (* b (- a)))) (pow (* x-scale y-scale) 2.0)))
(* (* b a) (* b (- a))))
(-
(+
(/
(/
(+
(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))
(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)))))))
(/ (* 4.0 (* (* b a) (* b (- a)))) (pow (* x-scale y-scale) 2.0))))(FPCore (a b angle x-scale y-scale)
:precision binary64
(if (<= y-scale -1.755251793459713e-50)
(*
(*
(/
(* a (* (sqrt 8.0) (cos (* 0.005555555555555556 (* angle PI)))))
(* (pow y-scale 2.0) x-scale))
(pow (* y-scale x-scale) 2.0))
0.25)
(if (<= y-scale 6.879951835079887e-110)
(*
-0.25
(log
(pow
(exp
(/
(sqrt
(*
(/
(* (* a a) (pow (cos (* 0.005555555555555556 (* angle PI))) 2.0))
(* y-scale y-scale))
(* (/ 8.0 (pow (* y-scale x-scale) 2.0)) (pow (* a b) 4.0))))
(* (* a b) (* b (- a)))))
(pow (* y-scale x-scale) 2.0))))
(*
(*
a
(* (sqrt 8.0) (* (cos (* 0.005555555555555556 (* angle PI))) x-scale)))
0.25))))double code(double a, double b, double angle, double x_45_scale, double y_45_scale) {
return -sqrt(((2.0 * ((4.0 * ((b * a) * (b * -a))) / pow((x_45_scale * y_45_scale), 2.0))) * ((b * a) * (b * -a))) * (((((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)) - 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)))) / ((4.0 * ((b * a) * (b * -a))) / pow((x_45_scale * y_45_scale), 2.0));
}
double code(double a, double b, double angle, double x_45_scale, double y_45_scale) {
double tmp;
if (y_45_scale <= -1.755251793459713e-50) {
tmp = (((a * (sqrt(8.0) * cos(0.005555555555555556 * (angle * ((double) M_PI))))) / (pow(y_45_scale, 2.0) * x_45_scale)) * pow((y_45_scale * x_45_scale), 2.0)) * 0.25;
} else if (y_45_scale <= 6.879951835079887e-110) {
tmp = -0.25 * log(pow(exp(sqrt((((a * a) * pow(cos(0.005555555555555556 * (angle * ((double) M_PI))), 2.0)) / (y_45_scale * y_45_scale)) * ((8.0 / pow((y_45_scale * x_45_scale), 2.0)) * pow((a * b), 4.0))) / ((a * b) * (b * -a))), pow((y_45_scale * x_45_scale), 2.0)));
} else {
tmp = (a * (sqrt(8.0) * (cos(0.005555555555555556 * (angle * ((double) M_PI))) * x_45_scale))) * 0.25;
}
return tmp;
}



Bits error versus a



Bits error versus b



Bits error versus angle



Bits error versus x-scale



Bits error versus y-scale
Results
if y-scale < -1.7552517934597131e-50Initial program 63.9
Simplified63.9
Taylor expanded around inf 63.1
Simplified63.1
rmApplied add-cube-cbrt_binary6463.2
Applied sqrt-prod_binary6463.2
Applied cancel-sign-sub-inv_binary6463.2
Simplified63.2
Taylor expanded around inf 61.2
Simplified61.2
Taylor expanded around 0 53.4
if -1.7552517934597131e-50 < y-scale < 6.87995183507988736e-110Initial program 64.0
Simplified64.0
Taylor expanded around inf 64.0
Simplified64.0
rmApplied add-cube-cbrt_binary6464.0
Applied sqrt-prod_binary6464.0
Applied cancel-sign-sub-inv_binary6464.0
Simplified64.0
Taylor expanded around inf 63.5
Simplified63.5
rmApplied add-log-exp_binary6463.5
Simplified43.6
if 6.87995183507988736e-110 < y-scale Initial program 63.9
Simplified63.9
Taylor expanded around inf 63.4
Simplified63.4
rmApplied add-cube-cbrt_binary6463.4
Applied sqrt-prod_binary6463.4
Applied cancel-sign-sub-inv_binary6463.4
Simplified63.4
Taylor expanded around inf 61.1
Simplified61.1
Taylor expanded around 0 53.0
Simplified53.0
Final simplification49.7
herbie shell --seed 2021175
(FPCore (a b angle x-scale y-scale)
:name "b from scale-rotated-ellipse"
:precision binary64
(/ (- (sqrt (* (* (* 2.0 (/ (* 4.0 (* (* b a) (* b (- a)))) (pow (* x-scale y-scale) 2.0))) (* (* b a) (* b (- a)))) (- (+ (/ (/ (+ (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)) (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))))))) (/ (* 4.0 (* (* b a) (* b (- a)))) (pow (* x-scale y-scale) 2.0))))