\frac{x \cdot \left(y - z\right)}{t - z}\begin{array}{l}
\mathbf{if}\;z \le -9.5127373023801481 \cdot 10^{-222} \lor \neg \left(z \le 5.43817376472661147 \cdot 10^{-107}\right):\\
\;\;\;\;x \cdot \frac{y - z}{t - z}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{t - z} \cdot \left(y - z\right)\\
\end{array}double f(double x, double y, double z, double t) {
double r550228 = x;
double r550229 = y;
double r550230 = z;
double r550231 = r550229 - r550230;
double r550232 = r550228 * r550231;
double r550233 = t;
double r550234 = r550233 - r550230;
double r550235 = r550232 / r550234;
return r550235;
}
double f(double x, double y, double z, double t) {
double r550236 = z;
double r550237 = -9.512737302380148e-222;
bool r550238 = r550236 <= r550237;
double r550239 = 5.4381737647266115e-107;
bool r550240 = r550236 <= r550239;
double r550241 = !r550240;
bool r550242 = r550238 || r550241;
double r550243 = x;
double r550244 = y;
double r550245 = r550244 - r550236;
double r550246 = t;
double r550247 = r550246 - r550236;
double r550248 = r550245 / r550247;
double r550249 = r550243 * r550248;
double r550250 = r550243 / r550247;
double r550251 = r550250 * r550245;
double r550252 = r550242 ? r550249 : r550251;
return r550252;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 11.9 |
|---|---|
| Target | 2.2 |
| Herbie | 2.4 |
if z < -9.512737302380148e-222 or 5.4381737647266115e-107 < z Initial program 13.3
rmApplied *-un-lft-identity13.3
Applied times-frac1.2
Simplified1.2
if -9.512737302380148e-222 < z < 5.4381737647266115e-107Initial program 6.7
rmApplied associate-/l*5.6
rmApplied associate-/r/6.9
Final simplification2.4
herbie shell --seed 2020047 +o rules:numerics
(FPCore (x y z t)
:name "Graphics.Rendering.Chart.Plot.AreaSpots:renderAreaSpots4D from Chart-1.5.3"
:precision binary64
:herbie-target
(/ x (/ (- t z) (- y z)))
(/ (* x (- y z)) (- t z)))