#include #include #include #include #include char *name = "Octave 3.8, jcobi/3"; double f_if(float alpha, float beta) { float r26141 = alpha; float r26142 = beta; float r26143 = r26141 + r26142; float r26144 = r26142 * r26141; float r26145 = r26143 + r26144; float r26146 = 1.0; float r26147 = r26145 + r26146; float r26148 = 2; float r26149 = 1; float r26150 = r26148 * r26149; float r26151 = r26143 + r26150; float r26152 = r26147 / r26151; float r26153 = r26152 / r26151; float r26154 = r26151 + r26146; float r26155 = r26153 / r26154; return r26155; } double f_id(double alpha, double beta) { double r26156 = alpha; double r26157 = beta; double r26158 = r26156 + r26157; double r26159 = r26157 * r26156; double r26160 = r26158 + r26159; double r26161 = 1.0; double r26162 = r26160 + r26161; double r26163 = 2; double r26164 = 1; double r26165 = r26163 * r26164; double r26166 = r26158 + r26165; double r26167 = r26162 / r26166; double r26168 = r26167 / r26166; double r26169 = r26166 + r26161; double r26170 = r26168 / r26169; return r26170; } double f_of(float alpha, float beta) { float r26171 = beta; float r26172 = alpha; float r26173 = r26171 * r26172; float r26174 = r26171 + r26172; float r26175 = r26173 + r26174; float r26176 = 1.0; float r26177 = r26175 + r26176; float r26178 = 2; float r26179 = r26174 + r26178; float r26180 = r26177 / r26179; float r26181 = r26180 / r26179; float r26182 = r26176 + r26179; float r26183 = r26181 / r26182; float r26184 = 1.7835983687668721e+308; bool r26185 = r26183 <= r26184; float r26186 = 0.25; float r26187 = r26186 * r26174; float r26188 = 0.5; float r26189 = r26187 + r26188; float r26190 = r26176 + r26178; float r26191 = r26174 + r26190; float r26192 = r26179 * r26191; float r26193 = r26189 / r26192; float r26194 = r26185 ? r26183 : r26193; return r26194; } double f_od(double alpha, double beta) { double r26195 = beta; double r26196 = alpha; double r26197 = r26195 * r26196; double r26198 = r26195 + r26196; double r26199 = r26197 + r26198; double r26200 = 1.0; double r26201 = r26199 + r26200; double r26202 = 2; double r26203 = r26198 + r26202; double r26204 = r26201 / r26203; double r26205 = r26204 / r26203; double r26206 = r26200 + r26203; double r26207 = r26205 / r26206; double r26208 = 1.7835983687668721e+308; bool r26209 = r26207 <= r26208; double r26210 = 0.25; double r26211 = r26210 * r26198; double r26212 = 0.5; double r26213 = r26211 + r26212; double r26214 = r26200 + r26202; double r26215 = r26198 + r26214; double r26216 = r26203 * r26215; double r26217 = r26213 / r26216; double r26218 = r26209 ? r26207 : r26217; return r26218; } 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 r26219, r26220, r26221, r26222, r26223, r26224, r26225, r26226, r26227, r26228, r26229, r26230, r26231, r26232, r26233; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init(r26219); mpfr_init(r26220); mpfr_init(r26221); mpfr_init(r26222); mpfr_init(r26223); mpfr_init_set_str(r26224, "1.0", 10, MPFR_RNDN); mpfr_init(r26225); mpfr_init_set_str(r26226, "2", 10, MPFR_RNDN); mpfr_init_set_str(r26227, "1", 10, MPFR_RNDN); mpfr_init(r26228); mpfr_init(r26229); mpfr_init(r26230); mpfr_init(r26231); mpfr_init(r26232); mpfr_init(r26233); } double f_im(double alpha, double beta) { mpfr_set_d(r26219, alpha, MPFR_RNDN); mpfr_set_d(r26220, beta, MPFR_RNDN); mpfr_add(r26221, r26219, r26220, MPFR_RNDN); mpfr_mul(r26222, r26220, r26219, MPFR_RNDN); mpfr_add(r26223, r26221, r26222, MPFR_RNDN); ; mpfr_add(r26225, r26223, r26224, MPFR_RNDN); ; ; mpfr_mul(r26228, r26226, r26227, MPFR_RNDN); mpfr_add(r26229, r26221, r26228, MPFR_RNDN); mpfr_div(r26230, r26225, r26229, MPFR_RNDN); mpfr_div(r26231, r26230, r26229, MPFR_RNDN); mpfr_add(r26232, r26229, r26224, MPFR_RNDN); mpfr_div(r26233, r26231, r26232, MPFR_RNDN); return mpfr_get_d(r26233, MPFR_RNDN); } static mpfr_t r26234, r26235, r26236, r26237, r26238, r26239, r26240, r26241, r26242, r26243, r26244, r26245, r26246, r26247, r26248, r26249, r26250, r26251, r26252, r26253, r26254, r26255, r26256, r26257; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init(r26234); mpfr_init(r26235); mpfr_init(r26236); mpfr_init(r26237); mpfr_init(r26238); mpfr_init_set_str(r26239, "1.0", 10, MPFR_RNDN); mpfr_init(r26240); mpfr_init_set_str(r26241, "2", 10, MPFR_RNDN); mpfr_init(r26242); mpfr_init(r26243); mpfr_init(r26244); mpfr_init(r26245); mpfr_init(r26246); mpfr_init_set_str(r26247, "1.7835983687668721e+308", 10, MPFR_RNDN); mpfr_init(r26248); mpfr_init_set_str(r26249, "0.25", 10, MPFR_RNDN); mpfr_init(r26250); mpfr_init_set_str(r26251, "0.5", 10, MPFR_RNDN); mpfr_init(r26252); mpfr_init(r26253); mpfr_init(r26254); mpfr_init(r26255); mpfr_init(r26256); mpfr_init(r26257); } double f_fm(double alpha, double beta) { mpfr_set_d(r26234, beta, MPFR_RNDN); mpfr_set_d(r26235, alpha, MPFR_RNDN); mpfr_mul(r26236, r26234, r26235, MPFR_RNDN); mpfr_add(r26237, r26234, r26235, MPFR_RNDN); mpfr_add(r26238, r26236, r26237, MPFR_RNDN); ; mpfr_add(r26240, r26238, r26239, MPFR_RNDN); ; mpfr_add(r26242, r26237, r26241, MPFR_RNDN); mpfr_div(r26243, r26240, r26242, MPFR_RNDN); mpfr_div(r26244, r26243, r26242, MPFR_RNDN); mpfr_add(r26245, r26239, r26242, MPFR_RNDN); mpfr_div(r26246, r26244, r26245, MPFR_RNDN); ; mpfr_set_si(r26248, mpfr_cmp(r26246, r26247) <= 0, MPFR_RNDN); ; mpfr_mul(r26250, r26249, r26237, MPFR_RNDN); ; mpfr_add(r26252, r26250, r26251, MPFR_RNDN); mpfr_add(r26253, r26239, r26241, MPFR_RNDN); mpfr_add(r26254, r26237, r26253, MPFR_RNDN); mpfr_mul(r26255, r26242, r26254, MPFR_RNDN); mpfr_div(r26256, r26252, r26255, MPFR_RNDN); if (mpfr_get_si(r26248, MPFR_RNDN)) { mpfr_set(r26257, r26246, MPFR_RNDN); } else { mpfr_set(r26257, r26256, MPFR_RNDN); }; return mpfr_get_d(r26257, MPFR_RNDN); } static mpfr_t r26258, r26259, r26260, r26261, r26262, r26263, r26264, r26265, r26266, r26267, r26268, r26269, r26270, r26271, r26272, r26273, r26274, r26275, r26276, r26277, r26278, r26279, r26280, r26281; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init(r26258); mpfr_init(r26259); mpfr_init(r26260); mpfr_init(r26261); mpfr_init(r26262); mpfr_init_set_str(r26263, "1.0", 10, MPFR_RNDN); mpfr_init(r26264); mpfr_init_set_str(r26265, "2", 10, MPFR_RNDN); mpfr_init(r26266); mpfr_init(r26267); mpfr_init(r26268); mpfr_init(r26269); mpfr_init(r26270); mpfr_init_set_str(r26271, "1.7835983687668721e+308", 10, MPFR_RNDN); mpfr_init(r26272); mpfr_init_set_str(r26273, "0.25", 10, MPFR_RNDN); mpfr_init(r26274); mpfr_init_set_str(r26275, "0.5", 10, MPFR_RNDN); mpfr_init(r26276); mpfr_init(r26277); mpfr_init(r26278); mpfr_init(r26279); mpfr_init(r26280); mpfr_init(r26281); } double f_dm(double alpha, double beta) { mpfr_set_d(r26258, beta, MPFR_RNDN); mpfr_set_d(r26259, alpha, MPFR_RNDN); mpfr_mul(r26260, r26258, r26259, MPFR_RNDN); mpfr_add(r26261, r26258, r26259, MPFR_RNDN); mpfr_add(r26262, r26260, r26261, MPFR_RNDN); ; mpfr_add(r26264, r26262, r26263, MPFR_RNDN); ; mpfr_add(r26266, r26261, r26265, MPFR_RNDN); mpfr_div(r26267, r26264, r26266, MPFR_RNDN); mpfr_div(r26268, r26267, r26266, MPFR_RNDN); mpfr_add(r26269, r26263, r26266, MPFR_RNDN); mpfr_div(r26270, r26268, r26269, MPFR_RNDN); ; mpfr_set_si(r26272, mpfr_cmp(r26270, r26271) <= 0, MPFR_RNDN); ; mpfr_mul(r26274, r26273, r26261, MPFR_RNDN); ; mpfr_add(r26276, r26274, r26275, MPFR_RNDN); mpfr_add(r26277, r26263, r26265, MPFR_RNDN); mpfr_add(r26278, r26261, r26277, MPFR_RNDN); mpfr_mul(r26279, r26266, r26278, MPFR_RNDN); mpfr_div(r26280, r26276, r26279, MPFR_RNDN); if (mpfr_get_si(r26272, MPFR_RNDN)) { mpfr_set(r26281, r26270, MPFR_RNDN); } else { mpfr_set(r26281, r26280, MPFR_RNDN); }; return mpfr_get_d(r26281, MPFR_RNDN); }