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

↓

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

double f(double d1, double d2, double d3) {
        double r2107010 = d1;
        double r2107011 = d2;
        double r2107012 = r2107010 * r2107011;
        double r2107013 = d3;
        double r2107014 = r2107010 * r2107013;
        double r2107015 = r2107012 + r2107014;
        return r2107015;
}

↓

double f(double d1, double d2, double d3) {
        double r2107016 = d3;
        double r2107017 = d2;
        double r2107018 = r2107016 + r2107017;
        double r2107019 = d1;
        double r2107020 = r2107018 * r2107019;
        return r2107020;
}

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.3
\[\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 2019101 
(FPCore (d1 d2 d3)
  :name "FastMath dist"
  (+.p16 (*.p16 d1 d2) (*.p16 d1 d3)))