Average Error: 3.6 → 2.8
Time: 13.5s
Precision: 64
\[-14 \le a \le -13 \land -3 \le b \le -2 \land 3 \le c \le 3.5 \land 12.5 \le d \le 13.5\]
\[\left(a + \left(b + \left(c + d\right)\right)\right) \cdot 2\]
↓
\[\mathsf{log1p}\left(\mathsf{expm1}\left(a + \left(\left(b + c\right) + d\right)\right)\right) \cdot 2\]
\left(a + \left(b + \left(c + d\right)\right)\right) \cdot 2
↓
\mathsf{log1p}\left(\mathsf{expm1}\left(a + \left(\left(b + c\right) + d\right)\right)\right) \cdot 2
double f(double a, double b, double c, double d) {
        double r64384 = a;
        double r64385 = b;
        double r64386 = c;
        double r64387 = d;
        double r64388 = r64386 + r64387;
        double r64389 = r64385 + r64388;
        double r64390 = r64384 + r64389;
        double r64391 = 2.0;
        double r64392 = r64390 * r64391;
        return r64392;
}
↓
double f(double a, double b, double c, double d) {
        double r64393 = a;
        double r64394 = b;
        double r64395 = c;
        double r64396 = r64394 + r64395;
        double r64397 = d;
        double r64398 = r64396 + r64397;
        double r64399 = r64393 + r64398;
        double r64400 = expm1(r64399);
        double r64401 = log1p(r64400);
        double r64402 = 2.0;
        double r64403 = r64401 * r64402;
        return r64403;
}

Error

Bits error versus a

Bits error versus b

Bits error versus c

Bits error versus d

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Target

Original3.6
Target3.8
Herbie2.8
\[\left(a + b\right) \cdot 2 + \left(c + d\right) \cdot 2\]

Derivation

  1. Initial program 3.6

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

  3. Applied associate-+r+2.8

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

  5. Applied log1p-expm1-u2.8

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

  6. Final simplification2.8

    \[\leadsto \mathsf{log1p}\left(\mathsf{expm1}\left(a + \left(\left(b + c\right) + d\right)\right)\right) \cdot 2\]

Reproduce

herbie shell --seed 2019303 +o rules:numerics
(FPCore (a b c d)
  :name "Expression, p6"
  :precision binary64
  :pre (and (<= -14 a -13) (<= -3 b -2) (<= 3 c 3.5) (<= 12.5 d 13.5))

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

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