\frac{x \cdot y}{z}\frac{\sqrt[3]{y} \cdot \sqrt[3]{y}}{\sqrt[3]{z} \cdot \sqrt[3]{z}} \cdot \frac{x}{\frac{\sqrt[3]{z}}{\sqrt[3]{y}}}double f(double x, double y, double z) {
double r733990 = x;
double r733991 = y;
double r733992 = r733990 * r733991;
double r733993 = z;
double r733994 = r733992 / r733993;
return r733994;
}
double f(double x, double y, double z) {
double r733995 = y;
double r733996 = cbrt(r733995);
double r733997 = r733996 * r733996;
double r733998 = z;
double r733999 = cbrt(r733998);
double r734000 = r733999 * r733999;
double r734001 = r733997 / r734000;
double r734002 = x;
double r734003 = r733999 / r733996;
double r734004 = r734002 / r734003;
double r734005 = r734001 * r734004;
return r734005;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 6.5 |
|---|---|
| Target | 6.2 |
| Herbie | 1.9 |
Initial program 6.5
rmApplied associate-/l*6.0
rmApplied add-cube-cbrt6.8
Applied add-cube-cbrt7.0
Applied times-frac7.0
Applied *-un-lft-identity7.0
Applied times-frac2.0
Simplified1.9
Final simplification1.9
herbie shell --seed 2019350
(FPCore (x y z)
:name "Diagrams.Solve.Tridiagonal:solveCyclicTriDiagonal from diagrams-solve-0.1, A"
:precision binary64
:herbie-target
(if (< z -4.262230790519429e-138) (/ (* x y) z) (if (< z 1.7042130660650472e-164) (/ x (/ z y)) (* (/ x z) y)))
(/ (* x y) z))