\left|\left(ew \cdot \sin t\right) \cdot \cos \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right|
\begin{array}{l}
t_1 := \frac{eh}{ew \cdot \tan t}\\
\left|\left(ew \cdot \sin t\right) \cdot \frac{1}{\sqrt{1 + t_1 \cdot t_1}} + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{\frac{eh}{ew}}{\tan t}\right)\right|
\end{array}
(FPCore (eh ew t) :precision binary64 (fabs (+ (* (* ew (sin t)) (cos (atan (/ (/ eh ew) (tan t))))) (* (* eh (cos t)) (sin (atan (/ (/ eh ew) (tan t))))))))
(FPCore (eh ew t)
:precision binary64
(let* ((t_1 (/ eh (* ew (tan t)))))
(fabs
(+
(* (* ew (sin t)) (/ 1.0 (sqrt (+ 1.0 (* t_1 t_1)))))
(* (* eh (cos t)) (sin (atan (/ (/ eh ew) (tan t)))))))))double code(double eh, double ew, double t) {
return fabs(((ew * sin(t)) * cos(atan((eh / ew) / tan(t)))) + ((eh * cos(t)) * sin(atan((eh / ew) / tan(t)))));
}
double code(double eh, double ew, double t) {
double t_1 = eh / (ew * tan(t));
return fabs(((ew * sin(t)) * (1.0 / sqrt(1.0 + (t_1 * t_1)))) + ((eh * cos(t)) * sin(atan((eh / ew) / tan(t)))));
}



Bits error versus eh



Bits error versus ew



Bits error versus t
Results
Initial program 0.1
rmApplied cos-atan_binary640.1
Simplified0.1
Final simplification0.1
herbie shell --seed 2021205
(FPCore (eh ew t)
:name "Example from Robby"
:precision binary64
(fabs (+ (* (* ew (sin t)) (cos (atan (/ (/ eh ew) (tan t))))) (* (* eh (cos t)) (sin (atan (/ (/ eh ew) (tan t))))))))