\frac{x \cdot \left(y - z\right)}{t - z}\begin{array}{l}
\mathbf{if}\;z \le -3.5433744152148751 \cdot 10^{-69} \lor \neg \left(z \le 1.20407222144680094 \cdot 10^{-212}\right):\\
\;\;\;\;\frac{x}{\frac{t - z}{y - z}}\\
\mathbf{else}:\\
\;\;\;\;\frac{x \cdot y + x \cdot \left(-z\right)}{t - z}\\
\end{array}double f(double x, double y, double z, double t) {
double r554677 = x;
double r554678 = y;
double r554679 = z;
double r554680 = r554678 - r554679;
double r554681 = r554677 * r554680;
double r554682 = t;
double r554683 = r554682 - r554679;
double r554684 = r554681 / r554683;
return r554684;
}
double f(double x, double y, double z, double t) {
double r554685 = z;
double r554686 = -3.543374415214875e-69;
bool r554687 = r554685 <= r554686;
double r554688 = 1.204072221446801e-212;
bool r554689 = r554685 <= r554688;
double r554690 = !r554689;
bool r554691 = r554687 || r554690;
double r554692 = x;
double r554693 = t;
double r554694 = r554693 - r554685;
double r554695 = y;
double r554696 = r554695 - r554685;
double r554697 = r554694 / r554696;
double r554698 = r554692 / r554697;
double r554699 = r554692 * r554695;
double r554700 = -r554685;
double r554701 = r554692 * r554700;
double r554702 = r554699 + r554701;
double r554703 = r554702 / r554694;
double r554704 = r554691 ? r554698 : r554703;
return r554704;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 11.6 |
|---|---|
| Target | 2.0 |
| Herbie | 2.2 |
if z < -3.543374415214875e-69 or 1.204072221446801e-212 < z Initial program 13.4
rmApplied associate-/l*1.0
if -3.543374415214875e-69 < z < 1.204072221446801e-212Initial program 5.9
rmApplied sub-neg5.9
Applied distribute-lft-in5.9
Final simplification2.2
herbie shell --seed 2020046 +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)))