#include #include #include #include #include char *name = "Octave 3.8, jcobi/4"; double f_if(float alpha, float beta, float i) { float r23648 = i; float r23649 = alpha; float r23650 = beta; float r23651 = r23649 + r23650; float r23652 = r23651 + r23648; float r23653 = r23648 * r23652; float r23654 = r23650 * r23649; float r23655 = r23654 + r23653; float r23656 = r23653 * r23655; float r23657 = 2; float r23658 = r23657 * r23648; float r23659 = r23651 + r23658; float r23660 = r23659 * r23659; float r23661 = r23656 / r23660; float r23662 = 1.0; float r23663 = r23660 - r23662; float r23664 = r23661 / r23663; return r23664; } double f_id(double alpha, double beta, double i) { double r23665 = i; double r23666 = alpha; double r23667 = beta; double r23668 = r23666 + r23667; double r23669 = r23668 + r23665; double r23670 = r23665 * r23669; double r23671 = r23667 * r23666; double r23672 = r23671 + r23670; double r23673 = r23670 * r23672; double r23674 = 2; double r23675 = r23674 * r23665; double r23676 = r23668 + r23675; double r23677 = r23676 * r23676; double r23678 = r23673 / r23677; double r23679 = 1.0; double r23680 = r23677 - r23679; double r23681 = r23678 / r23680; return r23681; } double f_of(float alpha, float beta, float i) { float r23682 = alpha; float r23683 = 1.3406027096911442e+154; bool r23684 = r23682 <= r23683; float r23685 = i; float r23686 = beta; float r23687 = r23682 + r23686; float r23688 = 2; float r23689 = r23685 * r23688; float r23690 = r23687 + r23689; float r23691 = r23685 / r23690; float r23692 = r23685 + r23686; float r23693 = r23682 + r23692; float r23694 = r23691 * r23693; float r23695 = r23686 * r23682; float r23696 = r23687 + r23685; float r23697 = r23685 * r23696; float r23698 = r23695 + r23697; float r23699 = r23688 * r23685; float r23700 = r23687 + r23699; float r23701 = r23698 / r23700; float r23702 = r23700 * r23700; float r23703 = 1.0; float r23704 = r23702 - r23703; float r23705 = r23701 / r23704; float r23706 = r23694 * r23705; float r23707 = 0; float r23708 = r23684 ? r23706 : r23707; return r23708; } double f_od(double alpha, double beta, double i) { double r23709 = alpha; double r23710 = 1.3406027096911442e+154; bool r23711 = r23709 <= r23710; double r23712 = i; double r23713 = beta; double r23714 = r23709 + r23713; double r23715 = 2; double r23716 = r23712 * r23715; double r23717 = r23714 + r23716; double r23718 = r23712 / r23717; double r23719 = r23712 + r23713; double r23720 = r23709 + r23719; double r23721 = r23718 * r23720; double r23722 = r23713 * r23709; double r23723 = r23714 + r23712; double r23724 = r23712 * r23723; double r23725 = r23722 + r23724; double r23726 = r23715 * r23712; double r23727 = r23714 + r23726; double r23728 = r23725 / r23727; double r23729 = r23727 * r23727; double r23730 = 1.0; double r23731 = r23729 - r23730; double r23732 = r23728 / r23731; double r23733 = r23721 * r23732; double r23734 = 0; double r23735 = r23711 ? r23733 : r23734; return r23735; } 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 r23736, r23737, r23738, r23739, r23740, r23741, r23742, r23743, r23744, r23745, r23746, r23747, r23748, r23749, r23750, r23751, r23752; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r23736); mpfr_init(r23737); mpfr_init(r23738); mpfr_init(r23739); mpfr_init(r23740); mpfr_init(r23741); mpfr_init(r23742); mpfr_init(r23743); mpfr_init(r23744); mpfr_init_set_str(r23745, "2", 10, MPFR_RNDN); mpfr_init(r23746); mpfr_init(r23747); mpfr_init(r23748); mpfr_init(r23749); mpfr_init_set_str(r23750, "1.0", 10, MPFR_RNDN); mpfr_init(r23751); mpfr_init(r23752); } double f_im(double alpha, double beta, double i) { mpfr_set_d(r23736, i, MPFR_RNDN); mpfr_set_d(r23737, alpha, MPFR_RNDN); mpfr_set_d(r23738, beta, MPFR_RNDN); mpfr_add(r23739, r23737, r23738, MPFR_RNDN); mpfr_add(r23740, r23739, r23736, MPFR_RNDN); mpfr_mul(r23741, r23736, r23740, MPFR_RNDN); mpfr_mul(r23742, r23738, r23737, MPFR_RNDN); mpfr_add(r23743, r23742, r23741, MPFR_RNDN); mpfr_mul(r23744, r23741, r23743, MPFR_RNDN); ; mpfr_mul(r23746, r23745, r23736, MPFR_RNDN); mpfr_add(r23747, r23739, r23746, MPFR_RNDN); mpfr_mul(r23748, r23747, r23747, MPFR_RNDN); mpfr_div(r23749, r23744, r23748, MPFR_RNDN); ; mpfr_sub(r23751, r23748, r23750, MPFR_RNDN); mpfr_div(r23752, r23749, r23751, MPFR_RNDN); return mpfr_get_d(r23752, MPFR_RNDN); } static mpfr_t 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, r23779; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r23753); mpfr_init_set_str(r23754, "1.3406027096911442e+154", 10, MPFR_RNDN); mpfr_init(r23755); mpfr_init(r23756); mpfr_init(r23757); mpfr_init(r23758); mpfr_init_set_str(r23759, "2", 10, MPFR_RNDN); mpfr_init(r23760); mpfr_init(r23761); mpfr_init(r23762); mpfr_init(r23763); mpfr_init(r23764); mpfr_init(r23765); mpfr_init(r23766); mpfr_init(r23767); mpfr_init(r23768); mpfr_init(r23769); mpfr_init(r23770); mpfr_init(r23771); mpfr_init(r23772); mpfr_init(r23773); mpfr_init_set_str(r23774, "1.0", 10, MPFR_RNDN); mpfr_init(r23775); mpfr_init(r23776); mpfr_init(r23777); mpfr_init_set_str(r23778, "0", 10, MPFR_RNDN); mpfr_init(r23779); } double f_fm(double alpha, double beta, double i) { mpfr_set_d(r23753, alpha, MPFR_RNDN); ; mpfr_set_si(r23755, mpfr_cmp(r23753, r23754) <= 0, MPFR_RNDN); mpfr_set_d(r23756, i, MPFR_RNDN); mpfr_set_d(r23757, beta, MPFR_RNDN); mpfr_add(r23758, r23753, r23757, MPFR_RNDN); ; mpfr_mul(r23760, r23756, r23759, MPFR_RNDN); mpfr_add(r23761, r23758, r23760, MPFR_RNDN); mpfr_div(r23762, r23756, r23761, MPFR_RNDN); mpfr_add(r23763, r23756, r23757, MPFR_RNDN); mpfr_add(r23764, r23753, r23763, MPFR_RNDN); mpfr_mul(r23765, r23762, r23764, MPFR_RNDN); mpfr_mul(r23766, r23757, r23753, MPFR_RNDN); mpfr_add(r23767, r23758, r23756, MPFR_RNDN); mpfr_mul(r23768, r23756, r23767, MPFR_RNDN); mpfr_add(r23769, r23766, r23768, MPFR_RNDN); mpfr_mul(r23770, r23759, r23756, MPFR_RNDN); mpfr_add(r23771, r23758, r23770, MPFR_RNDN); mpfr_div(r23772, r23769, r23771, MPFR_RNDN); mpfr_mul(r23773, r23771, r23771, MPFR_RNDN); ; mpfr_sub(r23775, r23773, r23774, MPFR_RNDN); mpfr_div(r23776, r23772, r23775, MPFR_RNDN); mpfr_mul(r23777, r23765, r23776, MPFR_RNDN); ; if (mpfr_get_si(r23755, MPFR_RNDN)) { mpfr_set(r23779, r23777, MPFR_RNDN); } else { mpfr_set(r23779, r23778, MPFR_RNDN); }; return mpfr_get_d(r23779, MPFR_RNDN); } static mpfr_t r23780, r23781, r23782, r23783, r23784, r23785, r23786, r23787, r23788, r23789, r23790, r23791, r23792, r23793, r23794, r23795, r23796, r23797, r23798, r23799, r23800, r23801, r23802, r23803, r23804, r23805, r23806; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r23780); mpfr_init_set_str(r23781, "1.3406027096911442e+154", 10, MPFR_RNDN); mpfr_init(r23782); mpfr_init(r23783); mpfr_init(r23784); mpfr_init(r23785); mpfr_init_set_str(r23786, "2", 10, MPFR_RNDN); mpfr_init(r23787); mpfr_init(r23788); mpfr_init(r23789); mpfr_init(r23790); mpfr_init(r23791); mpfr_init(r23792); mpfr_init(r23793); mpfr_init(r23794); mpfr_init(r23795); mpfr_init(r23796); mpfr_init(r23797); mpfr_init(r23798); mpfr_init(r23799); mpfr_init(r23800); mpfr_init_set_str(r23801, "1.0", 10, MPFR_RNDN); mpfr_init(r23802); mpfr_init(r23803); mpfr_init(r23804); mpfr_init_set_str(r23805, "0", 10, MPFR_RNDN); mpfr_init(r23806); } double f_dm(double alpha, double beta, double i) { mpfr_set_d(r23780, alpha, MPFR_RNDN); ; mpfr_set_si(r23782, mpfr_cmp(r23780, r23781) <= 0, MPFR_RNDN); mpfr_set_d(r23783, i, MPFR_RNDN); mpfr_set_d(r23784, beta, MPFR_RNDN); mpfr_add(r23785, r23780, r23784, MPFR_RNDN); ; mpfr_mul(r23787, r23783, r23786, MPFR_RNDN); mpfr_add(r23788, r23785, r23787, MPFR_RNDN); mpfr_div(r23789, r23783, r23788, MPFR_RNDN); mpfr_add(r23790, r23783, r23784, MPFR_RNDN); mpfr_add(r23791, r23780, r23790, MPFR_RNDN); mpfr_mul(r23792, r23789, r23791, MPFR_RNDN); mpfr_mul(r23793, r23784, r23780, MPFR_RNDN); mpfr_add(r23794, r23785, r23783, MPFR_RNDN); mpfr_mul(r23795, r23783, r23794, MPFR_RNDN); mpfr_add(r23796, r23793, r23795, MPFR_RNDN); mpfr_mul(r23797, r23786, r23783, MPFR_RNDN); mpfr_add(r23798, r23785, r23797, MPFR_RNDN); mpfr_div(r23799, r23796, r23798, MPFR_RNDN); mpfr_mul(r23800, r23798, r23798, MPFR_RNDN); ; mpfr_sub(r23802, r23800, r23801, MPFR_RNDN); mpfr_div(r23803, r23799, r23802, MPFR_RNDN); mpfr_mul(r23804, r23792, r23803, MPFR_RNDN); ; if (mpfr_get_si(r23782, MPFR_RNDN)) { mpfr_set(r23806, r23804, MPFR_RNDN); } else { mpfr_set(r23806, r23805, MPFR_RNDN); }; return mpfr_get_d(r23806, MPFR_RNDN); }