#include #include #include #include #include char *name = "Octave 3.8, jcobi/2"; double f_if(float alpha, float beta, float i) { float r22462 = alpha; float r22463 = beta; float r22464 = r22462 + r22463; float r22465 = r22463 - r22462; float r22466 = r22464 * r22465; float r22467 = 2; float r22468 = i; float r22469 = r22467 * r22468; float r22470 = r22464 + r22469; float r22471 = r22466 / r22470; float r22472 = 2.0; float r22473 = r22470 + r22472; float r22474 = r22471 / r22473; float r22475 = 1.0; float r22476 = r22474 + r22475; float r22477 = r22476 / r22472; return r22477; } double f_id(double alpha, double beta, double i) { double r22478 = alpha; double r22479 = beta; double r22480 = r22478 + r22479; double r22481 = r22479 - r22478; double r22482 = r22480 * r22481; double r22483 = 2; double r22484 = i; double r22485 = r22483 * r22484; double r22486 = r22480 + r22485; double r22487 = r22482 / r22486; double r22488 = 2.0; double r22489 = r22486 + r22488; double r22490 = r22487 / r22489; double r22491 = 1.0; double r22492 = r22490 + r22491; double r22493 = r22492 / r22488; return r22493; } double f_of(float alpha, float beta, float i) { float r22494 = alpha; float r22495 = beta; float r22496 = r22494 + r22495; float r22497 = 2; float r22498 = i; float r22499 = 2.0; float r22500 = fma(r22497, r22498, r22499); float r22501 = r22496 + r22500; float r22502 = r22496 / r22501; float r22503 = r22495 - r22494; float r22504 = fma(r22498, r22497, r22494); float r22505 = r22504 + r22495; float r22506 = cbrt(r22505); float r22507 = r22503 / r22506; float r22508 = r22506 * r22506; float r22509 = r22507 / r22508; float r22510 = 1.0; float r22511 = fma(r22502, r22509, r22510); float r22512 = 3; float r22513 = pow(r22511, r22512); float r22514 = cbrt(r22513); float r22515 = r22514 / r22499; return r22515; } double f_od(double alpha, double beta, double i) { double r22516 = alpha; double r22517 = beta; double r22518 = r22516 + r22517; double r22519 = 2; double r22520 = i; double r22521 = 2.0; double r22522 = fma(r22519, r22520, r22521); double r22523 = r22518 + r22522; double r22524 = r22518 / r22523; double r22525 = r22517 - r22516; double r22526 = fma(r22520, r22519, r22516); double r22527 = r22526 + r22517; double r22528 = cbrt(r22527); double r22529 = r22525 / r22528; double r22530 = r22528 * r22528; double r22531 = r22529 / r22530; double r22532 = 1.0; double r22533 = fma(r22524, r22531, r22532); double r22534 = 3; double r22535 = pow(r22533, r22534); double r22536 = cbrt(r22535); double r22537 = r22536 / r22521; return r22537; } 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 r22538, r22539, r22540, r22541, r22542, r22543, r22544, r22545, r22546, r22547, r22548, r22549, r22550, r22551, r22552, r22553; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init(r22538); mpfr_init(r22539); mpfr_init(r22540); mpfr_init(r22541); mpfr_init(r22542); mpfr_init_set_str(r22543, "2", 10, MPFR_RNDN); mpfr_init(r22544); mpfr_init(r22545); mpfr_init(r22546); mpfr_init(r22547); mpfr_init_set_str(r22548, "2.0", 10, MPFR_RNDN); mpfr_init(r22549); mpfr_init(r22550); mpfr_init_set_str(r22551, "1.0", 10, MPFR_RNDN); mpfr_init(r22552); mpfr_init(r22553); } double f_im(double alpha, double beta, double i) { mpfr_set_d(r22538, alpha, MPFR_RNDN); mpfr_set_d(r22539, beta, MPFR_RNDN); mpfr_add(r22540, r22538, r22539, MPFR_RNDN); mpfr_sub(r22541, r22539, r22538, MPFR_RNDN); mpfr_mul(r22542, r22540, r22541, MPFR_RNDN); ; mpfr_set_d(r22544, i, MPFR_RNDN); mpfr_mul(r22545, r22543, r22544, MPFR_RNDN); mpfr_add(r22546, r22540, r22545, MPFR_RNDN); mpfr_div(r22547, r22542, r22546, MPFR_RNDN); ; mpfr_add(r22549, r22546, r22548, MPFR_RNDN); mpfr_div(r22550, r22547, r22549, MPFR_RNDN); ; mpfr_add(r22552, r22550, r22551, MPFR_RNDN); mpfr_div(r22553, r22552, r22548, MPFR_RNDN); return mpfr_get_d(r22553, MPFR_RNDN); } static mpfr_t r22554, r22555, r22556, r22557, r22558, r22559, r22560, r22561, r22562, r22563, r22564, r22565, r22566, r22567, r22568, r22569, r22570, r22571, r22572, r22573, r22574, r22575; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r22554); mpfr_init(r22555); mpfr_init(r22556); mpfr_init_set_str(r22557, "2", 10, MPFR_RNDN); mpfr_init(r22558); mpfr_init_set_str(r22559, "2.0", 10, MPFR_RNDN); mpfr_init(r22560); mpfr_init(r22561); mpfr_init(r22562); mpfr_init(r22563); mpfr_init(r22564); mpfr_init(r22565); mpfr_init(r22566); mpfr_init(r22567); mpfr_init(r22568); mpfr_init(r22569); mpfr_init_set_str(r22570, "1.0", 10, MPFR_RNDN); mpfr_init(r22571); mpfr_init_set_str(r22572, "3", 10, MPFR_RNDN); mpfr_init(r22573); mpfr_init(r22574); mpfr_init(r22575); } double f_fm(double alpha, double beta, double i) { mpfr_set_d(r22554, alpha, MPFR_RNDN); mpfr_set_d(r22555, beta, MPFR_RNDN); mpfr_add(r22556, r22554, r22555, MPFR_RNDN); ; mpfr_set_d(r22558, i, MPFR_RNDN); ; mpfr_fma(r22560, r22557, r22558, r22559, MPFR_RNDN); mpfr_add(r22561, r22556, r22560, MPFR_RNDN); mpfr_div(r22562, r22556, r22561, MPFR_RNDN); mpfr_sub(r22563, r22555, r22554, MPFR_RNDN); mpfr_fma(r22564, r22558, r22557, r22554, MPFR_RNDN); mpfr_add(r22565, r22564, r22555, MPFR_RNDN); mpfr_cbrt(r22566, r22565, MPFR_RNDN); mpfr_div(r22567, r22563, r22566, MPFR_RNDN); mpfr_mul(r22568, r22566, r22566, MPFR_RNDN); mpfr_div(r22569, r22567, r22568, MPFR_RNDN); ; mpfr_fma(r22571, r22562, r22569, r22570, MPFR_RNDN); ; mpfr_pow(r22573, r22571, r22572, MPFR_RNDN); mpfr_cbrt(r22574, r22573, MPFR_RNDN); mpfr_div(r22575, r22574, r22559, MPFR_RNDN); return mpfr_get_d(r22575, MPFR_RNDN); } static mpfr_t r22576, r22577, r22578, r22579, r22580, r22581, r22582, r22583, r22584, r22585, r22586, r22587, r22588, r22589, r22590, r22591, r22592, r22593, r22594, r22595, r22596, r22597; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r22576); mpfr_init(r22577); mpfr_init(r22578); mpfr_init_set_str(r22579, "2", 10, MPFR_RNDN); mpfr_init(r22580); mpfr_init_set_str(r22581, "2.0", 10, MPFR_RNDN); mpfr_init(r22582); mpfr_init(r22583); mpfr_init(r22584); mpfr_init(r22585); mpfr_init(r22586); mpfr_init(r22587); mpfr_init(r22588); mpfr_init(r22589); mpfr_init(r22590); mpfr_init(r22591); mpfr_init_set_str(r22592, "1.0", 10, MPFR_RNDN); mpfr_init(r22593); mpfr_init_set_str(r22594, "3", 10, MPFR_RNDN); mpfr_init(r22595); mpfr_init(r22596); mpfr_init(r22597); } double f_dm(double alpha, double beta, double i) { mpfr_set_d(r22576, alpha, MPFR_RNDN); mpfr_set_d(r22577, beta, MPFR_RNDN); mpfr_add(r22578, r22576, r22577, MPFR_RNDN); ; mpfr_set_d(r22580, i, MPFR_RNDN); ; mpfr_fma(r22582, r22579, r22580, r22581, MPFR_RNDN); mpfr_add(r22583, r22578, r22582, MPFR_RNDN); mpfr_div(r22584, r22578, r22583, MPFR_RNDN); mpfr_sub(r22585, r22577, r22576, MPFR_RNDN); mpfr_fma(r22586, r22580, r22579, r22576, MPFR_RNDN); mpfr_add(r22587, r22586, r22577, MPFR_RNDN); mpfr_cbrt(r22588, r22587, MPFR_RNDN); mpfr_div(r22589, r22585, r22588, MPFR_RNDN); mpfr_mul(r22590, r22588, r22588, MPFR_RNDN); mpfr_div(r22591, r22589, r22590, MPFR_RNDN); ; mpfr_fma(r22593, r22584, r22591, r22592, MPFR_RNDN); ; mpfr_pow(r22595, r22593, r22594, MPFR_RNDN); mpfr_cbrt(r22596, r22595, MPFR_RNDN); mpfr_div(r22597, r22596, r22581, MPFR_RNDN); return mpfr_get_d(r22597, MPFR_RNDN); }