\frac{x \cdot y}{z}\frac{\sqrt[3]{x} \cdot \sqrt[3]{x}}{\frac{\sqrt[3]{z} \cdot \sqrt[3]{z}}{\sqrt[3]{y} \cdot \sqrt[3]{y}}} \cdot \frac{\sqrt[3]{x}}{\frac{\sqrt[3]{z}}{\sqrt[3]{y}}}double f(double x, double y, double z) {
double r838739 = x;
double r838740 = y;
double r838741 = r838739 * r838740;
double r838742 = z;
double r838743 = r838741 / r838742;
return r838743;
}
double f(double x, double y, double z) {
double r838744 = x;
double r838745 = cbrt(r838744);
double r838746 = r838745 * r838745;
double r838747 = z;
double r838748 = cbrt(r838747);
double r838749 = r838748 * r838748;
double r838750 = y;
double r838751 = cbrt(r838750);
double r838752 = r838751 * r838751;
double r838753 = r838749 / r838752;
double r838754 = r838746 / r838753;
double r838755 = r838748 / r838751;
double r838756 = r838745 / r838755;
double r838757 = r838754 * r838756;
return r838757;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 6.5 |
|---|---|
| Target | 6.0 |
| Herbie | 2.1 |
Initial program 6.5
rmApplied associate-/l*6.0
rmApplied add-cube-cbrt6.8
Applied add-cube-cbrt7.0
Applied times-frac7.0
Applied add-cube-cbrt7.1
Applied times-frac2.1
Final simplification2.1
herbie shell --seed 2020047 +o rules:numerics
(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))