Average Error: 31.1 → 0.6

Time: 38.6s

Precision: 64

Internal Precision: 384

\[\frac{\log \left(\sqrt{re \cdot re + im \cdot im}\right)}{\log 10}\]

↓

\[\frac{\log \left(\left(\sqrt[3]{\sqrt{re^2 + im^2}^*} \cdot \sqrt[3]{\sqrt{re^2 + im^2}^*}\right) \cdot \sqrt[3]{\sqrt{re^2 + im^2}^*}\right)}{\log 10}\]

Error

The X axis uses an exponential scale — Bits error versus `re`

Derivation

Initial program 31.1
\[\frac{\log \left(\sqrt{re \cdot re + im \cdot im}\right)}{\log 10}\]
Applied simplify0.6
\[\leadsto \color{blue}{\frac{\log \left(\sqrt{re^2 + im^2}^*\right)}{\log 10}}\]
Using strategy rm
Applied add-cube-cbrt0.6
\[\leadsto \frac{\log \color{blue}{\left(\left(\sqrt[3]{\sqrt{re^2 + im^2}^*} \cdot \sqrt[3]{\sqrt{re^2 + im^2}^*}\right) \cdot \sqrt[3]{\sqrt{re^2 + im^2}^*}\right)}}{\log 10}\]

Runtime

herbie shell --seed '#(1070960995 739739648 2531964651 3069671617 351857262 3877178482)' +o rules:numerics
(FPCore (re im)
  :name "math.log10 on complex, real part"
  (/ (log (sqrt (+ (* re re) (* im im)))) (log 10)))