x - \frac{y \cdot \left(z - t\right)}{a}\begin{array}{l}
\mathbf{if}\;y \cdot \left(z - t\right) \leq -3.327688885597853 \cdot 10^{+256}:\\
\;\;\;\;x - \left(z - t\right) \cdot \frac{y}{a}\\
\mathbf{elif}\;y \cdot \left(z - t\right) \leq 1.1031091379008978 \cdot 10^{+71}:\\
\;\;\;\;x - \frac{y \cdot z - y \cdot t}{a}\\
\mathbf{else}:\\
\;\;\;\;x - \frac{y}{\frac{a}{z - t}}\\
\end{array}(FPCore (x y z t a) :precision binary64 (- x (/ (* y (- z t)) a)))
(FPCore (x y z t a)
:precision binary64
(if (<= (* y (- z t)) -3.327688885597853e+256)
(- x (* (- z t) (/ y a)))
(if (<= (* y (- z t)) 1.1031091379008978e+71)
(- x (/ (- (* y z) (* y t)) a))
(- x (/ y (/ a (- z t)))))))double code(double x, double y, double z, double t, double a) {
return x - ((y * (z - t)) / a);
}
double code(double x, double y, double z, double t, double a) {
double tmp;
if ((y * (z - t)) <= -3.327688885597853e+256) {
tmp = x - ((z - t) * (y / a));
} else if ((y * (z - t)) <= 1.1031091379008978e+71) {
tmp = x - (((y * z) - (y * t)) / a);
} else {
tmp = x - (y / (a / (z - t)));
}
return tmp;
}














Bits error versus x














Bits error versus y














Bits error versus z














Bits error versus t














Bits error versus a
Results
| Original | 6.2 |
|---|---|
| Target | 0.8 |
| Herbie | 1.1 |
| Alternative 1 | |
|---|---|
| Error | 1.2 |
| Cost | 1730 |
| Alternative 2 | |
|---|---|
| Error | 2.5 |
| Cost | 904 |
| Alternative 3 | |
|---|---|
| Error | 5.6 |
| Cost | 904 |
| Alternative 4 | |
|---|---|
| Error | 9.8 |
| Cost | 840 |
| Alternative 5 | |
|---|---|
| Error | 9.4 |
| Cost | 776 |
| Alternative 6 | |
|---|---|
| Error | 10.1 |
| Cost | 776 |
| Alternative 7 | |
|---|---|
| Error | 12.6 |
| Cost | 776 |
| Alternative 8 | |
|---|---|
| Error | 19.3 |
| Cost | 448 |
| Alternative 9 | |
|---|---|
| Error | 61.9 |
| Cost | 64 |
| Alternative 10 | |
|---|---|
| Error | 61.9 |
| Cost | 64 |

if (*.f64 y (-.f64 z t)) < -3.32768888559785305e256Initial program 44.2
rmApplied associate-/l*_binary64_119580.3
rmApplied associate-/r/_binary64_119590.2
Simplified0.2
if -3.32768888559785305e256 < (*.f64 y (-.f64 z t)) < 1.10310913790089778e71Initial program 0.5
rmApplied sub-neg_binary64_120060.5
Applied distribute-rgt-in_binary64_119630.5
Simplified0.5
Simplified0.5
Simplified0.5
if 1.10310913790089778e71 < (*.f64 y (-.f64 z t)) Initial program 14.4
rmApplied associate-/l*_binary64_119583.5
Simplified3.5
Final simplification1.1
herbie shell --seed 2021044
(FPCore (x y z t a)
:name "Optimisation.CirclePacking:place from circle-packing-0.1.0.4, F"
:precision binary64
:herbie-target
(if (< y -1.0761266216389975e-10) (- x (/ 1.0 (/ (/ a (- z t)) y))) (if (< y 2.894426862792089e-49) (- x (/ (* y (- z t)) a)) (- x (/ y (/ a (- z t))))))
(- x (/ (* y (- z t)) a)))