\frac{x \cdot \left(\left(y - z\right) + 1\right)}{z}\frac{\sqrt[3]{\left(y - z\right) + 1} \cdot \sqrt[3]{\left(y - z\right) + 1}}{\sqrt[3]{z} \cdot \sqrt[3]{z}} \cdot \frac{x}{\frac{\sqrt[3]{z}}{\sqrt[3]{\left(y - z\right) + 1}}}double f(double x, double y, double z) {
double r667136 = x;
double r667137 = y;
double r667138 = z;
double r667139 = r667137 - r667138;
double r667140 = 1.0;
double r667141 = r667139 + r667140;
double r667142 = r667136 * r667141;
double r667143 = r667142 / r667138;
return r667143;
}
double f(double x, double y, double z) {
double r667144 = y;
double r667145 = z;
double r667146 = r667144 - r667145;
double r667147 = 1.0;
double r667148 = r667146 + r667147;
double r667149 = cbrt(r667148);
double r667150 = r667149 * r667149;
double r667151 = cbrt(r667145);
double r667152 = r667151 * r667151;
double r667153 = r667150 / r667152;
double r667154 = x;
double r667155 = r667151 / r667149;
double r667156 = r667154 / r667155;
double r667157 = r667153 * r667156;
return r667157;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 10.2 |
|---|---|
| Target | 0.6 |
| Herbie | 1.4 |
Initial program 10.2
rmApplied associate-/l*3.1
rmApplied add-cube-cbrt3.9
Applied add-cube-cbrt3.8
Applied times-frac3.8
Applied *-un-lft-identity3.8
Applied times-frac1.4
Simplified1.4
Final simplification1.4
herbie shell --seed 2020039
(FPCore (x y z)
:name "Diagrams.TwoD.Segment.Bernstein:evaluateBernstein from diagrams-lib-1.3.0.3"
:precision binary64
:herbie-target
(if (< x -2.71483106713436e-162) (- (* (+ 1 y) (/ x z)) x) (if (< x 3.874108816439546e-197) (* (* x (+ (- y z) 1)) (/ 1 z)) (- (* (+ 1 y) (/ x z)) x)))
(/ (* x (+ (- y z) 1)) z))