Average Error: 12.5 → 9.0
Time: 32.2s
Precision: 64
Internal Precision: 128
\[\frac{\sin ky}{\sqrt{{\left(\sin kx\right)}^{2} + {\left(\sin ky\right)}^{2}}} \cdot \sin th\]
\[\sin ky \cdot \frac{\sin th}{\sqrt{\left(\sin kx\right)^2 + \left(\sin ky\right)^2}^*}\]

Error

Bits error versus kx

Bits error versus ky

Bits error versus th

Derivation

  1. Initial program 12.5

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

    \[\leadsto \color{blue}{\frac{\sin th \cdot \sin ky}{\sqrt{\left(\sin kx\right)^2 + \left(\sin ky\right)^2}^*}}\]
  3. Using strategy rm
  4. Applied *-un-lft-identity11.4

    \[\leadsto \frac{\sin th \cdot \sin ky}{\color{blue}{1 \cdot \sqrt{\left(\sin kx\right)^2 + \left(\sin ky\right)^2}^*}}\]
  5. Applied times-frac9.0

    \[\leadsto \color{blue}{\frac{\sin th}{1} \cdot \frac{\sin ky}{\sqrt{\left(\sin kx\right)^2 + \left(\sin ky\right)^2}^*}}\]
  6. Simplified9.0

    \[\leadsto \color{blue}{\sin th} \cdot \frac{\sin ky}{\sqrt{\left(\sin kx\right)^2 + \left(\sin ky\right)^2}^*}\]
  7. Using strategy rm
  8. Applied expm1-log1p-u9.0

    \[\leadsto \sin th \cdot \color{blue}{(e^{\log_* (1 + \frac{\sin ky}{\sqrt{\left(\sin kx\right)^2 + \left(\sin ky\right)^2}^*})} - 1)^*}\]
  9. Using strategy rm
  10. Applied pow19.0

    \[\leadsto \sin th \cdot \color{blue}{{\left((e^{\log_* (1 + \frac{\sin ky}{\sqrt{\left(\sin kx\right)^2 + \left(\sin ky\right)^2}^*})} - 1)^*\right)}^{1}}\]
  11. Applied pow19.0

    \[\leadsto \color{blue}{{\left(\sin th\right)}^{1}} \cdot {\left((e^{\log_* (1 + \frac{\sin ky}{\sqrt{\left(\sin kx\right)^2 + \left(\sin ky\right)^2}^*})} - 1)^*\right)}^{1}\]
  12. Applied pow-prod-down9.0

    \[\leadsto \color{blue}{{\left(\sin th \cdot (e^{\log_* (1 + \frac{\sin ky}{\sqrt{\left(\sin kx\right)^2 + \left(\sin ky\right)^2}^*})} - 1)^*\right)}^{1}}\]
  13. Simplified9.0

    \[\leadsto {\color{blue}{\left(\sin ky \cdot \frac{\sin th}{\sqrt{\left(\sin kx\right)^2 + \left(\sin ky\right)^2}^*}\right)}}^{1}\]
  14. Final simplification9.0

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

Reproduce

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