#include #include #include #include #include char *name = "2cbrt (problem 3.3.4)"; double f_if(float x) { float r9421 = x; float r9422 = 1.0f; float r9423 = r9421 + r9422; float r9424 = cbrt(r9423); float r9425 = cbrt(r9421); float r9426 = r9424 - r9425; return r9426; } double f_id(double x) { double r9427 = x; double r9428 = 1.0; double r9429 = r9427 + r9428; double r9430 = cbrt(r9429); double r9431 = cbrt(r9427); double r9432 = r9430 - r9431; return r9432; } double f_of(float x) { float r9433 = 1.0f; float r9434 = x; float r9435 = r9433 + r9434; float r9436 = cbrt(r9435); float r9437 = cbrt(r9434); float r9438 = r9436 - r9437; float r9439 = 0.0003583561724968806f; bool r9440 = r9438 <= r9439; float r9441 = r9434 * r9434; float r9442 = r9437 / r9441; float r9443 = 0.06172839506172839f; float r9444 = r9443 / r9434; float r9445 = 0.1111111111111111f; float r9446 = r9444 - r9445; float r9447 = 0.3333333333333333f; float r9448 = r9447 / r9434; float r9449 = r9437 * r9448; float r9450 = fma(r9442, r9446, r9449); float r9451 = cbrt(r9437); float r9452 = r9451 * r9451; float r9453 = r9451 * r9452; float r9454 = r9436 - r9453; float r9455 = r9440 ? r9450 : r9454; return r9455; } double f_od(double x) { double r9456 = 1.0; double r9457 = x; double r9458 = r9456 + r9457; double r9459 = cbrt(r9458); double r9460 = cbrt(r9457); double r9461 = r9459 - r9460; double r9462 = 0.0003583561724968806; bool r9463 = r9461 <= r9462; double r9464 = r9457 * r9457; double r9465 = r9460 / r9464; double r9466 = 0.06172839506172839; double r9467 = r9466 / r9457; double r9468 = 0.1111111111111111; double r9469 = r9467 - r9468; double r9470 = 0.3333333333333333; double r9471 = r9470 / r9457; double r9472 = r9460 * r9471; double r9473 = fma(r9465, r9469, r9472); double r9474 = cbrt(r9460); double r9475 = r9474 * r9474; double r9476 = r9474 * r9475; double r9477 = r9459 - r9476; double r9478 = r9463 ? r9473 : r9477; return r9478; } void mpfr_fmod2(mpfr_t r, mpfr_t n, mpfr_t d, mpfr_rnd_t rmd) { mpfr_fmod(r, n, d, rmd); if (mpfr_cmp_ui(r, 0) < 0) mpfr_add(r, r, d, rmd); } static mpfr_t r9479, r9480, r9481, r9482, r9483, r9484; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r9479); mpfr_init_set_str(r9480, "1", 10, MPFR_RNDN); mpfr_init(r9481); mpfr_init(r9482); mpfr_init(r9483); mpfr_init(r9484); } double f_im(double x) { mpfr_set_d(r9479, x, MPFR_RNDN); ; mpfr_add(r9481, r9479, r9480, MPFR_RNDN); mpfr_cbrt(r9482, r9481, MPFR_RNDN); mpfr_cbrt(r9483, r9479, MPFR_RNDN); mpfr_sub(r9484, r9482, r9483, MPFR_RNDN); return mpfr_get_d(r9484, MPFR_RNDN); } static mpfr_t r9485, r9486, r9487, r9488, r9489, r9490, r9491, r9492, r9493, r9494, r9495, r9496, r9497, r9498, r9499, r9500, r9501, r9502, r9503, r9504, r9505, r9506, r9507; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init_set_str(r9485, "1", 10, MPFR_RNDN); mpfr_init(r9486); mpfr_init(r9487); mpfr_init(r9488); mpfr_init(r9489); mpfr_init(r9490); mpfr_init_set_str(r9491, "0.0003583561724968806", 10, MPFR_RNDN); mpfr_init(r9492); mpfr_init(r9493); mpfr_init(r9494); mpfr_init_set_str(r9495, "5/81", 10, MPFR_RNDN); mpfr_init(r9496); mpfr_init_set_str(r9497, "1/9", 10, MPFR_RNDN); mpfr_init(r9498); mpfr_init_set_str(r9499, "1/3", 10, MPFR_RNDN); mpfr_init(r9500); mpfr_init(r9501); mpfr_init(r9502); mpfr_init(r9503); mpfr_init(r9504); mpfr_init(r9505); mpfr_init(r9506); mpfr_init(r9507); } double f_fm(double x) { ; mpfr_set_d(r9486, x, MPFR_RNDN); mpfr_add(r9487, r9485, r9486, MPFR_RNDN); mpfr_cbrt(r9488, r9487, MPFR_RNDN); mpfr_cbrt(r9489, r9486, MPFR_RNDN); mpfr_sub(r9490, r9488, r9489, MPFR_RNDN); ; mpfr_set_si(r9492, mpfr_cmp(r9490, r9491) <= 0, MPFR_RNDN); mpfr_mul(r9493, r9486, r9486, MPFR_RNDN); mpfr_div(r9494, r9489, r9493, MPFR_RNDN); ; mpfr_div(r9496, r9495, r9486, MPFR_RNDN); ; mpfr_sub(r9498, r9496, r9497, MPFR_RNDN); ; mpfr_div(r9500, r9499, r9486, MPFR_RNDN); mpfr_mul(r9501, r9489, r9500, MPFR_RNDN); mpfr_fma(r9502, r9494, r9498, r9501, MPFR_RNDN); mpfr_cbrt(r9503, r9489, MPFR_RNDN); mpfr_mul(r9504, r9503, r9503, MPFR_RNDN); mpfr_mul(r9505, r9503, r9504, MPFR_RNDN); mpfr_sub(r9506, r9488, r9505, MPFR_RNDN); if (mpfr_get_si(r9492, MPFR_RNDN)) { mpfr_set(r9507, r9502, MPFR_RNDN); } else { mpfr_set(r9507, r9506, MPFR_RNDN); }; return mpfr_get_d(r9507, MPFR_RNDN); } static mpfr_t r9508, r9509, r9510, r9511, r9512, r9513, r9514, r9515, r9516, r9517, r9518, r9519, r9520, r9521, r9522, r9523, r9524, r9525, r9526, r9527, r9528, r9529, r9530; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init_set_str(r9508, "1", 10, MPFR_RNDN); mpfr_init(r9509); mpfr_init(r9510); mpfr_init(r9511); mpfr_init(r9512); mpfr_init(r9513); mpfr_init_set_str(r9514, "0.0003583561724968806", 10, MPFR_RNDN); mpfr_init(r9515); mpfr_init(r9516); mpfr_init(r9517); mpfr_init_set_str(r9518, "5/81", 10, MPFR_RNDN); mpfr_init(r9519); mpfr_init_set_str(r9520, "1/9", 10, MPFR_RNDN); mpfr_init(r9521); mpfr_init_set_str(r9522, "1/3", 10, MPFR_RNDN); mpfr_init(r9523); mpfr_init(r9524); mpfr_init(r9525); mpfr_init(r9526); mpfr_init(r9527); mpfr_init(r9528); mpfr_init(r9529); mpfr_init(r9530); } double f_dm(double x) { ; mpfr_set_d(r9509, x, MPFR_RNDN); mpfr_add(r9510, r9508, r9509, MPFR_RNDN); mpfr_cbrt(r9511, r9510, MPFR_RNDN); mpfr_cbrt(r9512, r9509, MPFR_RNDN); mpfr_sub(r9513, r9511, r9512, MPFR_RNDN); ; mpfr_set_si(r9515, mpfr_cmp(r9513, r9514) <= 0, MPFR_RNDN); mpfr_mul(r9516, r9509, r9509, MPFR_RNDN); mpfr_div(r9517, r9512, r9516, MPFR_RNDN); ; mpfr_div(r9519, r9518, r9509, MPFR_RNDN); ; mpfr_sub(r9521, r9519, r9520, MPFR_RNDN); ; mpfr_div(r9523, r9522, r9509, MPFR_RNDN); mpfr_mul(r9524, r9512, r9523, MPFR_RNDN); mpfr_fma(r9525, r9517, r9521, r9524, MPFR_RNDN); mpfr_cbrt(r9526, r9512, MPFR_RNDN); mpfr_mul(r9527, r9526, r9526, MPFR_RNDN); mpfr_mul(r9528, r9526, r9527, MPFR_RNDN); mpfr_sub(r9529, r9511, r9528, MPFR_RNDN); if (mpfr_get_si(r9515, MPFR_RNDN)) { mpfr_set(r9530, r9525, MPFR_RNDN); } else { mpfr_set(r9530, r9529, MPFR_RNDN); }; return mpfr_get_d(r9530, MPFR_RNDN); }