#include #include #include #include #include char *name = "Octave 3.8, jcobi/4, as called"; double f_if(float i) { float r23716 = i; float r23717 = r23716 * r23716; float r23718 = r23717 * r23717; float r23719 = 2; float r23720 = r23719 * r23716; float r23721 = r23720 * r23720; float r23722 = r23718 / r23721; float r23723 = 1.0; float r23724 = r23721 - r23723; float r23725 = r23722 / r23724; return r23725; } double f_id(double i) { double r23726 = i; double r23727 = r23726 * r23726; double r23728 = r23727 * r23727; double r23729 = 2; double r23730 = r23729 * r23726; double r23731 = r23730 * r23730; double r23732 = r23728 / r23731; double r23733 = 1.0; double r23734 = r23731 - r23733; double r23735 = r23732 / r23734; return r23735; } double f_of(float i) { float r23736 = i; float r23737 = 2845.110039665851; bool r23738 = r23736 <= r23737; float r23739 = 1; float r23740 = 2; float r23741 = r23739 / r23740; float r23742 = r23741 * r23736; float r23743 = r23742 * r23742; float r23744 = r23740 * r23736; float r23745 = r23744 * r23744; float r23746 = 1.0; float r23747 = r23745 - r23746; float r23748 = r23743 / r23747; float r23749 = 0.015625; float r23750 = r23749 / r23736; float r23751 = r23750 / r23736; float r23752 = 0.00390625; float r23753 = 4; float r23754 = pow(r23736, r23753); float r23755 = r23752 / r23754; float r23756 = 1/16; float r23757 = r23755 + r23756; float r23758 = r23751 + r23757; float r23759 = r23738 ? r23748 : r23758; return r23759; } double f_od(double i) { double r23760 = i; double r23761 = 2845.110039665851; bool r23762 = r23760 <= r23761; double r23763 = 1; double r23764 = 2; double r23765 = r23763 / r23764; double r23766 = r23765 * r23760; double r23767 = r23766 * r23766; double r23768 = r23764 * r23760; double r23769 = r23768 * r23768; double r23770 = 1.0; double r23771 = r23769 - r23770; double r23772 = r23767 / r23771; double r23773 = 0.015625; double r23774 = r23773 / r23760; double r23775 = r23774 / r23760; double r23776 = 0.00390625; double r23777 = 4; double r23778 = pow(r23760, r23777); double r23779 = r23776 / r23778; double r23780 = 1/16; double r23781 = r23779 + r23780; double r23782 = r23775 + r23781; double r23783 = r23762 ? r23772 : r23782; return r23783; } 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 r23784, r23785, r23786, r23787, r23788, r23789, r23790, r23791, r23792, r23793; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r23784); mpfr_init(r23785); mpfr_init(r23786); mpfr_init_set_str(r23787, "2", 10, MPFR_RNDN); mpfr_init(r23788); mpfr_init(r23789); mpfr_init(r23790); mpfr_init_set_str(r23791, "1.0", 10, MPFR_RNDN); mpfr_init(r23792); mpfr_init(r23793); } double f_im(double i) { mpfr_set_d(r23784, i, MPFR_RNDN); mpfr_mul(r23785, r23784, r23784, MPFR_RNDN); mpfr_mul(r23786, r23785, r23785, MPFR_RNDN); ; mpfr_mul(r23788, r23787, r23784, MPFR_RNDN); mpfr_mul(r23789, r23788, r23788, MPFR_RNDN); mpfr_div(r23790, r23786, r23789, MPFR_RNDN); ; mpfr_sub(r23792, r23789, r23791, MPFR_RNDN); mpfr_div(r23793, r23790, r23792, MPFR_RNDN); return mpfr_get_d(r23793, MPFR_RNDN); } static mpfr_t r23794, r23795, r23796, r23797, r23798, r23799, r23800, r23801, r23802, r23803, r23804, r23805, r23806, r23807, r23808, r23809, r23810, r23811, r23812, r23813, r23814, r23815, r23816, r23817; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r23794); mpfr_init_set_str(r23795, "2845.110039665851", 10, MPFR_RNDN); mpfr_init(r23796); mpfr_init_set_str(r23797, "1", 10, MPFR_RNDN); mpfr_init_set_str(r23798, "2", 10, MPFR_RNDN); mpfr_init(r23799); mpfr_init(r23800); mpfr_init(r23801); mpfr_init(r23802); mpfr_init(r23803); mpfr_init_set_str(r23804, "1.0", 10, MPFR_RNDN); mpfr_init(r23805); mpfr_init(r23806); mpfr_init_set_str(r23807, "0.015625", 10, MPFR_RNDN); mpfr_init(r23808); mpfr_init(r23809); mpfr_init_set_str(r23810, "0.00390625", 10, MPFR_RNDN); mpfr_init_set_str(r23811, "4", 10, MPFR_RNDN); mpfr_init(r23812); mpfr_init(r23813); mpfr_init_set_str(r23814, "1/16", 10, MPFR_RNDN); mpfr_init(r23815); mpfr_init(r23816); mpfr_init(r23817); } double f_fm(double i) { mpfr_set_d(r23794, i, MPFR_RNDN); ; mpfr_set_si(r23796, mpfr_cmp(r23794, r23795) <= 0, MPFR_RNDN); ; ; mpfr_div(r23799, r23797, r23798, MPFR_RNDN); mpfr_mul(r23800, r23799, r23794, MPFR_RNDN); mpfr_mul(r23801, r23800, r23800, MPFR_RNDN); mpfr_mul(r23802, r23798, r23794, MPFR_RNDN); mpfr_mul(r23803, r23802, r23802, MPFR_RNDN); ; mpfr_sub(r23805, r23803, r23804, MPFR_RNDN); mpfr_div(r23806, r23801, r23805, MPFR_RNDN); ; mpfr_div(r23808, r23807, r23794, MPFR_RNDN); mpfr_div(r23809, r23808, r23794, MPFR_RNDN); ; ; mpfr_pow(r23812, r23794, r23811, MPFR_RNDN); mpfr_div(r23813, r23810, r23812, MPFR_RNDN); ; mpfr_add(r23815, r23813, r23814, MPFR_RNDN); mpfr_add(r23816, r23809, r23815, MPFR_RNDN); if (mpfr_get_si(r23796, MPFR_RNDN)) { mpfr_set(r23817, r23806, MPFR_RNDN); } else { mpfr_set(r23817, r23816, MPFR_RNDN); }; return mpfr_get_d(r23817, MPFR_RNDN); } static mpfr_t r23818, r23819, r23820, r23821, r23822, r23823, r23824, r23825, r23826, r23827, r23828, r23829, r23830, r23831, r23832, r23833, r23834, r23835, r23836, r23837, r23838, r23839, r23840, r23841; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r23818); mpfr_init_set_str(r23819, "2845.110039665851", 10, MPFR_RNDN); mpfr_init(r23820); mpfr_init_set_str(r23821, "1", 10, MPFR_RNDN); mpfr_init_set_str(r23822, "2", 10, MPFR_RNDN); mpfr_init(r23823); mpfr_init(r23824); mpfr_init(r23825); mpfr_init(r23826); mpfr_init(r23827); mpfr_init_set_str(r23828, "1.0", 10, MPFR_RNDN); mpfr_init(r23829); mpfr_init(r23830); mpfr_init_set_str(r23831, "0.015625", 10, MPFR_RNDN); mpfr_init(r23832); mpfr_init(r23833); mpfr_init_set_str(r23834, "0.00390625", 10, MPFR_RNDN); mpfr_init_set_str(r23835, "4", 10, MPFR_RNDN); mpfr_init(r23836); mpfr_init(r23837); mpfr_init_set_str(r23838, "1/16", 10, MPFR_RNDN); mpfr_init(r23839); mpfr_init(r23840); mpfr_init(r23841); } double f_dm(double i) { mpfr_set_d(r23818, i, MPFR_RNDN); ; mpfr_set_si(r23820, mpfr_cmp(r23818, r23819) <= 0, MPFR_RNDN); ; ; mpfr_div(r23823, r23821, r23822, MPFR_RNDN); mpfr_mul(r23824, r23823, r23818, MPFR_RNDN); mpfr_mul(r23825, r23824, r23824, MPFR_RNDN); mpfr_mul(r23826, r23822, r23818, MPFR_RNDN); mpfr_mul(r23827, r23826, r23826, MPFR_RNDN); ; mpfr_sub(r23829, r23827, r23828, MPFR_RNDN); mpfr_div(r23830, r23825, r23829, MPFR_RNDN); ; mpfr_div(r23832, r23831, r23818, MPFR_RNDN); mpfr_div(r23833, r23832, r23818, MPFR_RNDN); ; ; mpfr_pow(r23836, r23818, r23835, MPFR_RNDN); mpfr_div(r23837, r23834, r23836, MPFR_RNDN); ; mpfr_add(r23839, r23837, r23838, MPFR_RNDN); mpfr_add(r23840, r23833, r23839, MPFR_RNDN); if (mpfr_get_si(r23820, MPFR_RNDN)) { mpfr_set(r23841, r23830, MPFR_RNDN); } else { mpfr_set(r23841, r23840, MPFR_RNDN); }; return mpfr_get_d(r23841, MPFR_RNDN); }