Average Error: 28.8 → 0.1
Time: 4.0s
Precision: 64
\[\frac{\left(x \cdot x + y \cdot y\right) - z \cdot z}{y \cdot 2}\]
↓
\[0.5 \cdot \left(\left(y + x \cdot \frac{x}{y}\right) - \left|z\right| \cdot \frac{\left|z\right|}{y}\right)\]
\frac{\left(x \cdot x + y \cdot y\right) - z \cdot z}{y \cdot 2}
↓
0.5 \cdot \left(\left(y + x \cdot \frac{x}{y}\right) - \left|z\right| \cdot \frac{\left|z\right|}{y}\right)
double f(double x, double y, double z) {
        double r760888 = x;
        double r760889 = r760888 * r760888;
        double r760890 = y;
        double r760891 = r760890 * r760890;
        double r760892 = r760889 + r760891;
        double r760893 = z;
        double r760894 = r760893 * r760893;
        double r760895 = r760892 - r760894;
        double r760896 = 2.0;
        double r760897 = r760890 * r760896;
        double r760898 = r760895 / r760897;
        return r760898;
}
↓
double f(double x, double y, double z) {
        double r760899 = 0.5;
        double r760900 = y;
        double r760901 = x;
        double r760902 = r760901 / r760900;
        double r760903 = r760901 * r760902;
        double r760904 = r760900 + r760903;
        double r760905 = z;
        double r760906 = fabs(r760905);
        double r760907 = r760906 / r760900;
        double r760908 = r760906 * r760907;
        double r760909 = r760904 - r760908;
        double r760910 = r760899 * r760909;
        return r760910;
}

Error

Bits error versus x

Bits error versus y

Bits error versus z

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Target

Original28.8
Target0.1
Herbie0.1
\[y \cdot 0.5 - \left(\frac{0.5}{y} \cdot \left(z + x\right)\right) \cdot \left(z - x\right)\]

Derivation

  1. Initial program 28.8

    \[\frac{\left(x \cdot x + y \cdot y\right) - z \cdot z}{y \cdot 2}\]

  2. Taylor expanded around 0 12.4

    \[\leadsto \color{blue}{\left(0.5 \cdot y + 0.5 \cdot \frac{{x}^{2}}{y}\right) - 0.5 \cdot \frac{{z}^{2}}{y}}\]

  3. Simplified12.4

    \[\leadsto \color{blue}{0.5 \cdot \left(\left(y + \frac{{x}^{2}}{y}\right) - \frac{{z}^{2}}{y}\right)}\]
  4. Using strategy rm

  5. Applied *-un-lft-identity12.4

    \[\leadsto 0.5 \cdot \left(\left(y + \frac{{x}^{2}}{\color{blue}{1 \cdot y}}\right) - \frac{{z}^{2}}{y}\right)\]

  6. Applied add-sqr-sqrt38.5

    \[\leadsto 0.5 \cdot \left(\left(y + \frac{{\color{blue}{\left(\sqrt{x} \cdot \sqrt{x}\right)}}^{2}}{1 \cdot y}\right) - \frac{{z}^{2}}{y}\right)\]

  7. Applied unpow-prod-down38.5

    \[\leadsto 0.5 \cdot \left(\left(y + \frac{\color{blue}{{\left(\sqrt{x}\right)}^{2} \cdot {\left(\sqrt{x}\right)}^{2}}}{1 \cdot y}\right) - \frac{{z}^{2}}{y}\right)\]

  8. Applied times-frac35.7

    \[\leadsto 0.5 \cdot \left(\left(y + \color{blue}{\frac{{\left(\sqrt{x}\right)}^{2}}{1} \cdot \frac{{\left(\sqrt{x}\right)}^{2}}{y}}\right) - \frac{{z}^{2}}{y}\right)\]

  9. Simplified35.7

    \[\leadsto 0.5 \cdot \left(\left(y + \color{blue}{x} \cdot \frac{{\left(\sqrt{x}\right)}^{2}}{y}\right) - \frac{{z}^{2}}{y}\right)\]

  10. Simplified6.9

    \[\leadsto 0.5 \cdot \left(\left(y + x \cdot \color{blue}{\frac{x}{y}}\right) - \frac{{z}^{2}}{y}\right)\]
  11. Using strategy rm

  12. Applied *-un-lft-identity6.9

    \[\leadsto 0.5 \cdot \left(\left(y + x \cdot \frac{x}{y}\right) - \frac{{z}^{2}}{\color{blue}{1 \cdot y}}\right)\]

  13. Applied add-sqr-sqrt6.9

    \[\leadsto 0.5 \cdot \left(\left(y + x \cdot \frac{x}{y}\right) - \frac{\color{blue}{\sqrt{{z}^{2}} \cdot \sqrt{{z}^{2}}}}{1 \cdot y}\right)\]

  14. Applied times-frac6.9

    \[\leadsto 0.5 \cdot \left(\left(y + x \cdot \frac{x}{y}\right) - \color{blue}{\frac{\sqrt{{z}^{2}}}{1} \cdot \frac{\sqrt{{z}^{2}}}{y}}\right)\]

  15. Simplified6.9

    \[\leadsto 0.5 \cdot \left(\left(y + x \cdot \frac{x}{y}\right) - \color{blue}{\left|z\right|} \cdot \frac{\sqrt{{z}^{2}}}{y}\right)\]

  16. Simplified0.1

    \[\leadsto 0.5 \cdot \left(\left(y + x \cdot \frac{x}{y}\right) - \left|z\right| \cdot \color{blue}{\frac{\left|z\right|}{y}}\right)\]

  17. Final simplification0.1

    \[\leadsto 0.5 \cdot \left(\left(y + x \cdot \frac{x}{y}\right) - \left|z\right| \cdot \frac{\left|z\right|}{y}\right)\]

Reproduce

herbie shell --seed 2020018 +o rules:numerics
(FPCore (x y z)
  :name "Diagrams.TwoD.Apollonian:initialConfig from diagrams-contrib-1.3.0.5, A"
  :precision binary64

  :herbie-target
  (- (* y 0.5) (* (* (/ 0.5 y) (+ z x)) (- z x)))

  (/ (- (+ (* x x) (* y y)) (* z z)) (* y 2)))