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

↓

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

double f(double d1, double d2, double d3) {
        double r2214568 = d1;
        double r2214569 = d2;
        double r2214570 = r2214568 * r2214569;
        double r2214571 = d3;
        double r2214572 = r2214568 * r2214571;
        double r2214573 = r2214570 + r2214572;
        return r2214573;
}

↓

double f(double d1, double d2, double d3) {
        double r2214574 = d3;
        double r2214575 = d2;
        double r2214576 = r2214574 + r2214575;
        double r2214577 = d1;
        double r2214578 = r2214576 * r2214577;
        return r2214578;
}

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