Diagrams.Solve.Polynomial:quartForm from diagrams-solve-0.1, B

Average Error: 0.0 → 0.0

Time: 4.4s

Precision: 64

Math

\[\left(\frac{1}{8} \cdot x - \frac{y \cdot z}{2}\right) + t\]

↓

\[t + \left(\frac{1}{8} \cdot x - \frac{y \cdot z}{2}\right)\]

C

double f(double x, double y, double z, double t) {
        double r31890596 = 1.0;
        double r31890597 = 8.0;
        double r31890598 = r31890596 / r31890597;
        double r31890599 = x;
        double r31890600 = r31890598 * r31890599;
        double r31890601 = y;
        double r31890602 = z;
        double r31890603 = r31890601 * r31890602;
        double r31890604 = 2.0;
        double r31890605 = r31890603 / r31890604;
        double r31890606 = r31890600 - r31890605;
        double r31890607 = t;
        double r31890608 = r31890606 + r31890607;
        return r31890608;
}

↓

double f(double x, double y, double z, double t) {
        double r31890609 = t;
        double r31890610 = 1.0;
        double r31890611 = 8.0;
        double r31890612 = r31890610 / r31890611;
        double r31890613 = x;
        double r31890614 = r31890612 * r31890613;
        double r31890615 = y;
        double r31890616 = z;
        double r31890617 = r31890615 * r31890616;
        double r31890618 = 2.0;
        double r31890619 = r31890617 / r31890618;
        double r31890620 = r31890614 - r31890619;
        double r31890621 = r31890609 + r31890620;
        return r31890621;
}

Error

Bits error versus x

Bits error versus y

Bits error versus z

Bits error versus t

Target

Original	0.0
Target	0.0
Herbie	0.0

\[\left(\frac{x}{8} + t\right) - \frac{z}{2} \cdot y\]

Derivation

1. Initial program 0.0

   \[\left(\frac{1}{8} \cdot x - \frac{y \cdot z}{2}\right) + t\]

2. Final simplification 0.0

   \[\leadsto t + \left(\frac{1}{8} \cdot x - \frac{y \cdot z}{2}\right)\]

Reproduce

herbie shell --seed 2019192 
(FPCore (x y z t)
  :name "Diagrams.Solve.Polynomial:quartForm  from diagrams-solve-0.1, B"

  :herbie-target
  (- (+ (/ x 8.0) t) (* (/ z 2.0) y))

  (+ (- (* (/ 1.0 8.0) x) (/ (* y z) 2.0)) t))