Diagrams.TwoD.Path.Metafont.Internal:hobbyF from diagrams-contrib-1.3.0.5

Average Error: 0.5 → 0.4

Time: 13.1s

Precision: binary64

Cost: 124416

Math

\[\frac{2 + \left(\left(\sqrt{2} \cdot \left(\sin x - \frac{\sin y}{16}\right)\right) \cdot \left(\sin y - \frac{\sin x}{16}\right)\right) \cdot \left(\cos x - \cos y\right)}{3 \cdot \left(\left(1 + \frac{\sqrt{5} - 1}{2} \cdot \cos x\right) + \frac{3 - \sqrt{5}}{2} \cdot \cos y\right)}\]

↓

\[\frac{2 + \left(\cos x - \cos y\right) \cdot \left(\left(\sin y - \frac{\sin x}{16}\right) \cdot \left(\sin x \cdot \sqrt{2} - \sqrt{2} \cdot {\left(\frac{\sqrt[3]{\sin y}}{\sqrt[3]{\sqrt[3]{16} \cdot \sqrt[3]{16}} \cdot \sqrt[3]{\sqrt[3]{16}}}\right)}^{3}\right)\right)}{3 + \left(\cos y \cdot \left(1.5 \cdot \left(3 - \sqrt{5}\right)\right) + 1.5 \cdot \left(\cos x \cdot \left(\sqrt{5} - 1\right)\right)\right)}\]

FPCore

(FPCore (x y)
 :precision binary64
 (/
  (+
   2.0
   (*
    (*
     (* (sqrt 2.0) (- (sin x) (/ (sin y) 16.0)))
     (- (sin y) (/ (sin x) 16.0)))
    (- (cos x) (cos y))))
  (*
   3.0
   (+
    (+ 1.0 (* (/ (- (sqrt 5.0) 1.0) 2.0) (cos x)))
    (* (/ (- 3.0 (sqrt 5.0)) 2.0) (cos y))))))

↓

(FPCore (x y)
 :precision binary64
 (/
  (+
   2.0
   (*
    (- (cos x) (cos y))
    (*
     (- (sin y) (/ (sin x) 16.0))
     (-
      (* (sin x) (sqrt 2.0))
      (*
       (sqrt 2.0)
       (pow
        (/
         (cbrt (sin y))
         (* (cbrt (* (cbrt 16.0) (cbrt 16.0))) (cbrt (cbrt 16.0))))
        3.0))))))
  (+
   3.0
   (+
    (* (cos y) (* 1.5 (- 3.0 (sqrt 5.0))))
    (* 1.5 (* (cos x) (- (sqrt 5.0) 1.0)))))))

C

double code(double x, double y) {
	return (2.0 + (((sqrt(2.0) * (sin(x) - (sin(y) / 16.0))) * (sin(y) - (sin(x) / 16.0))) * (cos(x) - cos(y)))) / (3.0 * ((1.0 + (((sqrt(5.0) - 1.0) / 2.0) * cos(x))) + (((3.0 - sqrt(5.0)) / 2.0) * cos(y))));
}

↓

double code(double x, double y) {
	return (2.0 + ((cos(x) - cos(y)) * ((sin(y) - (sin(x) / 16.0)) * ((sin(x) * sqrt(2.0)) - (sqrt(2.0) * pow((cbrt(sin(y)) / (cbrt(cbrt(16.0) * cbrt(16.0)) * cbrt(cbrt(16.0)))), 3.0)))))) / (3.0 + ((cos(y) * (1.5 * (3.0 - sqrt(5.0)))) + (1.5 * (cos(x) * (sqrt(5.0) - 1.0)))));
}

Error

Bits error versus x

Bits error versus y

Alternatives

Alternative 1
Error	0.4
Cost	72768

\[\frac{2 + \left(\cos x - \cos y\right) \cdot \left(\left(\sin y - \frac{\sin x}{16}\right) \cdot \left(\sqrt{2} \cdot \left(\sin x - \frac{\sin y}{16}\right)\right)\right)}{3 + \left(\cos y \cdot \left(1.5 \cdot \left(3 - \sqrt{5}\right)\right) + \frac{\cos x \cdot 6}{\sqrt{5} + 1}\right)}\]

Alternative 2
Error	0.4
Cost	72768

\[\frac{2 + \left(\cos x - \cos y\right) \cdot \left(\left(\sin y - \frac{\sin x}{16}\right) \cdot \left(\sqrt{2} \cdot \left(\sin x - \frac{\sin y}{16}\right)\right)\right)}{3 + \left(\cos y \cdot \left(1.5 \cdot \left(3 - \sqrt{5}\right)\right) + \left(\sqrt{5} - 1\right) \cdot \left(\cos x \cdot 1.5\right)\right)}\]

