math.exp on complex, imaginary part

Average Error: 0.0 → 0.0

Time: 12.3s

Precision: 64

• Falling back on racket egraph because egg-herbie package not installed

Math

\[e^{re} \cdot \sin im\]

↓

\[e^{re} \cdot \sin im\]

C

double f(double re, double im) {
        double r76319 = re;
        double r76320 = exp(r76319);
        double r76321 = im;
        double r76322 = sin(r76321);
        double r76323 = r76320 * r76322;
        return r76323;
}

↓

double f(double re, double im) {
        double r76324 = re;
        double r76325 = exp(r76324);
        double r76326 = im;
        double r76327 = sin(r76326);
        double r76328 = r76325 * r76327;
        return r76328;
}

Error

Bits error versus re

Bits error versus im

Derivation

1. Initial program 0.0

   \[e^{re} \cdot \sin im\]
2. Using strategy rm

3. Applied add-sqr-sqrt 0.0

   \[\leadsto \color{blue}{\left(\sqrt{e^{re}} \cdot \sqrt{e^{re}}\right)} \cdot \sin im\]

4. Applied associate-*l* 0.0

   \[\leadsto \color{blue}{\sqrt{e^{re}} \cdot \left(\sqrt{e^{re}} \cdot \sin im\right)}\]
5. Using strategy rm

6. Applied pow1 0.0

   \[\leadsto \sqrt{e^{re}} \cdot \left(\sqrt{e^{re}} \cdot \color{blue}{{\left(\sin im\right)}^{1}}\right)\]

7. Applied pow1 0.0

   \[\leadsto \sqrt{e^{re}} \cdot \left(\color{blue}{{\left(\sqrt{e^{re}}\right)}^{1}} \cdot {\left(\sin im\right)}^{1}\right)\]

8. Applied pow-prod-down 0.0

   \[\leadsto \sqrt{e^{re}} \cdot \color{blue}{{\left(\sqrt{e^{re}} \cdot \sin im\right)}^{1}}\]

9. Applied pow1 0.0

   \[\leadsto \color{blue}{{\left(\sqrt{e^{re}}\right)}^{1}} \cdot {\left(\sqrt{e^{re}} \cdot \sin im\right)}^{1}\]

10. Applied pow-prod-down 0.0

    \[\leadsto \color{blue}{{\left(\sqrt{e^{re}} \cdot \left(\sqrt{e^{re}} \cdot \sin im\right)\right)}^{1}}\]

11. Simplified 0.0

    \[\leadsto {\color{blue}{\left(\sin im \cdot e^{re}\right)}}^{1}\]

12. Final simplification 0.0

    \[\leadsto e^{re} \cdot \sin im\]

Reproduce

herbie shell --seed 2020047 +o rules:numerics
(FPCore (re im)
  :name "math.exp on complex, imaginary part"
  :precision binary64
  (* (exp re) (sin im)))