#include #include #include #include #include char *name = "Octave 3.8, jcobi/1"; double f_if(float alpha, float beta) { float r17241 = beta; float r17242 = alpha; float r17243 = r17241 - r17242; float r17244 = r17242 + r17241; float r17245 = 2.0f; float r17246 = r17244 + r17245; float r17247 = r17243 / r17246; float r17248 = 1.0f; float r17249 = r17247 + r17248; float r17250 = r17249 / r17245; return r17250; } double f_id(double alpha, double beta) { double r17251 = beta; double r17252 = alpha; double r17253 = r17251 - r17252; double r17254 = r17252 + r17251; double r17255 = 2.0; double r17256 = r17254 + r17255; double r17257 = r17253 / r17256; double r17258 = 1.0; double r17259 = r17257 + r17258; double r17260 = r17259 / r17255; return r17260; } double f_of(float alpha, float beta) { float r17261 = beta; float r17262 = 2.0f; float r17263 = r17261 / r17262; float r17264 = cbrt(r17263); float r17265 = alpha; float r17266 = r17265 + r17262; float r17267 = r17266 + r17261; float r17268 = cbrt(r17267); float r17269 = r17264 / r17268; float r17270 = r17269 * (r17269 * r17269); float r17271 = r17265 + r17261; float r17272 = r17271 + r17262; float r17273 = 1.0f; float r17274 = pow(r17272, r17273); float r17275 = r17265 / r17274; float r17276 = 1.0f; float r17277 = r17275 - r17276; float r17278 = r17277 / r17262; float r17279 = r17270 - r17278; return r17279; } double f_od(double alpha, double beta) { double r17280 = beta; double r17281 = 2.0; double r17282 = r17280 / r17281; double r17283 = cbrt(r17282); double r17284 = alpha; double r17285 = r17284 + r17281; double r17286 = r17285 + r17280; double r17287 = cbrt(r17286); double r17288 = r17283 / r17287; double r17289 = r17288 * (r17288 * r17288); double r17290 = r17284 + r17280; double r17291 = r17290 + r17281; double r17292 = 1.0; double r17293 = pow(r17291, r17292); double r17294 = r17284 / r17293; double r17295 = 1.0; double r17296 = r17294 - r17295; double r17297 = r17296 / r17281; double r17298 = r17289 - r17297; return r17298; } 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 r17299, r17300, r17301, r17302, r17303, r17304, r17305, r17306, r17307, r17308; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r17299); mpfr_init(r17300); mpfr_init(r17301); mpfr_init(r17302); mpfr_init_set_str(r17303, "2.0", 10, MPFR_RNDN); mpfr_init(r17304); mpfr_init(r17305); mpfr_init_set_str(r17306, "1.0", 10, MPFR_RNDN); mpfr_init(r17307); mpfr_init(r17308); } double f_im(double alpha, double beta) { mpfr_set_d(r17299, beta, MPFR_RNDN); mpfr_set_d(r17300, alpha, MPFR_RNDN); mpfr_sub(r17301, r17299, r17300, MPFR_RNDN); mpfr_add(r17302, r17300, r17299, MPFR_RNDN); ; mpfr_add(r17304, r17302, r17303, MPFR_RNDN); mpfr_div(r17305, r17301, r17304, MPFR_RNDN); ; mpfr_add(r17307, r17305, r17306, MPFR_RNDN); mpfr_div(r17308, r17307, r17303, MPFR_RNDN); return mpfr_get_d(r17308, MPFR_RNDN); } static mpfr_t r17309, r17310, r17311, r17312, r17313, r17314, r17315, r17316, r17317, r17318, r17319, r17320, r17321, r17322, r17323, r17324, r17325, r17326, r17327; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r17309); mpfr_init_set_str(r17310, "2.0", 10, MPFR_RNDN); mpfr_init(r17311); mpfr_init(r17312); mpfr_init(r17313); mpfr_init(r17314); mpfr_init(r17315); mpfr_init(r17316); mpfr_init(r17317); mpfr_init(r17318); mpfr_init(r17319); mpfr_init(r17320); mpfr_init_set_str(r17321, "1", 10, MPFR_RNDN); mpfr_init(r17322); mpfr_init(r17323); mpfr_init_set_str(r17324, "1.0", 10, MPFR_RNDN); mpfr_init(r17325); mpfr_init(r17326); mpfr_init(r17327); } double f_fm(double alpha, double beta) { mpfr_set_d(r17309, beta, MPFR_RNDN); ; mpfr_div(r17311, r17309, r17310, MPFR_RNDN); mpfr_cbrt(r17312, r17311, MPFR_RNDN); mpfr_set_d(r17313, alpha, MPFR_RNDN); mpfr_add(r17314, r17313, r17310, MPFR_RNDN); mpfr_add(r17315, r17314, r17309, MPFR_RNDN); mpfr_cbrt(r17316, r17315, MPFR_RNDN); mpfr_div(r17317, r17312, r17316, MPFR_RNDN); mpfr_mul(r17318, r17317, r17317, MPFR_RNDN); mpfr_mul(r17318, r17318, r17317, MPFR_RNDN); mpfr_add(r17319, r17313, r17309, MPFR_RNDN); mpfr_add(r17320, r17319, r17310, MPFR_RNDN); ; mpfr_pow(r17322, r17320, r17321, MPFR_RNDN); mpfr_div(r17323, r17313, r17322, MPFR_RNDN); ; mpfr_sub(r17325, r17323, r17324, MPFR_RNDN); mpfr_div(r17326, r17325, r17310, MPFR_RNDN); mpfr_sub(r17327, r17318, r17326, MPFR_RNDN); return mpfr_get_d(r17327, MPFR_RNDN); } static mpfr_t r17328, r17329, r17330, r17331, r17332, r17333, r17334, r17335, r17336, r17337, r17338, r17339, r17340, r17341, r17342, r17343, r17344, r17345, r17346; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r17328); mpfr_init_set_str(r17329, "2.0", 10, MPFR_RNDN); mpfr_init(r17330); mpfr_init(r17331); mpfr_init(r17332); mpfr_init(r17333); mpfr_init(r17334); mpfr_init(r17335); mpfr_init(r17336); mpfr_init(r17337); mpfr_init(r17338); mpfr_init(r17339); mpfr_init_set_str(r17340, "1", 10, MPFR_RNDN); mpfr_init(r17341); mpfr_init(r17342); mpfr_init_set_str(r17343, "1.0", 10, MPFR_RNDN); mpfr_init(r17344); mpfr_init(r17345); mpfr_init(r17346); } double f_dm(double alpha, double beta) { mpfr_set_d(r17328, beta, MPFR_RNDN); ; mpfr_div(r17330, r17328, r17329, MPFR_RNDN); mpfr_cbrt(r17331, r17330, MPFR_RNDN); mpfr_set_d(r17332, alpha, MPFR_RNDN); mpfr_add(r17333, r17332, r17329, MPFR_RNDN); mpfr_add(r17334, r17333, r17328, MPFR_RNDN); mpfr_cbrt(r17335, r17334, MPFR_RNDN); mpfr_div(r17336, r17331, r17335, MPFR_RNDN); mpfr_mul(r17337, r17336, r17336, MPFR_RNDN); mpfr_mul(r17337, r17337, r17336, MPFR_RNDN); mpfr_add(r17338, r17332, r17328, MPFR_RNDN); mpfr_add(r17339, r17338, r17329, MPFR_RNDN); ; mpfr_pow(r17341, r17339, r17340, MPFR_RNDN); mpfr_div(r17342, r17332, r17341, MPFR_RNDN); ; mpfr_sub(r17344, r17342, r17343, MPFR_RNDN); mpfr_div(r17345, r17344, r17329, MPFR_RNDN); mpfr_sub(r17346, r17337, r17345, MPFR_RNDN); return mpfr_get_d(r17346, MPFR_RNDN); }