\frac{x}{\left(y - z\right) \cdot \left(t - z\right)}\left(\frac{\sqrt[3]{x}}{\sqrt[3]{t - z}} \cdot \frac{\sqrt[3]{x}}{\sqrt[3]{t - z}}\right) \cdot \left(\frac{\frac{1}{\sqrt[3]{y - z}}}{\sqrt[3]{y - z}} \cdot \frac{\frac{\sqrt[3]{x}}{\sqrt[3]{y - z}}}{\sqrt[3]{t - z}}\right)double f(double x, double y, double z, double t) {
double r31088115 = x;
double r31088116 = y;
double r31088117 = z;
double r31088118 = r31088116 - r31088117;
double r31088119 = t;
double r31088120 = r31088119 - r31088117;
double r31088121 = r31088118 * r31088120;
double r31088122 = r31088115 / r31088121;
return r31088122;
}
double f(double x, double y, double z, double t) {
double r31088123 = x;
double r31088124 = cbrt(r31088123);
double r31088125 = t;
double r31088126 = z;
double r31088127 = r31088125 - r31088126;
double r31088128 = cbrt(r31088127);
double r31088129 = r31088124 / r31088128;
double r31088130 = r31088129 * r31088129;
double r31088131 = 1.0;
double r31088132 = y;
double r31088133 = r31088132 - r31088126;
double r31088134 = cbrt(r31088133);
double r31088135 = r31088131 / r31088134;
double r31088136 = r31088135 / r31088134;
double r31088137 = r31088124 / r31088134;
double r31088138 = r31088137 / r31088128;
double r31088139 = r31088136 * r31088138;
double r31088140 = r31088130 * r31088139;
return r31088140;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 7.6 |
|---|---|
| Target | 8.3 |
| Herbie | 1.2 |
Initial program 7.6
rmApplied associate-/r*2.1
rmApplied add-cube-cbrt2.7
Applied *-un-lft-identity2.7
Applied add-cube-cbrt2.8
Applied times-frac2.8
Applied times-frac1.2
Simplified1.3
rmApplied *-un-lft-identity1.3
Applied cbrt-prod1.3
Applied add-cube-cbrt1.4
Applied *-un-lft-identity1.4
Applied cbrt-prod1.4
Applied times-frac1.4
Applied times-frac1.2
Simplified1.2
Final simplification1.2
herbie shell --seed 2019171 +o rules:numerics
(FPCore (x y z t)
:name "Data.Random.Distribution.Triangular:triangularCDF from random-fu-0.2.6.2, B"
:herbie-target
(if (< (/ x (* (- y z) (- t z))) 0.0) (/ (/ x (- y z)) (- t z)) (* x (/ 1.0 (* (- y z) (- t z)))))
(/ x (* (- y z) (- t z))))