\[d1 \cdot d2 + d1 \cdot d3\]

↓

\[\left(d3 + d2\right) \cdot d1\]

double f(double d1, double d2, double d3) {
        double r2310753 = d1;
        double r2310754 = d2;
        double r2310755 = r2310753 * r2310754;
        double r2310756 = d3;
        double r2310757 = r2310753 * r2310756;
        double r2310758 = r2310755 + r2310757;
        return r2310758;
}

↓

double f(double d1, double d2, double d3) {
        double r2310759 = d3;
        double r2310760 = d2;
        double r2310761 = r2310759 + r2310760;
        double r2310762 = d1;
        double r2310763 = r2310761 * r2310762;
        return r2310763;
}

d1 \cdot d2 + d1 \cdot d3

↓

\left(d3 + d2\right) \cdot d1

Error

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

Derivation

Initial program 0.4
\[\frac{\left(d1 \cdot d2\right)}{\left(d1 \cdot d3\right)}\]
Simplified0.2
\[\leadsto \color{blue}{\left(\frac{d3}{d2}\right) \cdot d1}\]
Final simplification0.2
\[\leadsto \left(d3 + d2\right) \cdot d1\]

Reproduce

herbie shell --seed 2019102 
(FPCore (d1 d2 d3)
  :name "FastMath dist"
  (+.p16 (*.p16 d1 d2) (*.p16 d1 d3)))