\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}\frac{\left(\left(\mathsf{fma}\left(5.5, {33096}^{8}, \mathsf{fma}\left(77617 \cdot \mathsf{fma}\left({33096}^{4}, -121, \left(11 \cdot \left(\left(33096 \cdot 77617\right) \cdot \left(33096 \cdot 77617\right)\right) - {33096}^{6}\right) + -2\right), 77617, {33096}^{6} \cdot 333.75\right)\right)\right)\right) \cdot \left(\left(\mathsf{fma}\left(5.5, {33096}^{8}, \mathsf{fma}\left(77617 \cdot \mathsf{fma}\left({33096}^{4}, -121, \left(11 \cdot \left(\left(33096 \cdot 77617\right) \cdot \left(33096 \cdot 77617\right)\right) - {33096}^{6}\right) + -2\right), 77617, {33096}^{6} \cdot 333.75\right)\right)\right)\right) - \frac{77617}{2 \cdot 33096} \cdot \frac{77617}{2 \cdot 33096}}{\left(\left(\mathsf{fma}\left(5.5, {33096}^{8}, \mathsf{fma}\left(77617 \cdot \mathsf{fma}\left({33096}^{4}, -121, \left(11 \cdot \left(\left(33096 \cdot 77617\right) \cdot \left(33096 \cdot 77617\right)\right) - {33096}^{6}\right) + -2\right), 77617, {33096}^{6} \cdot 333.75\right)\right)\right)\right) - \frac{77617}{2 \cdot 33096}}double f() {
double r3251429 = 333.75;
double r3251430 = 33096.0;
double r3251431 = 6.0;
double r3251432 = pow(r3251430, r3251431);
double r3251433 = r3251429 * r3251432;
double r3251434 = 77617.0;
double r3251435 = r3251434 * r3251434;
double r3251436 = 11.0;
double r3251437 = r3251436 * r3251435;
double r3251438 = r3251430 * r3251430;
double r3251439 = r3251437 * r3251438;
double r3251440 = -r3251432;
double r3251441 = r3251439 + r3251440;
double r3251442 = -121.0;
double r3251443 = 4.0;
double r3251444 = pow(r3251430, r3251443);
double r3251445 = r3251442 * r3251444;
double r3251446 = r3251441 + r3251445;
double r3251447 = -2.0;
double r3251448 = r3251446 + r3251447;
double r3251449 = r3251435 * r3251448;
double r3251450 = r3251433 + r3251449;
double r3251451 = 5.5;
double r3251452 = 8.0;
double r3251453 = pow(r3251430, r3251452);
double r3251454 = r3251451 * r3251453;
double r3251455 = r3251450 + r3251454;
double r3251456 = 2.0;
double r3251457 = r3251456 * r3251430;
double r3251458 = r3251434 / r3251457;
double r3251459 = r3251455 + r3251458;
return r3251459;
}
double f() {
double r3251460 = 5.5;
double r3251461 = 33096.0;
double r3251462 = 8.0;
double r3251463 = pow(r3251461, r3251462);
double r3251464 = 77617.0;
double r3251465 = 4.0;
double r3251466 = pow(r3251461, r3251465);
double r3251467 = -121.0;
double r3251468 = 11.0;
double r3251469 = r3251461 * r3251464;
double r3251470 = r3251469 * r3251469;
double r3251471 = r3251468 * r3251470;
double r3251472 = 6.0;
double r3251473 = pow(r3251461, r3251472);
double r3251474 = r3251471 - r3251473;
double r3251475 = -2.0;
double r3251476 = r3251474 + r3251475;
double r3251477 = fma(r3251466, r3251467, r3251476);
double r3251478 = r3251464 * r3251477;
double r3251479 = 333.75;
double r3251480 = r3251473 * r3251479;
double r3251481 = fma(r3251478, r3251464, r3251480);
double r3251482 = fma(r3251460, r3251463, r3251481);
double r3251483 = /* ERROR: no posit support in C */;
double r3251484 = /* ERROR: no posit support in C */;
double r3251485 = r3251484 * r3251484;
double r3251486 = 2.0;
double r3251487 = r3251486 * r3251461;
double r3251488 = r3251464 / r3251487;
double r3251489 = r3251488 * r3251488;
double r3251490 = r3251485 - r3251489;
double r3251491 = r3251484 - r3251488;
double r3251492 = r3251490 / r3251491;
return r3251492;
}
Initial program 58.1
rmApplied insert-posit1656.8
Simplified56.8
rmApplied flip-+56.8
Final simplification56.8
herbie shell --seed 2019168 +o rules:numerics
(FPCore ()
:name "From Warwick Tucker's Validated Numerics"
(+ (+ (+ (* 333.75 (pow 33096.0 6.0)) (* (* 77617.0 77617.0) (+ (+ (+ (* (* 11.0 (* 77617.0 77617.0)) (* 33096.0 33096.0)) (- (pow 33096.0 6.0))) (* -121.0 (pow 33096.0 4.0))) -2.0))) (* 5.5 (pow 33096.0 8.0))) (/ 77617.0 (* 2.0 33096.0))))