#include #include #include #include #include char *name = "Numeric.SpecFunctions:invIncompleteBetaWorker from math-functions-0.1.5.2"; double f_if(float x, float y, float z, float t, float a, float b, float c) { float r1764600 = x; float r1764601 = y; float r1764602 = 2.0; float r1764603 = z; float r1764604 = t; float r1764605 = a; float r1764606 = r1764604 + r1764605; float r1764607 = sqrt(r1764606); float r1764608 = r1764603 * r1764607; float r1764609 = r1764608 / r1764604; float r1764610 = b; float r1764611 = c; float r1764612 = r1764610 - r1764611; float r1764613 = 5.0; float r1764614 = 6.0; float r1764615 = r1764613 / r1764614; float r1764616 = r1764605 + r1764615; float r1764617 = 3.0; float r1764618 = r1764604 * r1764617; float r1764619 = r1764602 / r1764618; float r1764620 = r1764616 - r1764619; float r1764621 = r1764612 * r1764620; float r1764622 = r1764609 - r1764621; float r1764623 = r1764602 * r1764622; float r1764624 = exp(r1764623); float r1764625 = r1764601 * r1764624; float r1764626 = r1764600 + r1764625; float r1764627 = r1764600 / r1764626; return r1764627; } double f_id(double x, double y, double z, double t, double a, double b, double c) { double r1764628 = x; double r1764629 = y; double r1764630 = 2.0; double r1764631 = z; double r1764632 = t; double r1764633 = a; double r1764634 = r1764632 + r1764633; double r1764635 = sqrt(r1764634); double r1764636 = r1764631 * r1764635; double r1764637 = r1764636 / r1764632; double r1764638 = b; double r1764639 = c; double r1764640 = r1764638 - r1764639; double r1764641 = 5.0; double r1764642 = 6.0; double r1764643 = r1764641 / r1764642; double r1764644 = r1764633 + r1764643; double r1764645 = 3.0; double r1764646 = r1764632 * r1764645; double r1764647 = r1764630 / r1764646; double r1764648 = r1764644 - r1764647; double r1764649 = r1764640 * r1764648; double r1764650 = r1764637 - r1764649; double r1764651 = r1764630 * r1764650; double r1764652 = exp(r1764651); double r1764653 = r1764629 * r1764652; double r1764654 = r1764628 + r1764653; double r1764655 = r1764628 / r1764654; return r1764655; } double f_of(float x, float y, float z, float t, float a, float b, float c) { float r1764656 = x; float r1764657 = y; float r1764658 = 2.0; float r1764659 = z; float r1764660 = t; float r1764661 = a; float r1764662 = r1764660 + r1764661; float r1764663 = sqrt(r1764662); float r1764664 = r1764659 * r1764663; float r1764665 = r1764664 / r1764660; float r1764666 = b; float r1764667 = c; float r1764668 = r1764666 - r1764667; float r1764669 = 5.0; float r1764670 = 6.0; float r1764671 = r1764669 / r1764670; float r1764672 = r1764661 + r1764671; float r1764673 = 3.0; float r1764674 = r1764660 * r1764673; float r1764675 = r1764658 / r1764674; float r1764676 = r1764672 - r1764675; float r1764677 = r1764668 * r1764676; float r1764678 = r1764665 - r1764677; float r1764679 = r1764658 * r1764678; float r1764680 = exp(r1764679); float r1764681 = r1764657 * r1764680; float r1764682 = r1764656 + r1764681; float r1764683 = r1764656 / r1764682; float r1764684 = r1764683 != r1764683; float r1764685 = 1; float r1764686 = exp(r1764658); float r1764687 = r1764667 - r1764666; float r1764688 = 0.8333333333333334; float r1764689 = r1764667 * r1764688; float r1764690 = fma(r1764687, r1764661, r1764689); float r1764691 = pow(r1764686, r1764690); float r1764692 = fma(r1764657, r1764691, r1764656); float r1764693 = r1764692 / r1764656; float r1764694 = r1764685 / r1764693; float r1764695 = +inf.0; bool r1764696 = r1764683 <= r1764695; float r1764697 = r1764659 / r1764660; float r1764698 = r1764697 * r1764663; float r1764699 = r1764671 + r1764661; float r1764700 = r1764658 / r1764660; float r1764701 = r1764700 / r1764673; float r1764702 = r1764699 - r1764701; float r1764703 = r1764702 * r1764668; float r1764704 = r1764698 - r1764703; float r1764705 = pow(r1764686, r1764704); float r1764706 = fma(r1764657, r1764705, r1764656); float r1764707 = r1764656 / r1764706; float r1764708 = r1764696 ? r1764707 : r1764707; float r1764709 = r1764684 ? r1764694 : r1764708; return r1764709; } double f_od(double x, double y, double z, double t, double a, double b, double c) { double r1764710 = x; double r1764711 = y; double r1764712 = 2.0; double r1764713 = z; double r1764714 = t; double r1764715 = a; double r1764716 = r1764714 + r1764715; double r1764717 = sqrt(r1764716); double r1764718 = r1764713 * r1764717; double r1764719 = r1764718 / r1764714; double r1764720 = b; double r1764721 = c; double r1764722 = r1764720 - r1764721; double r1764723 = 5.0; double r1764724 = 6.0; double r1764725 = r1764723 / r1764724; double r1764726 = r1764715 + r1764725; double r1764727 = 3.0; double r1764728 = r1764714 * r1764727; double r1764729 = r1764712 / r1764728; double r1764730 = r1764726 - r1764729; double r1764731 = r1764722 * r1764730; double r1764732 = r1764719 - r1764731; double r1764733 = r1764712 * r1764732; double r1764734 = exp(r1764733); double r1764735 = r1764711 * r1764734; double r1764736 = r1764710 + r1764735; double r1764737 = r1764710 / r1764736; double r1764738 = r1764737 != r1764737; double r1764739 = 1; double r1764740 = exp(r1764712); double r1764741 = r1764721 - r1764720; double r1764742 = 0.8333333333333334; double r1764743 = r1764721 * r1764742; double r1764744 = fma(r1764741, r1764715, r1764743); double r1764745 = pow(r1764740, r1764744); double r1764746 = fma(r1764711, r1764745, r1764710); double r1764747 = r1764746 / r1764710; double r1764748 = r1764739 / r1764747; double r1764749 = +inf.0; bool r1764750 = r1764737 <= r1764749; double r1764751 = r1764713 / r1764714; double r1764752 = r1764751 * r1764717; double r1764753 = r1764725 + r1764715; double r1764754 = r1764712 / r1764714; double r1764755 = r1764754 / r1764727; double r1764756 = r1764753 - r1764755; double r1764757 = r1764756 * r1764722; double r1764758 = r1764752 - r1764757; double r1764759 = pow(r1764740, r1764758); double r1764760 = fma(r1764711, r1764759, r1764710); double r1764761 = r1764710 / r1764760; double r1764762 = r1764750 ? r1764761 : r1764761; double r1764763 = r1764738 ? r1764748 : r1764762; return r1764763; } 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 r1764764, r1764765, r1764766, r1764767, r1764768, r1764769, r1764770, r1764771, r1764772, r1764773, r1764774, r1764775, r1764776, r1764777, r1764778, r1764779, r1764780, r1764781, r1764782, r1764783, r1764784, r1764785, r1764786, r1764787, r1764788, r1764789, r1764790, r1764791; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r1764764); mpfr_init(r1764765); mpfr_init_set_str(r1764766, "2.0", 10, MPFR_RNDN); mpfr_init(r1764767); mpfr_init(r1764768); mpfr_init(r1764769); mpfr_init(r1764770); mpfr_init(r1764771); mpfr_init(r1764772); mpfr_init(r1764773); mpfr_init(r1764774); mpfr_init(r1764775); mpfr_init(r1764776); mpfr_init_set_str(r1764777, "5.0", 10, MPFR_RNDN); mpfr_init_set_str(r1764778, "6.0", 10, MPFR_RNDN); mpfr_init(r1764779); mpfr_init(r1764780); mpfr_init_set_str(r1764781, "3.0", 10, MPFR_RNDN); mpfr_init(r1764782); mpfr_init(r1764783); mpfr_init(r1764784); mpfr_init(r1764785); mpfr_init(r1764786); mpfr_init(r1764787); mpfr_init(r1764788); mpfr_init(r1764789); mpfr_init(r1764790); mpfr_init(r1764791); } double f_im(double x, double y, double z, double t, double a, double b, double c) { mpfr_set_d(r1764764, x, MPFR_RNDN); mpfr_set_d(r1764765, y, MPFR_RNDN); ; mpfr_set_d(r1764767, z, MPFR_RNDN); mpfr_set_d(r1764768, t, MPFR_RNDN); mpfr_set_d(r1764769, a, MPFR_RNDN); mpfr_add(r1764770, r1764768, r1764769, MPFR_RNDN); mpfr_sqrt(r1764771, r1764770, MPFR_RNDN); mpfr_mul(r1764772, r1764767, r1764771, MPFR_RNDN); mpfr_div(r1764773, r1764772, r1764768, MPFR_RNDN); mpfr_set_d(r1764774, b, MPFR_RNDN); mpfr_set_d(r1764775, c, MPFR_RNDN); mpfr_sub(r1764776, r1764774, r1764775, MPFR_RNDN); ; ; mpfr_div(r1764779, r1764777, r1764778, MPFR_RNDN); mpfr_add(r1764780, r1764769, r1764779, MPFR_RNDN); ; mpfr_mul(r1764782, r1764768, r1764781, MPFR_RNDN); mpfr_div(r1764783, r1764766, r1764782, MPFR_RNDN); mpfr_sub(r1764784, r1764780, r1764783, MPFR_RNDN); mpfr_mul(r1764785, r1764776, r1764784, MPFR_RNDN); mpfr_sub(r1764786, r1764773, r1764785, MPFR_RNDN); mpfr_mul(r1764787, r1764766, r1764786, MPFR_RNDN); mpfr_exp(r1764788, r1764787, MPFR_RNDN); mpfr_mul(r1764789, r1764765, r1764788, MPFR_RNDN); mpfr_add(r1764790, r1764764, r1764789, MPFR_RNDN); mpfr_div(r1764791, r1764764, r1764790, MPFR_RNDN); return mpfr_get_d(r1764791, MPFR_RNDN); } static mpfr_t r1764792, r1764793, r1764794, r1764795, r1764796, r1764797, r1764798, r1764799, r1764800, r1764801, r1764802, r1764803, r1764804, r1764805, r1764806, r1764807, r1764808, r1764809, r1764810, r1764811, r1764812, r1764813, r1764814, r1764815, r1764816, r1764817, r1764818, r1764819, r1764820, r1764821, r1764822, r1764823, r1764824, r1764825, r1764826, r1764827, r1764828, r1764829, r1764830, r1764831, r1764832, r1764833, r1764834, r1764835, r1764836, r1764837, r1764838, r1764839, r1764840, r1764841, r1764842, r1764843, r1764844, r1764845; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r1764792); mpfr_init(r1764793); mpfr_init_set_str(r1764794, "2.0", 10, MPFR_RNDN); mpfr_init(r1764795); mpfr_init(r1764796); mpfr_init(r1764797); mpfr_init(r1764798); mpfr_init(r1764799); mpfr_init(r1764800); mpfr_init(r1764801); mpfr_init(r1764802); mpfr_init(r1764803); mpfr_init(r1764804); mpfr_init_set_str(r1764805, "5.0", 10, MPFR_RNDN); mpfr_init_set_str(r1764806, "6.0", 10, MPFR_RNDN); mpfr_init(r1764807); mpfr_init(r1764808); mpfr_init_set_str(r1764809, "3.0", 10, MPFR_RNDN); mpfr_init(r1764810); mpfr_init(r1764811); mpfr_init(r1764812); mpfr_init(r1764813); mpfr_init(r1764814); mpfr_init(r1764815); mpfr_init(r1764816); mpfr_init(r1764817); mpfr_init(r1764818); mpfr_init(r1764819); mpfr_init(r1764820); mpfr_init_set_str(r1764821, "1", 10, MPFR_RNDN); mpfr_init(r1764822); mpfr_init(r1764823); mpfr_init_set_str(r1764824, "0.8333333333333334", 10, MPFR_RNDN); mpfr_init(r1764825); mpfr_init(r1764826); mpfr_init(r1764827); mpfr_init(r1764828); mpfr_init(r1764829); mpfr_init(r1764830); mpfr_init_set_str(r1764831, "+inf.0", 10, MPFR_RNDN); mpfr_init(r1764832); mpfr_init(r1764833); mpfr_init(r1764834); mpfr_init(r1764835); mpfr_init(r1764836); mpfr_init(r1764837); mpfr_init(r1764838); mpfr_init(r1764839); mpfr_init(r1764840); mpfr_init(r1764841); mpfr_init(r1764842); mpfr_init(r1764843); mpfr_init(r1764844); mpfr_init(r1764845); } double f_fm(double x, double y, double z, double t, double a, double b, double c) { mpfr_set_d(r1764792, x, MPFR_RNDN); mpfr_set_d(r1764793, y, MPFR_RNDN); ; mpfr_set_d(r1764795, z, MPFR_RNDN); mpfr_set_d(r1764796, t, MPFR_RNDN); mpfr_set_d(r1764797, a, MPFR_RNDN); mpfr_add(r1764798, r1764796, r1764797, MPFR_RNDN); mpfr_sqrt(r1764799, r1764798, MPFR_RNDN); mpfr_mul(r1764800, r1764795, r1764799, MPFR_RNDN); mpfr_div(r1764801, r1764800, r1764796, MPFR_RNDN); mpfr_set_d(r1764802, b, MPFR_RNDN); mpfr_set_d(r1764803, c, MPFR_RNDN); mpfr_sub(r1764804, r1764802, r1764803, MPFR_RNDN); ; ; mpfr_div(r1764807, r1764805, r1764806, MPFR_RNDN); mpfr_add(r1764808, r1764797, r1764807, MPFR_RNDN); ; mpfr_mul(r1764810, r1764796, r1764809, MPFR_RNDN); mpfr_div(r1764811, r1764794, r1764810, MPFR_RNDN); mpfr_sub(r1764812, r1764808, r1764811, MPFR_RNDN); mpfr_mul(r1764813, r1764804, r1764812, MPFR_RNDN); mpfr_sub(r1764814, r1764801, r1764813, MPFR_RNDN); mpfr_mul(r1764815, r1764794, r1764814, MPFR_RNDN); mpfr_exp(r1764816, r1764815, MPFR_RNDN); mpfr_mul(r1764817, r1764793, r1764816, MPFR_RNDN); mpfr_add(r1764818, r1764792, r1764817, MPFR_RNDN); mpfr_div(r1764819, r1764792, r1764818, MPFR_RNDN); mpfr_set_si(r1764820, mpfr_cmp(r1764819, r1764819) != 0, MPFR_RNDN); ; mpfr_exp(r1764822, r1764794, MPFR_RNDN); mpfr_sub(r1764823, r1764803, r1764802, MPFR_RNDN); ; mpfr_mul(r1764825, r1764803, r1764824, MPFR_RNDN); mpfr_fma(r1764826, r1764823, r1764797, r1764825, MPFR_RNDN); mpfr_pow(r1764827, r1764822, r1764826, MPFR_RNDN); mpfr_fma(r1764828, r1764793, r1764827, r1764792, MPFR_RNDN); mpfr_div(r1764829, r1764828, r1764792, MPFR_RNDN); mpfr_div(r1764830, r1764821, r1764829, MPFR_RNDN); ; mpfr_set_si(r1764832, mpfr_cmp(r1764819, r1764831) <= 0, MPFR_RNDN); mpfr_div(r1764833, r1764795, r1764796, MPFR_RNDN); mpfr_mul(r1764834, r1764833, r1764799, MPFR_RNDN); mpfr_add(r1764835, r1764807, r1764797, MPFR_RNDN); mpfr_div(r1764836, r1764794, r1764796, MPFR_RNDN); mpfr_div(r1764837, r1764836, r1764809, MPFR_RNDN); mpfr_sub(r1764838, r1764835, r1764837, MPFR_RNDN); mpfr_mul(r1764839, r1764838, r1764804, MPFR_RNDN); mpfr_sub(r1764840, r1764834, r1764839, MPFR_RNDN); mpfr_pow(r1764841, r1764822, r1764840, MPFR_RNDN); mpfr_fma(r1764842, r1764793, r1764841, r1764792, MPFR_RNDN); mpfr_div(r1764843, r1764792, r1764842, MPFR_RNDN); if (mpfr_get_si(r1764832, MPFR_RNDN)) { mpfr_set(r1764844, r1764843, MPFR_RNDN); } else { mpfr_set(r1764844, r1764843, MPFR_RNDN); }; if (mpfr_get_si(r1764820, MPFR_RNDN)) { mpfr_set(r1764845, r1764830, MPFR_RNDN); } else { mpfr_set(r1764845, r1764844, MPFR_RNDN); }; return mpfr_get_d(r1764845, MPFR_RNDN); } static mpfr_t r1764846, r1764847, r1764848, r1764849, r1764850, r1764851, r1764852, r1764853, r1764854, r1764855, r1764856, r1764857, r1764858, r1764859, r1764860, r1764861, r1764862, r1764863, r1764864, r1764865, r1764866, r1764867, r1764868, r1764869, r1764870, r1764871, r1764872, r1764873, r1764874, r1764875, r1764876, r1764877, r1764878, r1764879, r1764880, r1764881, r1764882, r1764883, r1764884, r1764885, r1764886, r1764887, r1764888, r1764889, r1764890, r1764891, r1764892, r1764893, r1764894, r1764895, r1764896, r1764897, r1764898, r1764899; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r1764846); mpfr_init(r1764847); mpfr_init_set_str(r1764848, "2.0", 10, MPFR_RNDN); mpfr_init(r1764849); mpfr_init(r1764850); mpfr_init(r1764851); mpfr_init(r1764852); mpfr_init(r1764853); mpfr_init(r1764854); mpfr_init(r1764855); mpfr_init(r1764856); mpfr_init(r1764857); mpfr_init(r1764858); mpfr_init_set_str(r1764859, "5.0", 10, MPFR_RNDN); mpfr_init_set_str(r1764860, "6.0", 10, MPFR_RNDN); mpfr_init(r1764861); mpfr_init(r1764862); mpfr_init_set_str(r1764863, "3.0", 10, MPFR_RNDN); mpfr_init(r1764864); mpfr_init(r1764865); mpfr_init(r1764866); mpfr_init(r1764867); mpfr_init(r1764868); mpfr_init(r1764869); mpfr_init(r1764870); mpfr_init(r1764871); mpfr_init(r1764872); mpfr_init(r1764873); mpfr_init(r1764874); mpfr_init_set_str(r1764875, "1", 10, MPFR_RNDN); mpfr_init(r1764876); mpfr_init(r1764877); mpfr_init_set_str(r1764878, "0.8333333333333334", 10, MPFR_RNDN); mpfr_init(r1764879); mpfr_init(r1764880); mpfr_init(r1764881); mpfr_init(r1764882); mpfr_init(r1764883); mpfr_init(r1764884); mpfr_init_set_str(r1764885, "+inf.0", 10, MPFR_RNDN); mpfr_init(r1764886); mpfr_init(r1764887); mpfr_init(r1764888); mpfr_init(r1764889); mpfr_init(r1764890); mpfr_init(r1764891); mpfr_init(r1764892); mpfr_init(r1764893); mpfr_init(r1764894); mpfr_init(r1764895); mpfr_init(r1764896); mpfr_init(r1764897); mpfr_init(r1764898); mpfr_init(r1764899); } double f_dm(double x, double y, double z, double t, double a, double b, double c) { mpfr_set_d(r1764846, x, MPFR_RNDN); mpfr_set_d(r1764847, y, MPFR_RNDN); ; mpfr_set_d(r1764849, z, MPFR_RNDN); mpfr_set_d(r1764850, t, MPFR_RNDN); mpfr_set_d(r1764851, a, MPFR_RNDN); mpfr_add(r1764852, r1764850, r1764851, MPFR_RNDN); mpfr_sqrt(r1764853, r1764852, MPFR_RNDN); mpfr_mul(r1764854, r1764849, r1764853, MPFR_RNDN); mpfr_div(r1764855, r1764854, r1764850, MPFR_RNDN); mpfr_set_d(r1764856, b, MPFR_RNDN); mpfr_set_d(r1764857, c, MPFR_RNDN); mpfr_sub(r1764858, r1764856, r1764857, MPFR_RNDN); ; ; mpfr_div(r1764861, r1764859, r1764860, MPFR_RNDN); mpfr_add(r1764862, r1764851, r1764861, MPFR_RNDN); ; mpfr_mul(r1764864, r1764850, r1764863, MPFR_RNDN); mpfr_div(r1764865, r1764848, r1764864, MPFR_RNDN); mpfr_sub(r1764866, r1764862, r1764865, MPFR_RNDN); mpfr_mul(r1764867, r1764858, r1764866, MPFR_RNDN); mpfr_sub(r1764868, r1764855, r1764867, MPFR_RNDN); mpfr_mul(r1764869, r1764848, r1764868, MPFR_RNDN); mpfr_exp(r1764870, r1764869, MPFR_RNDN); mpfr_mul(r1764871, r1764847, r1764870, MPFR_RNDN); mpfr_add(r1764872, r1764846, r1764871, MPFR_RNDN); mpfr_div(r1764873, r1764846, r1764872, MPFR_RNDN); mpfr_set_si(r1764874, mpfr_cmp(r1764873, r1764873) != 0, MPFR_RNDN); ; mpfr_exp(r1764876, r1764848, MPFR_RNDN); mpfr_sub(r1764877, r1764857, r1764856, MPFR_RNDN); ; mpfr_mul(r1764879, r1764857, r1764878, MPFR_RNDN); mpfr_fma(r1764880, r1764877, r1764851, r1764879, MPFR_RNDN); mpfr_pow(r1764881, r1764876, r1764880, MPFR_RNDN); mpfr_fma(r1764882, r1764847, r1764881, r1764846, MPFR_RNDN); mpfr_div(r1764883, r1764882, r1764846, MPFR_RNDN); mpfr_div(r1764884, r1764875, r1764883, MPFR_RNDN); ; mpfr_set_si(r1764886, mpfr_cmp(r1764873, r1764885) <= 0, MPFR_RNDN); mpfr_div(r1764887, r1764849, r1764850, MPFR_RNDN); mpfr_mul(r1764888, r1764887, r1764853, MPFR_RNDN); mpfr_add(r1764889, r1764861, r1764851, MPFR_RNDN); mpfr_div(r1764890, r1764848, r1764850, MPFR_RNDN); mpfr_div(r1764891, r1764890, r1764863, MPFR_RNDN); mpfr_sub(r1764892, r1764889, r1764891, MPFR_RNDN); mpfr_mul(r1764893, r1764892, r1764858, MPFR_RNDN); mpfr_sub(r1764894, r1764888, r1764893, MPFR_RNDN); mpfr_pow(r1764895, r1764876, r1764894, MPFR_RNDN); mpfr_fma(r1764896, r1764847, r1764895, r1764846, MPFR_RNDN); mpfr_div(r1764897, r1764846, r1764896, MPFR_RNDN); if (mpfr_get_si(r1764886, MPFR_RNDN)) { mpfr_set(r1764898, r1764897, MPFR_RNDN); } else { mpfr_set(r1764898, r1764897, MPFR_RNDN); }; if (mpfr_get_si(r1764874, MPFR_RNDN)) { mpfr_set(r1764899, r1764884, MPFR_RNDN); } else { mpfr_set(r1764899, r1764898, MPFR_RNDN); }; return mpfr_get_d(r1764899, MPFR_RNDN); }