#include #include #include #include #include char *name = "Octave 3.8, jcobi/4"; double f_if(float alpha, float beta, float i) { float r25994 = i; float r25995 = alpha; float r25996 = beta; float r25997 = r25995 + r25996; float r25998 = r25997 + r25994; float r25999 = r25994 * r25998; float r26000 = r25996 * r25995; float r26001 = r26000 + r25999; float r26002 = r25999 * r26001; float r26003 = 2; float r26004 = r26003 * r25994; float r26005 = r25997 + r26004; float r26006 = r26005 * r26005; float r26007 = r26002 / r26006; float r26008 = 1.0; float r26009 = r26006 - r26008; float r26010 = r26007 / r26009; return r26010; } double f_id(double alpha, double beta, double i) { double r26011 = i; double r26012 = alpha; double r26013 = beta; double r26014 = r26012 + r26013; double r26015 = r26014 + r26011; double r26016 = r26011 * r26015; double r26017 = r26013 * r26012; double r26018 = r26017 + r26016; double r26019 = r26016 * r26018; double r26020 = 2; double r26021 = r26020 * r26011; double r26022 = r26014 + r26021; double r26023 = r26022 * r26022; double r26024 = r26019 / r26023; double r26025 = 1.0; double r26026 = r26023 - r26025; double r26027 = r26024 / r26026; return r26027; } double f_of(float alpha, float beta, float i) { float r26028 = alpha; float r26029 = beta; float r26030 = r26028 + r26029; float r26031 = 7.041788089365064e+146; bool r26032 = r26030 <= r26031; float r26033 = i; float r26034 = r26030 + r26033; float r26035 = r26033 * r26034; float r26036 = 2; float r26037 = r26036 * r26033; float r26038 = r26030 + r26037; float r26039 = r26035 / r26038; float r26040 = r26033 + r26028; float r26041 = r26029 + r26033; float r26042 = r26040 * r26041; float r26043 = r26037 + r26030; float r26044 = r26042 / r26043; float r26045 = r26039 * r26044; float r26046 = r26038 * r26038; float r26047 = 1.0; float r26048 = r26046 - r26047; float r26049 = r26045 / r26048; float r26050 = 0; float r26051 = r26032 ? r26049 : r26050; return r26051; } double f_od(double alpha, double beta, double i) { double r26052 = alpha; double r26053 = beta; double r26054 = r26052 + r26053; double r26055 = 7.041788089365064e+146; bool r26056 = r26054 <= r26055; double r26057 = i; double r26058 = r26054 + r26057; double r26059 = r26057 * r26058; double r26060 = 2; double r26061 = r26060 * r26057; double r26062 = r26054 + r26061; double r26063 = r26059 / r26062; double r26064 = r26057 + r26052; double r26065 = r26053 + r26057; double r26066 = r26064 * r26065; double r26067 = r26061 + r26054; double r26068 = r26066 / r26067; double r26069 = r26063 * r26068; double r26070 = r26062 * r26062; double r26071 = 1.0; double r26072 = r26070 - r26071; double r26073 = r26069 / r26072; double r26074 = 0; double r26075 = r26056 ? r26073 : r26074; return r26075; } 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 r26076, r26077, r26078, r26079, r26080, r26081, r26082, r26083, r26084, r26085, r26086, r26087, r26088, r26089, r26090, r26091, r26092; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init(r26076); mpfr_init(r26077); mpfr_init(r26078); mpfr_init(r26079); mpfr_init(r26080); mpfr_init(r26081); mpfr_init(r26082); mpfr_init(r26083); mpfr_init(r26084); mpfr_init_set_str(r26085, "2", 10, MPFR_RNDN); mpfr_init(r26086); mpfr_init(r26087); mpfr_init(r26088); mpfr_init(r26089); mpfr_init_set_str(r26090, "1.0", 10, MPFR_RNDN); mpfr_init(r26091); mpfr_init(r26092); } double f_im(double alpha, double beta, double i) { mpfr_set_d(r26076, i, MPFR_RNDN); mpfr_set_d(r26077, alpha, MPFR_RNDN); mpfr_set_d(r26078, beta, MPFR_RNDN); mpfr_add(r26079, r26077, r26078, MPFR_RNDN); mpfr_add(r26080, r26079, r26076, MPFR_RNDN); mpfr_mul(r26081, r26076, r26080, MPFR_RNDN); mpfr_mul(r26082, r26078, r26077, MPFR_RNDN); mpfr_add(r26083, r26082, r26081, MPFR_RNDN); mpfr_mul(r26084, r26081, r26083, MPFR_RNDN); ; mpfr_mul(r26086, r26085, r26076, MPFR_RNDN); mpfr_add(r26087, r26079, r26086, MPFR_RNDN); mpfr_mul(r26088, r26087, r26087, MPFR_RNDN); mpfr_div(r26089, r26084, r26088, MPFR_RNDN); ; mpfr_sub(r26091, r26088, r26090, MPFR_RNDN); mpfr_div(r26092, r26089, r26091, MPFR_RNDN); return mpfr_get_d(r26092, MPFR_RNDN); } static mpfr_t r26093, r26094, r26095, r26096, r26097, r26098, r26099, r26100, r26101, r26102, r26103, r26104, r26105, r26106, r26107, r26108, r26109, r26110, r26111, r26112, r26113, r26114, r26115, r26116; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init(r26093); mpfr_init(r26094); mpfr_init(r26095); mpfr_init_set_str(r26096, "7.041788089365064e+146", 10, MPFR_RNDN); mpfr_init(r26097); mpfr_init(r26098); mpfr_init(r26099); mpfr_init(r26100); mpfr_init_set_str(r26101, "2", 10, MPFR_RNDN); mpfr_init(r26102); mpfr_init(r26103); mpfr_init(r26104); mpfr_init(r26105); mpfr_init(r26106); mpfr_init(r26107); mpfr_init(r26108); mpfr_init(r26109); mpfr_init(r26110); mpfr_init(r26111); mpfr_init_set_str(r26112, "1.0", 10, MPFR_RNDN); mpfr_init(r26113); mpfr_init(r26114); mpfr_init_set_str(r26115, "0", 10, MPFR_RNDN); mpfr_init(r26116); } double f_fm(double alpha, double beta, double i) { mpfr_set_d(r26093, alpha, MPFR_RNDN); mpfr_set_d(r26094, beta, MPFR_RNDN); mpfr_add(r26095, r26093, r26094, MPFR_RNDN); ; mpfr_set_si(r26097, mpfr_cmp(r26095, r26096) <= 0, MPFR_RNDN); mpfr_set_d(r26098, i, MPFR_RNDN); mpfr_add(r26099, r26095, r26098, MPFR_RNDN); mpfr_mul(r26100, r26098, r26099, MPFR_RNDN); ; mpfr_mul(r26102, r26101, r26098, MPFR_RNDN); mpfr_add(r26103, r26095, r26102, MPFR_RNDN); mpfr_div(r26104, r26100, r26103, MPFR_RNDN); mpfr_add(r26105, r26098, r26093, MPFR_RNDN); mpfr_add(r26106, r26094, r26098, MPFR_RNDN); mpfr_mul(r26107, r26105, r26106, MPFR_RNDN); mpfr_add(r26108, r26102, r26095, MPFR_RNDN); mpfr_div(r26109, r26107, r26108, MPFR_RNDN); mpfr_mul(r26110, r26104, r26109, MPFR_RNDN); mpfr_mul(r26111, r26103, r26103, MPFR_RNDN); ; mpfr_sub(r26113, r26111, r26112, MPFR_RNDN); mpfr_div(r26114, r26110, r26113, MPFR_RNDN); ; if (mpfr_get_si(r26097, MPFR_RNDN)) { mpfr_set(r26116, r26114, MPFR_RNDN); } else { mpfr_set(r26116, r26115, MPFR_RNDN); }; return mpfr_get_d(r26116, MPFR_RNDN); } static mpfr_t r26117, r26118, r26119, r26120, r26121, r26122, r26123, r26124, r26125, r26126, r26127, r26128, r26129, r26130, r26131, r26132, r26133, r26134, r26135, r26136, r26137, r26138, r26139, r26140; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init(r26117); mpfr_init(r26118); mpfr_init(r26119); mpfr_init_set_str(r26120, "7.041788089365064e+146", 10, MPFR_RNDN); mpfr_init(r26121); mpfr_init(r26122); mpfr_init(r26123); mpfr_init(r26124); mpfr_init_set_str(r26125, "2", 10, MPFR_RNDN); mpfr_init(r26126); mpfr_init(r26127); mpfr_init(r26128); mpfr_init(r26129); mpfr_init(r26130); mpfr_init(r26131); mpfr_init(r26132); mpfr_init(r26133); mpfr_init(r26134); mpfr_init(r26135); mpfr_init_set_str(r26136, "1.0", 10, MPFR_RNDN); mpfr_init(r26137); mpfr_init(r26138); mpfr_init_set_str(r26139, "0", 10, MPFR_RNDN); mpfr_init(r26140); } double f_dm(double alpha, double beta, double i) { mpfr_set_d(r26117, alpha, MPFR_RNDN); mpfr_set_d(r26118, beta, MPFR_RNDN); mpfr_add(r26119, r26117, r26118, MPFR_RNDN); ; mpfr_set_si(r26121, mpfr_cmp(r26119, r26120) <= 0, MPFR_RNDN); mpfr_set_d(r26122, i, MPFR_RNDN); mpfr_add(r26123, r26119, r26122, MPFR_RNDN); mpfr_mul(r26124, r26122, r26123, MPFR_RNDN); ; mpfr_mul(r26126, r26125, r26122, MPFR_RNDN); mpfr_add(r26127, r26119, r26126, MPFR_RNDN); mpfr_div(r26128, r26124, r26127, MPFR_RNDN); mpfr_add(r26129, r26122, r26117, MPFR_RNDN); mpfr_add(r26130, r26118, r26122, MPFR_RNDN); mpfr_mul(r26131, r26129, r26130, MPFR_RNDN); mpfr_add(r26132, r26126, r26119, MPFR_RNDN); mpfr_div(r26133, r26131, r26132, MPFR_RNDN); mpfr_mul(r26134, r26128, r26133, MPFR_RNDN); mpfr_mul(r26135, r26127, r26127, MPFR_RNDN); ; mpfr_sub(r26137, r26135, r26136, MPFR_RNDN); mpfr_div(r26138, r26134, r26137, MPFR_RNDN); ; if (mpfr_get_si(r26121, MPFR_RNDN)) { mpfr_set(r26140, r26138, MPFR_RNDN); } else { mpfr_set(r26140, r26139, MPFR_RNDN); }; return mpfr_get_d(r26140, MPFR_RNDN); }