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

Average Error: 0.0 → 0.0

Time: 7.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 r33774481 = 1.0;
        double r33774482 = 8.0;
        double r33774483 = r33774481 / r33774482;
        double r33774484 = x;
        double r33774485 = r33774483 * r33774484;
        double r33774486 = y;
        double r33774487 = z;
        double r33774488 = r33774486 * r33774487;
        double r33774489 = 2.0;
        double r33774490 = r33774488 / r33774489;
        double r33774491 = r33774485 - r33774490;
        double r33774492 = t;
        double r33774493 = r33774491 + r33774492;
        return r33774493;
}

↓

double f(double x, double y, double z, double t) {
        double r33774494 = t;
        double r33774495 = 1.0;
        double r33774496 = 8.0;
        double r33774497 = r33774495 / r33774496;
        double r33774498 = x;
        double r33774499 = r33774497 * r33774498;
        double r33774500 = y;
        double r33774501 = z;
        double r33774502 = r33774500 * r33774501;
        double r33774503 = 2.0;
        double r33774504 = r33774502 / r33774503;
        double r33774505 = r33774499 - r33774504;
        double r33774506 = r33774494 + r33774505;
        return r33774506;
}

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 2019169 
(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))