Alternative 3
Error	0.4
Cost	72768

\[\frac{2 + \left(\cos x - \cos y\right) \cdot \left(\left(\sin y - \frac{\sin x}{16}\right) \cdot \left(\sqrt{2} \cdot \left(\sin x - \frac{\sin y}{16}\right)\right)\right)}{3 + \left(\cos y \cdot \left(1.5 \cdot \left(3 - \sqrt{5}\right)\right) + 1.5 \cdot \left(\cos x \cdot \left(\sqrt{5} - 1\right)\right)\right)}\]

Alternative 4
Error	11.6
Cost	66696

\[\begin{array}{l} \mathbf{if}\;x \leq -0.028121830058273272 \lor \neg \left(x \leq 0.023668285587148744\right):\\ \;\;\;\;\frac{2 + \left(\cos x - \cos y\right) \cdot \left(\left(\sin y - \frac{\sin x}{16}\right) \cdot \left(\sin x \cdot \sqrt{2}\right)\right)}{3 + \left(\cos y \cdot \left(1.5 \cdot \left(3 - \sqrt{5}\right)\right) + 1.5 \cdot \left(\cos x \cdot \left(\sqrt{5} - 1\right)\right)\right)}\\ \mathbf{else}:\\ \;\;\;\;\frac{2 + \left(\cos x - \cos y\right) \cdot \left(\left(\sin y - \frac{\sin x}{16}\right) \cdot \left(\sqrt{2} \cdot \left(x - \sin y \cdot 0.0625\right)\right)\right)}{3 + \left(\cos y \cdot \left(1.5 \cdot \left(3 - \sqrt{5}\right)\right) + 1.5 \cdot \left(\cos x \cdot \left(\sqrt{5} - 1\right)\right)\right)}\\ \end{array}\]

Alternative 5
Error	11.7
Cost	66568

\[\begin{array}{l} \mathbf{if}\;x \leq -0.0001083584555664739 \lor \neg \left(x \leq 1.7349909248656208 \cdot 10^{-05}\right):\\ \;\;\;\;\frac{2 + \left(\cos x - \cos y\right) \cdot \left(\left(\sin y - \frac{\sin x}{16}\right) \cdot \left(\sin x \cdot \sqrt{2}\right)\right)}{3 + \left(\cos y \cdot \left(1.5 \cdot \left(3 - \sqrt{5}\right)\right) + 1.5 \cdot \left(\cos x \cdot \left(\sqrt{5} - 1\right)\right)\right)}\\ \mathbf{else}:\\ \;\;\;\;\frac{2 + \left(\cos x - \cos y\right) \cdot \left(\left(\sin y - \frac{\sin x}{16}\right) \cdot \left(\sqrt{2} \cdot \left(\sin x - \frac{\sin y}{16}\right)\right)\right)}{3 + \left(\cos y \cdot \left(1.5 \cdot \left(3 - \sqrt{5}\right)\right) + 1.5 \cdot \left(\sqrt{5} - 1\right)\right)}\\ \end{array}\]

Alternative 6
Error	23.2
Cost	66112

\[\frac{2 + \left(\cos x - \cos y\right) \cdot \left(\left(\sin y - \frac{\sin x}{16}\right) \cdot \left(\sin x \cdot \sqrt{2}\right)\right)}{3 + \left(\cos y \cdot \left(1.5 \cdot \left(3 - \sqrt{5}\right)\right) + 1.5 \cdot \left(\cos x \cdot \left(\sqrt{5} - 1\right)\right)\right)}\]

Alternative 7
Error	52.1
Cost	64

\[1\]

Derivation

1. Initial program 0.5

   \[\frac{2 + \left(\left(\sqrt{2} \cdot \left(\sin x - \frac{\sin y}{16}\right)\right) \cdot \left(\sin y - \frac{\sin x}{16}\right)\right) \cdot \left(\cos x - \cos y\right)}{3 \cdot \left(\left(1 + \frac{\sqrt{5} - 1}{2} \cdot \cos x\right) + \frac{3 - \sqrt{5}}{2} \cdot \cos y\right)}\]

2. Simplified 0.4

   \[\leadsto \color{blue}{\frac{2 + \left(\left(\sqrt{2} \cdot \left(\sin x - \frac{\sin y}{16}\right)\right) \cdot \left(\sin y - \frac{\sin x}{16}\right)\right) \cdot \left(\cos x - \cos y\right)}{3 + \left(\cos y \cdot \left(1.5 \cdot \left(3 - \sqrt{5}\right)\right) + 1.5 \cdot \left(\cos x \cdot \left(\sqrt{5} - 1\right)\right)\right)}}\]
3. Using strategy rm

