Average Error: 12.8 → 13.0
Time: 45.6s
Precision: 64
Internal Precision: 128
\[\frac{\sin ky}{\sqrt{{\left(\sin kx\right)}^{2} + {\left(\sin ky\right)}^{2}}} \cdot \sin th\]
\[\sin th \cdot \frac{\sin ky}{\sqrt{\sin ky \cdot \sin ky + \left(\sqrt[3]{\sin kx} \cdot \sqrt[3]{\sin kx}\right) \cdot \left(\left(\sqrt[3]{\sin kx} \cdot \sqrt[3]{\sin kx}\right) \cdot \left(\sqrt[3]{\sin kx} \cdot \sqrt[3]{\sin kx}\right)\right)}}\]

Error

Bits error versus kx

Bits error versus ky

Bits error versus th

Derivation

  1. Initial program 12.8

    \[\frac{\sin ky}{\sqrt{{\left(\sin kx\right)}^{2} + {\left(\sin ky\right)}^{2}}} \cdot \sin th\]
  2. Simplified12.8

    \[\leadsto \color{blue}{\sin th \cdot \frac{\sin ky}{\sqrt{\sin kx \cdot \sin kx + \sin ky \cdot \sin ky}}}\]
  3. Using strategy rm
  4. Applied add-cube-cbrt12.9

    \[\leadsto \sin th \cdot \frac{\sin ky}{\sqrt{\sin kx \cdot \color{blue}{\left(\left(\sqrt[3]{\sin kx} \cdot \sqrt[3]{\sin kx}\right) \cdot \sqrt[3]{\sin kx}\right)} + \sin ky \cdot \sin ky}}\]
  5. Applied add-cube-cbrt13.0

    \[\leadsto \sin th \cdot \frac{\sin ky}{\sqrt{\color{blue}{\left(\left(\sqrt[3]{\sin kx} \cdot \sqrt[3]{\sin kx}\right) \cdot \sqrt[3]{\sin kx}\right)} \cdot \left(\left(\sqrt[3]{\sin kx} \cdot \sqrt[3]{\sin kx}\right) \cdot \sqrt[3]{\sin kx}\right) + \sin ky \cdot \sin ky}}\]
  6. Applied swap-sqr13.0

    \[\leadsto \sin th \cdot \frac{\sin ky}{\sqrt{\color{blue}{\left(\left(\sqrt[3]{\sin kx} \cdot \sqrt[3]{\sin kx}\right) \cdot \left(\sqrt[3]{\sin kx} \cdot \sqrt[3]{\sin kx}\right)\right) \cdot \left(\sqrt[3]{\sin kx} \cdot \sqrt[3]{\sin kx}\right)} + \sin ky \cdot \sin ky}}\]
  7. Final simplification13.0

    \[\leadsto \sin th \cdot \frac{\sin ky}{\sqrt{\sin ky \cdot \sin ky + \left(\sqrt[3]{\sin kx} \cdot \sqrt[3]{\sin kx}\right) \cdot \left(\left(\sqrt[3]{\sin kx} \cdot \sqrt[3]{\sin kx}\right) \cdot \left(\sqrt[3]{\sin kx} \cdot \sqrt[3]{\sin kx}\right)\right)}}\]

Reproduce

herbie shell --seed 2019094 
(FPCore (kx ky th)
  :name "Toniolo and Linder, Equation (3b), real"
  (* (/ (sin ky) (sqrt (+ (pow (sin kx) 2) (pow (sin ky) 2)))) (sin th)))