#include #include #include #include #include char *name = "Octave 3.8, jcobi/3"; double f_if(float alpha, float beta) { float r23944 = alpha; float r23945 = beta; float r23946 = r23944 + r23945; float r23947 = r23945 * r23944; float r23948 = r23946 + r23947; float r23949 = 1.0; float r23950 = r23948 + r23949; float r23951 = 2; float r23952 = 1; float r23953 = r23951 * r23952; float r23954 = r23946 + r23953; float r23955 = r23950 / r23954; float r23956 = r23955 / r23954; float r23957 = r23954 + r23949; float r23958 = r23956 / r23957; return r23958; } double f_id(double alpha, double beta) { double r23959 = alpha; double r23960 = beta; double r23961 = r23959 + r23960; double r23962 = r23960 * r23959; double r23963 = r23961 + r23962; double r23964 = 1.0; double r23965 = r23963 + r23964; double r23966 = 2; double r23967 = 1; double r23968 = r23966 * r23967; double r23969 = r23961 + r23968; double r23970 = r23965 / r23969; double r23971 = r23970 / r23969; double r23972 = r23969 + r23964; double r23973 = r23971 / r23972; return r23973; } double f_of(float alpha, float beta) { float r23974 = beta; float r23975 = alpha; float r23976 = r23974 + r23975; float r23977 = 9.029168664560661e+206; bool r23978 = r23976 <= r23977; float r23979 = r23974 * r23975; float r23980 = r23976 + r23979; float r23981 = 1.0; float r23982 = r23980 + r23981; float r23983 = 2; float r23984 = r23976 + r23983; float r23985 = r23982 / r23984; float r23986 = r23985 / r23984; float r23987 = r23984 + r23981; float r23988 = r23986 / r23987; float r23989 = 0.25; float r23990 = r23989 * r23976; float r23991 = 0.5; float r23992 = r23990 + r23991; float r23993 = r23983 + r23981; float r23994 = r23976 + r23993; float r23995 = r23984 * r23994; float r23996 = r23992 / r23995; float r23997 = r23978 ? r23988 : r23996; return r23997; } double f_od(double alpha, double beta) { double r23998 = beta; double r23999 = alpha; double r24000 = r23998 + r23999; double r24001 = 9.029168664560661e+206; bool r24002 = r24000 <= r24001; double r24003 = r23998 * r23999; double r24004 = r24000 + r24003; double r24005 = 1.0; double r24006 = r24004 + r24005; double r24007 = 2; double r24008 = r24000 + r24007; double r24009 = r24006 / r24008; double r24010 = r24009 / r24008; double r24011 = r24008 + r24005; double r24012 = r24010 / r24011; double r24013 = 0.25; double r24014 = r24013 * r24000; double r24015 = 0.5; double r24016 = r24014 + r24015; double r24017 = r24007 + r24005; double r24018 = r24000 + r24017; double r24019 = r24008 * r24018; double r24020 = r24016 / r24019; double r24021 = r24002 ? r24012 : r24020; return r24021; } 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 r24022, r24023, r24024, r24025, r24026, r24027, r24028, r24029, r24030, r24031, r24032, r24033, r24034, r24035, r24036; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r24022); mpfr_init(r24023); mpfr_init(r24024); mpfr_init(r24025); mpfr_init(r24026); mpfr_init_set_str(r24027, "1.0", 10, MPFR_RNDN); mpfr_init(r24028); mpfr_init_set_str(r24029, "2", 10, MPFR_RNDN); mpfr_init_set_str(r24030, "1", 10, MPFR_RNDN); mpfr_init(r24031); mpfr_init(r24032); mpfr_init(r24033); mpfr_init(r24034); mpfr_init(r24035); mpfr_init(r24036); } double f_im(double alpha, double beta) { mpfr_set_d(r24022, alpha, MPFR_RNDN); mpfr_set_d(r24023, beta, MPFR_RNDN); mpfr_add(r24024, r24022, r24023, MPFR_RNDN); mpfr_mul(r24025, r24023, r24022, MPFR_RNDN); mpfr_add(r24026, r24024, r24025, MPFR_RNDN); ; mpfr_add(r24028, r24026, r24027, MPFR_RNDN); ; ; mpfr_mul(r24031, r24029, r24030, MPFR_RNDN); mpfr_add(r24032, r24024, r24031, MPFR_RNDN); mpfr_div(r24033, r24028, r24032, MPFR_RNDN); mpfr_div(r24034, r24033, r24032, MPFR_RNDN); mpfr_add(r24035, r24032, r24027, MPFR_RNDN); mpfr_div(r24036, r24034, r24035, MPFR_RNDN); return mpfr_get_d(r24036, MPFR_RNDN); } static mpfr_t r24037, r24038, r24039, r24040, r24041, r24042, r24043, r24044, r24045, r24046, r24047, r24048, r24049, r24050, r24051, r24052, r24053, r24054, r24055, r24056, r24057, r24058, r24059, r24060; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r24037); mpfr_init(r24038); mpfr_init(r24039); mpfr_init_set_str(r24040, "9.029168664560661e+206", 10, MPFR_RNDN); mpfr_init(r24041); mpfr_init(r24042); mpfr_init(r24043); mpfr_init_set_str(r24044, "1.0", 10, MPFR_RNDN); mpfr_init(r24045); mpfr_init_set_str(r24046, "2", 10, MPFR_RNDN); mpfr_init(r24047); mpfr_init(r24048); mpfr_init(r24049); mpfr_init(r24050); mpfr_init(r24051); mpfr_init_set_str(r24052, "0.25", 10, MPFR_RNDN); mpfr_init(r24053); mpfr_init_set_str(r24054, "0.5", 10, MPFR_RNDN); mpfr_init(r24055); mpfr_init(r24056); mpfr_init(r24057); mpfr_init(r24058); mpfr_init(r24059); mpfr_init(r24060); } double f_fm(double alpha, double beta) { mpfr_set_d(r24037, beta, MPFR_RNDN); mpfr_set_d(r24038, alpha, MPFR_RNDN); mpfr_add(r24039, r24037, r24038, MPFR_RNDN); ; mpfr_set_si(r24041, mpfr_cmp(r24039, r24040) <= 0, MPFR_RNDN); mpfr_mul(r24042, r24037, r24038, MPFR_RNDN); mpfr_add(r24043, r24039, r24042, MPFR_RNDN); ; mpfr_add(r24045, r24043, r24044, MPFR_RNDN); ; mpfr_add(r24047, r24039, r24046, MPFR_RNDN); mpfr_div(r24048, r24045, r24047, MPFR_RNDN); mpfr_div(r24049, r24048, r24047, MPFR_RNDN); mpfr_add(r24050, r24047, r24044, MPFR_RNDN); mpfr_div(r24051, r24049, r24050, MPFR_RNDN); ; mpfr_mul(r24053, r24052, r24039, MPFR_RNDN); ; mpfr_add(r24055, r24053, r24054, MPFR_RNDN); mpfr_add(r24056, r24046, r24044, MPFR_RNDN); mpfr_add(r24057, r24039, r24056, MPFR_RNDN); mpfr_mul(r24058, r24047, r24057, MPFR_RNDN); mpfr_div(r24059, r24055, r24058, MPFR_RNDN); if (mpfr_get_si(r24041, MPFR_RNDN)) { mpfr_set(r24060, r24051, MPFR_RNDN); } else { mpfr_set(r24060, r24059, MPFR_RNDN); }; return mpfr_get_d(r24060, MPFR_RNDN); } static mpfr_t r24061, r24062, r24063, r24064, r24065, r24066, r24067, r24068, r24069, r24070, r24071, r24072, r24073, r24074, r24075, r24076, r24077, r24078, r24079, r24080, r24081, r24082, r24083, r24084; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r24061); mpfr_init(r24062); mpfr_init(r24063); mpfr_init_set_str(r24064, "9.029168664560661e+206", 10, MPFR_RNDN); mpfr_init(r24065); mpfr_init(r24066); mpfr_init(r24067); mpfr_init_set_str(r24068, "1.0", 10, MPFR_RNDN); mpfr_init(r24069); mpfr_init_set_str(r24070, "2", 10, MPFR_RNDN); mpfr_init(r24071); mpfr_init(r24072); mpfr_init(r24073); mpfr_init(r24074); mpfr_init(r24075); mpfr_init_set_str(r24076, "0.25", 10, MPFR_RNDN); mpfr_init(r24077); mpfr_init_set_str(r24078, "0.5", 10, MPFR_RNDN); mpfr_init(r24079); mpfr_init(r24080); mpfr_init(r24081); mpfr_init(r24082); mpfr_init(r24083); mpfr_init(r24084); } double f_dm(double alpha, double beta) { mpfr_set_d(r24061, beta, MPFR_RNDN); mpfr_set_d(r24062, alpha, MPFR_RNDN); mpfr_add(r24063, r24061, r24062, MPFR_RNDN); ; mpfr_set_si(r24065, mpfr_cmp(r24063, r24064) <= 0, MPFR_RNDN); mpfr_mul(r24066, r24061, r24062, MPFR_RNDN); mpfr_add(r24067, r24063, r24066, MPFR_RNDN); ; mpfr_add(r24069, r24067, r24068, MPFR_RNDN); ; mpfr_add(r24071, r24063, r24070, MPFR_RNDN); mpfr_div(r24072, r24069, r24071, MPFR_RNDN); mpfr_div(r24073, r24072, r24071, MPFR_RNDN); mpfr_add(r24074, r24071, r24068, MPFR_RNDN); mpfr_div(r24075, r24073, r24074, MPFR_RNDN); ; mpfr_mul(r24077, r24076, r24063, MPFR_RNDN); ; mpfr_add(r24079, r24077, r24078, MPFR_RNDN); mpfr_add(r24080, r24070, r24068, MPFR_RNDN); mpfr_add(r24081, r24063, r24080, MPFR_RNDN); mpfr_mul(r24082, r24071, r24081, MPFR_RNDN); mpfr_div(r24083, r24079, r24082, MPFR_RNDN); if (mpfr_get_si(r24065, MPFR_RNDN)) { mpfr_set(r24084, r24075, MPFR_RNDN); } else { mpfr_set(r24084, r24083, MPFR_RNDN); }; return mpfr_get_d(r24084, MPFR_RNDN); }