Average Error: 58.1 → 58.1
Time: 3.5s
Precision: 64
\[\left(\left(333.75 \cdot {33096}^{6} + \left(77617 \cdot 77617\right) \cdot \left(\left(\left(\left(11 \cdot \left(77617 \cdot 77617\right)\right) \cdot \left(33096 \cdot 33096\right) + \left(-{33096}^{6}\right)\right) + -121 \cdot {33096}^{4}\right) + -2\right)\right) + 5.5 \cdot {33096}^{8}\right) + \frac{77617}{2 \cdot 33096}\]
\[\sqrt[3]{{\left(\mathsf{fma}\left(77617, 77617 \cdot \left(\left(11 \cdot \left(77617 \cdot 77617\right)\right) \cdot \left(33096 \cdot 33096\right) - \left({33096}^{6} - \mathsf{fma}\left({33096}^{4}, -121, -2\right)\right)\right), \mathsf{fma}\left(333.75, {33096}^{6}, \mathsf{fma}\left({33096}^{8}, 5.5, \frac{77617}{2 \cdot 33096}\right)\right)\right)\right)}^{3}}\]
\left(\left(333.75 \cdot {33096}^{6} + \left(77617 \cdot 77617\right) \cdot \left(\left(\left(\left(11 \cdot \left(77617 \cdot 77617\right)\right) \cdot \left(33096 \cdot 33096\right) + \left(-{33096}^{6}\right)\right) + -121 \cdot {33096}^{4}\right) + -2\right)\right) + 5.5 \cdot {33096}^{8}\right) + \frac{77617}{2 \cdot 33096}
\sqrt[3]{{\left(\mathsf{fma}\left(77617, 77617 \cdot \left(\left(11 \cdot \left(77617 \cdot 77617\right)\right) \cdot \left(33096 \cdot 33096\right) - \left({33096}^{6} - \mathsf{fma}\left({33096}^{4}, -121, -2\right)\right)\right), \mathsf{fma}\left(333.75, {33096}^{6}, \mathsf{fma}\left({33096}^{8}, 5.5, \frac{77617}{2 \cdot 33096}\right)\right)\right)\right)}^{3}}
double f() {
        double r70580 = 333.75;
        double r70581 = 33096.0;
        double r70582 = 6.0;
        double r70583 = pow(r70581, r70582);
        double r70584 = r70580 * r70583;
        double r70585 = 77617.0;
        double r70586 = r70585 * r70585;
        double r70587 = 11.0;
        double r70588 = r70587 * r70586;
        double r70589 = r70581 * r70581;
        double r70590 = r70588 * r70589;
        double r70591 = -r70583;
        double r70592 = r70590 + r70591;
        double r70593 = -121.0;
        double r70594 = 4.0;
        double r70595 = pow(r70581, r70594);
        double r70596 = r70593 * r70595;
        double r70597 = r70592 + r70596;
        double r70598 = -2.0;
        double r70599 = r70597 + r70598;
        double r70600 = r70586 * r70599;
        double r70601 = r70584 + r70600;
        double r70602 = 5.5;
        double r70603 = 8.0;
        double r70604 = pow(r70581, r70603);
        double r70605 = r70602 * r70604;
        double r70606 = r70601 + r70605;
        double r70607 = 2.0;
        double r70608 = r70607 * r70581;
        double r70609 = r70585 / r70608;
        double r70610 = r70606 + r70609;
        return r70610;
}

double f() {
        double r70611 = 77617.0;
        double r70612 = 11.0;
        double r70613 = r70611 * r70611;
        double r70614 = r70612 * r70613;
        double r70615 = 33096.0;
        double r70616 = r70615 * r70615;
        double r70617 = r70614 * r70616;
        double r70618 = 6.0;
        double r70619 = pow(r70615, r70618);
        double r70620 = 4.0;
        double r70621 = pow(r70615, r70620);
        double r70622 = -121.0;
        double r70623 = -2.0;
        double r70624 = fma(r70621, r70622, r70623);
        double r70625 = r70619 - r70624;
        double r70626 = r70617 - r70625;
        double r70627 = r70611 * r70626;
        double r70628 = 333.75;
        double r70629 = 8.0;
        double r70630 = pow(r70615, r70629);
        double r70631 = 5.5;
        double r70632 = 2.0;
        double r70633 = r70632 * r70615;
        double r70634 = r70611 / r70633;
        double r70635 = fma(r70630, r70631, r70634);
        double r70636 = fma(r70628, r70619, r70635);
        double r70637 = fma(r70611, r70627, r70636);
        double r70638 = 3.0;
        double r70639 = pow(r70637, r70638);
        double r70640 = cbrt(r70639);
        return r70640;
}

