#include #include #include #include #include char *name = "Octave 3.8, jcobi/4"; double f_if(float alpha, float beta, float i) { float r26107 = i; float r26108 = alpha; float r26109 = beta; float r26110 = r26108 + r26109; float r26111 = r26110 + r26107; float r26112 = r26107 * r26111; float r26113 = r26109 * r26108; float r26114 = r26113 + r26112; float r26115 = r26112 * r26114; float r26116 = 2; float r26117 = r26116 * r26107; float r26118 = r26110 + r26117; float r26119 = r26118 * r26118; float r26120 = r26115 / r26119; float r26121 = 1.0; float r26122 = r26119 - r26121; float r26123 = r26120 / r26122; return r26123; } double f_id(double alpha, double beta, double i) { double r26124 = i; double r26125 = alpha; double r26126 = beta; double r26127 = r26125 + r26126; double r26128 = r26127 + r26124; double r26129 = r26124 * r26128; double r26130 = r26126 * r26125; double r26131 = r26130 + r26129; double r26132 = r26129 * r26131; double r26133 = 2; double r26134 = r26133 * r26124; double r26135 = r26127 + r26134; double r26136 = r26135 * r26135; double r26137 = r26132 / r26136; double r26138 = 1.0; double r26139 = r26136 - r26138; double r26140 = r26137 / r26139; return r26140; } double f_of(float alpha, float beta, float i) { float r26141 = alpha; float r26142 = beta; float r26143 = r26141 + r26142; float r26144 = 5.893473026189091e+153; bool r26145 = r26143 <= r26144; float r26146 = i; float r26147 = r26143 + r26146; float r26148 = r26146 * r26147; float r26149 = 2; float r26150 = r26149 * r26146; float r26151 = r26143 + r26150; float r26152 = r26148 / r26151; float r26153 = r26146 + r26141; float r26154 = r26142 + r26146; float r26155 = r26153 * r26154; float r26156 = r26150 + r26143; float r26157 = r26155 / r26156; float r26158 = r26152 * r26157; float r26159 = r26151 * r26151; float r26160 = 1.0; float r26161 = r26159 - r26160; float r26162 = r26158 / r26161; float r26163 = 0; float r26164 = r26145 ? r26162 : r26163; return r26164; } double f_od(double alpha, double beta, double i) { double r26165 = alpha; double r26166 = beta; double r26167 = r26165 + r26166; double r26168 = 5.893473026189091e+153; bool r26169 = r26167 <= r26168; double r26170 = i; double r26171 = r26167 + r26170; double r26172 = r26170 * r26171; double r26173 = 2; double r26174 = r26173 * r26170; double r26175 = r26167 + r26174; double r26176 = r26172 / r26175; double r26177 = r26170 + r26165; double r26178 = r26166 + r26170; double r26179 = r26177 * r26178; double r26180 = r26174 + r26167; double r26181 = r26179 / r26180; double r26182 = r26176 * r26181; double r26183 = r26175 * r26175; double r26184 = 1.0; double r26185 = r26183 - r26184; double r26186 = r26182 / r26185; double r26187 = 0; double r26188 = r26169 ? r26186 : r26187; return r26188; } 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 r26189, r26190, r26191, r26192, r26193, r26194, r26195, r26196, r26197, r26198, r26199, r26200, r26201, r26202, r26203, r26204, r26205; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r26189); mpfr_init(r26190); mpfr_init(r26191); mpfr_init(r26192); mpfr_init(r26193); mpfr_init(r26194); mpfr_init(r26195); mpfr_init(r26196); mpfr_init(r26197); mpfr_init_set_str(r26198, "2", 10, MPFR_RNDN); mpfr_init(r26199); mpfr_init(r26200); mpfr_init(r26201); mpfr_init(r26202); mpfr_init_set_str(r26203, "1.0", 10, MPFR_RNDN); mpfr_init(r26204); mpfr_init(r26205); } double f_im(double alpha, double beta, double i) { mpfr_set_d(r26189, i, MPFR_RNDN); mpfr_set_d(r26190, alpha, MPFR_RNDN); mpfr_set_d(r26191, beta, MPFR_RNDN); mpfr_add(r26192, r26190, r26191, MPFR_RNDN); mpfr_add(r26193, r26192, r26189, MPFR_RNDN); mpfr_mul(r26194, r26189, r26193, MPFR_RNDN); mpfr_mul(r26195, r26191, r26190, MPFR_RNDN); mpfr_add(r26196, r26195, r26194, MPFR_RNDN); mpfr_mul(r26197, r26194, r26196, MPFR_RNDN); ; mpfr_mul(r26199, r26198, r26189, MPFR_RNDN); mpfr_add(r26200, r26192, r26199, MPFR_RNDN); mpfr_mul(r26201, r26200, r26200, MPFR_RNDN); mpfr_div(r26202, r26197, r26201, MPFR_RNDN); ; mpfr_sub(r26204, r26201, r26203, MPFR_RNDN); mpfr_div(r26205, r26202, r26204, MPFR_RNDN); return mpfr_get_d(r26205, MPFR_RNDN); } static mpfr_t r26206, r26207, r26208, r26209, r26210, r26211, r26212, r26213, r26214, r26215, r26216, r26217, r26218, r26219, r26220, r26221, r26222, r26223, r26224, r26225, r26226, r26227, r26228, r26229; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r26206); mpfr_init(r26207); mpfr_init(r26208); mpfr_init_set_str(r26209, "5.893473026189091e+153", 10, MPFR_RNDN); mpfr_init(r26210); mpfr_init(r26211); mpfr_init(r26212); mpfr_init(r26213); mpfr_init_set_str(r26214, "2", 10, MPFR_RNDN); mpfr_init(r26215); mpfr_init(r26216); mpfr_init(r26217); mpfr_init(r26218); mpfr_init(r26219); mpfr_init(r26220); mpfr_init(r26221); mpfr_init(r26222); mpfr_init(r26223); mpfr_init(r26224); mpfr_init_set_str(r26225, "1.0", 10, MPFR_RNDN); mpfr_init(r26226); mpfr_init(r26227); mpfr_init_set_str(r26228, "0", 10, MPFR_RNDN); mpfr_init(r26229); } double f_fm(double alpha, double beta, double i) { mpfr_set_d(r26206, alpha, MPFR_RNDN); mpfr_set_d(r26207, beta, MPFR_RNDN); mpfr_add(r26208, r26206, r26207, MPFR_RNDN); ; mpfr_set_si(r26210, mpfr_cmp(r26208, r26209) <= 0, MPFR_RNDN); mpfr_set_d(r26211, i, MPFR_RNDN); mpfr_add(r26212, r26208, r26211, MPFR_RNDN); mpfr_mul(r26213, r26211, r26212, MPFR_RNDN); ; mpfr_mul(r26215, r26214, r26211, MPFR_RNDN); mpfr_add(r26216, r26208, r26215, MPFR_RNDN); mpfr_div(r26217, r26213, r26216, MPFR_RNDN); mpfr_add(r26218, r26211, r26206, MPFR_RNDN); mpfr_add(r26219, r26207, r26211, MPFR_RNDN); mpfr_mul(r26220, r26218, r26219, MPFR_RNDN); mpfr_add(r26221, r26215, r26208, MPFR_RNDN); mpfr_div(r26222, r26220, r26221, MPFR_RNDN); mpfr_mul(r26223, r26217, r26222, MPFR_RNDN); mpfr_mul(r26224, r26216, r26216, MPFR_RNDN); ; mpfr_sub(r26226, r26224, r26225, MPFR_RNDN); mpfr_div(r26227, r26223, r26226, MPFR_RNDN); ; if (mpfr_get_si(r26210, MPFR_RNDN)) { mpfr_set(r26229, r26227, MPFR_RNDN); } else { mpfr_set(r26229, r26228, MPFR_RNDN); }; return mpfr_get_d(r26229, MPFR_RNDN); } static mpfr_t r26230, r26231, r26232, r26233, r26234, r26235, r26236, r26237, r26238, r26239, r26240, r26241, r26242, r26243, r26244, r26245, r26246, r26247, r26248, r26249, r26250, r26251, r26252, r26253; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r26230); mpfr_init(r26231); mpfr_init(r26232); mpfr_init_set_str(r26233, "5.893473026189091e+153", 10, MPFR_RNDN); mpfr_init(r26234); mpfr_init(r26235); mpfr_init(r26236); mpfr_init(r26237); mpfr_init_set_str(r26238, "2", 10, MPFR_RNDN); mpfr_init(r26239); mpfr_init(r26240); mpfr_init(r26241); mpfr_init(r26242); mpfr_init(r26243); mpfr_init(r26244); mpfr_init(r26245); mpfr_init(r26246); mpfr_init(r26247); mpfr_init(r26248); mpfr_init_set_str(r26249, "1.0", 10, MPFR_RNDN); mpfr_init(r26250); mpfr_init(r26251); mpfr_init_set_str(r26252, "0", 10, MPFR_RNDN); mpfr_init(r26253); } double f_dm(double alpha, double beta, double i) { mpfr_set_d(r26230, alpha, MPFR_RNDN); mpfr_set_d(r26231, beta, MPFR_RNDN); mpfr_add(r26232, r26230, r26231, MPFR_RNDN); ; mpfr_set_si(r26234, mpfr_cmp(r26232, r26233) <= 0, MPFR_RNDN); mpfr_set_d(r26235, i, MPFR_RNDN); mpfr_add(r26236, r26232, r26235, MPFR_RNDN); mpfr_mul(r26237, r26235, r26236, MPFR_RNDN); ; mpfr_mul(r26239, r26238, r26235, MPFR_RNDN); mpfr_add(r26240, r26232, r26239, MPFR_RNDN); mpfr_div(r26241, r26237, r26240, MPFR_RNDN); mpfr_add(r26242, r26235, r26230, MPFR_RNDN); mpfr_add(r26243, r26231, r26235, MPFR_RNDN); mpfr_mul(r26244, r26242, r26243, MPFR_RNDN); mpfr_add(r26245, r26239, r26232, MPFR_RNDN); mpfr_div(r26246, r26244, r26245, MPFR_RNDN); mpfr_mul(r26247, r26241, r26246, MPFR_RNDN); mpfr_mul(r26248, r26240, r26240, MPFR_RNDN); ; mpfr_sub(r26250, r26248, r26249, MPFR_RNDN); mpfr_div(r26251, r26247, r26250, MPFR_RNDN); ; if (mpfr_get_si(r26234, MPFR_RNDN)) { mpfr_set(r26253, r26251, MPFR_RNDN); } else { mpfr_set(r26253, r26252, MPFR_RNDN); }; return mpfr_get_d(r26253, MPFR_RNDN); }