#include #include #include #include #include char *name = "Octave 3.8, jcobi/3"; double f_if(float alpha, float beta) { float r23459 = alpha; float r23460 = beta; float r23461 = r23459 + r23460; float r23462 = r23460 * r23459; float r23463 = r23461 + r23462; float r23464 = 1.0; float r23465 = r23463 + r23464; float r23466 = 2; float r23467 = 1; float r23468 = r23466 * r23467; float r23469 = r23461 + r23468; float r23470 = r23465 / r23469; float r23471 = r23470 / r23469; float r23472 = r23469 + r23464; float r23473 = r23471 / r23472; return r23473; } double f_id(double alpha, double beta) { double r23474 = alpha; double r23475 = beta; double r23476 = r23474 + r23475; double r23477 = r23475 * r23474; double r23478 = r23476 + r23477; double r23479 = 1.0; double r23480 = r23478 + r23479; double r23481 = 2; double r23482 = 1; double r23483 = r23481 * r23482; double r23484 = r23476 + r23483; double r23485 = r23480 / r23484; double r23486 = r23485 / r23484; double r23487 = r23484 + r23479; double r23488 = r23486 / r23487; return r23488; } double f_of(float alpha, float beta) { float r23489 = beta; float r23490 = alpha; float r23491 = r23489 + r23490; float r23492 = 5.479649207112688e+198; bool r23493 = r23491 <= r23492; float r23494 = r23489 * r23490; float r23495 = r23491 + r23494; float r23496 = 1.0; float r23497 = r23495 + r23496; float r23498 = 2; float r23499 = r23491 + r23498; float r23500 = r23497 / r23499; float r23501 = r23500 / r23499; float r23502 = r23499 + r23496; float r23503 = r23501 / r23502; float r23504 = 0.25; float r23505 = r23504 * r23491; float r23506 = 0.5; float r23507 = r23505 + r23506; float r23508 = r23498 + r23496; float r23509 = r23491 + r23508; float r23510 = r23499 * r23509; float r23511 = r23507 / r23510; float r23512 = r23493 ? r23503 : r23511; return r23512; } double f_od(double alpha, double beta) { double r23513 = beta; double r23514 = alpha; double r23515 = r23513 + r23514; double r23516 = 5.479649207112688e+198; bool r23517 = r23515 <= r23516; double r23518 = r23513 * r23514; double r23519 = r23515 + r23518; double r23520 = 1.0; double r23521 = r23519 + r23520; double r23522 = 2; double r23523 = r23515 + r23522; double r23524 = r23521 / r23523; double r23525 = r23524 / r23523; double r23526 = r23523 + r23520; double r23527 = r23525 / r23526; double r23528 = 0.25; double r23529 = r23528 * r23515; double r23530 = 0.5; double r23531 = r23529 + r23530; double r23532 = r23522 + r23520; double r23533 = r23515 + r23532; double r23534 = r23523 * r23533; double r23535 = r23531 / r23534; double r23536 = r23517 ? r23527 : r23535; return r23536; } 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 r23537, r23538, r23539, r23540, r23541, r23542, r23543, r23544, r23545, r23546, r23547, r23548, r23549, r23550, r23551; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r23537); mpfr_init(r23538); mpfr_init(r23539); mpfr_init(r23540); mpfr_init(r23541); mpfr_init_set_str(r23542, "1.0", 10, MPFR_RNDN); mpfr_init(r23543); mpfr_init_set_str(r23544, "2", 10, MPFR_RNDN); mpfr_init_set_str(r23545, "1", 10, MPFR_RNDN); mpfr_init(r23546); mpfr_init(r23547); mpfr_init(r23548); mpfr_init(r23549); mpfr_init(r23550); mpfr_init(r23551); } double f_im(double alpha, double beta) { mpfr_set_d(r23537, alpha, MPFR_RNDN); mpfr_set_d(r23538, beta, MPFR_RNDN); mpfr_add(r23539, r23537, r23538, MPFR_RNDN); mpfr_mul(r23540, r23538, r23537, MPFR_RNDN); mpfr_add(r23541, r23539, r23540, MPFR_RNDN); ; mpfr_add(r23543, r23541, r23542, MPFR_RNDN); ; ; mpfr_mul(r23546, r23544, r23545, MPFR_RNDN); mpfr_add(r23547, r23539, r23546, MPFR_RNDN); mpfr_div(r23548, r23543, r23547, MPFR_RNDN); mpfr_div(r23549, r23548, r23547, MPFR_RNDN); mpfr_add(r23550, r23547, r23542, MPFR_RNDN); mpfr_div(r23551, r23549, r23550, MPFR_RNDN); return mpfr_get_d(r23551, MPFR_RNDN); } static mpfr_t r23552, r23553, r23554, r23555, r23556, r23557, r23558, r23559, r23560, r23561, r23562, r23563, r23564, r23565, r23566, r23567, r23568, r23569, r23570, r23571, r23572, r23573, r23574, r23575; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r23552); mpfr_init(r23553); mpfr_init(r23554); mpfr_init_set_str(r23555, "5.479649207112688e+198", 10, MPFR_RNDN); mpfr_init(r23556); mpfr_init(r23557); mpfr_init(r23558); mpfr_init_set_str(r23559, "1.0", 10, MPFR_RNDN); mpfr_init(r23560); mpfr_init_set_str(r23561, "2", 10, MPFR_RNDN); mpfr_init(r23562); mpfr_init(r23563); mpfr_init(r23564); mpfr_init(r23565); mpfr_init(r23566); mpfr_init_set_str(r23567, "0.25", 10, MPFR_RNDN); mpfr_init(r23568); mpfr_init_set_str(r23569, "0.5", 10, MPFR_RNDN); mpfr_init(r23570); mpfr_init(r23571); mpfr_init(r23572); mpfr_init(r23573); mpfr_init(r23574); mpfr_init(r23575); } double f_fm(double alpha, double beta) { mpfr_set_d(r23552, beta, MPFR_RNDN); mpfr_set_d(r23553, alpha, MPFR_RNDN); mpfr_add(r23554, r23552, r23553, MPFR_RNDN); ; mpfr_set_si(r23556, mpfr_cmp(r23554, r23555) <= 0, MPFR_RNDN); mpfr_mul(r23557, r23552, r23553, MPFR_RNDN); mpfr_add(r23558, r23554, r23557, MPFR_RNDN); ; mpfr_add(r23560, r23558, r23559, MPFR_RNDN); ; mpfr_add(r23562, r23554, r23561, MPFR_RNDN); mpfr_div(r23563, r23560, r23562, MPFR_RNDN); mpfr_div(r23564, r23563, r23562, MPFR_RNDN); mpfr_add(r23565, r23562, r23559, MPFR_RNDN); mpfr_div(r23566, r23564, r23565, MPFR_RNDN); ; mpfr_mul(r23568, r23567, r23554, MPFR_RNDN); ; mpfr_add(r23570, r23568, r23569, MPFR_RNDN); mpfr_add(r23571, r23561, r23559, MPFR_RNDN); mpfr_add(r23572, r23554, r23571, MPFR_RNDN); mpfr_mul(r23573, r23562, r23572, MPFR_RNDN); mpfr_div(r23574, r23570, r23573, MPFR_RNDN); if (mpfr_get_si(r23556, MPFR_RNDN)) { mpfr_set(r23575, r23566, MPFR_RNDN); } else { mpfr_set(r23575, r23574, MPFR_RNDN); }; return mpfr_get_d(r23575, MPFR_RNDN); } static mpfr_t r23576, r23577, r23578, r23579, r23580, r23581, r23582, r23583, r23584, r23585, r23586, r23587, r23588, r23589, r23590, r23591, r23592, r23593, r23594, r23595, r23596, r23597, r23598, r23599; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r23576); mpfr_init(r23577); mpfr_init(r23578); mpfr_init_set_str(r23579, "5.479649207112688e+198", 10, MPFR_RNDN); mpfr_init(r23580); mpfr_init(r23581); mpfr_init(r23582); mpfr_init_set_str(r23583, "1.0", 10, MPFR_RNDN); mpfr_init(r23584); mpfr_init_set_str(r23585, "2", 10, MPFR_RNDN); mpfr_init(r23586); mpfr_init(r23587); mpfr_init(r23588); mpfr_init(r23589); mpfr_init(r23590); mpfr_init_set_str(r23591, "0.25", 10, MPFR_RNDN); mpfr_init(r23592); mpfr_init_set_str(r23593, "0.5", 10, MPFR_RNDN); mpfr_init(r23594); mpfr_init(r23595); mpfr_init(r23596); mpfr_init(r23597); mpfr_init(r23598); mpfr_init(r23599); } double f_dm(double alpha, double beta) { mpfr_set_d(r23576, beta, MPFR_RNDN); mpfr_set_d(r23577, alpha, MPFR_RNDN); mpfr_add(r23578, r23576, r23577, MPFR_RNDN); ; mpfr_set_si(r23580, mpfr_cmp(r23578, r23579) <= 0, MPFR_RNDN); mpfr_mul(r23581, r23576, r23577, MPFR_RNDN); mpfr_add(r23582, r23578, r23581, MPFR_RNDN); ; mpfr_add(r23584, r23582, r23583, MPFR_RNDN); ; mpfr_add(r23586, r23578, r23585, MPFR_RNDN); mpfr_div(r23587, r23584, r23586, MPFR_RNDN); mpfr_div(r23588, r23587, r23586, MPFR_RNDN); mpfr_add(r23589, r23586, r23583, MPFR_RNDN); mpfr_div(r23590, r23588, r23589, MPFR_RNDN); ; mpfr_mul(r23592, r23591, r23578, MPFR_RNDN); ; mpfr_add(r23594, r23592, r23593, MPFR_RNDN); mpfr_add(r23595, r23585, r23583, MPFR_RNDN); mpfr_add(r23596, r23578, r23595, MPFR_RNDN); mpfr_mul(r23597, r23586, r23596, MPFR_RNDN); mpfr_div(r23598, r23594, r23597, MPFR_RNDN); if (mpfr_get_si(r23580, MPFR_RNDN)) { mpfr_set(r23599, r23590, MPFR_RNDN); } else { mpfr_set(r23599, r23598, MPFR_RNDN); }; return mpfr_get_d(r23599, MPFR_RNDN); }