#include #include #include #include #include char *name = "jeff quadratic root 2"; double f_if(float a, float b, float c) { float r20772 = b; float r20773 = 0; bool r20774 = r20772 >= r20773; float r20775 = 2; float r20776 = c; float r20777 = r20775 * r20776; float r20778 = -r20772; float r20779 = r20772 * r20772; float r20780 = 4; float r20781 = a; float r20782 = r20780 * r20781; float r20783 = r20782 * r20776; float r20784 = r20779 - r20783; float r20785 = sqrt(r20784); float r20786 = r20778 - r20785; float r20787 = r20777 / r20786; float r20788 = r20778 + r20785; float r20789 = r20775 * r20781; float r20790 = r20788 / r20789; float r20791 = r20774 ? r20787 : r20790; return r20791; } double f_id(double a, double b, double c) { double r20792 = b; double r20793 = 0; bool r20794 = r20792 >= r20793; double r20795 = 2; double r20796 = c; double r20797 = r20795 * r20796; double r20798 = -r20792; double r20799 = r20792 * r20792; double r20800 = 4; double r20801 = a; double r20802 = r20800 * r20801; double r20803 = r20802 * r20796; double r20804 = r20799 - r20803; double r20805 = sqrt(r20804); double r20806 = r20798 - r20805; double r20807 = r20797 / r20806; double r20808 = r20798 + r20805; double r20809 = r20795 * r20801; double r20810 = r20808 / r20809; double r20811 = r20794 ? r20807 : r20810; return r20811; } double f_of(float a, float b, float c) { float r20812 = b; float r20813 = -5.4941265907654905e+109; bool r20814 = r20812 <= r20813; float r20815 = 0; bool r20816 = r20812 >= r20815; float r20817 = c; float r20818 = r20817 + r20817; float r20819 = -r20812; float r20820 = r20812 * r20812; float r20821 = a; float r20822 = r20821 * r20817; float r20823 = 4; float r20824 = r20822 * r20823; float r20825 = r20820 - r20824; float r20826 = sqrt(r20825); float r20827 = r20819 - r20826; float r20828 = r20818 / r20827; float r20829 = r20817 / r20812; float r20830 = 1; float r20831 = r20829 / r20830; float r20832 = r20812 - r20819; float r20833 = r20821 + r20821; float r20834 = r20832 / r20833; float r20835 = r20831 - r20834; float r20836 = r20816 ? r20828 : r20835; float r20837 = 6.886487059856606e+121; bool r20838 = r20812 <= r20837; float r20839 = 2; float r20840 = r20839 * r20817; float r20841 = r20823 * r20821; float r20842 = r20841 * r20817; float r20843 = r20820 - r20842; float r20844 = sqrt(r20843); float r20845 = r20819 - r20844; float r20846 = r20840 / r20845; float r20847 = r20823 * r20817; float r20848 = r20821 * r20847; float r20849 = r20820 - r20848; float r20850 = cbrt(r20849); float r20851 = fabs(r20850); float r20852 = cbrt(r20843); float r20853 = sqrt(r20852); float r20854 = r20851 * r20853; float r20855 = r20819 + r20854; float r20856 = cbrt(r20855); float r20857 = r20856 * r20856; float r20858 = r20857 * r20856; float r20859 = r20839 * r20821; float r20860 = r20858 / r20859; float r20861 = r20816 ? r20846 : r20860; float r20862 = +inf.0; bool r20863 = r20812 <= r20862; float r20864 = r20821 / r20812; float r20865 = r20817 * r20864; float r20866 = r20865 - r20812; float r20867 = r20817 / r20866; float r20868 = r20819 + r20844; float r20869 = r20868 / r20833; float r20870 = r20816 ? r20867 : r20869; float r20871 = r20863 ? r20870 : r20870; float r20872 = r20838 ? r20861 : r20871; float r20873 = r20814 ? r20836 : r20872; return r20873; } double f_od(double a, double b, double c) { double r20874 = b; double r20875 = -5.4941265907654905e+109; bool r20876 = r20874 <= r20875; double r20877 = 0; bool r20878 = r20874 >= r20877; double r20879 = c; double r20880 = r20879 + r20879; double r20881 = -r20874; double r20882 = r20874 * r20874; double r20883 = a; double r20884 = r20883 * r20879; double r20885 = 4; double r20886 = r20884 * r20885; double r20887 = r20882 - r20886; double r20888 = sqrt(r20887); double r20889 = r20881 - r20888; double r20890 = r20880 / r20889; double r20891 = r20879 / r20874; double r20892 = 1; double r20893 = r20891 / r20892; double r20894 = r20874 - r20881; double r20895 = r20883 + r20883; double r20896 = r20894 / r20895; double r20897 = r20893 - r20896; double r20898 = r20878 ? r20890 : r20897; double r20899 = 6.886487059856606e+121; bool r20900 = r20874 <= r20899; double r20901 = 2; double r20902 = r20901 * r20879; double r20903 = r20885 * r20883; double r20904 = r20903 * r20879; double r20905 = r20882 - r20904; double r20906 = sqrt(r20905); double r20907 = r20881 - r20906; double r20908 = r20902 / r20907; double r20909 = r20885 * r20879; double r20910 = r20883 * r20909; double r20911 = r20882 - r20910; double r20912 = cbrt(r20911); double r20913 = fabs(r20912); double r20914 = cbrt(r20905); double r20915 = sqrt(r20914); double r20916 = r20913 * r20915; double r20917 = r20881 + r20916; double r20918 = cbrt(r20917); double r20919 = r20918 * r20918; double r20920 = r20919 * r20918; double r20921 = r20901 * r20883; double r20922 = r20920 / r20921; double r20923 = r20878 ? r20908 : r20922; double r20924 = +inf.0; bool r20925 = r20874 <= r20924; double r20926 = r20883 / r20874; double r20927 = r20879 * r20926; double r20928 = r20927 - r20874; double r20929 = r20879 / r20928; double r20930 = r20881 + r20906; double r20931 = r20930 / r20895; double r20932 = r20878 ? r20929 : r20931; double r20933 = r20925 ? r20932 : r20932; double r20934 = r20900 ? r20923 : r20933; double r20935 = r20876 ? r20898 : r20934; return r20935; } 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 r20936, r20937, r20938, r20939, r20940, r20941, r20942, r20943, r20944, r20945, r20946, r20947, r20948, r20949, r20950, r20951, r20952, r20953, r20954, r20955; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r20936); mpfr_init_set_str(r20937, "0", 10, MPFR_RNDN); mpfr_init(r20938); mpfr_init_set_str(r20939, "2", 10, MPFR_RNDN); mpfr_init(r20940); mpfr_init(r20941); mpfr_init(r20942); mpfr_init(r20943); mpfr_init_set_str(r20944, "4", 10, MPFR_RNDN); mpfr_init(r20945); mpfr_init(r20946); mpfr_init(r20947); mpfr_init(r20948); mpfr_init(r20949); mpfr_init(r20950); mpfr_init(r20951); mpfr_init(r20952); mpfr_init(r20953); mpfr_init(r20954); mpfr_init(r20955); } double f_im(double a, double b, double c) { mpfr_set_d(r20936, b, MPFR_RNDN); ; mpfr_set_si(r20938, mpfr_cmp(r20936, r20937) >= 0, MPFR_RNDN); ; mpfr_set_d(r20940, c, MPFR_RNDN); mpfr_mul(r20941, r20939, r20940, MPFR_RNDN); mpfr_neg(r20942, r20936, MPFR_RNDN); mpfr_mul(r20943, r20936, r20936, MPFR_RNDN); ; mpfr_set_d(r20945, a, MPFR_RNDN); mpfr_mul(r20946, r20944, r20945, MPFR_RNDN); mpfr_mul(r20947, r20946, r20940, MPFR_RNDN); mpfr_sub(r20948, r20943, r20947, MPFR_RNDN); mpfr_sqrt(r20949, r20948, MPFR_RNDN); mpfr_sub(r20950, r20942, r20949, MPFR_RNDN); mpfr_div(r20951, r20941, r20950, MPFR_RNDN); mpfr_add(r20952, r20942, r20949, MPFR_RNDN); mpfr_mul(r20953, r20939, r20945, MPFR_RNDN); mpfr_div(r20954, r20952, r20953, MPFR_RNDN); if (mpfr_get_si(r20938, MPFR_RNDN)) { mpfr_set(r20955, r20951, MPFR_RNDN); } else { mpfr_set(r20955, r20954, MPFR_RNDN); }; return mpfr_get_d(r20955, MPFR_RNDN); } static mpfr_t r20956, r20957, r20958, r20959, r20960, r20961, r20962, r20963, r20964, r20965, r20966, r20967, r20968, r20969, r20970, r20971, r20972, r20973, r20974, r20975, r20976, r20977, r20978, r20979, r20980, r20981, r20982, r20983, r20984, r20985, r20986, r20987, r20988, r20989, r20990, r20991, r20992, r20993, r20994, r20995, r20996, r20997, r20998, r20999, r21000, r21001, r21002, r21003, r21004, r21005, r21006, r21007, r21008, r21009, r21010, r21011, r21012, r21013, r21014, r21015, r21016, r21017; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r20956); mpfr_init_set_str(r20957, "-5.4941265907654905e+109", 10, MPFR_RNDN); mpfr_init(r20958); mpfr_init_set_str(r20959, "0", 10, MPFR_RNDN); mpfr_init(r20960); mpfr_init(r20961); mpfr_init(r20962); mpfr_init(r20963); mpfr_init(r20964); mpfr_init(r20965); mpfr_init(r20966); mpfr_init_set_str(r20967, "4", 10, MPFR_RNDN); mpfr_init(r20968); mpfr_init(r20969); mpfr_init(r20970); mpfr_init(r20971); mpfr_init(r20972); mpfr_init(r20973); mpfr_init_set_str(r20974, "1", 10, MPFR_RNDN); mpfr_init(r20975); mpfr_init(r20976); mpfr_init(r20977); mpfr_init(r20978); mpfr_init(r20979); mpfr_init(r20980); mpfr_init_set_str(r20981, "6.886487059856606e+121", 10, MPFR_RNDN); mpfr_init(r20982); mpfr_init_set_str(r20983, "2", 10, MPFR_RNDN); mpfr_init(r20984); mpfr_init(r20985); mpfr_init(r20986); mpfr_init(r20987); mpfr_init(r20988); mpfr_init(r20989); mpfr_init(r20990); mpfr_init(r20991); mpfr_init(r20992); mpfr_init(r20993); mpfr_init(r20994); mpfr_init(r20995); mpfr_init(r20996); mpfr_init(r20997); mpfr_init(r20998); mpfr_init(r20999); mpfr_init(r21000); mpfr_init(r21001); mpfr_init(r21002); mpfr_init(r21003); mpfr_init(r21004); mpfr_init(r21005); mpfr_init_set_str(r21006, "+inf.0", 10, MPFR_RNDN); mpfr_init(r21007); mpfr_init(r21008); mpfr_init(r21009); mpfr_init(r21010); mpfr_init(r21011); mpfr_init(r21012); mpfr_init(r21013); mpfr_init(r21014); mpfr_init(r21015); mpfr_init(r21016); mpfr_init(r21017); } double f_fm(double a, double b, double c) { mpfr_set_d(r20956, b, MPFR_RNDN); ; mpfr_set_si(r20958, mpfr_cmp(r20956, r20957) <= 0, MPFR_RNDN); ; mpfr_set_si(r20960, mpfr_cmp(r20956, r20959) >= 0, MPFR_RNDN); mpfr_set_d(r20961, c, MPFR_RNDN); mpfr_add(r20962, r20961, r20961, MPFR_RNDN); mpfr_neg(r20963, r20956, MPFR_RNDN); mpfr_mul(r20964, r20956, r20956, MPFR_RNDN); mpfr_set_d(r20965, a, MPFR_RNDN); mpfr_mul(r20966, r20965, r20961, MPFR_RNDN); ; mpfr_mul(r20968, r20966, r20967, MPFR_RNDN); mpfr_sub(r20969, r20964, r20968, MPFR_RNDN); mpfr_sqrt(r20970, r20969, MPFR_RNDN); mpfr_sub(r20971, r20963, r20970, MPFR_RNDN); mpfr_div(r20972, r20962, r20971, MPFR_RNDN); mpfr_div(r20973, r20961, r20956, MPFR_RNDN); ; mpfr_div(r20975, r20973, r20974, MPFR_RNDN); mpfr_sub(r20976, r20956, r20963, MPFR_RNDN); mpfr_add(r20977, r20965, r20965, MPFR_RNDN); mpfr_div(r20978, r20976, r20977, MPFR_RNDN); mpfr_sub(r20979, r20975, r20978, MPFR_RNDN); if (mpfr_get_si(r20960, MPFR_RNDN)) { mpfr_set(r20980, r20972, MPFR_RNDN); } else { mpfr_set(r20980, r20979, MPFR_RNDN); }; ; mpfr_set_si(r20982, mpfr_cmp(r20956, r20981) <= 0, MPFR_RNDN); ; mpfr_mul(r20984, r20983, r20961, MPFR_RNDN); mpfr_mul(r20985, r20967, r20965, MPFR_RNDN); mpfr_mul(r20986, r20985, r20961, MPFR_RNDN); mpfr_sub(r20987, r20964, r20986, MPFR_RNDN); mpfr_sqrt(r20988, r20987, MPFR_RNDN); mpfr_sub(r20989, r20963, r20988, MPFR_RNDN); mpfr_div(r20990, r20984, r20989, MPFR_RNDN); mpfr_mul(r20991, r20967, r20961, MPFR_RNDN); mpfr_mul(r20992, r20965, r20991, MPFR_RNDN); mpfr_sub(r20993, r20964, r20992, MPFR_RNDN); mpfr_cbrt(r20994, r20993, MPFR_RNDN); mpfr_abs(r20995, r20994, MPFR_RNDN); mpfr_cbrt(r20996, r20987, MPFR_RNDN); mpfr_sqrt(r20997, r20996, MPFR_RNDN); mpfr_mul(r20998, r20995, r20997, MPFR_RNDN); mpfr_add(r20999, r20963, r20998, MPFR_RNDN); mpfr_cbrt(r21000, r20999, MPFR_RNDN); mpfr_mul(r21001, r21000, r21000, MPFR_RNDN); mpfr_mul(r21002, r21001, r21000, MPFR_RNDN); mpfr_mul(r21003, r20983, r20965, MPFR_RNDN); mpfr_div(r21004, r21002, r21003, MPFR_RNDN); if (mpfr_get_si(r20960, MPFR_RNDN)) { mpfr_set(r21005, r20990, MPFR_RNDN); } else { mpfr_set(r21005, r21004, MPFR_RNDN); }; ; mpfr_set_si(r21007, mpfr_cmp(r20956, r21006) <= 0, MPFR_RNDN); mpfr_div(r21008, r20965, r20956, MPFR_RNDN); mpfr_mul(r21009, r20961, r21008, MPFR_RNDN); mpfr_sub(r21010, r21009, r20956, MPFR_RNDN); mpfr_div(r21011, r20961, r21010, MPFR_RNDN); mpfr_add(r21012, r20963, r20988, MPFR_RNDN); mpfr_div(r21013, r21012, r20977, MPFR_RNDN); if (mpfr_get_si(r20960, MPFR_RNDN)) { mpfr_set(r21014, r21011, MPFR_RNDN); } else { mpfr_set(r21014, r21013, MPFR_RNDN); }; if (mpfr_get_si(r21007, MPFR_RNDN)) { mpfr_set(r21015, r21014, MPFR_RNDN); } else { mpfr_set(r21015, r21014, MPFR_RNDN); }; if (mpfr_get_si(r20982, MPFR_RNDN)) { mpfr_set(r21016, r21005, MPFR_RNDN); } else { mpfr_set(r21016, r21015, MPFR_RNDN); }; if (mpfr_get_si(r20958, MPFR_RNDN)) { mpfr_set(r21017, r20980, MPFR_RNDN); } else { mpfr_set(r21017, r21016, MPFR_RNDN); }; return mpfr_get_d(r21017, MPFR_RNDN); } static mpfr_t r21018, r21019, r21020, r21021, r21022, r21023, r21024, r21025, r21026, r21027, r21028, r21029, r21030, r21031, r21032, r21033, r21034, r21035, r21036, r21037, r21038, r21039, r21040, r21041, r21042, r21043, r21044, r21045, r21046, r21047, r21048, r21049, r21050, r21051, r21052, r21053, r21054, r21055, r21056, r21057, r21058, r21059, r21060, r21061, r21062, r21063, r21064, r21065, r21066, r21067, r21068, r21069, r21070, r21071, r21072, r21073, r21074, r21075, r21076, r21077, r21078, r21079; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r21018); mpfr_init_set_str(r21019, "-5.4941265907654905e+109", 10, MPFR_RNDN); mpfr_init(r21020); mpfr_init_set_str(r21021, "0", 10, MPFR_RNDN); mpfr_init(r21022); mpfr_init(r21023); mpfr_init(r21024); mpfr_init(r21025); mpfr_init(r21026); mpfr_init(r21027); mpfr_init(r21028); mpfr_init_set_str(r21029, "4", 10, MPFR_RNDN); mpfr_init(r21030); mpfr_init(r21031); mpfr_init(r21032); mpfr_init(r21033); mpfr_init(r21034); mpfr_init(r21035); mpfr_init_set_str(r21036, "1", 10, MPFR_RNDN); mpfr_init(r21037); mpfr_init(r21038); mpfr_init(r21039); mpfr_init(r21040); mpfr_init(r21041); mpfr_init(r21042); mpfr_init_set_str(r21043, "6.886487059856606e+121", 10, MPFR_RNDN); mpfr_init(r21044); mpfr_init_set_str(r21045, "2", 10, MPFR_RNDN); mpfr_init(r21046); mpfr_init(r21047); mpfr_init(r21048); mpfr_init(r21049); mpfr_init(r21050); mpfr_init(r21051); mpfr_init(r21052); mpfr_init(r21053); mpfr_init(r21054); mpfr_init(r21055); mpfr_init(r21056); mpfr_init(r21057); mpfr_init(r21058); mpfr_init(r21059); mpfr_init(r21060); mpfr_init(r21061); mpfr_init(r21062); mpfr_init(r21063); mpfr_init(r21064); mpfr_init(r21065); mpfr_init(r21066); mpfr_init(r21067); mpfr_init_set_str(r21068, "+inf.0", 10, MPFR_RNDN); mpfr_init(r21069); mpfr_init(r21070); mpfr_init(r21071); mpfr_init(r21072); mpfr_init(r21073); mpfr_init(r21074); mpfr_init(r21075); mpfr_init(r21076); mpfr_init(r21077); mpfr_init(r21078); mpfr_init(r21079); } double f_dm(double a, double b, double c) { mpfr_set_d(r21018, b, MPFR_RNDN); ; mpfr_set_si(r21020, mpfr_cmp(r21018, r21019) <= 0, MPFR_RNDN); ; mpfr_set_si(r21022, mpfr_cmp(r21018, r21021) >= 0, MPFR_RNDN); mpfr_set_d(r21023, c, MPFR_RNDN); mpfr_add(r21024, r21023, r21023, MPFR_RNDN); mpfr_neg(r21025, r21018, MPFR_RNDN); mpfr_mul(r21026, r21018, r21018, MPFR_RNDN); mpfr_set_d(r21027, a, MPFR_RNDN); mpfr_mul(r21028, r21027, r21023, MPFR_RNDN); ; mpfr_mul(r21030, r21028, r21029, MPFR_RNDN); mpfr_sub(r21031, r21026, r21030, MPFR_RNDN); mpfr_sqrt(r21032, r21031, MPFR_RNDN); mpfr_sub(r21033, r21025, r21032, MPFR_RNDN); mpfr_div(r21034, r21024, r21033, MPFR_RNDN); mpfr_div(r21035, r21023, r21018, MPFR_RNDN); ; mpfr_div(r21037, r21035, r21036, MPFR_RNDN); mpfr_sub(r21038, r21018, r21025, MPFR_RNDN); mpfr_add(r21039, r21027, r21027, MPFR_RNDN); mpfr_div(r21040, r21038, r21039, MPFR_RNDN); mpfr_sub(r21041, r21037, r21040, MPFR_RNDN); if (mpfr_get_si(r21022, MPFR_RNDN)) { mpfr_set(r21042, r21034, MPFR_RNDN); } else { mpfr_set(r21042, r21041, MPFR_RNDN); }; ; mpfr_set_si(r21044, mpfr_cmp(r21018, r21043) <= 0, MPFR_RNDN); ; mpfr_mul(r21046, r21045, r21023, MPFR_RNDN); mpfr_mul(r21047, r21029, r21027, MPFR_RNDN); mpfr_mul(r21048, r21047, r21023, MPFR_RNDN); mpfr_sub(r21049, r21026, r21048, MPFR_RNDN); mpfr_sqrt(r21050, r21049, MPFR_RNDN); mpfr_sub(r21051, r21025, r21050, MPFR_RNDN); mpfr_div(r21052, r21046, r21051, MPFR_RNDN); mpfr_mul(r21053, r21029, r21023, MPFR_RNDN); mpfr_mul(r21054, r21027, r21053, MPFR_RNDN); mpfr_sub(r21055, r21026, r21054, MPFR_RNDN); mpfr_cbrt(r21056, r21055, MPFR_RNDN); mpfr_abs(r21057, r21056, MPFR_RNDN); mpfr_cbrt(r21058, r21049, MPFR_RNDN); mpfr_sqrt(r21059, r21058, MPFR_RNDN); mpfr_mul(r21060, r21057, r21059, MPFR_RNDN); mpfr_add(r21061, r21025, r21060, MPFR_RNDN); mpfr_cbrt(r21062, r21061, MPFR_RNDN); mpfr_mul(r21063, r21062, r21062, MPFR_RNDN); mpfr_mul(r21064, r21063, r21062, MPFR_RNDN); mpfr_mul(r21065, r21045, r21027, MPFR_RNDN); mpfr_div(r21066, r21064, r21065, MPFR_RNDN); if (mpfr_get_si(r21022, MPFR_RNDN)) { mpfr_set(r21067, r21052, MPFR_RNDN); } else { mpfr_set(r21067, r21066, MPFR_RNDN); }; ; mpfr_set_si(r21069, mpfr_cmp(r21018, r21068) <= 0, MPFR_RNDN); mpfr_div(r21070, r21027, r21018, MPFR_RNDN); mpfr_mul(r21071, r21023, r21070, MPFR_RNDN); mpfr_sub(r21072, r21071, r21018, MPFR_RNDN); mpfr_div(r21073, r21023, r21072, MPFR_RNDN); mpfr_add(r21074, r21025, r21050, MPFR_RNDN); mpfr_div(r21075, r21074, r21039, MPFR_RNDN); if (mpfr_get_si(r21022, MPFR_RNDN)) { mpfr_set(r21076, r21073, MPFR_RNDN); } else { mpfr_set(r21076, r21075, MPFR_RNDN); }; if (mpfr_get_si(r21069, MPFR_RNDN)) { mpfr_set(r21077, r21076, MPFR_RNDN); } else { mpfr_set(r21077, r21076, MPFR_RNDN); }; if (mpfr_get_si(r21044, MPFR_RNDN)) { mpfr_set(r21078, r21067, MPFR_RNDN); } else { mpfr_set(r21078, r21077, MPFR_RNDN); }; if (mpfr_get_si(r21020, MPFR_RNDN)) { mpfr_set(r21079, r21042, MPFR_RNDN); } else { mpfr_set(r21079, r21078, MPFR_RNDN); }; return mpfr_get_d(r21079, MPFR_RNDN); }