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

Average Error: 0.0 → 0.0

Time: 19.7s

Precision: 64

Math

\[\left(\frac{1.0}{8.0} \cdot x - \frac{y \cdot z}{2.0}\right) + t\]

↓

\[t + \left(\frac{1.0}{8.0} \cdot x - \frac{y \cdot z}{2.0}\right)\]

C

double f(double x, double y, double z, double t) {
        double r36455513 = 1.0;
        double r36455514 = 8.0;
        double r36455515 = r36455513 / r36455514;
        double r36455516 = x;
        double r36455517 = r36455515 * r36455516;
        double r36455518 = y;
        double r36455519 = z;
        double r36455520 = r36455518 * r36455519;
        double r36455521 = 2.0;
        double r36455522 = r36455520 / r36455521;
        double r36455523 = r36455517 - r36455522;
        double r36455524 = t;
        double r36455525 = r36455523 + r36455524;
        return r36455525;
}

↓

double f(double x, double y, double z, double t) {
        double r36455526 = t;
        double r36455527 = 1.0;
        double r36455528 = 8.0;
        double r36455529 = r36455527 / r36455528;
        double r36455530 = x;
        double r36455531 = r36455529 * r36455530;
        double r36455532 = y;
        double r36455533 = z;
        double r36455534 = r36455532 * r36455533;
        double r36455535 = 2.0;
        double r36455536 = r36455534 / r36455535;
        double r36455537 = r36455531 - r36455536;
        double r36455538 = r36455526 + r36455537;
        return r36455538;
}

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.0} + t\right) - \frac{z}{2.0} \cdot y\]

Derivation

1. Initial program 0.0

   \[\left(\frac{1.0}{8.0} \cdot x - \frac{y \cdot z}{2.0}\right) + t\]

2. Final simplification 0.0

   \[\leadsto t + \left(\frac{1.0}{8.0} \cdot x - \frac{y \cdot z}{2.0}\right)\]

Reproduce

herbie shell --seed 2019162 
(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))