\frac{x}{\left(y - z\right) \cdot \left(t - z\right)}\frac{\sqrt[3]{x} \cdot \sqrt[3]{x}}{y - z} \cdot \frac{\sqrt[3]{x}}{t - z}double f(double x, double y, double z, double t) {
double r896530 = x;
double r896531 = y;
double r896532 = z;
double r896533 = r896531 - r896532;
double r896534 = t;
double r896535 = r896534 - r896532;
double r896536 = r896533 * r896535;
double r896537 = r896530 / r896536;
return r896537;
}
double f(double x, double y, double z, double t) {
double r896538 = x;
double r896539 = cbrt(r896538);
double r896540 = r896539 * r896539;
double r896541 = y;
double r896542 = z;
double r896543 = r896541 - r896542;
double r896544 = r896540 / r896543;
double r896545 = t;
double r896546 = r896545 - r896542;
double r896547 = r896539 / r896546;
double r896548 = r896544 * r896547;
return r896548;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 7.9 |
|---|---|
| Target | 8.6 |
| Herbie | 1.8 |
Initial program 7.9
rmApplied add-cube-cbrt8.4
Applied times-frac1.8
Final simplification1.8
herbie shell --seed 2019353 +o rules:numerics
(FPCore (x y z t)
:name "Data.Random.Distribution.Triangular:triangularCDF from random-fu-0.2.6.2, B"
:precision binary64
:herbie-target
(if (< (/ x (* (- y z) (- t z))) 0.0) (/ (/ x (- y z)) (- t z)) (* x (/ 1 (* (- y z) (- t z)))))
(/ x (* (- y z) (- t z))))