\frac{x \cdot x}{y \cdot y} + \frac{z \cdot z}{t \cdot t}\mathsf{fma}\left(\frac{z}{t}, \frac{z}{t}, \sqrt[3]{{\left(\left|\frac{x}{y}\right|\right)}^{\frac{3}{2}}} \cdot {\left(\left|\frac{x}{y}\right|\right)}^{\frac{3}{2}}\right)double code(double x, double y, double z, double t) {
return (((x * x) / (y * y)) + ((z * z) / (t * t)));
}
double code(double x, double y, double z, double t) {
return fma((z / t), (z / t), (cbrt(pow(fabs((x / y)), 1.5)) * pow(fabs((x / y)), 1.5)));
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 33.7 |
|---|---|
| Target | 0.4 |
| Herbie | 0.6 |
Initial program 33.7
Simplified19.3
rmApplied add-sqr-sqrt19.3
Simplified19.3
Simplified0.4
rmApplied add-sqr-sqrt0.5
Applied associate-*l*0.5
Simplified0.5
rmApplied add-cbrt-cube0.6
Simplified0.6
Final simplification0.6
herbie shell --seed 2020049 +o rules:numerics
(FPCore (x y z t)
:name "Graphics.Rasterific.Svg.PathConverter:arcToSegments from rasterific-svg-0.2.3.1"
:precision binary64
:herbie-target
(+ (pow (/ x y) 2) (pow (/ z t) 2))
(+ (/ (* x x) (* y y)) (/ (* z z) (* t t))))