#include #include #include #include #include char *name = "Octave 3.8, jcobi/4, as called"; double f_if(float i) { float r25599 = i; float r25600 = r25599 * r25599; float r25601 = r25600 * r25600; float r25602 = 2; float r25603 = r25602 * r25599; float r25604 = r25603 * r25603; float r25605 = r25601 / r25604; float r25606 = 1.0; float r25607 = r25604 - r25606; float r25608 = r25605 / r25607; return r25608; } double f_id(double i) { double r25609 = i; double r25610 = r25609 * r25609; double r25611 = r25610 * r25610; double r25612 = 2; double r25613 = r25612 * r25609; double r25614 = r25613 * r25613; double r25615 = r25611 / r25614; double r25616 = 1.0; double r25617 = r25614 - r25616; double r25618 = r25615 / r25617; return r25618; } double f_of(float i) { float r25619 = i; float r25620 = 226.15481187430868; bool r25621 = r25619 <= r25620; float r25622 = 2; float r25623 = r25619 / r25622; float r25624 = r25623 * r25623; float r25625 = r25622 * r25619; float r25626 = 1.0; float r25627 = -r25626; float r25628 = fma(r25625, r25625, r25627); float r25629 = r25624 / r25628; float r25630 = expm1(r25629); float r25631 = log1p(r25630); float r25632 = 0.00390625; float r25633 = 4; float r25634 = pow(r25619, r25633); float r25635 = r25632 / r25634; float r25636 = 0.015625; float r25637 = r25636 / r25619; float r25638 = r25637 / r25619; float r25639 = r25635 + r25638; float r25640 = 1/16; float r25641 = r25639 + r25640; float r25642 = r25621 ? r25631 : r25641; return r25642; } double f_od(double i) { double r25643 = i; double r25644 = 226.15481187430868; bool r25645 = r25643 <= r25644; double r25646 = 2; double r25647 = r25643 / r25646; double r25648 = r25647 * r25647; double r25649 = r25646 * r25643; double r25650 = 1.0; double r25651 = -r25650; double r25652 = fma(r25649, r25649, r25651); double r25653 = r25648 / r25652; double r25654 = expm1(r25653); double r25655 = log1p(r25654); double r25656 = 0.00390625; double r25657 = 4; double r25658 = pow(r25643, r25657); double r25659 = r25656 / r25658; double r25660 = 0.015625; double r25661 = r25660 / r25643; double r25662 = r25661 / r25643; double r25663 = r25659 + r25662; double r25664 = 1/16; double r25665 = r25663 + r25664; double r25666 = r25645 ? r25655 : r25665; return r25666; } 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 r25667, r25668, r25669, r25670, r25671, r25672, r25673, r25674, r25675, r25676; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init(r25667); mpfr_init(r25668); mpfr_init(r25669); mpfr_init_set_str(r25670, "2", 10, MPFR_RNDN); mpfr_init(r25671); mpfr_init(r25672); mpfr_init(r25673); mpfr_init_set_str(r25674, "1.0", 10, MPFR_RNDN); mpfr_init(r25675); mpfr_init(r25676); } double f_im(double i) { mpfr_set_d(r25667, i, MPFR_RNDN); mpfr_mul(r25668, r25667, r25667, MPFR_RNDN); mpfr_mul(r25669, r25668, r25668, MPFR_RNDN); ; mpfr_mul(r25671, r25670, r25667, MPFR_RNDN); mpfr_mul(r25672, r25671, r25671, MPFR_RNDN); mpfr_div(r25673, r25669, r25672, MPFR_RNDN); ; mpfr_sub(r25675, r25672, r25674, MPFR_RNDN); mpfr_div(r25676, r25673, r25675, MPFR_RNDN); return mpfr_get_d(r25676, MPFR_RNDN); } static mpfr_t r25677, r25678, r25679, r25680, r25681, r25682, r25683, r25684, r25685, r25686, r25687, r25688, r25689, r25690, r25691, r25692, r25693, r25694, r25695, r25696, r25697, r25698, r25699, r25700; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init(r25677); mpfr_init_set_str(r25678, "226.15481187430868", 10, MPFR_RNDN); mpfr_init(r25679); mpfr_init_set_str(r25680, "2", 10, MPFR_RNDN); mpfr_init(r25681); mpfr_init(r25682); mpfr_init(r25683); mpfr_init_set_str(r25684, "1.0", 10, MPFR_RNDN); mpfr_init(r25685); mpfr_init(r25686); mpfr_init(r25687); mpfr_init(r25688); mpfr_init(r25689); mpfr_init_set_str(r25690, "0.00390625", 10, MPFR_RNDN); mpfr_init_set_str(r25691, "4", 10, MPFR_RNDN); mpfr_init(r25692); mpfr_init(r25693); mpfr_init_set_str(r25694, "0.015625", 10, MPFR_RNDN); mpfr_init(r25695); mpfr_init(r25696); mpfr_init(r25697); mpfr_init_set_str(r25698, "1/16", 10, MPFR_RNDN); mpfr_init(r25699); mpfr_init(r25700); } double f_fm(double i) { mpfr_set_d(r25677, i, MPFR_RNDN); ; mpfr_set_si(r25679, mpfr_cmp(r25677, r25678) <= 0, MPFR_RNDN); ; mpfr_div(r25681, r25677, r25680, MPFR_RNDN); mpfr_mul(r25682, r25681, r25681, MPFR_RNDN); mpfr_mul(r25683, r25680, r25677, MPFR_RNDN); ; mpfr_neg(r25685, r25684, MPFR_RNDN); mpfr_fma(r25686, r25683, r25683, r25685, MPFR_RNDN); mpfr_div(r25687, r25682, r25686, MPFR_RNDN); mpfr_expm1(r25688, r25687, MPFR_RNDN); mpfr_log1p(r25689, r25688, MPFR_RNDN); ; ; mpfr_pow(r25692, r25677, r25691, MPFR_RNDN); mpfr_div(r25693, r25690, r25692, MPFR_RNDN); ; mpfr_div(r25695, r25694, r25677, MPFR_RNDN); mpfr_div(r25696, r25695, r25677, MPFR_RNDN); mpfr_add(r25697, r25693, r25696, MPFR_RNDN); ; mpfr_add(r25699, r25697, r25698, MPFR_RNDN); if (mpfr_get_si(r25679, MPFR_RNDN)) { mpfr_set(r25700, r25689, MPFR_RNDN); } else { mpfr_set(r25700, r25699, MPFR_RNDN); }; return mpfr_get_d(r25700, MPFR_RNDN); } static mpfr_t r25701, r25702, r25703, r25704, r25705, r25706, r25707, r25708, r25709, r25710, r25711, r25712, r25713, r25714, r25715, r25716, r25717, r25718, r25719, r25720, r25721, r25722, r25723, r25724; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init(r25701); mpfr_init_set_str(r25702, "226.15481187430868", 10, MPFR_RNDN); mpfr_init(r25703); mpfr_init_set_str(r25704, "2", 10, MPFR_RNDN); mpfr_init(r25705); mpfr_init(r25706); mpfr_init(r25707); mpfr_init_set_str(r25708, "1.0", 10, MPFR_RNDN); mpfr_init(r25709); mpfr_init(r25710); mpfr_init(r25711); mpfr_init(r25712); mpfr_init(r25713); mpfr_init_set_str(r25714, "0.00390625", 10, MPFR_RNDN); mpfr_init_set_str(r25715, "4", 10, MPFR_RNDN); mpfr_init(r25716); mpfr_init(r25717); mpfr_init_set_str(r25718, "0.015625", 10, MPFR_RNDN); mpfr_init(r25719); mpfr_init(r25720); mpfr_init(r25721); mpfr_init_set_str(r25722, "1/16", 10, MPFR_RNDN); mpfr_init(r25723); mpfr_init(r25724); } double f_dm(double i) { mpfr_set_d(r25701, i, MPFR_RNDN); ; mpfr_set_si(r25703, mpfr_cmp(r25701, r25702) <= 0, MPFR_RNDN); ; mpfr_div(r25705, r25701, r25704, MPFR_RNDN); mpfr_mul(r25706, r25705, r25705, MPFR_RNDN); mpfr_mul(r25707, r25704, r25701, MPFR_RNDN); ; mpfr_neg(r25709, r25708, MPFR_RNDN); mpfr_fma(r25710, r25707, r25707, r25709, MPFR_RNDN); mpfr_div(r25711, r25706, r25710, MPFR_RNDN); mpfr_expm1(r25712, r25711, MPFR_RNDN); mpfr_log1p(r25713, r25712, MPFR_RNDN); ; ; mpfr_pow(r25716, r25701, r25715, MPFR_RNDN); mpfr_div(r25717, r25714, r25716, MPFR_RNDN); ; mpfr_div(r25719, r25718, r25701, MPFR_RNDN); mpfr_div(r25720, r25719, r25701, MPFR_RNDN); mpfr_add(r25721, r25717, r25720, MPFR_RNDN); ; mpfr_add(r25723, r25721, r25722, MPFR_RNDN); if (mpfr_get_si(r25703, MPFR_RNDN)) { mpfr_set(r25724, r25713, MPFR_RNDN); } else { mpfr_set(r25724, r25723, MPFR_RNDN); }; return mpfr_get_d(r25724, MPFR_RNDN); }