Average Error: 0.0 → 0.0
Time: 2.2s
Precision: 64
\[56789 \le a \le 98765 \land 0.0 \le b \le 1 \land 0.0 \le c \le 0.0016773000000000001 \land 0.0 \le d \le 0.0016773000000000001\]
\[a \cdot \left(\left(b + c\right) + d\right)\]
↓
\[\mathsf{fma}\left(d, a, \mathsf{fma}\left(a, b, a \cdot c\right)\right)\]
a \cdot \left(\left(b + c\right) + d\right)
↓
\mathsf{fma}\left(d, a, \mathsf{fma}\left(a, b, a \cdot c\right)\right)
double f(double a, double b, double c, double d) {
        double r92922 = a;
        double r92923 = b;
        double r92924 = c;
        double r92925 = r92923 + r92924;
        double r92926 = d;
        double r92927 = r92925 + r92926;
        double r92928 = r92922 * r92927;
        return r92928;
}
↓
double f(double a, double b, double c, double d) {
        double r92929 = d;
        double r92930 = a;
        double r92931 = b;
        double r92932 = c;
        double r92933 = r92930 * r92932;
        double r92934 = fma(r92930, r92931, r92933);
        double r92935 = fma(r92929, r92930, r92934);
        return r92935;
}

Error

Bits error versus a

Bits error versus b

Bits error versus c

Bits error versus d

Target

Original0.0
Target0.0
Herbie0.0
\[a \cdot b + a \cdot \left(c + d\right)\]

Derivation

  1. Initial program 0.0

    \[a \cdot \left(\left(b + c\right) + d\right)\]
  2. Using strategy rm

  3. Applied pow10.0

    \[\leadsto a \cdot \color{blue}{{\left(\left(b + c\right) + d\right)}^{1}}\]

  4. Applied pow10.0

    \[\leadsto \color{blue}{{a}^{1}} \cdot {\left(\left(b + c\right) + d\right)}^{1}\]

  5. Applied pow-prod-down0.0

    \[\leadsto \color{blue}{{\left(a \cdot \left(\left(b + c\right) + d\right)\right)}^{1}}\]

  6. Simplified0.0

    \[\leadsto {\color{blue}{\left(\mathsf{fma}\left(d, a, \mathsf{fma}\left(a, b, a \cdot c\right)\right)\right)}}^{1}\]

  7. Final simplification0.0

    \[\leadsto \mathsf{fma}\left(d, a, \mathsf{fma}\left(a, b, a \cdot c\right)\right)\]

Reproduce

herbie shell --seed 2020035 +o rules:numerics
(FPCore (a b c d)
  :name "Expression, p14"
  :precision binary64
  :pre (and (<= 56789 a 98765) (<= 0.0 b 1) (<= 0.0 c 0.0016773) (<= 0.0 d 0.0016773))

  :herbie-target
  (+ (* a b) (* a (+ c d)))

  (* a (+ (+ b c) d)))