\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|\left|ew \cdot \left(\cos \tan^{-1} \left(\frac{eh}{ew \cdot \tan t}\right) \cdot \sin t\right) + \left(eh \cdot \cos t\right) \cdot \sin \tan^{-1} \left(\frac{1}{\sqrt[3]{\tan t} \cdot \sqrt[3]{\tan t}} \cdot \frac{\frac{eh}{ew}}{\sqrt[3]{\tan t}}\right)\right|(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
(fabs
(+
(* ew (* (cos (atan (/ eh (* ew (tan t))))) (sin t)))
(*
(* eh (cos t))
(sin
(atan
(*
(/ 1.0 (* (cbrt (tan t)) (cbrt (tan t))))
(/ (/ eh ew) (cbrt (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) {
return fabs((ew * (cos(atan(eh / (ew * tan(t)))) * sin(t))) + ((eh * cos(t)) * sin(atan((1.0 / (cbrt(tan(t)) * cbrt(tan(t)))) * ((eh / ew) / cbrt(tan(t)))))));
}



Bits error versus eh



Bits error versus ew



Bits error versus t
Results
Initial program 0.1
rmApplied associate-*l*_binary640.1
Simplified0.1
rmApplied add-cube-cbrt_binary640.1
Applied *-un-lft-identity_binary640.1
Applied times-frac_binary640.1
Final simplification0.1
herbie shell --seed 2021118
(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))))))))