4. Applied add-cube-cbrt_binary64_6592 0.4

   \[\leadsto \frac{2 + \left(\left(\sqrt{2} \cdot \left(\sin x - \frac{\sin y}{\color{blue}{\left(\sqrt[3]{16} \cdot \sqrt[3]{16}\right) \cdot \sqrt[3]{16}}}\right)\right) \cdot \left(\sin y - \frac{\sin x}{16}\right)\right) \cdot \left(\cos x - \cos y\right)}{3 + \left(\cos y \cdot \left(1.5 \cdot \left(3 - \sqrt{5}\right)\right) + 1.5 \cdot \left(\cos x \cdot \left(\sqrt{5} - 1\right)\right)\right)}\]

5. Applied add-cube-cbrt_binary64_6592 0.4

   \[\leadsto \frac{2 + \left(\left(\sqrt{2} \cdot \left(\sin x - \frac{\color{blue}{\left(\sqrt[3]{\sin y} \cdot \sqrt[3]{\sin y}\right) \cdot \sqrt[3]{\sin y}}}{\left(\sqrt[3]{16} \cdot \sqrt[3]{16}\right) \cdot \sqrt[3]{16}}\right)\right) \cdot \left(\sin y - \frac{\sin x}{16}\right)\right) \cdot \left(\cos x - \cos y\right)}{3 + \left(\cos y \cdot \left(1.5 \cdot \left(3 - \sqrt{5}\right)\right) + 1.5 \cdot \left(\cos x \cdot \left(\sqrt{5} - 1\right)\right)\right)}\]

6. Applied times-frac_binary64_6563 0.4

   \[\leadsto \frac{2 + \left(\left(\sqrt{2} \cdot \left(\sin x - \color{blue}{\frac{\sqrt[3]{\sin y} \cdot \sqrt[3]{\sin y}}{\sqrt[3]{16} \cdot \sqrt[3]{16}} \cdot \frac{\sqrt[3]{\sin y}}{\sqrt[3]{16}}}\right)\right) \cdot \left(\sin y - \frac{\sin x}{16}\right)\right) \cdot \left(\cos x - \cos y\right)}{3 + \left(\cos y \cdot \left(1.5 \cdot \left(3 - \sqrt{5}\right)\right) + 1.5 \cdot \left(\cos x \cdot \left(\sqrt{5} - 1\right)\right)\right)}\]

7. Applied cancel-sign-sub-inv_binary64_6523 0.4

   \[\leadsto \frac{2 + \left(\left(\sqrt{2} \cdot \color{blue}{\left(\sin x + \left(-\frac{\sqrt[3]{\sin y} \cdot \sqrt[3]{\sin y}}{\sqrt[3]{16} \cdot \sqrt[3]{16}}\right) \cdot \frac{\sqrt[3]{\sin y}}{\sqrt[3]{16}}\right)}\right) \cdot \left(\sin y - \frac{\sin x}{16}\right)\right) \cdot \left(\cos x - \cos y\right)}{3 + \left(\cos y \cdot \left(1.5 \cdot \left(3 - \sqrt{5}\right)\right) + 1.5 \cdot \left(\cos x \cdot \left(\sqrt{5} - 1\right)\right)\right)}\]

8. Applied distribute-rgt-in_binary64_6507 0.4

   \[\leadsto \frac{2 + \left(\color{blue}{\left(\sin x \cdot \sqrt{2} + \left(\left(-\frac{\sqrt[3]{\sin y} \cdot \sqrt[3]{\sin y}}{\sqrt[3]{16} \cdot \sqrt[3]{16}}\right) \cdot \frac{\sqrt[3]{\sin y}}{\sqrt[3]{16}}\right) \cdot \sqrt{2}\right)} \cdot \left(\sin y - \frac{\sin x}{16}\right)\right) \cdot \left(\cos x - \cos y\right)}{3 + \left(\cos y \cdot \left(1.5 \cdot \left(3 - \sqrt{5}\right)\right) + 1.5 \cdot \left(\cos x \cdot \left(\sqrt{5} - 1\right)\right)\right)}\]

9. Simplified 0.4

   \[\leadsto \frac{2 + \left(\left(\color{blue}{\sqrt{2} \cdot \sin x} + \left(\left(-\frac{\sqrt[3]{\sin y} \cdot \sqrt[3]{\sin y}}{\sqrt[3]{16} \cdot \sqrt[3]{16}}\right) \cdot \frac{\sqrt[3]{\sin y}}{\sqrt[3]{16}}\right) \cdot \sqrt{2}\right) \cdot \left(\sin y - \frac{\sin x}{16}\right)\right) \cdot \left(\cos x - \cos y\right)}{3 + \left(\cos y \cdot \left(1.5 \cdot \left(3 - \sqrt{5}\right)\right) + 1.5 \cdot \left(\cos x \cdot \left(\sqrt{5} - 1\right)\right)\right)}\]