Error

Derivation

  1. Initial program 58.1

    \[\left(\left(333.75 \cdot {33096}^{6} + \left(77617 \cdot 77617\right) \cdot \left(\left(\left(\left(11 \cdot \left(77617 \cdot 77617\right)\right) \cdot \left(33096 \cdot 33096\right) + \left(-{33096}^{6}\right)\right) + -121 \cdot {33096}^{4}\right) + -2\right)\right) + 5.5 \cdot {33096}^{8}\right) + \frac{77617}{2 \cdot 33096}\]
  2. Simplified58.1

    \[\leadsto \color{blue}{\mathsf{fma}\left(77617 \cdot 77617, \left(11 \cdot \left(77617 \cdot 77617\right)\right) \cdot \left(33096 \cdot 33096\right) - \left({33096}^{6} - \mathsf{fma}\left({33096}^{4}, -121, -2\right)\right), 333.75 \cdot {33096}^{6} + \mathsf{fma}\left({33096}^{8}, 5.5, \frac{77617}{2 \cdot 33096}\right)\right)}\]
  3. Using strategy rm
  4. Applied add-cbrt-cube58.1

    \[\leadsto \color{blue}{\sqrt[3]{\left(\mathsf{fma}\left(77617 \cdot 77617, \left(11 \cdot \left(77617 \cdot 77617\right)\right) \cdot \left(33096 \cdot 33096\right) - \left({33096}^{6} - \mathsf{fma}\left({33096}^{4}, -121, -2\right)\right), 333.75 \cdot {33096}^{6} + \mathsf{fma}\left({33096}^{8}, 5.5, \frac{77617}{2 \cdot 33096}\right)\right) \cdot \mathsf{fma}\left(77617 \cdot 77617, \left(11 \cdot \left(77617 \cdot 77617\right)\right) \cdot \left(33096 \cdot 33096\right) - \left({33096}^{6} - \mathsf{fma}\left({33096}^{4}, -121, -2\right)\right), 333.75 \cdot {33096}^{6} + \mathsf{fma}\left({33096}^{8}, 5.5, \frac{77617}{2 \cdot 33096}\right)\right)\right) \cdot \mathsf{fma}\left(77617 \cdot 77617, \left(11 \cdot \left(77617 \cdot 77617\right)\right) \cdot \left(33096 \cdot 33096\right) - \left({33096}^{6} - \mathsf{fma}\left({33096}^{4}, -121, -2\right)\right), 333.75 \cdot {33096}^{6} + \mathsf{fma}\left({33096}^{8}, 5.5, \frac{77617}{2 \cdot 33096}\right)\right)}}\]
  5. Simplified58.1

    \[\leadsto \sqrt[3]{\color{blue}{{\left(\mathsf{fma}\left(77617, 77617 \cdot \left(\left(11 \cdot \left(77617 \cdot 77617\right)\right) \cdot \left(33096 \cdot 33096\right) - \left({33096}^{6} - \mathsf{fma}\left({33096}^{4}, -121, -2\right)\right)\right), \mathsf{fma}\left(333.75, {33096}^{6}, \mathsf{fma}\left({33096}^{8}, 5.5, \frac{77617}{2 \cdot 33096}\right)\right)\right)\right)}^{3}}}\]
  6. Final simplification58.1

    \[\leadsto \sqrt[3]{{\left(\mathsf{fma}\left(77617, 77617 \cdot \left(\left(11 \cdot \left(77617 \cdot 77617\right)\right) \cdot \left(33096 \cdot 33096\right) - \left({33096}^{6} - \mathsf{fma}\left({33096}^{4}, -121, -2\right)\right)\right), \mathsf{fma}\left(333.75, {33096}^{6}, \mathsf{fma}\left({33096}^{8}, 5.5, \frac{77617}{2 \cdot 33096}\right)\right)\right)\right)}^{3}}\]

Reproduce

herbie shell --seed 2020083 +o rules:numerics
(FPCore ()
  :name "From Warwick Tucker's Validated Numerics"
  :precision binary64
  (+ (+ (+ (* 333.75 (pow 33096 6)) (* (* 77617 77617) (+ (+ (+ (* (* 11 (* 77617 77617)) (* 33096 33096)) (- (pow 33096 6))) (* -121 (pow 33096 4))) -2))) (* 5.5 (pow 33096 8))) (/ 77617 (* 2 33096))))