#include #include #include #include #include char *name = "Octave 3.8, jcobi/4"; double f_if(float alpha, float beta, float i) { float r17001 = i; float r17002 = alpha; float r17003 = beta; float r17004 = r17002 + r17003; float r17005 = r17004 + r17001; float r17006 = r17001 * r17005; float r17007 = r17003 * r17002; float r17008 = r17007 + r17006; float r17009 = r17006 * r17008; float r17010 = 2.0f; float r17011 = r17010 * r17001; float r17012 = r17004 + r17011; float r17013 = r17012 * r17012; float r17014 = r17009 / r17013; float r17015 = 1.0f; float r17016 = r17013 - r17015; float r17017 = r17014 / r17016; return r17017; } double f_id(double alpha, double beta, double i) { double r17018 = i; double r17019 = alpha; double r17020 = beta; double r17021 = r17019 + r17020; double r17022 = r17021 + r17018; double r17023 = r17018 * r17022; double r17024 = r17020 * r17019; double r17025 = r17024 + r17023; double r17026 = r17023 * r17025; double r17027 = 2.0; double r17028 = r17027 * r17018; double r17029 = r17021 + r17028; double r17030 = r17029 * r17029; double r17031 = r17026 / r17030; double r17032 = 1.0; double r17033 = r17030 - r17032; double r17034 = r17031 / r17033; return r17034; } double f_of(float alpha, float beta, float i) { float r17035 = alpha; float r17036 = 2.069742248392879e+108f; bool r17037 = r17035 <= r17036; float r17038 = i; float r17039 = beta; float r17040 = r17038 + r17035; float r17041 = r17039 + r17040; float r17042 = r17038 * r17041; float r17043 = r17039 + r17035; float r17044 = 2.0f; float r17045 = r17044 * r17038; float r17046 = r17043 + r17045; float r17047 = r17046 * r17046; float r17048 = r17035 * r17039; float r17049 = r17048 + r17042; float r17050 = r17047 / r17049; float r17051 = r17042 / r17050; float r17052 = 1.0f; float r17053 = r17047 - r17052; float r17054 = r17051 / r17053; float r17055 = 0.0f; float r17056 = r17037 ? r17054 : r17055; return r17056; } double f_od(double alpha, double beta, double i) { double r17057 = alpha; double r17058 = 2.069742248392879e+108; bool r17059 = r17057 <= r17058; double r17060 = i; double r17061 = beta; double r17062 = r17060 + r17057; double r17063 = r17061 + r17062; double r17064 = r17060 * r17063; double r17065 = r17061 + r17057; double r17066 = 2.0; double r17067 = r17066 * r17060; double r17068 = r17065 + r17067; double r17069 = r17068 * r17068; double r17070 = r17057 * r17061; double r17071 = r17070 + r17064; double r17072 = r17069 / r17071; double r17073 = r17064 / r17072; double r17074 = 1.0; double r17075 = r17069 - r17074; double r17076 = r17073 / r17075; double r17077 = 0.0; double r17078 = r17059 ? r17076 : r17077; return r17078; } 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 r17079, r17080, r17081, r17082, r17083, r17084, r17085, r17086, r17087, r17088, r17089, r17090, r17091, r17092, r17093, r17094, r17095; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r17079); mpfr_init(r17080); mpfr_init(r17081); mpfr_init(r17082); mpfr_init(r17083); mpfr_init(r17084); mpfr_init(r17085); mpfr_init(r17086); mpfr_init(r17087); mpfr_init_set_str(r17088, "2", 10, MPFR_RNDN); mpfr_init(r17089); mpfr_init(r17090); mpfr_init(r17091); mpfr_init(r17092); mpfr_init_set_str(r17093, "1.0", 10, MPFR_RNDN); mpfr_init(r17094); mpfr_init(r17095); } double f_im(double alpha, double beta, double i) { mpfr_set_d(r17079, i, MPFR_RNDN); mpfr_set_d(r17080, alpha, MPFR_RNDN); mpfr_set_d(r17081, beta, MPFR_RNDN); mpfr_add(r17082, r17080, r17081, MPFR_RNDN); mpfr_add(r17083, r17082, r17079, MPFR_RNDN); mpfr_mul(r17084, r17079, r17083, MPFR_RNDN); mpfr_mul(r17085, r17081, r17080, MPFR_RNDN); mpfr_add(r17086, r17085, r17084, MPFR_RNDN); mpfr_mul(r17087, r17084, r17086, MPFR_RNDN); ; mpfr_mul(r17089, r17088, r17079, MPFR_RNDN); mpfr_add(r17090, r17082, r17089, MPFR_RNDN); mpfr_mul(r17091, r17090, r17090, MPFR_RNDN); mpfr_div(r17092, r17087, r17091, MPFR_RNDN); ; mpfr_sub(r17094, r17091, r17093, MPFR_RNDN); mpfr_div(r17095, r17092, r17094, MPFR_RNDN); return mpfr_get_d(r17095, MPFR_RNDN); } static mpfr_t r17096, r17097, r17098, r17099, r17100, r17101, r17102, r17103, r17104, r17105, r17106, r17107, r17108, r17109, r17110, r17111, r17112, r17113, r17114, r17115, r17116, r17117; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r17096); mpfr_init_set_str(r17097, "2.069742248392879e+108", 10, MPFR_RNDN); mpfr_init(r17098); mpfr_init(r17099); mpfr_init(r17100); mpfr_init(r17101); mpfr_init(r17102); mpfr_init(r17103); mpfr_init(r17104); mpfr_init_set_str(r17105, "2", 10, MPFR_RNDN); mpfr_init(r17106); mpfr_init(r17107); mpfr_init(r17108); mpfr_init(r17109); mpfr_init(r17110); mpfr_init(r17111); mpfr_init(r17112); mpfr_init_set_str(r17113, "1.0", 10, MPFR_RNDN); mpfr_init(r17114); mpfr_init(r17115); mpfr_init_set_str(r17116, "0", 10, MPFR_RNDN); mpfr_init(r17117); } double f_fm(double alpha, double beta, double i) { mpfr_set_d(r17096, alpha, MPFR_RNDN); ; mpfr_set_si(r17098, mpfr_cmp(r17096, r17097) <= 0, MPFR_RNDN); mpfr_set_d(r17099, i, MPFR_RNDN); mpfr_set_d(r17100, beta, MPFR_RNDN); mpfr_add(r17101, r17099, r17096, MPFR_RNDN); mpfr_add(r17102, r17100, r17101, MPFR_RNDN); mpfr_mul(r17103, r17099, r17102, MPFR_RNDN); mpfr_add(r17104, r17100, r17096, MPFR_RNDN); ; mpfr_mul(r17106, r17105, r17099, MPFR_RNDN); mpfr_add(r17107, r17104, r17106, MPFR_RNDN); mpfr_sqr(r17108, r17107, MPFR_RNDN); mpfr_mul(r17109, r17096, r17100, MPFR_RNDN); mpfr_add(r17110, r17109, r17103, MPFR_RNDN); mpfr_div(r17111, r17108, r17110, MPFR_RNDN); mpfr_div(r17112, r17103, r17111, MPFR_RNDN); ; mpfr_sub(r17114, r17108, r17113, MPFR_RNDN); mpfr_div(r17115, r17112, r17114, MPFR_RNDN); ; if (mpfr_get_si(r17098, MPFR_RNDN)) { mpfr_set(r17117, r17115, MPFR_RNDN); } else { mpfr_set(r17117, r17116, MPFR_RNDN); }; return mpfr_get_d(r17117, MPFR_RNDN); } static mpfr_t r17118, r17119, r17120, r17121, r17122, r17123, r17124, r17125, r17126, r17127, r17128, r17129, r17130, r17131, r17132, r17133, r17134, r17135, r17136, r17137, r17138, r17139; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r17118); mpfr_init_set_str(r17119, "2.069742248392879e+108", 10, MPFR_RNDN); mpfr_init(r17120); mpfr_init(r17121); mpfr_init(r17122); mpfr_init(r17123); mpfr_init(r17124); mpfr_init(r17125); mpfr_init(r17126); mpfr_init_set_str(r17127, "2", 10, MPFR_RNDN); mpfr_init(r17128); mpfr_init(r17129); mpfr_init(r17130); mpfr_init(r17131); mpfr_init(r17132); mpfr_init(r17133); mpfr_init(r17134); mpfr_init_set_str(r17135, "1.0", 10, MPFR_RNDN); mpfr_init(r17136); mpfr_init(r17137); mpfr_init_set_str(r17138, "0", 10, MPFR_RNDN); mpfr_init(r17139); } double f_dm(double alpha, double beta, double i) { mpfr_set_d(r17118, alpha, MPFR_RNDN); ; mpfr_set_si(r17120, mpfr_cmp(r17118, r17119) <= 0, MPFR_RNDN); mpfr_set_d(r17121, i, MPFR_RNDN); mpfr_set_d(r17122, beta, MPFR_RNDN); mpfr_add(r17123, r17121, r17118, MPFR_RNDN); mpfr_add(r17124, r17122, r17123, MPFR_RNDN); mpfr_mul(r17125, r17121, r17124, MPFR_RNDN); mpfr_add(r17126, r17122, r17118, MPFR_RNDN); ; mpfr_mul(r17128, r17127, r17121, MPFR_RNDN); mpfr_add(r17129, r17126, r17128, MPFR_RNDN); mpfr_sqr(r17130, r17129, MPFR_RNDN); mpfr_mul(r17131, r17118, r17122, MPFR_RNDN); mpfr_add(r17132, r17131, r17125, MPFR_RNDN); mpfr_div(r17133, r17130, r17132, MPFR_RNDN); mpfr_div(r17134, r17125, r17133, MPFR_RNDN); ; mpfr_sub(r17136, r17130, r17135, MPFR_RNDN); mpfr_div(r17137, r17134, r17136, MPFR_RNDN); ; if (mpfr_get_si(r17120, MPFR_RNDN)) { mpfr_set(r17139, r17137, MPFR_RNDN); } else { mpfr_set(r17139, r17138, MPFR_RNDN); }; return mpfr_get_d(r17139, MPFR_RNDN); }