10. Simplified 0.4

    \[\leadsto \frac{2 + \left(\left(\sqrt{2} \cdot \sin x + \color{blue}{\sqrt{2} \cdot \left(-{\left(\frac{\sqrt[3]{\sin y}}{\sqrt[3]{16}}\right)}^{3}\right)}\right) \cdot \left(\sin y - \frac{\sin x}{16}\right)\right) \cdot \left(\cos x - \cos y\right)}{3 + \left(\cos y \cdot \left(1.5 \cdot \left(3 - \sqrt{5}\right)\right) + 1.5 \cdot \left(\cos x \cdot \left(\sqrt{5} - 1\right)\right)\right)}\]
11. Using strategy rm

12. Applied add-cube-cbrt_binary64_6592 0.4

    \[\leadsto \frac{2 + \left(\left(\sqrt{2} \cdot \sin x + \sqrt{2} \cdot \left(-{\left(\frac{\sqrt[3]{\sin y}}{\color{blue}{\left(\sqrt[3]{\sqrt[3]{16}} \cdot \sqrt[3]{\sqrt[3]{16}}\right) \cdot \sqrt[3]{\sqrt[3]{16}}}}\right)}^{3}\right)\right) \cdot \left(\sin y - \frac{\sin x}{16}\right)\right) \cdot \left(\cos x - \cos y\right)}{3 + \left(\cos y \cdot \left(1.5 \cdot \left(3 - \sqrt{5}\right)\right) + 1.5 \cdot \left(\cos x \cdot \left(\sqrt{5} - 1\right)\right)\right)}\]

13. Simplified 0.4

    \[\leadsto \frac{2 + \left(\left(\sqrt{2} \cdot \sin x + \sqrt{2} \cdot \left(-{\left(\frac{\sqrt[3]{\sin y}}{\color{blue}{\sqrt[3]{\sqrt[3]{16} \cdot \sqrt[3]{16}}} \cdot \sqrt[3]{\sqrt[3]{16}}}\right)}^{3}\right)\right) \cdot \left(\sin y - \frac{\sin x}{16}\right)\right) \cdot \left(\cos x - \cos y\right)}{3 + \left(\cos y \cdot \left(1.5 \cdot \left(3 - \sqrt{5}\right)\right) + 1.5 \cdot \left(\cos x \cdot \left(\sqrt{5} - 1\right)\right)\right)}\]

14. Simplified 0.4

    \[\leadsto \color{blue}{\frac{2 + \left(\cos x - \cos y\right) \cdot \left(\left(\sin y - \frac{\sin x}{16}\right) \cdot \left(\sqrt{2} \cdot \sin x - \sqrt{2} \cdot {\left(\frac{\sqrt[3]{\sin y}}{\sqrt[3]{\sqrt[3]{16} \cdot \sqrt[3]{16}} \cdot \sqrt[3]{\sqrt[3]{16}}}\right)}^{3}\right)\right)}{3 + \left(\cos y \cdot \left(1.5 \cdot \left(3 - \sqrt{5}\right)\right) + 1.5 \cdot \left(\cos x \cdot \left(\sqrt{5} - 1\right)\right)\right)}}\]

15. Final simplification 0.4

    \[\leadsto \frac{2 + \left(\cos x - \cos y\right) \cdot \left(\left(\sin y - \frac{\sin x}{16}\right) \cdot \left(\sin x \cdot \sqrt{2} - \sqrt{2} \cdot {\left(\frac{\sqrt[3]{\sin y}}{\sqrt[3]{\sqrt[3]{16} \cdot \sqrt[3]{16}} \cdot \sqrt[3]{\sqrt[3]{16}}}\right)}^{3}\right)\right)}{3 + \left(\cos y \cdot \left(1.5 \cdot \left(3 - \sqrt{5}\right)\right) + 1.5 \cdot \left(\cos x \cdot \left(\sqrt{5} - 1\right)\right)\right)}\]

Reproduce

herbie shell --seed 2021044 
(FPCore (x y)
  :name "Diagrams.TwoD.Path.Metafont.Internal:hobbyF from diagrams-contrib-1.3.0.5"
  :precision binary64
  (/ (+ 2.0 (* (* (* (sqrt 2.0) (- (sin x) (/ (sin y) 16.0))) (- (sin y) (/ (sin x) 16.0))) (- (cos x) (cos y)))) (* 3.0 (+ (+ 1.0 (* (/ (- (sqrt 5.0) 1.0) 2.0) (cos x))) (* (/ (- 3.0 (sqrt 5.0)) 2.0) (cos y))))))