#include #include #include #include #include char *name = "Octave 3.8, jcobi/3"; double f_if(float alpha, float beta) { float r17132 = alpha; float r17133 = beta; float r17134 = r17132 + r17133; float r17135 = r17133 * r17132; float r17136 = r17134 + r17135; float r17137 = 1.0f; float r17138 = r17136 + r17137; float r17139 = 2.0f; float r17140 = 1.0f; float r17141 = r17139 * r17140; float r17142 = r17134 + r17141; float r17143 = r17138 / r17142; float r17144 = r17143 / r17142; float r17145 = r17142 + r17137; float r17146 = r17144 / r17145; return r17146; } double f_id(double alpha, double beta) { double r17147 = alpha; double r17148 = beta; double r17149 = r17147 + r17148; double r17150 = r17148 * r17147; double r17151 = r17149 + r17150; double r17152 = 1.0; double r17153 = r17151 + r17152; double r17154 = 2.0; double r17155 = 1.0; double r17156 = r17154 * r17155; double r17157 = r17149 + r17156; double r17158 = r17153 / r17157; double r17159 = r17158 / r17157; double r17160 = r17157 + r17152; double r17161 = r17159 / r17160; return r17161; } double f_of(float alpha, float beta) { float r17162 = alpha; float r17163 = 1.2543987852557203e+178f; bool r17164 = r17162 <= r17163; float r17165 = 1.0f; float r17166 = r17162 + r17165; float r17167 = beta; float r17168 = r17167 * r17162; float r17169 = r17167 + r17168; float r17170 = r17166 + r17169; float r17171 = 2.0f; float r17172 = r17171 + r17167; float r17173 = r17162 + r17172; float r17174 = r17170 / r17173; float r17175 = r17166 + r17172; float r17176 = r17174 / r17175; float r17177 = r17176 / r17173; float r17178 = 0.25f; float r17179 = r17162 + r17167; float r17180 = r17178 * r17179; float r17181 = 0.5f; float r17182 = r17180 + r17181; float r17183 = r17179 + r17171; float r17184 = r17171 + r17165; float r17185 = r17179 + r17184; float r17186 = r17183 * r17185; float r17187 = r17182 / r17186; float r17188 = r17164 ? r17177 : r17187; return r17188; } double f_od(double alpha, double beta) { double r17189 = alpha; double r17190 = 1.2543987852557203e+178; bool r17191 = r17189 <= r17190; double r17192 = 1.0; double r17193 = r17189 + r17192; double r17194 = beta; double r17195 = r17194 * r17189; double r17196 = r17194 + r17195; double r17197 = r17193 + r17196; double r17198 = 2.0; double r17199 = r17198 + r17194; double r17200 = r17189 + r17199; double r17201 = r17197 / r17200; double r17202 = r17193 + r17199; double r17203 = r17201 / r17202; double r17204 = r17203 / r17200; double r17205 = 0.25; double r17206 = r17189 + r17194; double r17207 = r17205 * r17206; double r17208 = 0.5; double r17209 = r17207 + r17208; double r17210 = r17206 + r17198; double r17211 = r17198 + r17192; double r17212 = r17206 + r17211; double r17213 = r17210 * r17212; double r17214 = r17209 / r17213; double r17215 = r17191 ? r17204 : r17214; return r17215; } 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 r17216, r17217, r17218, r17219, r17220, r17221, r17222, r17223, r17224, r17225, r17226, r17227, r17228, r17229, r17230; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r17216); mpfr_init(r17217); mpfr_init(r17218); mpfr_init(r17219); mpfr_init(r17220); mpfr_init_set_str(r17221, "1.0", 10, MPFR_RNDN); mpfr_init(r17222); mpfr_init_set_str(r17223, "2", 10, MPFR_RNDN); mpfr_init_set_str(r17224, "1", 10, MPFR_RNDN); mpfr_init(r17225); mpfr_init(r17226); mpfr_init(r17227); mpfr_init(r17228); mpfr_init(r17229); mpfr_init(r17230); } double f_im(double alpha, double beta) { mpfr_set_d(r17216, alpha, MPFR_RNDN); mpfr_set_d(r17217, beta, MPFR_RNDN); mpfr_add(r17218, r17216, r17217, MPFR_RNDN); mpfr_mul(r17219, r17217, r17216, MPFR_RNDN); mpfr_add(r17220, r17218, r17219, MPFR_RNDN); ; mpfr_add(r17222, r17220, r17221, MPFR_RNDN); ; ; mpfr_mul(r17225, r17223, r17224, MPFR_RNDN); mpfr_add(r17226, r17218, r17225, MPFR_RNDN); mpfr_div(r17227, r17222, r17226, MPFR_RNDN); mpfr_div(r17228, r17227, r17226, MPFR_RNDN); mpfr_add(r17229, r17226, r17221, MPFR_RNDN); mpfr_div(r17230, r17228, r17229, MPFR_RNDN); return mpfr_get_d(r17230, MPFR_RNDN); } static mpfr_t r17231, r17232, r17233, r17234, r17235, r17236, r17237, r17238, r17239, r17240, r17241, r17242, r17243, r17244, r17245, r17246, r17247, r17248, r17249, r17250, r17251, r17252, r17253, r17254, r17255, r17256, r17257; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r17231); mpfr_init_set_str(r17232, "1.2543987852557203e+178", 10, MPFR_RNDN); mpfr_init(r17233); mpfr_init_set_str(r17234, "1.0", 10, MPFR_RNDN); mpfr_init(r17235); mpfr_init(r17236); mpfr_init(r17237); mpfr_init(r17238); mpfr_init(r17239); mpfr_init_set_str(r17240, "2", 10, MPFR_RNDN); mpfr_init(r17241); mpfr_init(r17242); mpfr_init(r17243); mpfr_init(r17244); mpfr_init(r17245); mpfr_init(r17246); mpfr_init_set_str(r17247, "0.25", 10, MPFR_RNDN); mpfr_init(r17248); mpfr_init(r17249); mpfr_init_set_str(r17250, "0.5", 10, MPFR_RNDN); mpfr_init(r17251); mpfr_init(r17252); mpfr_init(r17253); mpfr_init(r17254); mpfr_init(r17255); mpfr_init(r17256); mpfr_init(r17257); } double f_fm(double alpha, double beta) { mpfr_set_d(r17231, alpha, MPFR_RNDN); ; mpfr_set_si(r17233, mpfr_cmp(r17231, r17232) <= 0, MPFR_RNDN); ; mpfr_add(r17235, r17231, r17234, MPFR_RNDN); mpfr_set_d(r17236, beta, MPFR_RNDN); mpfr_mul(r17237, r17236, r17231, MPFR_RNDN); mpfr_add(r17238, r17236, r17237, MPFR_RNDN); mpfr_add(r17239, r17235, r17238, MPFR_RNDN); ; mpfr_add(r17241, r17240, r17236, MPFR_RNDN); mpfr_add(r17242, r17231, r17241, MPFR_RNDN); mpfr_div(r17243, r17239, r17242, MPFR_RNDN); mpfr_add(r17244, r17235, r17241, MPFR_RNDN); mpfr_div(r17245, r17243, r17244, MPFR_RNDN); mpfr_div(r17246, r17245, r17242, MPFR_RNDN); ; mpfr_add(r17248, r17231, r17236, MPFR_RNDN); mpfr_mul(r17249, r17247, r17248, MPFR_RNDN); ; mpfr_add(r17251, r17249, r17250, MPFR_RNDN); mpfr_add(r17252, r17248, r17240, MPFR_RNDN); mpfr_add(r17253, r17240, r17234, MPFR_RNDN); mpfr_add(r17254, r17248, r17253, MPFR_RNDN); mpfr_mul(r17255, r17252, r17254, MPFR_RNDN); mpfr_div(r17256, r17251, r17255, MPFR_RNDN); if (mpfr_get_si(r17233, MPFR_RNDN)) { mpfr_set(r17257, r17246, MPFR_RNDN); } else { mpfr_set(r17257, r17256, MPFR_RNDN); }; return mpfr_get_d(r17257, MPFR_RNDN); } static mpfr_t r17258, r17259, r17260, r17261, r17262, r17263, r17264, r17265, r17266, r17267, r17268, r17269, r17270, r17271, r17272, r17273, r17274, r17275, r17276, r17277, r17278, r17279, r17280, r17281, r17282, r17283, r17284; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r17258); mpfr_init_set_str(r17259, "1.2543987852557203e+178", 10, MPFR_RNDN); mpfr_init(r17260); mpfr_init_set_str(r17261, "1.0", 10, MPFR_RNDN); mpfr_init(r17262); mpfr_init(r17263); mpfr_init(r17264); mpfr_init(r17265); mpfr_init(r17266); mpfr_init_set_str(r17267, "2", 10, MPFR_RNDN); mpfr_init(r17268); mpfr_init(r17269); mpfr_init(r17270); mpfr_init(r17271); mpfr_init(r17272); mpfr_init(r17273); mpfr_init_set_str(r17274, "0.25", 10, MPFR_RNDN); mpfr_init(r17275); mpfr_init(r17276); mpfr_init_set_str(r17277, "0.5", 10, MPFR_RNDN); mpfr_init(r17278); mpfr_init(r17279); mpfr_init(r17280); mpfr_init(r17281); mpfr_init(r17282); mpfr_init(r17283); mpfr_init(r17284); } double f_dm(double alpha, double beta) { mpfr_set_d(r17258, alpha, MPFR_RNDN); ; mpfr_set_si(r17260, mpfr_cmp(r17258, r17259) <= 0, MPFR_RNDN); ; mpfr_add(r17262, r17258, r17261, MPFR_RNDN); mpfr_set_d(r17263, beta, MPFR_RNDN); mpfr_mul(r17264, r17263, r17258, MPFR_RNDN); mpfr_add(r17265, r17263, r17264, MPFR_RNDN); mpfr_add(r17266, r17262, r17265, MPFR_RNDN); ; mpfr_add(r17268, r17267, r17263, MPFR_RNDN); mpfr_add(r17269, r17258, r17268, MPFR_RNDN); mpfr_div(r17270, r17266, r17269, MPFR_RNDN); mpfr_add(r17271, r17262, r17268, MPFR_RNDN); mpfr_div(r17272, r17270, r17271, MPFR_RNDN); mpfr_div(r17273, r17272, r17269, MPFR_RNDN); ; mpfr_add(r17275, r17258, r17263, MPFR_RNDN); mpfr_mul(r17276, r17274, r17275, MPFR_RNDN); ; mpfr_add(r17278, r17276, r17277, MPFR_RNDN); mpfr_add(r17279, r17275, r17267, MPFR_RNDN); mpfr_add(r17280, r17267, r17261, MPFR_RNDN); mpfr_add(r17281, r17275, r17280, MPFR_RNDN); mpfr_mul(r17282, r17279, r17281, MPFR_RNDN); mpfr_div(r17283, r17278, r17282, MPFR_RNDN); if (mpfr_get_si(r17260, MPFR_RNDN)) { mpfr_set(r17284, r17273, MPFR_RNDN); } else { mpfr_set(r17284, r17283, MPFR_RNDN); }; return mpfr_get_d(r17284, MPFR_RNDN); }