Optimisation.CirclePacking:place from circle-packing-0.1.0.4, F

Average Error: 6.0 → 0.8

Time: 14.4s

Precision: 64

Math

\[x - \frac{y \cdot \left(z - t\right)}{a}\]

↓

\[\begin{array}{l} \mathbf{if}\;\frac{y \cdot \left(z - t\right)}{a} \le -8.963886535039347471937028690168546917254 \cdot 10^{299} \lor \neg \left(\frac{y \cdot \left(z - t\right)}{a} \le 1.947585218067920906632704854880810808419 \cdot 10^{229}\right):\\ \;\;\;\;x - \frac{y}{a} \cdot \left(z - t\right)\\ \mathbf{else}:\\ \;\;\;\;x - \frac{y \cdot \left(z - t\right)}{a}\\ \end{array}\]

C

double f(double x, double y, double z, double t, double a) {
        double r238891 = x;
        double r238892 = y;
        double r238893 = z;
        double r238894 = t;
        double r238895 = r238893 - r238894;
        double r238896 = r238892 * r238895;
        double r238897 = a;
        double r238898 = r238896 / r238897;
        double r238899 = r238891 - r238898;
        return r238899;
}

↓

double f(double x, double y, double z, double t, double a) {
        double r238900 = y;
        double r238901 = z;
        double r238902 = t;
        double r238903 = r238901 - r238902;
        double r238904 = r238900 * r238903;
        double r238905 = a;
        double r238906 = r238904 / r238905;
        double r238907 = -8.963886535039347e+299;
        bool r238908 = r238906 <= r238907;
        double r238909 = 1.947585218067921e+229;
        bool r238910 = r238906 <= r238909;
        double r238911 = !r238910;
        bool r238912 = r238908 || r238911;
        double r238913 = x;
        double r238914 = r238900 / r238905;
        double r238915 = r238914 * r238903;
        double r238916 = r238913 - r238915;
        double r238917 = r238913 - r238906;
        double r238918 = r238912 ? r238916 : r238917;
        return r238918;
}

Error

Bits error versus x

Bits error versus y

Bits error versus z

Bits error versus t

Bits error versus a

Target

Original	6.0
Target	0.7
Herbie	0.8

\[\begin{array}{l} \mathbf{if}\;y \lt -1.07612662163899753216593153715602325729 \cdot 10^{-10}:\\ \;\;\;\;x - \frac{1}{\frac{\frac{a}{z - t}}{y}}\\ \mathbf{elif}\;y \lt 2.894426862792089097262541964056085749132 \cdot 10^{-49}:\\ \;\;\;\;x - \frac{y \cdot \left(z - t\right)}{a}\\ \mathbf{else}:\\ \;\;\;\;x - \frac{y}{\frac{a}{z - t}}\\ \end{array}\]

Derivation

1. Split input into 2 regimes

2. if (/ (* y (- z t)) a) < -8.963886535039347e+299 or 1.947585218067921e+229 < (/ (* y (- z t)) a)

   1. Initial program 41.8

      \[x - \frac{y \cdot \left(z - t\right)}{a}\]
   2. Using strategy rm

   3. Applied associate-/l* 9.0

      \[\leadsto x - \color{blue}{\frac{y}{\frac{a}{z - t}}}\]
   4. Using strategy rm

   5. Applied associate-/r/ 3.6

      \[\leadsto x - \color{blue}{\frac{y}{a} \cdot \left(z - t\right)}\]

   if -8.963886535039347e+299 < (/ (* y (- z t)) a) < 1.947585218067921e+229

   1. Initial program 0.4

      \[x - \frac{y \cdot \left(z - t\right)}{a}\]
3. Recombined 2 regimes into one program.

4. Final simplification 0.8

   \[\leadsto \begin{array}{l} \mathbf{if}\;\frac{y \cdot \left(z - t\right)}{a} \le -8.963886535039347471937028690168546917254 \cdot 10^{299} \lor \neg \left(\frac{y \cdot \left(z - t\right)}{a} \le 1.947585218067920906632704854880810808419 \cdot 10^{229}\right):\\ \;\;\;\;x - \frac{y}{a} \cdot \left(z - t\right)\\ \mathbf{else}:\\ \;\;\;\;x - \frac{y \cdot \left(z - t\right)}{a}\\ \end{array}\]

Reproduce

herbie shell --seed 2019325 
(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 (/ (/ a (- z t)) y))) (if (< y 2.894426862792089e-49) (- x (/ (* y (- z t)) a)) (- x (/ y (/ a (- z t))))))

  (- x (/ (* y (- z t)) a)))