#include #include #include #include #include char *name = "Octave 3.8, jcobi/4"; double f_if(float alpha, float beta, float i) { float r17109 = i; float r17110 = alpha; float r17111 = beta; float r17112 = r17110 + r17111; float r17113 = r17112 + r17109; float r17114 = r17109 * r17113; float r17115 = r17111 * r17110; float r17116 = r17115 + r17114; float r17117 = r17114 * r17116; float r17118 = 2.0f; float r17119 = r17118 * r17109; float r17120 = r17112 + r17119; float r17121 = r17120 * r17120; float r17122 = r17117 / r17121; float r17123 = 1.0f; float r17124 = r17121 - r17123; float r17125 = r17122 / r17124; return r17125; } double f_id(double alpha, double beta, double i) { double r17126 = i; double r17127 = alpha; double r17128 = beta; double r17129 = r17127 + r17128; double r17130 = r17129 + r17126; double r17131 = r17126 * r17130; double r17132 = r17128 * r17127; double r17133 = r17132 + r17131; double r17134 = r17131 * r17133; double r17135 = 2.0; double r17136 = r17135 * r17126; double r17137 = r17129 + r17136; double r17138 = r17137 * r17137; double r17139 = r17134 / r17138; double r17140 = 1.0; double r17141 = r17138 - r17140; double r17142 = r17139 / r17141; return r17142; } double f_of(float alpha, float beta, float i) { float r17143 = beta; float r17144 = 1.8414905412296409e+137f; bool r17145 = r17143 <= r17144; float r17146 = i; float r17147 = alpha; float r17148 = r17146 + r17147; float r17149 = r17143 + r17148; float r17150 = r17146 * r17149; float r17151 = sqrt(r17150); float r17152 = r17143 + r17147; float r17153 = 2.0f; float r17154 = r17153 * r17146; float r17155 = r17152 + r17154; float r17156 = r17155 * r17155; float r17157 = r17147 * r17143; float r17158 = r17157 + r17150; float r17159 = r17156 / r17158; float r17160 = sqrt(r17159); float r17161 = r17151 / r17160; float r17162 = 1.0f; float r17163 = r17156 - r17162; float r17164 = sqrt(r17163); float r17165 = r17161 / r17164; float r17166 = r17165 * r17165; float r17167 = 0.0f; float r17168 = r17145 ? r17166 : r17167; return r17168; } double f_od(double alpha, double beta, double i) { double r17169 = beta; double r17170 = 1.8414905412296409e+137; bool r17171 = r17169 <= r17170; double r17172 = i; double r17173 = alpha; double r17174 = r17172 + r17173; double r17175 = r17169 + r17174; double r17176 = r17172 * r17175; double r17177 = sqrt(r17176); double r17178 = r17169 + r17173; double r17179 = 2.0; double r17180 = r17179 * r17172; double r17181 = r17178 + r17180; double r17182 = r17181 * r17181; double r17183 = r17173 * r17169; double r17184 = r17183 + r17176; double r17185 = r17182 / r17184; double r17186 = sqrt(r17185); double r17187 = r17177 / r17186; double r17188 = 1.0; double r17189 = r17182 - r17188; double r17190 = sqrt(r17189); double r17191 = r17187 / r17190; double r17192 = r17191 * r17191; double r17193 = 0.0; double r17194 = r17171 ? r17192 : r17193; return r17194; } 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 r17195, r17196, r17197, r17198, r17199, r17200, r17201, r17202, r17203, r17204, r17205, r17206, r17207, r17208, r17209, r17210, r17211; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r17195); mpfr_init(r17196); mpfr_init(r17197); mpfr_init(r17198); mpfr_init(r17199); mpfr_init(r17200); mpfr_init(r17201); mpfr_init(r17202); mpfr_init(r17203); mpfr_init_set_str(r17204, "2", 10, MPFR_RNDN); mpfr_init(r17205); mpfr_init(r17206); mpfr_init(r17207); mpfr_init(r17208); mpfr_init_set_str(r17209, "1.0", 10, MPFR_RNDN); mpfr_init(r17210); mpfr_init(r17211); } double f_im(double alpha, double beta, double i) { mpfr_set_d(r17195, i, MPFR_RNDN); mpfr_set_d(r17196, alpha, MPFR_RNDN); mpfr_set_d(r17197, beta, MPFR_RNDN); mpfr_add(r17198, r17196, r17197, MPFR_RNDN); mpfr_add(r17199, r17198, r17195, MPFR_RNDN); mpfr_mul(r17200, r17195, r17199, MPFR_RNDN); mpfr_mul(r17201, r17197, r17196, MPFR_RNDN); mpfr_add(r17202, r17201, r17200, MPFR_RNDN); mpfr_mul(r17203, r17200, r17202, MPFR_RNDN); ; mpfr_mul(r17205, r17204, r17195, MPFR_RNDN); mpfr_add(r17206, r17198, r17205, MPFR_RNDN); mpfr_mul(r17207, r17206, r17206, MPFR_RNDN); mpfr_div(r17208, r17203, r17207, MPFR_RNDN); ; mpfr_sub(r17210, r17207, r17209, MPFR_RNDN); mpfr_div(r17211, r17208, r17210, MPFR_RNDN); return mpfr_get_d(r17211, MPFR_RNDN); } static mpfr_t r17212, r17213, r17214, r17215, r17216, r17217, r17218, r17219, r17220, r17221, r17222, r17223, r17224, r17225, r17226, r17227, r17228, r17229, r17230, r17231, r17232, r17233, r17234, r17235, r17236, r17237; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r17212); mpfr_init_set_str(r17213, "1.8414905412296409e+137", 10, MPFR_RNDN); mpfr_init(r17214); mpfr_init(r17215); mpfr_init(r17216); mpfr_init(r17217); mpfr_init(r17218); mpfr_init(r17219); mpfr_init(r17220); mpfr_init(r17221); mpfr_init_set_str(r17222, "2", 10, MPFR_RNDN); mpfr_init(r17223); mpfr_init(r17224); mpfr_init(r17225); mpfr_init(r17226); mpfr_init(r17227); mpfr_init(r17228); mpfr_init(r17229); mpfr_init(r17230); mpfr_init_set_str(r17231, "1.0", 10, MPFR_RNDN); mpfr_init(r17232); mpfr_init(r17233); mpfr_init(r17234); mpfr_init(r17235); mpfr_init_set_str(r17236, "0", 10, MPFR_RNDN); mpfr_init(r17237); } double f_fm(double alpha, double beta, double i) { mpfr_set_d(r17212, beta, MPFR_RNDN); ; mpfr_set_si(r17214, mpfr_cmp(r17212, r17213) <= 0, MPFR_RNDN); mpfr_set_d(r17215, i, MPFR_RNDN); mpfr_set_d(r17216, alpha, MPFR_RNDN); mpfr_add(r17217, r17215, r17216, MPFR_RNDN); mpfr_add(r17218, r17212, r17217, MPFR_RNDN); mpfr_mul(r17219, r17215, r17218, MPFR_RNDN); mpfr_sqrt(r17220, r17219, MPFR_RNDN); mpfr_add(r17221, r17212, r17216, MPFR_RNDN); ; mpfr_mul(r17223, r17222, r17215, MPFR_RNDN); mpfr_add(r17224, r17221, r17223, MPFR_RNDN); mpfr_sqr(r17225, r17224, MPFR_RNDN); mpfr_mul(r17226, r17216, r17212, MPFR_RNDN); mpfr_add(r17227, r17226, r17219, MPFR_RNDN); mpfr_div(r17228, r17225, r17227, MPFR_RNDN); mpfr_sqrt(r17229, r17228, MPFR_RNDN); mpfr_div(r17230, r17220, r17229, MPFR_RNDN); ; mpfr_sub(r17232, r17225, r17231, MPFR_RNDN); mpfr_sqrt(r17233, r17232, MPFR_RNDN); mpfr_div(r17234, r17230, r17233, MPFR_RNDN); mpfr_sqr(r17235, r17234, MPFR_RNDN); ; if (mpfr_get_si(r17214, MPFR_RNDN)) { mpfr_set(r17237, r17235, MPFR_RNDN); } else { mpfr_set(r17237, r17236, MPFR_RNDN); }; return mpfr_get_d(r17237, MPFR_RNDN); } static mpfr_t r17238, r17239, r17240, r17241, r17242, r17243, r17244, r17245, r17246, r17247, r17248, r17249, r17250, r17251, r17252, r17253, r17254, r17255, r17256, r17257, r17258, r17259, r17260, r17261, r17262, r17263; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r17238); mpfr_init_set_str(r17239, "1.8414905412296409e+137", 10, MPFR_RNDN); mpfr_init(r17240); mpfr_init(r17241); mpfr_init(r17242); mpfr_init(r17243); mpfr_init(r17244); mpfr_init(r17245); mpfr_init(r17246); mpfr_init(r17247); mpfr_init_set_str(r17248, "2", 10, MPFR_RNDN); mpfr_init(r17249); mpfr_init(r17250); mpfr_init(r17251); mpfr_init(r17252); mpfr_init(r17253); mpfr_init(r17254); mpfr_init(r17255); mpfr_init(r17256); mpfr_init_set_str(r17257, "1.0", 10, MPFR_RNDN); mpfr_init(r17258); mpfr_init(r17259); mpfr_init(r17260); mpfr_init(r17261); mpfr_init_set_str(r17262, "0", 10, MPFR_RNDN); mpfr_init(r17263); } double f_dm(double alpha, double beta, double i) { mpfr_set_d(r17238, beta, MPFR_RNDN); ; mpfr_set_si(r17240, mpfr_cmp(r17238, r17239) <= 0, MPFR_RNDN); mpfr_set_d(r17241, i, MPFR_RNDN); mpfr_set_d(r17242, alpha, MPFR_RNDN); mpfr_add(r17243, r17241, r17242, MPFR_RNDN); mpfr_add(r17244, r17238, r17243, MPFR_RNDN); mpfr_mul(r17245, r17241, r17244, MPFR_RNDN); mpfr_sqrt(r17246, r17245, MPFR_RNDN); mpfr_add(r17247, r17238, r17242, MPFR_RNDN); ; mpfr_mul(r17249, r17248, r17241, MPFR_RNDN); mpfr_add(r17250, r17247, r17249, MPFR_RNDN); mpfr_sqr(r17251, r17250, MPFR_RNDN); mpfr_mul(r17252, r17242, r17238, MPFR_RNDN); mpfr_add(r17253, r17252, r17245, MPFR_RNDN); mpfr_div(r17254, r17251, r17253, MPFR_RNDN); mpfr_sqrt(r17255, r17254, MPFR_RNDN); mpfr_div(r17256, r17246, r17255, MPFR_RNDN); ; mpfr_sub(r17258, r17251, r17257, MPFR_RNDN); mpfr_sqrt(r17259, r17258, MPFR_RNDN); mpfr_div(r17260, r17256, r17259, MPFR_RNDN); mpfr_sqr(r17261, r17260, MPFR_RNDN); ; if (mpfr_get_si(r17240, MPFR_RNDN)) { mpfr_set(r17263, r17261, MPFR_RNDN); } else { mpfr_set(r17263, r17262, MPFR_RNDN); }; return mpfr_get_d(r17263, MPFR_RNDN); }