#include #include #include #include #include char *name = "Toniolo and Linder, Equation (7)"; double f_if(float x, float l, float t) { float r30461875 = 2; float r30461876 = sqrt(r30461875); float r30461877 = t; float r30461878 = r30461876 * r30461877; float r30461879 = x; float r30461880 = 1; float r30461881 = r30461879 + r30461880; float r30461882 = r30461879 - r30461880; float r30461883 = r30461881 / r30461882; float r30461884 = l; float r30461885 = r30461884 * r30461884; float r30461886 = r30461877 * r30461877; float r30461887 = r30461875 * r30461886; float r30461888 = r30461885 + r30461887; float r30461889 = r30461883 * r30461888; float r30461890 = r30461889 - r30461885; float r30461891 = sqrt(r30461890); float r30461892 = r30461878 / r30461891; return r30461892; } double f_id(double x, double l, double t) { double r30461893 = 2; double r30461894 = sqrt(r30461893); double r30461895 = t; double r30461896 = r30461894 * r30461895; double r30461897 = x; double r30461898 = 1; double r30461899 = r30461897 + r30461898; double r30461900 = r30461897 - r30461898; double r30461901 = r30461899 / r30461900; double r30461902 = l; double r30461903 = r30461902 * r30461902; double r30461904 = r30461895 * r30461895; double r30461905 = r30461893 * r30461904; double r30461906 = r30461903 + r30461905; double r30461907 = r30461901 * r30461906; double r30461908 = r30461907 - r30461903; double r30461909 = sqrt(r30461908); double r30461910 = r30461896 / r30461909; return r30461910; } double f_of(float x, float l, float t) { float r30461911 = t; float r30461912 = -1.2723937600129117e+124; bool r30461913 = r30461911 <= r30461912; float r30461914 = 2; float r30461915 = sqrt(r30461914); float r30461916 = r30461911 * r30461915; float r30461917 = r30461911 / r30461915; float r30461918 = x; float r30461919 = r30461918 * r30461918; float r30461920 = r30461917 / r30461919; float r30461921 = r30461911 / r30461918; float r30461922 = r30461914 / r30461918; float r30461923 = r30461922 / r30461915; float r30461924 = r30461921 * r30461923; float r30461925 = r30461920 - r30461924; float r30461926 = r30461923 + r30461915; float r30461927 = r30461911 * r30461926; float r30461928 = r30461925 - r30461927; float r30461929 = r30461916 / r30461928; float r30461930 = 2.639444841744209e-258; bool r30461931 = r30461911 <= r30461930; float r30461932 = pow(r30461911, r30461914); float r30461933 = r30461914 * r30461932; float r30461934 = 4; float r30461935 = r30461932 / r30461918; float r30461936 = r30461934 * r30461935; float r30461937 = l; float r30461938 = r30461918 / r30461937; float r30461939 = r30461937 / r30461938; float r30461940 = r30461914 * r30461939; float r30461941 = r30461936 + r30461940; float r30461942 = r30461933 + r30461941; float r30461943 = sqrt(r30461942); float r30461944 = sqrt(r30461943); float r30461945 = r30461915 / r30461944; float r30461946 = r30461911 / r30461944; float r30461947 = r30461945 * r30461946; float r30461948 = 1.0741008202836691e-159; bool r30461949 = r30461911 <= r30461948; float r30461950 = r30461915 + r30461923; float r30461951 = r30461911 * r30461950; float r30461952 = r30461911 + r30461911; float r30461953 = r30461919 * r30461915; float r30461954 = r30461952 / r30461953; float r30461955 = r30461921 / r30461918; float r30461956 = r30461955 / r30461915; float r30461957 = r30461954 - r30461956; float r30461958 = r30461951 + r30461957; float r30461959 = r30461916 / r30461958; float r30461960 = 1.0330302136080365e+41; bool r30461961 = r30461911 <= r30461960; float r30461962 = r30461961 ? r30461947 : r30461959; float r30461963 = r30461949 ? r30461959 : r30461962; float r30461964 = r30461931 ? r30461947 : r30461963; float r30461965 = r30461913 ? r30461929 : r30461964; return r30461965; } double f_od(double x, double l, double t) { double r30461966 = t; double r30461967 = -1.2723937600129117e+124; bool r30461968 = r30461966 <= r30461967; double r30461969 = 2; double r30461970 = sqrt(r30461969); double r30461971 = r30461966 * r30461970; double r30461972 = r30461966 / r30461970; double r30461973 = x; double r30461974 = r30461973 * r30461973; double r30461975 = r30461972 / r30461974; double r30461976 = r30461966 / r30461973; double r30461977 = r30461969 / r30461973; double r30461978 = r30461977 / r30461970; double r30461979 = r30461976 * r30461978; double r30461980 = r30461975 - r30461979; double r30461981 = r30461978 + r30461970; double r30461982 = r30461966 * r30461981; double r30461983 = r30461980 - r30461982; double r30461984 = r30461971 / r30461983; double r30461985 = 2.639444841744209e-258; bool r30461986 = r30461966 <= r30461985; double r30461987 = pow(r30461966, r30461969); double r30461988 = r30461969 * r30461987; double r30461989 = 4; double r30461990 = r30461987 / r30461973; double r30461991 = r30461989 * r30461990; double r30461992 = l; double r30461993 = r30461973 / r30461992; double r30461994 = r30461992 / r30461993; double r30461995 = r30461969 * r30461994; double r30461996 = r30461991 + r30461995; double r30461997 = r30461988 + r30461996; double r30461998 = sqrt(r30461997); double r30461999 = sqrt(r30461998); double r30462000 = r30461970 / r30461999; double r30462001 = r30461966 / r30461999; double r30462002 = r30462000 * r30462001; double r30462003 = 1.0741008202836691e-159; bool r30462004 = r30461966 <= r30462003; double r30462005 = r30461970 + r30461978; double r30462006 = r30461966 * r30462005; double r30462007 = r30461966 + r30461966; double r30462008 = r30461974 * r30461970; double r30462009 = r30462007 / r30462008; double r30462010 = r30461976 / r30461973; double r30462011 = r30462010 / r30461970; double r30462012 = r30462009 - r30462011; double r30462013 = r30462006 + r30462012; double r30462014 = r30461971 / r30462013; double r30462015 = 1.0330302136080365e+41; bool r30462016 = r30461966 <= r30462015; double r30462017 = r30462016 ? r30462002 : r30462014; double r30462018 = r30462004 ? r30462014 : r30462017; double r30462019 = r30461986 ? r30462002 : r30462018; double r30462020 = r30461968 ? r30461984 : r30462019; return r30462020; } 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 r30462021, r30462022, r30462023, r30462024, r30462025, r30462026, r30462027, r30462028, r30462029, r30462030, r30462031, r30462032, r30462033, r30462034, r30462035, r30462036, r30462037, r30462038; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init_set_str(r30462021, "2", 10, MPFR_RNDN); mpfr_init(r30462022); mpfr_init(r30462023); mpfr_init(r30462024); mpfr_init(r30462025); mpfr_init_set_str(r30462026, "1", 10, MPFR_RNDN); mpfr_init(r30462027); mpfr_init(r30462028); mpfr_init(r30462029); mpfr_init(r30462030); mpfr_init(r30462031); mpfr_init(r30462032); mpfr_init(r30462033); mpfr_init(r30462034); mpfr_init(r30462035); mpfr_init(r30462036); mpfr_init(r30462037); mpfr_init(r30462038); } double f_im(double x, double l, double t) { ; mpfr_sqrt(r30462022, r30462021, MPFR_RNDN); mpfr_set_d(r30462023, t, MPFR_RNDN); mpfr_mul(r30462024, r30462022, r30462023, MPFR_RNDN); mpfr_set_d(r30462025, x, MPFR_RNDN); ; mpfr_add(r30462027, r30462025, r30462026, MPFR_RNDN); mpfr_sub(r30462028, r30462025, r30462026, MPFR_RNDN); mpfr_div(r30462029, r30462027, r30462028, MPFR_RNDN); mpfr_set_d(r30462030, l, MPFR_RNDN); mpfr_mul(r30462031, r30462030, r30462030, MPFR_RNDN); mpfr_mul(r30462032, r30462023, r30462023, MPFR_RNDN); mpfr_mul(r30462033, r30462021, r30462032, MPFR_RNDN); mpfr_add(r30462034, r30462031, r30462033, MPFR_RNDN); mpfr_mul(r30462035, r30462029, r30462034, MPFR_RNDN); mpfr_sub(r30462036, r30462035, r30462031, MPFR_RNDN); mpfr_sqrt(r30462037, r30462036, MPFR_RNDN); mpfr_div(r30462038, r30462024, r30462037, MPFR_RNDN); return mpfr_get_d(r30462038, MPFR_RNDN); } static mpfr_t r30462039, r30462040, r30462041, r30462042, r30462043, r30462044, r30462045, r30462046, r30462047, r30462048, r30462049, r30462050, r30462051, r30462052, r30462053, r30462054, r30462055, r30462056, r30462057, r30462058, r30462059, r30462060, r30462061, r30462062, r30462063, r30462064, r30462065, r30462066, r30462067, r30462068, r30462069, r30462070, r30462071, r30462072, r30462073, r30462074, r30462075, r30462076, r30462077, r30462078, r30462079, r30462080, r30462081, r30462082, r30462083, r30462084, r30462085, r30462086, r30462087, r30462088, r30462089, r30462090, r30462091, r30462092, r30462093; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r30462039); mpfr_init_set_str(r30462040, "-1.2723937600129117e+124", 10, MPFR_RNDN); mpfr_init(r30462041); mpfr_init_set_str(r30462042, "2", 10, MPFR_RNDN); mpfr_init(r30462043); mpfr_init(r30462044); mpfr_init(r30462045); mpfr_init(r30462046); mpfr_init(r30462047); mpfr_init(r30462048); mpfr_init(r30462049); mpfr_init(r30462050); mpfr_init(r30462051); mpfr_init(r30462052); mpfr_init(r30462053); mpfr_init(r30462054); mpfr_init(r30462055); mpfr_init(r30462056); mpfr_init(r30462057); mpfr_init_set_str(r30462058, "2.639444841744209e-258", 10, MPFR_RNDN); mpfr_init(r30462059); mpfr_init(r30462060); mpfr_init(r30462061); mpfr_init_set_str(r30462062, "4", 10, MPFR_RNDN); mpfr_init(r30462063); mpfr_init(r30462064); mpfr_init(r30462065); mpfr_init(r30462066); mpfr_init(r30462067); mpfr_init(r30462068); mpfr_init(r30462069); mpfr_init(r30462070); mpfr_init(r30462071); mpfr_init(r30462072); mpfr_init(r30462073); mpfr_init(r30462074); mpfr_init(r30462075); mpfr_init_set_str(r30462076, "1.0741008202836691e-159", 10, MPFR_RNDN); mpfr_init(r30462077); mpfr_init(r30462078); mpfr_init(r30462079); mpfr_init(r30462080); mpfr_init(r30462081); mpfr_init(r30462082); mpfr_init(r30462083); mpfr_init(r30462084); mpfr_init(r30462085); mpfr_init(r30462086); mpfr_init(r30462087); mpfr_init_set_str(r30462088, "1.0330302136080365e+41", 10, MPFR_RNDN); mpfr_init(r30462089); mpfr_init(r30462090); mpfr_init(r30462091); mpfr_init(r30462092); mpfr_init(r30462093); } double f_fm(double x, double l, double t) { mpfr_set_d(r30462039, t, MPFR_RNDN); ; mpfr_set_si(r30462041, mpfr_cmp(r30462039, r30462040) <= 0, MPFR_RNDN); ; mpfr_sqrt(r30462043, r30462042, MPFR_RNDN); mpfr_mul(r30462044, r30462039, r30462043, MPFR_RNDN); mpfr_div(r30462045, r30462039, r30462043, MPFR_RNDN); mpfr_set_d(r30462046, x, MPFR_RNDN); mpfr_mul(r30462047, r30462046, r30462046, MPFR_RNDN); mpfr_div(r30462048, r30462045, r30462047, MPFR_RNDN); mpfr_div(r30462049, r30462039, r30462046, MPFR_RNDN); mpfr_div(r30462050, r30462042, r30462046, MPFR_RNDN); mpfr_div(r30462051, r30462050, r30462043, MPFR_RNDN); mpfr_mul(r30462052, r30462049, r30462051, MPFR_RNDN); mpfr_sub(r30462053, r30462048, r30462052, MPFR_RNDN); mpfr_add(r30462054, r30462051, r30462043, MPFR_RNDN); mpfr_mul(r30462055, r30462039, r30462054, MPFR_RNDN); mpfr_sub(r30462056, r30462053, r30462055, MPFR_RNDN); mpfr_div(r30462057, r30462044, r30462056, MPFR_RNDN); ; mpfr_set_si(r30462059, mpfr_cmp(r30462039, r30462058) <= 0, MPFR_RNDN); mpfr_pow(r30462060, r30462039, r30462042, MPFR_RNDN); mpfr_mul(r30462061, r30462042, r30462060, MPFR_RNDN); ; mpfr_div(r30462063, r30462060, r30462046, MPFR_RNDN); mpfr_mul(r30462064, r30462062, r30462063, MPFR_RNDN); mpfr_set_d(r30462065, l, MPFR_RNDN); mpfr_div(r30462066, r30462046, r30462065, MPFR_RNDN); mpfr_div(r30462067, r30462065, r30462066, MPFR_RNDN); mpfr_mul(r30462068, r30462042, r30462067, MPFR_RNDN); mpfr_add(r30462069, r30462064, r30462068, MPFR_RNDN); mpfr_add(r30462070, r30462061, r30462069, MPFR_RNDN); mpfr_sqrt(r30462071, r30462070, MPFR_RNDN); mpfr_sqrt(r30462072, r30462071, MPFR_RNDN); mpfr_div(r30462073, r30462043, r30462072, MPFR_RNDN); mpfr_div(r30462074, r30462039, r30462072, MPFR_RNDN); mpfr_mul(r30462075, r30462073, r30462074, MPFR_RNDN); ; mpfr_set_si(r30462077, mpfr_cmp(r30462039, r30462076) <= 0, MPFR_RNDN); mpfr_add(r30462078, r30462043, r30462051, MPFR_RNDN); mpfr_mul(r30462079, r30462039, r30462078, MPFR_RNDN); mpfr_add(r30462080, r30462039, r30462039, MPFR_RNDN); mpfr_mul(r30462081, r30462047, r30462043, MPFR_RNDN); mpfr_div(r30462082, r30462080, r30462081, MPFR_RNDN); mpfr_div(r30462083, r30462049, r30462046, MPFR_RNDN); mpfr_div(r30462084, r30462083, r30462043, MPFR_RNDN); mpfr_sub(r30462085, r30462082, r30462084, MPFR_RNDN); mpfr_add(r30462086, r30462079, r30462085, MPFR_RNDN); mpfr_div(r30462087, r30462044, r30462086, MPFR_RNDN); ; mpfr_set_si(r30462089, mpfr_cmp(r30462039, r30462088) <= 0, MPFR_RNDN); if (mpfr_get_si(r30462089, MPFR_RNDN)) { mpfr_set(r30462090, r30462075, MPFR_RNDN); } else { mpfr_set(r30462090, r30462087, MPFR_RNDN); }; if (mpfr_get_si(r30462077, MPFR_RNDN)) { mpfr_set(r30462091, r30462087, MPFR_RNDN); } else { mpfr_set(r30462091, r30462090, MPFR_RNDN); }; if (mpfr_get_si(r30462059, MPFR_RNDN)) { mpfr_set(r30462092, r30462075, MPFR_RNDN); } else { mpfr_set(r30462092, r30462091, MPFR_RNDN); }; if (mpfr_get_si(r30462041, MPFR_RNDN)) { mpfr_set(r30462093, r30462057, MPFR_RNDN); } else { mpfr_set(r30462093, r30462092, MPFR_RNDN); }; return mpfr_get_d(r30462093, MPFR_RNDN); } static mpfr_t r30462094, r30462095, r30462096, r30462097, r30462098, r30462099, r30462100, r30462101, r30462102, r30462103, r30462104, r30462105, r30462106, r30462107, r30462108, r30462109, r30462110, r30462111, r30462112, r30462113, r30462114, r30462115, r30462116, r30462117, r30462118, r30462119, r30462120, r30462121, r30462122, r30462123, r30462124, r30462125, r30462126, r30462127, r30462128, r30462129, r30462130, r30462131, r30462132, r30462133, r30462134, r30462135, r30462136, r30462137, r30462138, r30462139, r30462140, r30462141, r30462142, r30462143, r30462144, r30462145, r30462146, r30462147, r30462148; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r30462094); mpfr_init_set_str(r30462095, "-1.2723937600129117e+124", 10, MPFR_RNDN); mpfr_init(r30462096); mpfr_init_set_str(r30462097, "2", 10, MPFR_RNDN); mpfr_init(r30462098); mpfr_init(r30462099); mpfr_init(r30462100); mpfr_init(r30462101); mpfr_init(r30462102); mpfr_init(r30462103); mpfr_init(r30462104); mpfr_init(r30462105); mpfr_init(r30462106); mpfr_init(r30462107); mpfr_init(r30462108); mpfr_init(r30462109); mpfr_init(r30462110); mpfr_init(r30462111); mpfr_init(r30462112); mpfr_init_set_str(r30462113, "2.639444841744209e-258", 10, MPFR_RNDN); mpfr_init(r30462114); mpfr_init(r30462115); mpfr_init(r30462116); mpfr_init_set_str(r30462117, "4", 10, MPFR_RNDN); mpfr_init(r30462118); mpfr_init(r30462119); mpfr_init(r30462120); mpfr_init(r30462121); mpfr_init(r30462122); mpfr_init(r30462123); mpfr_init(r30462124); mpfr_init(r30462125); mpfr_init(r30462126); mpfr_init(r30462127); mpfr_init(r30462128); mpfr_init(r30462129); mpfr_init(r30462130); mpfr_init_set_str(r30462131, "1.0741008202836691e-159", 10, MPFR_RNDN); mpfr_init(r30462132); mpfr_init(r30462133); mpfr_init(r30462134); mpfr_init(r30462135); mpfr_init(r30462136); mpfr_init(r30462137); mpfr_init(r30462138); mpfr_init(r30462139); mpfr_init(r30462140); mpfr_init(r30462141); mpfr_init(r30462142); mpfr_init_set_str(r30462143, "1.0330302136080365e+41", 10, MPFR_RNDN); mpfr_init(r30462144); mpfr_init(r30462145); mpfr_init(r30462146); mpfr_init(r30462147); mpfr_init(r30462148); } double f_dm(double x, double l, double t) { mpfr_set_d(r30462094, t, MPFR_RNDN); ; mpfr_set_si(r30462096, mpfr_cmp(r30462094, r30462095) <= 0, MPFR_RNDN); ; mpfr_sqrt(r30462098, r30462097, MPFR_RNDN); mpfr_mul(r30462099, r30462094, r30462098, MPFR_RNDN); mpfr_div(r30462100, r30462094, r30462098, MPFR_RNDN); mpfr_set_d(r30462101, x, MPFR_RNDN); mpfr_mul(r30462102, r30462101, r30462101, MPFR_RNDN); mpfr_div(r30462103, r30462100, r30462102, MPFR_RNDN); mpfr_div(r30462104, r30462094, r30462101, MPFR_RNDN); mpfr_div(r30462105, r30462097, r30462101, MPFR_RNDN); mpfr_div(r30462106, r30462105, r30462098, MPFR_RNDN); mpfr_mul(r30462107, r30462104, r30462106, MPFR_RNDN); mpfr_sub(r30462108, r30462103, r30462107, MPFR_RNDN); mpfr_add(r30462109, r30462106, r30462098, MPFR_RNDN); mpfr_mul(r30462110, r30462094, r30462109, MPFR_RNDN); mpfr_sub(r30462111, r30462108, r30462110, MPFR_RNDN); mpfr_div(r30462112, r30462099, r30462111, MPFR_RNDN); ; mpfr_set_si(r30462114, mpfr_cmp(r30462094, r30462113) <= 0, MPFR_RNDN); mpfr_pow(r30462115, r30462094, r30462097, MPFR_RNDN); mpfr_mul(r30462116, r30462097, r30462115, MPFR_RNDN); ; mpfr_div(r30462118, r30462115, r30462101, MPFR_RNDN); mpfr_mul(r30462119, r30462117, r30462118, MPFR_RNDN); mpfr_set_d(r30462120, l, MPFR_RNDN); mpfr_div(r30462121, r30462101, r30462120, MPFR_RNDN); mpfr_div(r30462122, r30462120, r30462121, MPFR_RNDN); mpfr_mul(r30462123, r30462097, r30462122, MPFR_RNDN); mpfr_add(r30462124, r30462119, r30462123, MPFR_RNDN); mpfr_add(r30462125, r30462116, r30462124, MPFR_RNDN); mpfr_sqrt(r30462126, r30462125, MPFR_RNDN); mpfr_sqrt(r30462127, r30462126, MPFR_RNDN); mpfr_div(r30462128, r30462098, r30462127, MPFR_RNDN); mpfr_div(r30462129, r30462094, r30462127, MPFR_RNDN); mpfr_mul(r30462130, r30462128, r30462129, MPFR_RNDN); ; mpfr_set_si(r30462132, mpfr_cmp(r30462094, r30462131) <= 0, MPFR_RNDN); mpfr_add(r30462133, r30462098, r30462106, MPFR_RNDN); mpfr_mul(r30462134, r30462094, r30462133, MPFR_RNDN); mpfr_add(r30462135, r30462094, r30462094, MPFR_RNDN); mpfr_mul(r30462136, r30462102, r30462098, MPFR_RNDN); mpfr_div(r30462137, r30462135, r30462136, MPFR_RNDN); mpfr_div(r30462138, r30462104, r30462101, MPFR_RNDN); mpfr_div(r30462139, r30462138, r30462098, MPFR_RNDN); mpfr_sub(r30462140, r30462137, r30462139, MPFR_RNDN); mpfr_add(r30462141, r30462134, r30462140, MPFR_RNDN); mpfr_div(r30462142, r30462099, r30462141, MPFR_RNDN); ; mpfr_set_si(r30462144, mpfr_cmp(r30462094, r30462143) <= 0, MPFR_RNDN); if (mpfr_get_si(r30462144, MPFR_RNDN)) { mpfr_set(r30462145, r30462130, MPFR_RNDN); } else { mpfr_set(r30462145, r30462142, MPFR_RNDN); }; if (mpfr_get_si(r30462132, MPFR_RNDN)) { mpfr_set(r30462146, r30462142, MPFR_RNDN); } else { mpfr_set(r30462146, r30462145, MPFR_RNDN); }; if (mpfr_get_si(r30462114, MPFR_RNDN)) { mpfr_set(r30462147, r30462130, MPFR_RNDN); } else { mpfr_set(r30462147, r30462146, MPFR_RNDN); }; if (mpfr_get_si(r30462096, MPFR_RNDN)) { mpfr_set(r30462148, r30462112, MPFR_RNDN); } else { mpfr_set(r30462148, r30462147, MPFR_RNDN); }; return mpfr_get_d(r30462148, MPFR_RNDN); }