Complex division, real part

Average Error: 26.5 → 13.8

Time: 4.2s

Precision: 64

Math

\[\frac{a \cdot c + b \cdot d}{c \cdot c + d \cdot d}\]

↓

\[\begin{array}{l} \mathbf{if}\;c \le -4.1536238602847571 \cdot 10^{76}:\\ \;\;\;\;\frac{-1 \cdot a}{\mathsf{hypot}\left(c, d\right) \cdot 1}\\ \mathbf{elif}\;c \le 1.7187819587037017 \cdot 10^{178}:\\ \;\;\;\;\frac{\frac{\mathsf{fma}\left(a, c, b \cdot d\right)}{\mathsf{hypot}\left(c, d\right)}}{\mathsf{hypot}\left(c, d\right) \cdot 1}\\ \mathbf{else}:\\ \;\;\;\;\frac{a}{\mathsf{hypot}\left(c, d\right) \cdot 1}\\ \end{array}\]

C

double f(double a, double b, double c, double d) {
        double r108455 = a;
        double r108456 = c;
        double r108457 = r108455 * r108456;
        double r108458 = b;
        double r108459 = d;
        double r108460 = r108458 * r108459;
        double r108461 = r108457 + r108460;
        double r108462 = r108456 * r108456;
        double r108463 = r108459 * r108459;
        double r108464 = r108462 + r108463;
        double r108465 = r108461 / r108464;
        return r108465;
}

↓

double f(double a, double b, double c, double d) {
        double r108466 = c;
        double r108467 = -4.153623860284757e+76;
        bool r108468 = r108466 <= r108467;
        double r108469 = -1.0;
        double r108470 = a;
        double r108471 = r108469 * r108470;
        double r108472 = d;
        double r108473 = hypot(r108466, r108472);
        double r108474 = 1.0;
        double r108475 = r108473 * r108474;
        double r108476 = r108471 / r108475;
        double r108477 = 1.7187819587037017e+178;
        bool r108478 = r108466 <= r108477;
        double r108479 = b;
        double r108480 = r108479 * r108472;
        double r108481 = fma(r108470, r108466, r108480);
        double r108482 = r108481 / r108473;
        double r108483 = r108482 / r108475;
        double r108484 = r108470 / r108475;
        double r108485 = r108478 ? r108483 : r108484;
        double r108486 = r108468 ? r108476 : r108485;
        return r108486;
}

Error

Bits error versus a

Bits error versus b

Bits error versus c

Bits error versus d

Target

Original	26.5
Target	0.3
Herbie	13.8

\[\begin{array}{l} \mathbf{if}\;\left|d\right| \lt \left|c\right|:\\ \;\;\;\;\frac{a + b \cdot \frac{d}{c}}{c + d \cdot \frac{d}{c}}\\ \mathbf{else}:\\ \;\;\;\;\frac{b + a \cdot \frac{c}{d}}{d + c \cdot \frac{c}{d}}\\ \end{array}\]

Derivation

1. Split input into 3 regimes

