math.log10 on complex, imaginary part

Average Error: 0.8 → 0.1

Time: 3.6s

Precision: 64

Math

\[\frac{\tan^{-1}_* \frac{im}{re}}{\log 10}\]

↓

\[\sqrt{\sqrt{\frac{1}{\sqrt{\log 10}}}} \cdot \left(\sqrt{\sqrt{\frac{1}{\sqrt{\log 10}}}} \cdot \left(\sqrt{\frac{1}{\sqrt{\log 10}}} \cdot \left(\tan^{-1}_* \frac{im}{re} \cdot \sqrt{\frac{1}{\log 10}}\right)\right)\right)\]

C

double f(double re, double im) {
        double r29580 = im;
        double r29581 = re;
        double r29582 = atan2(r29580, r29581);
        double r29583 = 10.0;
        double r29584 = log(r29583);
        double r29585 = r29582 / r29584;
        return r29585;
}

↓

double f(double re, double im) {
        double r29586 = 1.0;
        double r29587 = 10.0;
        double r29588 = log(r29587);
        double r29589 = sqrt(r29588);
        double r29590 = r29586 / r29589;
        double r29591 = sqrt(r29590);
        double r29592 = sqrt(r29591);
        double r29593 = im;
        double r29594 = re;
        double r29595 = atan2(r29593, r29594);
        double r29596 = r29586 / r29588;
        double r29597 = sqrt(r29596);
        double r29598 = r29595 * r29597;
        double r29599 = r29591 * r29598;
        double r29600 = r29592 * r29599;
        double r29601 = r29592 * r29600;
        return r29601;
}

Error

Bits error versus re

Bits error versus im

Derivation

1. Initial program 0.8

   \[\frac{\tan^{-1}_* \frac{im}{re}}{\log 10}\]
2. Using strategy rm

3. Applied add-sqr-sqrt 0.8

   \[\leadsto \frac{\tan^{-1}_* \frac{im}{re}}{\color{blue}{\sqrt{\log 10} \cdot \sqrt{\log 10}}}\]

4. Applied *-un-lft-identity 0.8

   \[\leadsto \frac{\color{blue}{1 \cdot \tan^{-1}_* \frac{im}{re}}}{\sqrt{\log 10} \cdot \sqrt{\log 10}}\]

5. Applied times-frac 0.8

   \[\leadsto \color{blue}{\frac{1}{\sqrt{\log 10}} \cdot \frac{\tan^{-1}_* \frac{im}{re}}{\sqrt{\log 10}}}\]

6. Taylor expanded around 0 0.8

   \[\leadsto \frac{1}{\sqrt{\log 10}} \cdot \color{blue}{\left(\tan^{-1}_* \frac{im}{re} \cdot \sqrt{\frac{1}{\log 10}}\right)}\]
7. Using strategy rm

8. Applied add-sqr-sqrt 0.8

   \[\leadsto \color{blue}{\left(\sqrt{\frac{1}{\sqrt{\log 10}}} \cdot \sqrt{\frac{1}{\sqrt{\log 10}}}\right)} \cdot \left(\tan^{-1}_* \frac{im}{re} \cdot \sqrt{\frac{1}{\log 10}}\right)\]

9. Applied associate-*l* 0.8

   \[\leadsto \color{blue}{\sqrt{\frac{1}{\sqrt{\log 10}}} \cdot \left(\sqrt{\frac{1}{\sqrt{\log 10}}} \cdot \left(\tan^{-1}_* \frac{im}{re} \cdot \sqrt{\frac{1}{\log 10}}\right)\right)}\]
10. Using strategy rm

11. Applied add-sqr-sqrt 0.8

    \[\leadsto \sqrt{\color{blue}{\sqrt{\frac{1}{\sqrt{\log 10}}} \cdot \sqrt{\frac{1}{\sqrt{\log 10}}}}} \cdot \left(\sqrt{\frac{1}{\sqrt{\log 10}}} \cdot \left(\tan^{-1}_* \frac{im}{re} \cdot \sqrt{\frac{1}{\log 10}}\right)\right)\]

12. Applied sqrt-prod 0.1

    \[\leadsto \color{blue}{\left(\sqrt{\sqrt{\frac{1}{\sqrt{\log 10}}}} \cdot \sqrt{\sqrt{\frac{1}{\sqrt{\log 10}}}}\right)} \cdot \left(\sqrt{\frac{1}{\sqrt{\log 10}}} \cdot \left(\tan^{-1}_* \frac{im}{re} \cdot \sqrt{\frac{1}{\log 10}}\right)\right)\]

13. Applied associate-*l* 0.1

    \[\leadsto \color{blue}{\sqrt{\sqrt{\frac{1}{\sqrt{\log 10}}}} \cdot \left(\sqrt{\sqrt{\frac{1}{\sqrt{\log 10}}}} \cdot \left(\sqrt{\frac{1}{\sqrt{\log 10}}} \cdot \left(\tan^{-1}_* \frac{im}{re} \cdot \sqrt{\frac{1}{\log 10}}\right)\right)\right)}\]

14. Final simplification 0.1

    \[\leadsto \sqrt{\sqrt{\frac{1}{\sqrt{\log 10}}}} \cdot \left(\sqrt{\sqrt{\frac{1}{\sqrt{\log 10}}}} \cdot \left(\sqrt{\frac{1}{\sqrt{\log 10}}} \cdot \left(\tan^{-1}_* \frac{im}{re} \cdot \sqrt{\frac{1}{\log 10}}\right)\right)\right)\]

Reproduce

herbie shell --seed 2020002 +o rules:numerics
(FPCore (re im)
  :name "math.log10 on complex, imaginary part"
  :precision binary64
  (/ (atan2 im re) (log 10)))