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

Average Error: 0.2 → 0.2

Time: 1.3s

Precision: 64

Math

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

↓

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

C

double f(double x, double y) {
        double r653617 = x;
        double r653618 = y;
        double r653619 = 3.0;
        double r653620 = r653618 * r653619;
        double r653621 = r653617 / r653620;
        return r653621;
}

↓

double f(double x, double y) {
        double r653622 = x;
        double r653623 = 3.0;
        double r653624 = r653622 / r653623;
        double r653625 = y;
        double r653626 = r653624 / r653625;
        double r653627 = 1.0;
        double r653628 = sqrt(r653627);
        double r653629 = r653626 * r653628;
        return r653629;
}

Error

Bits error versus x

Bits error versus y

Target

Original	0.2
Target	0.3
Herbie	0.2

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

Derivation

1. Initial program 0.2

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

3. Applied *-un-lft-identity 0.2

   \[\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 2020024 
(FPCore (x y)
  :name "Diagrams.Solve.Polynomial:cubForm  from diagrams-solve-0.1, C"
  :precision binary64

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

  (/ x (* y 3)))