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

↓

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

double f(double d1, double d2, double d3) {
        double r785736 = d1;
        double r785737 = d2;
        double r785738 = r785736 * r785737;
        double r785739 = d3;
        double r785740 = r785736 * r785739;
        double r785741 = r785738 + r785740;
        return r785741;
}

↓

double f(double d1, double d2, double d3) {
        double r785742 = d3;
        double r785743 = d2;
        double r785744 = r785742 + r785743;
        double r785745 = d1;
        double r785746 = r785744 * r785745;
        return r785746;
}

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 +o rules:numerics
(FPCore (d1 d2 d3)
  :name "FastMath dist"
  (+.p16 (*.p16 d1 d2) (*.p16 d1 d3)))