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