Average Error: 33.6 → 0.4
Time: 20.6s
Precision: 64
\[\frac{x \cdot x}{y \cdot y} + \frac{z \cdot z}{t \cdot t}\]
\[\mathsf{fma}\left(\frac{x}{y}, \frac{x}{y}, \frac{\frac{z}{t}}{\frac{t}{z}}\right)\]
\frac{x \cdot x}{y \cdot y} + \frac{z \cdot z}{t \cdot t}
\mathsf{fma}\left(\frac{x}{y}, \frac{x}{y}, \frac{\frac{z}{t}}{\frac{t}{z}}\right)
double f(double x, double y, double z, double t) {
        double r435510 = x;
        double r435511 = r435510 * r435510;
        double r435512 = y;
        double r435513 = r435512 * r435512;
        double r435514 = r435511 / r435513;
        double r435515 = z;
        double r435516 = r435515 * r435515;
        double r435517 = t;
        double r435518 = r435517 * r435517;
        double r435519 = r435516 / r435518;
        double r435520 = r435514 + r435519;
        return r435520;
}

double f(double x, double y, double z, double t) {
        double r435521 = x;
        double r435522 = y;
        double r435523 = r435521 / r435522;
        double r435524 = z;
        double r435525 = t;
        double r435526 = r435524 / r435525;
        double r435527 = r435525 / r435524;
        double r435528 = r435526 / r435527;
        double r435529 = fma(r435523, r435523, r435528);
        return r435529;
}

Error

Bits error versus x

Bits error versus y

Bits error versus z

Bits error versus t

Target

Original33.6
Target0.4
Herbie0.4
\[{\left(\frac{x}{y}\right)}^{2} + {\left(\frac{z}{t}\right)}^{2}\]

Derivation

  1. Initial program 33.6

    \[\frac{x \cdot x}{y \cdot y} + \frac{z \cdot z}{t \cdot t}\]
  2. Simplified18.7

    \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{x}{y}, \frac{x}{y}, \frac{z \cdot z}{t \cdot t}\right)}\]
  3. Using strategy rm
  4. Applied associate-/l*12.9

    \[\leadsto \mathsf{fma}\left(\frac{x}{y}, \frac{x}{y}, \color{blue}{\frac{z}{\frac{t \cdot t}{z}}}\right)\]
  5. Simplified3.9

    \[\leadsto \mathsf{fma}\left(\frac{x}{y}, \frac{x}{y}, \frac{z}{\color{blue}{t \cdot \frac{t}{z}}}\right)\]
  6. Using strategy rm
  7. Applied associate-/r*0.4

    \[\leadsto \mathsf{fma}\left(\frac{x}{y}, \frac{x}{y}, \color{blue}{\frac{\frac{z}{t}}{\frac{t}{z}}}\right)\]
  8. Final simplification0.4

    \[\leadsto \mathsf{fma}\left(\frac{x}{y}, \frac{x}{y}, \frac{\frac{z}{t}}{\frac{t}{z}}\right)\]

Reproduce

herbie shell --seed 2019323 +o rules:numerics
(FPCore (x y z t)
  :name "Graphics.Rasterific.Svg.PathConverter:arcToSegments from rasterific-svg-0.2.3.1"
  :precision binary64

  :herbie-target
  (+ (pow (/ x y) 2) (pow (/ z t) 2))

  (+ (/ (* x x) (* y y)) (/ (* z z) (* t t))))