Numeric.Integration.TanhSinh:everywhere from integration-0.2.1

Average Error: 5.1 → 5.1

Time: 14.0s

Precision: 64

Math

\[x \cdot \left(1 + y \cdot y\right)\]

↓

\[\mathsf{fma}\left(y, y, 1\right) \cdot x\]

C

double f(double x, double y) {
        double r18952332 = x;
        double r18952333 = 1.0;
        double r18952334 = y;
        double r18952335 = r18952334 * r18952334;
        double r18952336 = r18952333 + r18952335;
        double r18952337 = r18952332 * r18952336;
        return r18952337;
}

↓

double f(double x, double y) {
        double r18952338 = y;
        double r18952339 = 1.0;
        double r18952340 = fma(r18952338, r18952338, r18952339);
        double r18952341 = x;
        double r18952342 = r18952340 * r18952341;
        return r18952342;
}

Error

Bits error versus x

Bits error versus y

Target

Original	5.1
Target	0.1
Herbie	5.1

\[x + \left(x \cdot y\right) \cdot y\]

Derivation

1. Initial program 5.1

   \[x \cdot \left(1 + y \cdot y\right)\]

2. Simplified 5.1

   \[\leadsto \color{blue}{\mathsf{fma}\left(y, y, 1\right) \cdot x}\]

3. Final simplification 5.1

   \[\leadsto \mathsf{fma}\left(y, y, 1\right) \cdot x\]

Reproduce

herbie shell --seed 2019172 +o rules:numerics
(FPCore (x y)
  :name "Numeric.Integration.TanhSinh:everywhere from integration-0.2.1"

  :herbie-target
  (+ x (* (* x y) y))

  (* x (+ 1.0 (* y y))))