Diagrams.Solve.Polynomial:cubForm from diagrams-solve-0.1, C

Average Error: 0.3 → 0.2

Time: 11.9s

Precision: 64

• Falling back on racket egraph because egg-herbie package not installed

Math

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

↓

\[\frac{\frac{x}{3}}{y} \cdot \sqrt{1}\]

C

double f(double x, double y) {
        double r804050 = x;
        double r804051 = y;
        double r804052 = 3.0;
        double r804053 = r804051 * r804052;
        double r804054 = r804050 / r804053;
        return r804054;
}

↓

double f(double x, double y) {
        double r804055 = x;
        double r804056 = 3.0;
        double r804057 = r804055 / r804056;
        double r804058 = y;
        double r804059 = r804057 / r804058;
        double r804060 = 1.0;
        double r804061 = sqrt(r804060);
        double r804062 = r804059 * r804061;
        return r804062;
}

Error

Bits error versus x

Bits error versus y

Target

Original	0.3
Target	0.3
Herbie	0.2

\[\frac{\frac{x}{y}}{3}\]

Derivation

1. Initial program 0.3

   \[\frac{x}{y \cdot 3}\]
2. Using strategy rm

3. Applied *-un-lft-identity 0.3

   \[\leadsto \frac{\color{blue}{1 \cdot x}}{y \cdot 3}\]

4. Applied times-frac 0.3

   \[\leadsto \color{blue}{\frac{1}{y} \cdot \frac{x}{3}}\]
5. Using strategy rm

6. Applied *-un-lft-identity 0.3

   \[\leadsto \frac{1}{\color{blue}{1 \cdot y}} \cdot \frac{x}{3}\]

7. Applied add-sqr-sqrt 0.3

   \[\leadsto \frac{\color{blue}{\sqrt{1} \cdot \sqrt{1}}}{1 \cdot y} \cdot \frac{x}{3}\]

8. Applied times-frac 0.3

   \[\leadsto \color{blue}{\left(\frac{\sqrt{1}}{1} \cdot \frac{\sqrt{1}}{y}\right)} \cdot \frac{x}{3}\]

9. Applied associate-*l* 0.3

   \[\leadsto \color{blue}{\frac{\sqrt{1}}{1} \cdot \left(\frac{\sqrt{1}}{y} \cdot \frac{x}{3}\right)}\]

10. Simplified 0.2

    \[\leadsto \frac{\sqrt{1}}{1} \cdot \color{blue}{\frac{\frac{x}{3}}{y}}\]

11. Final simplification 0.2

    \[\leadsto \frac{\frac{x}{3}}{y} \cdot \sqrt{1}\]

Reproduce

herbie shell --seed 2020042 
(FPCore (x y)
  :name "Diagrams.Solve.Polynomial:cubForm  from diagrams-solve-0.1, C"
  :precision binary64

  :herbie-target
  (/ (/ x y) 3)

  (/ x (* y 3)))