#include #include #include #include #include char *name = "Toniolo and Linder, Equation (10+)"; double f_if(float t, float l, float k) { float r20919 = 2; float r20920 = t; float r20921 = 3; float r20922 = pow(r20920, r20921); float r20923 = l; float r20924 = r20923 * r20923; float r20925 = r20922 / r20924; float r20926 = k; float r20927 = sin(r20926); float r20928 = r20925 * r20927; float r20929 = tan(r20926); float r20930 = r20928 * r20929; float r20931 = 1; float r20932 = r20926 / r20920; float r20933 = pow(r20932, r20919); float r20934 = r20931 + r20933; float r20935 = r20934 + r20931; float r20936 = r20930 * r20935; float r20937 = r20919 / r20936; return r20937; } double f_id(double t, double l, double k) { double r20938 = 2; double r20939 = t; double r20940 = 3; double r20941 = pow(r20939, r20940); double r20942 = l; double r20943 = r20942 * r20942; double r20944 = r20941 / r20943; double r20945 = k; double r20946 = sin(r20945); double r20947 = r20944 * r20946; double r20948 = tan(r20945); double r20949 = r20947 * r20948; double r20950 = 1; double r20951 = r20945 / r20939; double r20952 = pow(r20951, r20938); double r20953 = r20950 + r20952; double r20954 = r20953 + r20950; double r20955 = r20949 * r20954; double r20956 = r20938 / r20955; return r20956; } double f_of(float t, float l, float k) { float r20957 = 2; float r20958 = t; float r20959 = l; float r20960 = r20958 / r20959; float r20961 = k; float r20962 = sin(r20961); float r20963 = r20958 * r20962; float r20964 = r20960 * r20963; float r20965 = r20964 * r20962; float r20966 = 1; float r20967 = r20961 / r20958; float r20968 = pow(r20967, r20957); float r20969 = r20966 + r20968; float r20970 = r20969 + r20966; float r20971 = r20965 * r20970; float r20972 = r20959 / r20958; float r20973 = cos(r20961); float r20974 = r20972 * r20973; float r20975 = r20971 / r20974; float r20976 = r20957 / r20975; return r20976; } double f_od(double t, double l, double k) { double r20977 = 2; double r20978 = t; double r20979 = l; double r20980 = r20978 / r20979; double r20981 = k; double r20982 = sin(r20981); double r20983 = r20978 * r20982; double r20984 = r20980 * r20983; double r20985 = r20984 * r20982; double r20986 = 1; double r20987 = r20981 / r20978; double r20988 = pow(r20987, r20977); double r20989 = r20986 + r20988; double r20990 = r20989 + r20986; double r20991 = r20985 * r20990; double r20992 = r20979 / r20978; double r20993 = cos(r20981); double r20994 = r20992 * r20993; double r20995 = r20991 / r20994; double r20996 = r20977 / r20995; return r20996; } 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 r20997, r20998, r20999, r21000, r21001, r21002, r21003, r21004, r21005, r21006, r21007, r21008, r21009, r21010, r21011, r21012, r21013, r21014, r21015; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r20997, "2", 10, MPFR_RNDN); mpfr_init(r20998); mpfr_init_set_str(r20999, "3", 10, MPFR_RNDN); mpfr_init(r21000); mpfr_init(r21001); mpfr_init(r21002); mpfr_init(r21003); mpfr_init(r21004); mpfr_init(r21005); mpfr_init(r21006); mpfr_init(r21007); mpfr_init(r21008); mpfr_init_set_str(r21009, "1", 10, MPFR_RNDN); mpfr_init(r21010); mpfr_init(r21011); mpfr_init(r21012); mpfr_init(r21013); mpfr_init(r21014); mpfr_init(r21015); } double f_im(double t, double l, double k) { ; mpfr_set_d(r20998, t, MPFR_RNDN); ; mpfr_pow(r21000, r20998, r20999, MPFR_RNDN); mpfr_set_d(r21001, l, MPFR_RNDN); mpfr_mul(r21002, r21001, r21001, MPFR_RNDN); mpfr_div(r21003, r21000, r21002, MPFR_RNDN); mpfr_set_d(r21004, k, MPFR_RNDN); mpfr_sin(r21005, r21004, MPFR_RNDN); mpfr_mul(r21006, r21003, r21005, MPFR_RNDN); mpfr_tan(r21007, r21004, MPFR_RNDN); mpfr_mul(r21008, r21006, r21007, MPFR_RNDN); ; mpfr_div(r21010, r21004, r20998, MPFR_RNDN); mpfr_pow(r21011, r21010, r20997, MPFR_RNDN); mpfr_add(r21012, r21009, r21011, MPFR_RNDN); mpfr_add(r21013, r21012, r21009, MPFR_RNDN); mpfr_mul(r21014, r21008, r21013, MPFR_RNDN); mpfr_div(r21015, r20997, r21014, MPFR_RNDN); return mpfr_get_d(r21015, MPFR_RNDN); } static mpfr_t r21016, r21017, r21018, r21019, r21020, r21021, r21022, r21023, r21024, r21025, r21026, r21027, r21028, r21029, r21030, r21031, r21032, r21033, r21034, r21035; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r21016, "2", 10, MPFR_RNDN); mpfr_init(r21017); mpfr_init(r21018); mpfr_init(r21019); mpfr_init(r21020); mpfr_init(r21021); mpfr_init(r21022); mpfr_init(r21023); mpfr_init(r21024); mpfr_init_set_str(r21025, "1", 10, MPFR_RNDN); mpfr_init(r21026); mpfr_init(r21027); mpfr_init(r21028); mpfr_init(r21029); mpfr_init(r21030); mpfr_init(r21031); mpfr_init(r21032); mpfr_init(r21033); mpfr_init(r21034); mpfr_init(r21035); } double f_fm(double t, double l, double k) { ; mpfr_set_d(r21017, t, MPFR_RNDN); mpfr_set_d(r21018, l, MPFR_RNDN); mpfr_div(r21019, r21017, r21018, MPFR_RNDN); mpfr_set_d(r21020, k, MPFR_RNDN); mpfr_sin(r21021, r21020, MPFR_RNDN); mpfr_mul(r21022, r21017, r21021, MPFR_RNDN); mpfr_mul(r21023, r21019, r21022, MPFR_RNDN); mpfr_mul(r21024, r21023, r21021, MPFR_RNDN); ; mpfr_div(r21026, r21020, r21017, MPFR_RNDN); mpfr_pow(r21027, r21026, r21016, MPFR_RNDN); mpfr_add(r21028, r21025, r21027, MPFR_RNDN); mpfr_add(r21029, r21028, r21025, MPFR_RNDN); mpfr_mul(r21030, r21024, r21029, MPFR_RNDN); mpfr_div(r21031, r21018, r21017, MPFR_RNDN); mpfr_cos(r21032, r21020, MPFR_RNDN); mpfr_mul(r21033, r21031, r21032, MPFR_RNDN); mpfr_div(r21034, r21030, r21033, MPFR_RNDN); mpfr_div(r21035, r21016, r21034, MPFR_RNDN); return mpfr_get_d(r21035, MPFR_RNDN); } static mpfr_t r21036, r21037, r21038, r21039, r21040, r21041, r21042, r21043, r21044, r21045, r21046, r21047, r21048, r21049, r21050, r21051, r21052, r21053, r21054, r21055; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r21036, "2", 10, MPFR_RNDN); mpfr_init(r21037); mpfr_init(r21038); mpfr_init(r21039); mpfr_init(r21040); mpfr_init(r21041); mpfr_init(r21042); mpfr_init(r21043); mpfr_init(r21044); mpfr_init_set_str(r21045, "1", 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); } double f_dm(double t, double l, double k) { ; mpfr_set_d(r21037, t, MPFR_RNDN); mpfr_set_d(r21038, l, MPFR_RNDN); mpfr_div(r21039, r21037, r21038, MPFR_RNDN); mpfr_set_d(r21040, k, MPFR_RNDN); mpfr_sin(r21041, r21040, MPFR_RNDN); mpfr_mul(r21042, r21037, r21041, MPFR_RNDN); mpfr_mul(r21043, r21039, r21042, MPFR_RNDN); mpfr_mul(r21044, r21043, r21041, MPFR_RNDN); ; mpfr_div(r21046, r21040, r21037, MPFR_RNDN); mpfr_pow(r21047, r21046, r21036, MPFR_RNDN); mpfr_add(r21048, r21045, r21047, MPFR_RNDN); mpfr_add(r21049, r21048, r21045, MPFR_RNDN); mpfr_mul(r21050, r21044, r21049, MPFR_RNDN); mpfr_div(r21051, r21038, r21037, MPFR_RNDN); mpfr_cos(r21052, r21040, MPFR_RNDN); mpfr_mul(r21053, r21051, r21052, MPFR_RNDN); mpfr_div(r21054, r21050, r21053, MPFR_RNDN); mpfr_div(r21055, r21036, r21054, MPFR_RNDN); return mpfr_get_d(r21055, MPFR_RNDN); }