#include #include #include #include #include char *name = "Octave 3.8, jcobi/3"; double f_if(float alpha, float beta) { float r23622 = alpha; float r23623 = beta; float r23624 = r23622 + r23623; float r23625 = r23623 * r23622; float r23626 = r23624 + r23625; float r23627 = 1.0; float r23628 = r23626 + r23627; float r23629 = 2; float r23630 = 1; float r23631 = r23629 * r23630; float r23632 = r23624 + r23631; float r23633 = r23628 / r23632; float r23634 = r23633 / r23632; float r23635 = r23632 + r23627; float r23636 = r23634 / r23635; return r23636; } double f_id(double alpha, double beta) { double r23637 = alpha; double r23638 = beta; double r23639 = r23637 + r23638; double r23640 = r23638 * r23637; double r23641 = r23639 + r23640; double r23642 = 1.0; double r23643 = r23641 + r23642; double r23644 = 2; double r23645 = 1; double r23646 = r23644 * r23645; double r23647 = r23639 + r23646; double r23648 = r23643 / r23647; double r23649 = r23648 / r23647; double r23650 = r23647 + r23642; double r23651 = r23649 / r23650; return r23651; } double f_of(float alpha, float beta) { float r23652 = alpha; float r23653 = 1.7925079542235113e+163; bool r23654 = r23652 <= r23653; float r23655 = beta; float r23656 = r23652 + r23655; float r23657 = r23655 * r23652; float r23658 = r23656 + r23657; float r23659 = 1.0; float r23660 = r23658 + r23659; float r23661 = 2; float r23662 = r23656 + r23661; float r23663 = r23660 / r23662; float r23664 = r23663 / r23662; float r23665 = r23662 + r23659; float r23666 = r23664 / r23665; float r23667 = 1; float r23668 = r23667 / r23652; float r23669 = 2.0; float r23670 = r23669 / r23652; float r23671 = r23670 - r23659; float r23672 = fma(r23668, r23671, r23667); float r23673 = r23661 + r23659; float r23674 = r23656 + r23673; float r23675 = r23661 + r23652; float r23676 = r23655 + r23675; float r23677 = r23674 * r23676; float r23678 = r23672 / r23677; float r23679 = r23654 ? r23666 : r23678; return r23679; } double f_od(double alpha, double beta) { double r23680 = alpha; double r23681 = 1.7925079542235113e+163; bool r23682 = r23680 <= r23681; double r23683 = beta; double r23684 = r23680 + r23683; double r23685 = r23683 * r23680; double r23686 = r23684 + r23685; double r23687 = 1.0; double r23688 = r23686 + r23687; double r23689 = 2; double r23690 = r23684 + r23689; double r23691 = r23688 / r23690; double r23692 = r23691 / r23690; double r23693 = r23690 + r23687; double r23694 = r23692 / r23693; double r23695 = 1; double r23696 = r23695 / r23680; double r23697 = 2.0; double r23698 = r23697 / r23680; double r23699 = r23698 - r23687; double r23700 = fma(r23696, r23699, r23695); double r23701 = r23689 + r23687; double r23702 = r23684 + r23701; double r23703 = r23689 + r23680; double r23704 = r23683 + r23703; double r23705 = r23702 * r23704; double r23706 = r23700 / r23705; double r23707 = r23682 ? r23694 : r23706; return r23707; } 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 r23708, r23709, r23710, r23711, r23712, r23713, r23714, r23715, r23716, r23717, r23718, r23719, r23720, r23721, r23722; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init(r23708); mpfr_init(r23709); mpfr_init(r23710); mpfr_init(r23711); mpfr_init(r23712); mpfr_init_set_str(r23713, "1.0", 10, MPFR_RNDN); mpfr_init(r23714); mpfr_init_set_str(r23715, "2", 10, MPFR_RNDN); mpfr_init_set_str(r23716, "1", 10, MPFR_RNDN); mpfr_init(r23717); mpfr_init(r23718); mpfr_init(r23719); mpfr_init(r23720); mpfr_init(r23721); mpfr_init(r23722); } double f_im(double alpha, double beta) { mpfr_set_d(r23708, alpha, MPFR_RNDN); mpfr_set_d(r23709, beta, MPFR_RNDN); mpfr_add(r23710, r23708, r23709, MPFR_RNDN); mpfr_mul(r23711, r23709, r23708, MPFR_RNDN); mpfr_add(r23712, r23710, r23711, MPFR_RNDN); ; mpfr_add(r23714, r23712, r23713, MPFR_RNDN); ; ; mpfr_mul(r23717, r23715, r23716, MPFR_RNDN); mpfr_add(r23718, r23710, r23717, MPFR_RNDN); mpfr_div(r23719, r23714, r23718, MPFR_RNDN); mpfr_div(r23720, r23719, r23718, MPFR_RNDN); mpfr_add(r23721, r23718, r23713, MPFR_RNDN); mpfr_div(r23722, r23720, r23721, MPFR_RNDN); return mpfr_get_d(r23722, MPFR_RNDN); } static mpfr_t r23723, r23724, r23725, r23726, r23727, r23728, r23729, r23730, r23731, r23732, r23733, r23734, r23735, r23736, r23737, r23738, r23739, r23740, r23741, r23742, r23743, r23744, r23745, r23746, r23747, r23748, r23749, r23750; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init(r23723); mpfr_init_set_str(r23724, "1.7925079542235113e+163", 10, MPFR_RNDN); mpfr_init(r23725); mpfr_init(r23726); mpfr_init(r23727); mpfr_init(r23728); mpfr_init(r23729); mpfr_init_set_str(r23730, "1.0", 10, MPFR_RNDN); mpfr_init(r23731); mpfr_init_set_str(r23732, "2", 10, MPFR_RNDN); mpfr_init(r23733); mpfr_init(r23734); mpfr_init(r23735); mpfr_init(r23736); mpfr_init(r23737); mpfr_init_set_str(r23738, "1", 10, MPFR_RNDN); mpfr_init(r23739); mpfr_init_set_str(r23740, "2.0", 10, MPFR_RNDN); mpfr_init(r23741); mpfr_init(r23742); mpfr_init(r23743); mpfr_init(r23744); mpfr_init(r23745); mpfr_init(r23746); mpfr_init(r23747); mpfr_init(r23748); mpfr_init(r23749); mpfr_init(r23750); } double f_fm(double alpha, double beta) { mpfr_set_d(r23723, alpha, MPFR_RNDN); ; mpfr_set_si(r23725, mpfr_cmp(r23723, r23724) <= 0, MPFR_RNDN); mpfr_set_d(r23726, beta, MPFR_RNDN); mpfr_add(r23727, r23723, r23726, MPFR_RNDN); mpfr_mul(r23728, r23726, r23723, MPFR_RNDN); mpfr_add(r23729, r23727, r23728, MPFR_RNDN); ; mpfr_add(r23731, r23729, r23730, MPFR_RNDN); ; mpfr_add(r23733, r23727, r23732, MPFR_RNDN); mpfr_div(r23734, r23731, r23733, MPFR_RNDN); mpfr_div(r23735, r23734, r23733, MPFR_RNDN); mpfr_add(r23736, r23733, r23730, MPFR_RNDN); mpfr_div(r23737, r23735, r23736, MPFR_RNDN); ; mpfr_div(r23739, r23738, r23723, MPFR_RNDN); ; mpfr_div(r23741, r23740, r23723, MPFR_RNDN); mpfr_sub(r23742, r23741, r23730, MPFR_RNDN); mpfr_fma(r23743, r23739, r23742, r23738, MPFR_RNDN); mpfr_add(r23744, r23732, r23730, MPFR_RNDN); mpfr_add(r23745, r23727, r23744, MPFR_RNDN); mpfr_add(r23746, r23732, r23723, MPFR_RNDN); mpfr_add(r23747, r23726, r23746, MPFR_RNDN); mpfr_mul(r23748, r23745, r23747, MPFR_RNDN); mpfr_div(r23749, r23743, r23748, MPFR_RNDN); if (mpfr_get_si(r23725, MPFR_RNDN)) { mpfr_set(r23750, r23737, MPFR_RNDN); } else { mpfr_set(r23750, r23749, MPFR_RNDN); }; return mpfr_get_d(r23750, MPFR_RNDN); } static mpfr_t r23751, r23752, r23753, r23754, r23755, r23756, r23757, r23758, r23759, r23760, r23761, r23762, r23763, r23764, r23765, r23766, r23767, r23768, r23769, r23770, r23771, r23772, r23773, r23774, r23775, r23776, r23777, r23778; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init(r23751); mpfr_init_set_str(r23752, "1.7925079542235113e+163", 10, MPFR_RNDN); mpfr_init(r23753); mpfr_init(r23754); mpfr_init(r23755); mpfr_init(r23756); mpfr_init(r23757); mpfr_init_set_str(r23758, "1.0", 10, MPFR_RNDN); mpfr_init(r23759); mpfr_init_set_str(r23760, "2", 10, MPFR_RNDN); mpfr_init(r23761); mpfr_init(r23762); mpfr_init(r23763); mpfr_init(r23764); mpfr_init(r23765); mpfr_init_set_str(r23766, "1", 10, MPFR_RNDN); mpfr_init(r23767); mpfr_init_set_str(r23768, "2.0", 10, MPFR_RNDN); mpfr_init(r23769); mpfr_init(r23770); mpfr_init(r23771); mpfr_init(r23772); mpfr_init(r23773); mpfr_init(r23774); mpfr_init(r23775); mpfr_init(r23776); mpfr_init(r23777); mpfr_init(r23778); } double f_dm(double alpha, double beta) { mpfr_set_d(r23751, alpha, MPFR_RNDN); ; mpfr_set_si(r23753, mpfr_cmp(r23751, r23752) <= 0, MPFR_RNDN); mpfr_set_d(r23754, beta, MPFR_RNDN); mpfr_add(r23755, r23751, r23754, MPFR_RNDN); mpfr_mul(r23756, r23754, r23751, MPFR_RNDN); mpfr_add(r23757, r23755, r23756, MPFR_RNDN); ; mpfr_add(r23759, r23757, r23758, MPFR_RNDN); ; mpfr_add(r23761, r23755, r23760, MPFR_RNDN); mpfr_div(r23762, r23759, r23761, MPFR_RNDN); mpfr_div(r23763, r23762, r23761, MPFR_RNDN); mpfr_add(r23764, r23761, r23758, MPFR_RNDN); mpfr_div(r23765, r23763, r23764, MPFR_RNDN); ; mpfr_div(r23767, r23766, r23751, MPFR_RNDN); ; mpfr_div(r23769, r23768, r23751, MPFR_RNDN); mpfr_sub(r23770, r23769, r23758, MPFR_RNDN); mpfr_fma(r23771, r23767, r23770, r23766, MPFR_RNDN); mpfr_add(r23772, r23760, r23758, MPFR_RNDN); mpfr_add(r23773, r23755, r23772, MPFR_RNDN); mpfr_add(r23774, r23760, r23751, MPFR_RNDN); mpfr_add(r23775, r23754, r23774, MPFR_RNDN); mpfr_mul(r23776, r23773, r23775, MPFR_RNDN); mpfr_div(r23777, r23771, r23776, MPFR_RNDN); if (mpfr_get_si(r23753, MPFR_RNDN)) { mpfr_set(r23778, r23765, MPFR_RNDN); } else { mpfr_set(r23778, r23777, MPFR_RNDN); }; return mpfr_get_d(r23778, MPFR_RNDN); }