Average Error: 0.0 → 0.0
Time: 14.3s
Precision: 64
\[56789 \le a \le 98765 \land 0.0 \le b \le 1 \land 0.0 \le c \le 0.001677300000000000058247850986958837893326 \land 0.0 \le d \le 0.001677300000000000058247850986958837893326\]
\[a \cdot \left(\left(b + c\right) + d\right)\]
↓
\[\mathsf{fma}\left(b + c, a, d \cdot a\right)\]
a \cdot \left(\left(b + c\right) + d\right)
↓
\mathsf{fma}\left(b + c, a, d \cdot a\right)
double f(double a, double b, double c, double d) {
        double r98762 = a;
        double r98763 = b;
        double r98764 = c;
        double r98765 = r98763 + r98764;
        double r98766 = d;
        double r98767 = r98765 + r98766;
        double r98768 = r98762 * r98767;
        return r98768;
}
↓
double f(double a, double b, double c, double d) {
        double r98769 = b;
        double r98770 = c;
        double r98771 = r98769 + r98770;
        double r98772 = a;
        double r98773 = d;
        double r98774 = r98773 * r98772;
        double r98775 = fma(r98771, r98772, r98774);
        return r98775;
}

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 distribute-lft-in0.0

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

  4. Simplified0.0

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

  5. Simplified0.0

    \[\leadsto \left(b + c\right) \cdot a + \color{blue}{d \cdot a}\]
  6. Using strategy rm

  7. Applied fma-def0.0

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

  8. Final simplification0.0

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

Reproduce

herbie shell --seed 2019303 +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.0016773000000000001) (<= 0.0 d 0.0016773000000000001))

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

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