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




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 33.6 |
|---|---|
| Target | 0.4 |
| Herbie | 1.2 |
Initial program 33.6
Simplified19.2
rmApplied times-frac0.4
rmApplied add-cube-cbrt1.2
rmApplied add-sqr-sqrt1.2
Final simplification1.2
herbie shell --seed 2020121 +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))))