#include #include #include #include #include char *name = "Octave 3.8, jcobi/2"; double f_if(float alpha, float beta, float i) { float r23979 = alpha; float r23980 = beta; float r23981 = r23979 + r23980; float r23982 = r23980 - r23979; float r23983 = r23981 * r23982; float r23984 = 2; float r23985 = i; float r23986 = r23984 * r23985; float r23987 = r23981 + r23986; float r23988 = r23983 / r23987; float r23989 = 2.0; float r23990 = r23987 + r23989; float r23991 = r23988 / r23990; float r23992 = 1.0; float r23993 = r23991 + r23992; float r23994 = r23993 / r23989; return r23994; } double f_id(double alpha, double beta, double i) { double r23995 = alpha; double r23996 = beta; double r23997 = r23995 + r23996; double r23998 = r23996 - r23995; double r23999 = r23997 * r23998; double r24000 = 2; double r24001 = i; double r24002 = r24000 * r24001; double r24003 = r23997 + r24002; double r24004 = r23999 / r24003; double r24005 = 2.0; double r24006 = r24003 + r24005; double r24007 = r24004 / r24006; double r24008 = 1.0; double r24009 = r24007 + r24008; double r24010 = r24009 / r24005; return r24010; } double f_of(float alpha, float beta, float i) { float r24011 = beta; float r24012 = alpha; float r24013 = r24011 + r24012; float r24014 = 2; float r24015 = i; float r24016 = fma(r24014, r24015, r24011); float r24017 = 2.0; float r24018 = r24017 + r24012; float r24019 = r24016 + r24018; float r24020 = r24013 / r24019; float r24021 = fma(r24015, r24014, r24013); float r24022 = cbrt(r24021); float r24023 = r24020 / r24022; float r24024 = r24011 - r24012; float r24025 = cbrt(r24024); float r24026 = r24025 * r24025; float r24027 = r24026 / r24022; float r24028 = r24023 * r24027; float r24029 = fma(r24014, r24015, r24013); float r24030 = cbrt(r24029); float r24031 = r24025 / r24030; float r24032 = r24028 * r24031; float r24033 = 1.0; float r24034 = r24032 + r24033; float r24035 = r24034 / r24017; return r24035; } double f_od(double alpha, double beta, double i) { double r24036 = beta; double r24037 = alpha; double r24038 = r24036 + r24037; double r24039 = 2; double r24040 = i; double r24041 = fma(r24039, r24040, r24036); double r24042 = 2.0; double r24043 = r24042 + r24037; double r24044 = r24041 + r24043; double r24045 = r24038 / r24044; double r24046 = fma(r24040, r24039, r24038); double r24047 = cbrt(r24046); double r24048 = r24045 / r24047; double r24049 = r24036 - r24037; double r24050 = cbrt(r24049); double r24051 = r24050 * r24050; double r24052 = r24051 / r24047; double r24053 = r24048 * r24052; double r24054 = fma(r24039, r24040, r24038); double r24055 = cbrt(r24054); double r24056 = r24050 / r24055; double r24057 = r24053 * r24056; double r24058 = 1.0; double r24059 = r24057 + r24058; double r24060 = r24059 / r24042; return r24060; } 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 r24061, r24062, r24063, r24064, r24065, r24066, r24067, r24068, r24069, r24070, r24071, r24072, r24073, r24074, r24075, r24076; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r24061); mpfr_init(r24062); mpfr_init(r24063); mpfr_init(r24064); mpfr_init(r24065); mpfr_init_set_str(r24066, "2", 10, MPFR_RNDN); mpfr_init(r24067); mpfr_init(r24068); mpfr_init(r24069); mpfr_init(r24070); mpfr_init_set_str(r24071, "2.0", 10, MPFR_RNDN); mpfr_init(r24072); mpfr_init(r24073); mpfr_init_set_str(r24074, "1.0", 10, MPFR_RNDN); mpfr_init(r24075); mpfr_init(r24076); } double f_im(double alpha, double beta, double i) { mpfr_set_d(r24061, alpha, MPFR_RNDN); mpfr_set_d(r24062, beta, MPFR_RNDN); mpfr_add(r24063, r24061, r24062, MPFR_RNDN); mpfr_sub(r24064, r24062, r24061, MPFR_RNDN); mpfr_mul(r24065, r24063, r24064, MPFR_RNDN); ; mpfr_set_d(r24067, i, MPFR_RNDN); mpfr_mul(r24068, r24066, r24067, MPFR_RNDN); mpfr_add(r24069, r24063, r24068, MPFR_RNDN); mpfr_div(r24070, r24065, r24069, MPFR_RNDN); ; mpfr_add(r24072, r24069, r24071, MPFR_RNDN); mpfr_div(r24073, r24070, r24072, MPFR_RNDN); ; mpfr_add(r24075, r24073, r24074, MPFR_RNDN); mpfr_div(r24076, r24075, r24071, MPFR_RNDN); return mpfr_get_d(r24076, MPFR_RNDN); } static mpfr_t r24077, r24078, r24079, r24080, r24081, r24082, r24083, r24084, r24085, r24086, r24087, r24088, r24089, r24090, r24091, r24092, r24093, r24094, r24095, r24096, r24097, r24098, r24099, r24100, r24101; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r24077); mpfr_init(r24078); mpfr_init(r24079); mpfr_init_set_str(r24080, "2", 10, MPFR_RNDN); mpfr_init(r24081); mpfr_init(r24082); mpfr_init_set_str(r24083, "2.0", 10, MPFR_RNDN); mpfr_init(r24084); mpfr_init(r24085); mpfr_init(r24086); mpfr_init(r24087); mpfr_init(r24088); mpfr_init(r24089); mpfr_init(r24090); mpfr_init(r24091); mpfr_init(r24092); mpfr_init(r24093); mpfr_init(r24094); mpfr_init(r24095); mpfr_init(r24096); mpfr_init(r24097); mpfr_init(r24098); mpfr_init_set_str(r24099, "1.0", 10, MPFR_RNDN); mpfr_init(r24100); mpfr_init(r24101); } double f_fm(double alpha, double beta, double i) { mpfr_set_d(r24077, beta, MPFR_RNDN); mpfr_set_d(r24078, alpha, MPFR_RNDN); mpfr_add(r24079, r24077, r24078, MPFR_RNDN); ; mpfr_set_d(r24081, i, MPFR_RNDN); mpfr_fma(r24082, r24080, r24081, r24077, MPFR_RNDN); ; mpfr_add(r24084, r24083, r24078, MPFR_RNDN); mpfr_add(r24085, r24082, r24084, MPFR_RNDN); mpfr_div(r24086, r24079, r24085, MPFR_RNDN); mpfr_fma(r24087, r24081, r24080, r24079, MPFR_RNDN); mpfr_cbrt(r24088, r24087, MPFR_RNDN); mpfr_div(r24089, r24086, r24088, MPFR_RNDN); mpfr_sub(r24090, r24077, r24078, MPFR_RNDN); mpfr_cbrt(r24091, r24090, MPFR_RNDN); mpfr_mul(r24092, r24091, r24091, MPFR_RNDN); mpfr_div(r24093, r24092, r24088, MPFR_RNDN); mpfr_mul(r24094, r24089, r24093, MPFR_RNDN); mpfr_fma(r24095, r24080, r24081, r24079, MPFR_RNDN); mpfr_cbrt(r24096, r24095, MPFR_RNDN); mpfr_div(r24097, r24091, r24096, MPFR_RNDN); mpfr_mul(r24098, r24094, r24097, MPFR_RNDN); ; mpfr_add(r24100, r24098, r24099, MPFR_RNDN); mpfr_div(r24101, r24100, r24083, MPFR_RNDN); return mpfr_get_d(r24101, MPFR_RNDN); } static mpfr_t r24102, r24103, r24104, r24105, r24106, r24107, r24108, r24109, r24110, r24111, r24112, r24113, r24114, r24115, r24116, r24117, r24118, r24119, r24120, r24121, r24122, r24123, r24124, r24125, r24126; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r24102); mpfr_init(r24103); mpfr_init(r24104); mpfr_init_set_str(r24105, "2", 10, MPFR_RNDN); mpfr_init(r24106); mpfr_init(r24107); mpfr_init_set_str(r24108, "2.0", 10, MPFR_RNDN); mpfr_init(r24109); mpfr_init(r24110); mpfr_init(r24111); mpfr_init(r24112); mpfr_init(r24113); mpfr_init(r24114); mpfr_init(r24115); mpfr_init(r24116); mpfr_init(r24117); mpfr_init(r24118); mpfr_init(r24119); mpfr_init(r24120); mpfr_init(r24121); mpfr_init(r24122); mpfr_init(r24123); mpfr_init_set_str(r24124, "1.0", 10, MPFR_RNDN); mpfr_init(r24125); mpfr_init(r24126); } double f_dm(double alpha, double beta, double i) { mpfr_set_d(r24102, beta, MPFR_RNDN); mpfr_set_d(r24103, alpha, MPFR_RNDN); mpfr_add(r24104, r24102, r24103, MPFR_RNDN); ; mpfr_set_d(r24106, i, MPFR_RNDN); mpfr_fma(r24107, r24105, r24106, r24102, MPFR_RNDN); ; mpfr_add(r24109, r24108, r24103, MPFR_RNDN); mpfr_add(r24110, r24107, r24109, MPFR_RNDN); mpfr_div(r24111, r24104, r24110, MPFR_RNDN); mpfr_fma(r24112, r24106, r24105, r24104, MPFR_RNDN); mpfr_cbrt(r24113, r24112, MPFR_RNDN); mpfr_div(r24114, r24111, r24113, MPFR_RNDN); mpfr_sub(r24115, r24102, r24103, MPFR_RNDN); mpfr_cbrt(r24116, r24115, MPFR_RNDN); mpfr_mul(r24117, r24116, r24116, MPFR_RNDN); mpfr_div(r24118, r24117, r24113, MPFR_RNDN); mpfr_mul(r24119, r24114, r24118, MPFR_RNDN); mpfr_fma(r24120, r24105, r24106, r24104, MPFR_RNDN); mpfr_cbrt(r24121, r24120, MPFR_RNDN); mpfr_div(r24122, r24116, r24121, MPFR_RNDN); mpfr_mul(r24123, r24119, r24122, MPFR_RNDN); ; mpfr_add(r24125, r24123, r24124, MPFR_RNDN); mpfr_div(r24126, r24125, r24108, MPFR_RNDN); return mpfr_get_d(r24126, MPFR_RNDN); }