\frac{x \cdot x}{y \cdot y} + \frac{z \cdot z}{t \cdot t}\frac{x}{y} \cdot \frac{x}{y} + \left(\frac{z}{t} \cdot \left(\sqrt[3]{\frac{z}{t}} \cdot \left(\sqrt[3]{\frac{1}{\sqrt[3]{t} \cdot \sqrt[3]{t}}} \cdot \sqrt[3]{\frac{z}{\sqrt[3]{t}}}\right)\right)\right) \cdot \sqrt[3]{\frac{z}{t}}double f(double x, double y, double z, double t) {
double r820656 = x;
double r820657 = r820656 * r820656;
double r820658 = y;
double r820659 = r820658 * r820658;
double r820660 = r820657 / r820659;
double r820661 = z;
double r820662 = r820661 * r820661;
double r820663 = t;
double r820664 = r820663 * r820663;
double r820665 = r820662 / r820664;
double r820666 = r820660 + r820665;
return r820666;
}
double f(double x, double y, double z, double t) {
double r820667 = x;
double r820668 = y;
double r820669 = r820667 / r820668;
double r820670 = r820669 * r820669;
double r820671 = z;
double r820672 = t;
double r820673 = r820671 / r820672;
double r820674 = cbrt(r820673);
double r820675 = 1.0;
double r820676 = cbrt(r820672);
double r820677 = r820676 * r820676;
double r820678 = r820675 / r820677;
double r820679 = cbrt(r820678);
double r820680 = r820671 / r820676;
double r820681 = cbrt(r820680);
double r820682 = r820679 * r820681;
double r820683 = r820674 * r820682;
double r820684 = r820673 * r820683;
double r820685 = r820684 * r820674;
double r820686 = r820670 + r820685;
return r820686;
}




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.8 |
Initial program 33.7
rmApplied times-frac18.9
rmApplied times-frac0.4
rmApplied add-cube-cbrt0.8
Applied associate-*r*0.8
rmApplied add-cube-cbrt0.8
Applied *-un-lft-identity0.8
Applied times-frac0.8
Applied cbrt-prod0.8
Final simplification0.8
herbie shell --seed 2020046
(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))))