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

Average Error: 0.1 → 0.1

Time: 5.1s

Precision: 64

Math

\[x - \frac{3}{8} \cdot y\]

↓

\[x - \frac{3}{8} \cdot y\]

C

double f(double x, double y) {
        double r225737 = x;
        double r225738 = 3.0;
        double r225739 = 8.0;
        double r225740 = r225738 / r225739;
        double r225741 = y;
        double r225742 = r225740 * r225741;
        double r225743 = r225737 - r225742;
        return r225743;
}

↓

double f(double x, double y) {
        double r225744 = x;
        double r225745 = 3.0;
        double r225746 = 8.0;
        double r225747 = r225745 / r225746;
        double r225748 = y;
        double r225749 = r225747 * r225748;
        double r225750 = r225744 - r225749;
        return r225750;
}

Error

Bits error versus x

Bits error versus y

Derivation

1. Initial program 0.1

   \[x - \frac{3}{8} \cdot y\]

2. Final simplification 0.1

   \[\leadsto x - \frac{3}{8} \cdot y\]

Reproduce

herbie shell --seed 2020021 +o rules:numerics
(FPCore (x y)
  :name "Diagrams.Solve.Polynomial:quartForm  from diagrams-solve-0.1, A"
  :precision binary64
  (- x (* (/ 3 8) y)))