\frac{x}{\left(y - z\right) \cdot \left(t - z\right)}\frac{\frac{\sqrt[3]{x} \cdot \sqrt[3]{x}}{\sqrt[3]{y - z} \cdot \sqrt[3]{y - z}}}{\sqrt[3]{t - z} \cdot \sqrt[3]{t - z}} \cdot \frac{\frac{\sqrt[3]{x}}{\sqrt[3]{y - z}}}{\sqrt[3]{t - z}}double f(double x, double y, double z, double t) {
double r825707 = x;
double r825708 = y;
double r825709 = z;
double r825710 = r825708 - r825709;
double r825711 = t;
double r825712 = r825711 - r825709;
double r825713 = r825710 * r825712;
double r825714 = r825707 / r825713;
return r825714;
}
double f(double x, double y, double z, double t) {
double r825715 = x;
double r825716 = cbrt(r825715);
double r825717 = r825716 * r825716;
double r825718 = y;
double r825719 = z;
double r825720 = r825718 - r825719;
double r825721 = cbrt(r825720);
double r825722 = r825721 * r825721;
double r825723 = r825717 / r825722;
double r825724 = t;
double r825725 = r825724 - r825719;
double r825726 = cbrt(r825725);
double r825727 = r825726 * r825726;
double r825728 = r825723 / r825727;
double r825729 = r825716 / r825721;
double r825730 = r825729 / r825726;
double r825731 = r825728 * r825730;
return r825731;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 7.3 |
|---|---|
| Target | 8.1 |
| Herbie | 1.3 |
Initial program 7.3
rmApplied associate-/r*2.1
rmApplied add-cube-cbrt2.7
Applied add-cube-cbrt3.0
Applied add-cube-cbrt3.1
Applied times-frac3.1
Applied times-frac1.3
Final simplification1.3
herbie shell --seed 2020062 +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))))