2. if c < -4.153623860284757e+76

   1. Initial program 38.7

      \[\frac{a \cdot c + b \cdot d}{c \cdot c + d \cdot d}\]
   2. Using strategy rm

   3. Applied add-sqr-sqrt 38.7

      \[\leadsto \frac{a \cdot c + b \cdot d}{\color{blue}{\sqrt{c \cdot c + d \cdot d} \cdot \sqrt{c \cdot c + d \cdot d}}}\]

   4. Applied *-un-lft-identity 38.7

      \[\leadsto \frac{\color{blue}{1 \cdot \left(a \cdot c + b \cdot d\right)}}{\sqrt{c \cdot c + d \cdot d} \cdot \sqrt{c \cdot c + d \cdot d}}\]

   5. Applied times-frac 38.7

      \[\leadsto \color{blue}{\frac{1}{\sqrt{c \cdot c + d \cdot d}} \cdot \frac{a \cdot c + b \cdot d}{\sqrt{c \cdot c + d \cdot d}}}\]

   6. Simplified 38.7

      \[\leadsto \color{blue}{\frac{1}{\mathsf{hypot}\left(c, d\right) \cdot 1}} \cdot \frac{a \cdot c + b \cdot d}{\sqrt{c \cdot c + d \cdot d}}\]

   7. Simplified 26.7

      \[\leadsto \frac{1}{\mathsf{hypot}\left(c, d\right) \cdot 1} \cdot \color{blue}{\frac{\mathsf{fma}\left(a, c, b \cdot d\right)}{\mathsf{hypot}\left(c, d\right) \cdot 1}}\]
   8. Using strategy rm

   9. Applied associate-*r/ 26.7

      \[\leadsto \color{blue}{\frac{\frac{1}{\mathsf{hypot}\left(c, d\right) \cdot 1} \cdot \mathsf{fma}\left(a, c, b \cdot d\right)}{\mathsf{hypot}\left(c, d\right) \cdot 1}}\]

   10. Simplified 26.6

       \[\leadsto \frac{\color{blue}{\frac{\mathsf{fma}\left(a, c, b \cdot d\right)}{\mathsf{hypot}\left(c, d\right)}}}{\mathsf{hypot}\left(c, d\right) \cdot 1}\]

   11. Taylor expanded around -inf 18.0

       \[\leadsto \frac{\color{blue}{-1 \cdot a}}{\mathsf{hypot}\left(c, d\right) \cdot 1}\]

   if -4.153623860284757e+76 < c < 1.7187819587037017e+178

   1. Initial program 20.0

      \[\frac{a \cdot c + b \cdot d}{c \cdot c + d \cdot d}\]
   2. Using strategy rm

   3. Applied add-sqr-sqrt 20.0

      \[\leadsto \frac{a \cdot c + b \cdot d}{\color{blue}{\sqrt{c \cdot c + d \cdot d} \cdot \sqrt{c \cdot c + d \cdot d}}}\]

   4. Applied *-un-lft-identity 20.0

      \[\leadsto \frac{\color{blue}{1 \cdot \left(a \cdot c + b \cdot d\right)}}{\sqrt{c \cdot c + d \cdot d} \cdot \sqrt{c \cdot c + d \cdot d}}\]

   5. Applied times-frac 20.0

      \[\leadsto \color{blue}{\frac{1}{\sqrt{c \cdot c + d \cdot d}} \cdot \frac{a \cdot c + b \cdot d}{\sqrt{c \cdot c + d \cdot d}}}\]

   6. Simplified 20.0

      \[\leadsto \color{blue}{\frac{1}{\mathsf{hypot}\left(c, d\right) \cdot 1}} \cdot \frac{a \cdot c + b \cdot d}{\sqrt{c \cdot c + d \cdot d}}\]

   7. Simplified 12.6

      \[\leadsto \frac{1}{\mathsf{hypot}\left(c, d\right) \cdot 1} \cdot \color{blue}{\frac{\mathsf{fma}\left(a, c, b \cdot d\right)}{\mathsf{hypot}\left(c, d\right) \cdot 1}}\]
   8. Using strategy rm

   9. Applied associate-*r/ 12.6

      \[\leadsto \color{blue}{\frac{\frac{1}{\mathsf{hypot}\left(c, d\right) \cdot 1} \cdot \mathsf{fma}\left(a, c, b \cdot d\right)}{\mathsf{hypot}\left(c, d\right) \cdot 1}}\]

   10. Simplified 12.5

       \[\leadsto \frac{\color{blue}{\frac{\mathsf{fma}\left(a, c, b \cdot d\right)}{\mathsf{hypot}\left(c, d\right)}}}{\mathsf{hypot}\left(c, d\right) \cdot 1}\]

   if 1.7187819587037017e+178 < c

   1. Initial program 45.1

      \[\frac{a \cdot c + b \cdot d}{c \cdot c + d \cdot d}\]
   2. Using strategy rm

   3. Applied add-sqr-sqrt 45.1

      \[\leadsto \frac{a \cdot c + b \cdot d}{\color{blue}{\sqrt{c \cdot c + d \cdot d} \cdot \sqrt{c \cdot c + d \cdot d}}}\]

   4. Applied *-un-lft-identity 45.1

      \[\leadsto \frac{\color{blue}{1 \cdot \left(a \cdot c + b \cdot d\right)}}{\sqrt{c \cdot c + d \cdot d} \cdot \sqrt{c \cdot c + d \cdot d}}\]

   5. Applied times-frac 45.1

      \[\leadsto \color{blue}{\frac{1}{\sqrt{c \cdot c + d \cdot d}} \cdot \frac{a \cdot c + b \cdot d}{\sqrt{c \cdot c + d \cdot d}}}\]

   6. Simplified 45.1

      \[\leadsto \color{blue}{\frac{1}{\mathsf{hypot}\left(c, d\right) \cdot 1}} \cdot \frac{a \cdot c + b \cdot d}{\sqrt{c \cdot c + d \cdot d}}\]

   7. Simplified 31.0

      \[\leadsto \frac{1}{\mathsf{hypot}\left(c, d\right) \cdot 1} \cdot \color{blue}{\frac{\mathsf{fma}\left(a, c, b \cdot d\right)}{\mathsf{hypot}\left(c, d\right) \cdot 1}}\]
   8. Using strategy rm

   9. Applied associate-*r/ 31.0

      \[\leadsto \color{blue}{\frac{\frac{1}{\mathsf{hypot}\left(c, d\right) \cdot 1} \cdot \mathsf{fma}\left(a, c, b \cdot d\right)}{\mathsf{hypot}\left(c, d\right) \cdot 1}}\]

   10. Simplified 30.9

       \[\leadsto \frac{\color{blue}{\frac{\mathsf{fma}\left(a, c, b \cdot d\right)}{\mathsf{hypot}\left(c, d\right)}}}{\mathsf{hypot}\left(c, d\right) \cdot 1}\]

   11. Taylor expanded around inf 14.4

       \[\leadsto \frac{\color{blue}{a}}{\mathsf{hypot}\left(c, d\right) \cdot 1}\]
3. Recombined 3 regimes into one program.

4. Final simplification 13.8

   \[\leadsto \begin{array}{l} \mathbf{if}\;c \le -4.1536238602847571 \cdot 10^{76}:\\ \;\;\;\;\frac{-1 \cdot a}{\mathsf{hypot}\left(c, d\right) \cdot 1}\\ \mathbf{elif}\;c \le 1.7187819587037017 \cdot 10^{178}:\\ \;\;\;\;\frac{\frac{\mathsf{fma}\left(a, c, b \cdot d\right)}{\mathsf{hypot}\left(c, d\right)}}{\mathsf{hypot}\left(c, d\right) \cdot 1}\\ \mathbf{else}:\\ \;\;\;\;\frac{a}{\mathsf{hypot}\left(c, d\right) \cdot 1}\\ \end{array}\]

Reproduce

herbie shell --seed 2020027 +o rules:numerics
(FPCore (a b c d)
  :name "Complex division, real part"
  :precision binary64

  :herbie-target
  (if (< (fabs d) (fabs c)) (/ (+ a (* b (/ d c))) (+ c (* d (/ d c)))) (/ (+ b (* a (/ c d))) (+ d (* c (/ c d)))))

  (/ (+ (* a c) (* b d)) (+ (* c c) (* d d))))