Average Error: 0.1 → 0.1
Time: 21.1s
Precision: 64
\[0.0 \lt m \land 0.0 \lt v \land v \lt 0.25\]
\[\left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right) \cdot \left(1 - m\right)\]
↓
\[1 \cdot \left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right) + \left(\mathsf{fma}\left(1, m, m \cdot \frac{m \cdot m}{v}\right) - 1 \cdot \frac{{m}^{2}}{v}\right)\]
\left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right) \cdot \left(1 - m\right)
↓
1 \cdot \left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right) + \left(\mathsf{fma}\left(1, m, m \cdot \frac{m \cdot m}{v}\right) - 1 \cdot \frac{{m}^{2}}{v}\right)
double f(double m, double v) {
        double r23605 = m;
        double r23606 = 1.0;
        double r23607 = r23606 - r23605;
        double r23608 = r23605 * r23607;
        double r23609 = v;
        double r23610 = r23608 / r23609;
        double r23611 = r23610 - r23606;
        double r23612 = r23611 * r23607;
        return r23612;
}
↓
double f(double m, double v) {
        double r23613 = 1.0;
        double r23614 = m;
        double r23615 = r23613 - r23614;
        double r23616 = r23614 * r23615;
        double r23617 = v;
        double r23618 = r23616 / r23617;
        double r23619 = r23618 - r23613;
        double r23620 = r23613 * r23619;
        double r23621 = r23614 * r23614;
        double r23622 = r23621 / r23617;
        double r23623 = r23614 * r23622;
        double r23624 = fma(r23613, r23614, r23623);
        double r23625 = 2.0;
        double r23626 = pow(r23614, r23625);
        double r23627 = r23626 / r23617;
        double r23628 = r23613 * r23627;
        double r23629 = r23624 - r23628;
        double r23630 = r23620 + r23629;
        return r23630;
}

Error

Bits error versus m

Bits error versus v

Derivation

  1. Initial program 0.1

    \[\left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right) \cdot \left(1 - m\right)\]
  2. Using strategy rm

  3. Applied sub-neg0.1

    \[\leadsto \left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right) \cdot \color{blue}{\left(1 + \left(-m\right)\right)}\]

  4. Applied distribute-lft-in0.1

    \[\leadsto \color{blue}{\left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right) \cdot 1 + \left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right) \cdot \left(-m\right)}\]

  5. Simplified0.1

    \[\leadsto \color{blue}{1 \cdot \left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right)} + \left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right) \cdot \left(-m\right)\]

  6. Simplified0.1

    \[\leadsto 1 \cdot \left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right) + \color{blue}{\left(-m\right) \cdot \left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right)}\]

  7. Taylor expanded around 0 0.1

    \[\leadsto 1 \cdot \left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right) + \color{blue}{\left(\left(1 \cdot m + \frac{{m}^{3}}{v}\right) - 1 \cdot \frac{{m}^{2}}{v}\right)}\]

  8. Simplified0.1

    \[\leadsto 1 \cdot \left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right) + \color{blue}{\left(\mathsf{fma}\left(1, m, \frac{{m}^{3}}{v}\right) - 1 \cdot \frac{{m}^{2}}{v}\right)}\]
  9. Using strategy rm

  10. Applied *-un-lft-identity0.1

    \[\leadsto 1 \cdot \left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right) + \left(\mathsf{fma}\left(1, m, \frac{{m}^{3}}{\color{blue}{1 \cdot v}}\right) - 1 \cdot \frac{{m}^{2}}{v}\right)\]

  11. Applied cube-mult0.1

    \[\leadsto 1 \cdot \left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right) + \left(\mathsf{fma}\left(1, m, \frac{\color{blue}{m \cdot \left(m \cdot m\right)}}{1 \cdot v}\right) - 1 \cdot \frac{{m}^{2}}{v}\right)\]

  12. Applied times-frac0.1

    \[\leadsto 1 \cdot \left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right) + \left(\mathsf{fma}\left(1, m, \color{blue}{\frac{m}{1} \cdot \frac{m \cdot m}{v}}\right) - 1 \cdot \frac{{m}^{2}}{v}\right)\]

  13. Simplified0.1

    \[\leadsto 1 \cdot \left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right) + \left(\mathsf{fma}\left(1, m, \color{blue}{m} \cdot \frac{m \cdot m}{v}\right) - 1 \cdot \frac{{m}^{2}}{v}\right)\]

  14. Final simplification0.1

    \[\leadsto 1 \cdot \left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right) + \left(\mathsf{fma}\left(1, m, m \cdot \frac{m \cdot m}{v}\right) - 1 \cdot \frac{{m}^{2}}{v}\right)\]

Reproduce

herbie shell --seed 2019326 +o rules:numerics
(FPCore (m v)
  :name "b parameter of renormalized beta distribution"
  :precision binary64
  :pre (and (< 0.0 m) (< 0.0 v) (< v 0.25))
  (* (- (/ (* m (- 1 m)) v) 1) (- 1 m)))