Henrywood and Agarwal, Equation (13)

Average Error: 59.3 → 33.8

Time: 28.5s

Precision: 64

• could not determine a ground truth for program body

  1. c0 = -1.8426955311883228e+230
  2. w = 1.4316726427427919e-198
  3. h = 2.0944197576728813e-33
  4. D = -1.263279512918104e+22
  5. d = 5.382419207738136e+203
  6. M = 1.5330596284314128e-114

Math

\[\frac{c0}{2 \cdot w} \cdot \left(\frac{c0 \cdot \left(d \cdot d\right)}{\left(w \cdot h\right) \cdot \left(D \cdot D\right)} + \sqrt{\frac{c0 \cdot \left(d \cdot d\right)}{\left(w \cdot h\right) \cdot \left(D \cdot D\right)} \cdot \frac{c0 \cdot \left(d \cdot d\right)}{\left(w \cdot h\right) \cdot \left(D \cdot D\right)} - M \cdot M}\right)\]

↓

\[0\]

C

double f(double c0, double w, double h, double D, double d, double M) {
        double r145439 = c0;
        double r145440 = 2.0;
        double r145441 = w;
        double r145442 = r145440 * r145441;
        double r145443 = r145439 / r145442;
        double r145444 = d;
        double r145445 = r145444 * r145444;
        double r145446 = r145439 * r145445;
        double r145447 = h;
        double r145448 = r145441 * r145447;
        double r145449 = D;
        double r145450 = r145449 * r145449;
        double r145451 = r145448 * r145450;
        double r145452 = r145446 / r145451;
        double r145453 = r145452 * r145452;
        double r145454 = M;
        double r145455 = r145454 * r145454;
        double r145456 = r145453 - r145455;
        double r145457 = sqrt(r145456);
        double r145458 = r145452 + r145457;
        double r145459 = r145443 * r145458;
        return r145459;
}

↓

double f(double __attribute__((unused)) c0, double __attribute__((unused)) w, double __attribute__((unused)) h, double __attribute__((unused)) D, double __attribute__((unused)) d, double __attribute__((unused)) M) {
        double r145460 = 0.0;
        return r145460;
}

Error

Bits error versus c0

Bits error versus w

Bits error versus h

Bits error versus D

Bits error versus d

Bits error versus M

Derivation

1. Initial program 59.3

   \[\frac{c0}{2 \cdot w} \cdot \left(\frac{c0 \cdot \left(d \cdot d\right)}{\left(w \cdot h\right) \cdot \left(D \cdot D\right)} + \sqrt{\frac{c0 \cdot \left(d \cdot d\right)}{\left(w \cdot h\right) \cdot \left(D \cdot D\right)} \cdot \frac{c0 \cdot \left(d \cdot d\right)}{\left(w \cdot h\right) \cdot \left(D \cdot D\right)} - M \cdot M}\right)\]

2. Taylor expanded around inf 35.7

   \[\leadsto \frac{c0}{2 \cdot w} \cdot \color{blue}{0}\]
3. Using strategy rm

4. Applied add-cbrt-cube 35.7

   \[\leadsto \frac{c0}{2 \cdot w} \cdot \color{blue}{\sqrt[3]{\left(0 \cdot 0\right) \cdot 0}}\]

5. Applied add-cbrt-cube 42.3

   \[\leadsto \frac{c0}{2 \cdot \color{blue}{\sqrt[3]{\left(w \cdot w\right) \cdot w}}} \cdot \sqrt[3]{\left(0 \cdot 0\right) \cdot 0}\]

6. Applied add-cbrt-cube 42.3

   \[\leadsto \frac{c0}{\color{blue}{\sqrt[3]{\left(2 \cdot 2\right) \cdot 2}} \cdot \sqrt[3]{\left(w \cdot w\right) \cdot w}} \cdot \sqrt[3]{\left(0 \cdot 0\right) \cdot 0}\]

7. Applied cbrt-unprod 42.3

   \[\leadsto \frac{c0}{\color{blue}{\sqrt[3]{\left(\left(2 \cdot 2\right) \cdot 2\right) \cdot \left(\left(w \cdot w\right) \cdot w\right)}}} \cdot \sqrt[3]{\left(0 \cdot 0\right) \cdot 0}\]

8. Applied add-cbrt-cube 48.7

   \[\leadsto \frac{\color{blue}{\sqrt[3]{\left(c0 \cdot c0\right) \cdot c0}}}{\sqrt[3]{\left(\left(2 \cdot 2\right) \cdot 2\right) \cdot \left(\left(w \cdot w\right) \cdot w\right)}} \cdot \sqrt[3]{\left(0 \cdot 0\right) \cdot 0}\]

9. Applied cbrt-undiv 49.1

   \[\leadsto \color{blue}{\sqrt[3]{\frac{\left(c0 \cdot c0\right) \cdot c0}{\left(\left(2 \cdot 2\right) \cdot 2\right) \cdot \left(\left(w \cdot w\right) \cdot w\right)}}} \cdot \sqrt[3]{\left(0 \cdot 0\right) \cdot 0}\]

10. Applied cbrt-unprod 49.1

    \[\leadsto \color{blue}{\sqrt[3]{\frac{\left(c0 \cdot c0\right) \cdot c0}{\left(\left(2 \cdot 2\right) \cdot 2\right) \cdot \left(\left(w \cdot w\right) \cdot w\right)} \cdot \left(\left(0 \cdot 0\right) \cdot 0\right)}}\]

11. Simplified 33.8

    \[\leadsto \sqrt[3]{\color{blue}{0}}\]

12. Final simplification 33.8

    \[\leadsto 0\]

Reproduce

herbie shell --seed 2019326 
(FPCore (c0 w h D d M)
  :name "Henrywood and Agarwal, Equation (13)"
  :precision binary64
  (* (/ c0 (* 2 w)) (+ (/ (* c0 (* d d)) (* (* w h) (* D D))) (sqrt (- (* (/ (* c0 (* d d)) (* (* w h) (* D D))) (/ (* c0 (* d d)) (* (* w h) (* D D)))) (* M M))))))