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

Average Error: 0.1 → 0.1

Time: 12.2s

Precision: 64

Math

\[x - \left(y \cdot 4\right) \cdot z\]

↓

\[x - \left(4 \cdot y\right) \cdot z\]

C

double f(double x, double y, double z) {
        double r11562813 = x;
        double r11562814 = y;
        double r11562815 = 4.0;
        double r11562816 = r11562814 * r11562815;
        double r11562817 = z;
        double r11562818 = r11562816 * r11562817;
        double r11562819 = r11562813 - r11562818;
        return r11562819;
}

↓

double f(double x, double y, double z) {
        double r11562820 = x;
        double r11562821 = 4.0;
        double r11562822 = y;
        double r11562823 = r11562821 * r11562822;
        double r11562824 = z;
        double r11562825 = r11562823 * r11562824;
        double r11562826 = r11562820 - r11562825;
        return r11562826;
}

Error

Bits error versus x

Bits error versus y

Bits error versus z

Derivation

1. Initial program 0.1

   \[x - \left(y \cdot 4\right) \cdot z\]

2. Final simplification 0.1

   \[\leadsto x - \left(4 \cdot y\right) \cdot z\]

Reproduce

herbie shell --seed 2019170 
(FPCore (x y z)
  :name "Diagrams.Solve.Polynomial:quadForm from diagrams-solve-0.1, A"
  (- x (* (* y 4.0) z)))