\frac{x \cdot x}{y \cdot y} + \frac{z \cdot z}{t \cdot t}\left(\sqrt[3]{\frac{x}{y} \cdot \frac{x}{y}} \cdot \sqrt[3]{\frac{x}{y} \cdot \frac{x}{y}}\right) \cdot \sqrt[3]{\frac{x}{y} \cdot \frac{x}{y}} + \frac{z}{t} \cdot \frac{z}{t}double f(double x, double y, double z, double t) {
double r724909 = x;
double r724910 = r724909 * r724909;
double r724911 = y;
double r724912 = r724911 * r724911;
double r724913 = r724910 / r724912;
double r724914 = z;
double r724915 = r724914 * r724914;
double r724916 = t;
double r724917 = r724916 * r724916;
double r724918 = r724915 / r724917;
double r724919 = r724913 + r724918;
return r724919;
}
double f(double x, double y, double z, double t) {
double r724920 = x;
double r724921 = y;
double r724922 = r724920 / r724921;
double r724923 = r724922 * r724922;
double r724924 = cbrt(r724923);
double r724925 = r724924 * r724924;
double r724926 = r724925 * r724924;
double r724927 = z;
double r724928 = t;
double r724929 = r724927 / r724928;
double r724930 = r724929 * r724929;
double r724931 = r724926 + r724930;
return r724931;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 33.9 |
|---|---|
| Target | 0.4 |
| Herbie | 0.8 |
Initial program 33.9
rmApplied times-frac19.5
rmApplied times-frac0.4
rmApplied add-cube-cbrt0.8
Final simplification0.8
herbie shell --seed 2020042
(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))))