#include #include #include #include #include char *name = "kepler1"; double f_if(float x1, float x2, float x3, float x4) { float r24645 = x1; float r24646 = x4; float r24647 = r24645 * r24646; float r24648 = -r24645; float r24649 = x2; float r24650 = r24648 + r24649; float r24651 = x3; float r24652 = r24650 + r24651; float r24653 = r24652 - r24646; float r24654 = r24647 * r24653; float r24655 = r24645 - r24649; float r24656 = r24655 + r24651; float r24657 = r24656 + r24646; float r24658 = r24649 * r24657; float r24659 = r24654 + r24658; float r24660 = r24645 + r24649; float r24661 = r24660 - r24651; float r24662 = r24661 + r24646; float r24663 = r24651 * r24662; float r24664 = r24659 + r24663; float r24665 = r24649 * r24651; float r24666 = r24665 * r24646; float r24667 = r24664 - r24666; float r24668 = r24645 * r24651; float r24669 = r24667 - r24668; float r24670 = r24645 * r24649; float r24671 = r24669 - r24670; float r24672 = r24671 - r24646; return r24672; } double f_id(double x1, double x2, double x3, double x4) { double r24673 = x1; double r24674 = x4; double r24675 = r24673 * r24674; double r24676 = -r24673; double r24677 = x2; double r24678 = r24676 + r24677; double r24679 = x3; double r24680 = r24678 + r24679; double r24681 = r24680 - r24674; double r24682 = r24675 * r24681; double r24683 = r24673 - r24677; double r24684 = r24683 + r24679; double r24685 = r24684 + r24674; double r24686 = r24677 * r24685; double r24687 = r24682 + r24686; double r24688 = r24673 + r24677; double r24689 = r24688 - r24679; double r24690 = r24689 + r24674; double r24691 = r24679 * r24690; double r24692 = r24687 + r24691; double r24693 = r24677 * r24679; double r24694 = r24693 * r24674; double r24695 = r24692 - r24694; double r24696 = r24673 * r24679; double r24697 = r24695 - r24696; double r24698 = r24673 * r24677; double r24699 = r24697 - r24698; double r24700 = r24699 - r24674; return r24700; } double f_of(float x1, float x2, float x3, float x4) { float r24701 = x1; float r24702 = x4; float r24703 = r24701 * r24702; float r24704 = x2; float r24705 = x3; float r24706 = r24705 - r24702; float r24707 = r24704 + r24706; float r24708 = -r24701; float r24709 = r24707 + r24708; float r24710 = r24703 * r24709; float r24711 = r24704 * r24702; float r24712 = r24701 - r24711; float r24713 = r24706 - r24704; float r24714 = r24712 - r24713; float r24715 = r24714 * r24705; float r24716 = r24710 + r24715; float r24717 = r24705 + r24702; float r24718 = r24701 - r24704; float r24719 = r24717 + r24718; float r24720 = r24704 * r24719; float r24721 = r24716 + r24720; float r24722 = 3; float r24723 = pow(r24721, r24722); float r24724 = cbrt(r24723); float r24725 = r24704 + r24705; float r24726 = r24725 * r24701; float r24727 = r24726 + r24702; float r24728 = r24724 - r24727; return r24728; } double f_od(double x1, double x2, double x3, double x4) { double r24729 = x1; double r24730 = x4; double r24731 = r24729 * r24730; double r24732 = x2; double r24733 = x3; double r24734 = r24733 - r24730; double r24735 = r24732 + r24734; double r24736 = -r24729; double r24737 = r24735 + r24736; double r24738 = r24731 * r24737; double r24739 = r24732 * r24730; double r24740 = r24729 - r24739; double r24741 = r24734 - r24732; double r24742 = r24740 - r24741; double r24743 = r24742 * r24733; double r24744 = r24738 + r24743; double r24745 = r24733 + r24730; double r24746 = r24729 - r24732; double r24747 = r24745 + r24746; double r24748 = r24732 * r24747; double r24749 = r24744 + r24748; double r24750 = 3; double r24751 = pow(r24749, r24750); double r24752 = cbrt(r24751); double r24753 = r24732 + r24733; double r24754 = r24753 * r24729; double r24755 = r24754 + r24730; double r24756 = r24752 - r24755; return r24756; } 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 r24757, r24758, r24759, r24760, r24761, r24762, r24763, r24764, r24765, r24766, r24767, r24768, r24769, r24770, r24771, r24772, r24773, r24774, r24775, r24776, r24777, r24778, r24779, r24780, r24781, r24782, r24783, r24784; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r24757); mpfr_init(r24758); mpfr_init(r24759); mpfr_init(r24760); mpfr_init(r24761); mpfr_init(r24762); mpfr_init(r24763); mpfr_init(r24764); mpfr_init(r24765); mpfr_init(r24766); mpfr_init(r24767); mpfr_init(r24768); mpfr_init(r24769); mpfr_init(r24770); mpfr_init(r24771); mpfr_init(r24772); mpfr_init(r24773); mpfr_init(r24774); mpfr_init(r24775); mpfr_init(r24776); mpfr_init(r24777); mpfr_init(r24778); mpfr_init(r24779); mpfr_init(r24780); mpfr_init(r24781); mpfr_init(r24782); mpfr_init(r24783); mpfr_init(r24784); } double f_im(double x1, double x2, double x3, double x4) { mpfr_set_d(r24757, x1, MPFR_RNDN); mpfr_set_d(r24758, x4, MPFR_RNDN); mpfr_mul(r24759, r24757, r24758, MPFR_RNDN); mpfr_neg(r24760, r24757, MPFR_RNDN); mpfr_set_d(r24761, x2, MPFR_RNDN); mpfr_add(r24762, r24760, r24761, MPFR_RNDN); mpfr_set_d(r24763, x3, MPFR_RNDN); mpfr_add(r24764, r24762, r24763, MPFR_RNDN); mpfr_sub(r24765, r24764, r24758, MPFR_RNDN); mpfr_mul(r24766, r24759, r24765, MPFR_RNDN); mpfr_sub(r24767, r24757, r24761, MPFR_RNDN); mpfr_add(r24768, r24767, r24763, MPFR_RNDN); mpfr_add(r24769, r24768, r24758, MPFR_RNDN); mpfr_mul(r24770, r24761, r24769, MPFR_RNDN); mpfr_add(r24771, r24766, r24770, MPFR_RNDN); mpfr_add(r24772, r24757, r24761, MPFR_RNDN); mpfr_sub(r24773, r24772, r24763, MPFR_RNDN); mpfr_add(r24774, r24773, r24758, MPFR_RNDN); mpfr_mul(r24775, r24763, r24774, MPFR_RNDN); mpfr_add(r24776, r24771, r24775, MPFR_RNDN); mpfr_mul(r24777, r24761, r24763, MPFR_RNDN); mpfr_mul(r24778, r24777, r24758, MPFR_RNDN); mpfr_sub(r24779, r24776, r24778, MPFR_RNDN); mpfr_mul(r24780, r24757, r24763, MPFR_RNDN); mpfr_sub(r24781, r24779, r24780, MPFR_RNDN); mpfr_mul(r24782, r24757, r24761, MPFR_RNDN); mpfr_sub(r24783, r24781, r24782, MPFR_RNDN); mpfr_sub(r24784, r24783, r24758, MPFR_RNDN); return mpfr_get_d(r24784, MPFR_RNDN); } static mpfr_t r24785, r24786, r24787, r24788, r24789, r24790, r24791, r24792, r24793, r24794, r24795, r24796, r24797, r24798, r24799, r24800, r24801, r24802, r24803, r24804, r24805, r24806, r24807, r24808, r24809, r24810, r24811, r24812; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r24785); mpfr_init(r24786); mpfr_init(r24787); mpfr_init(r24788); mpfr_init(r24789); mpfr_init(r24790); mpfr_init(r24791); mpfr_init(r24792); mpfr_init(r24793); mpfr_init(r24794); mpfr_init(r24795); mpfr_init(r24796); mpfr_init(r24797); mpfr_init(r24798); mpfr_init(r24799); mpfr_init(r24800); mpfr_init(r24801); mpfr_init(r24802); mpfr_init(r24803); mpfr_init(r24804); mpfr_init(r24805); mpfr_init_set_str(r24806, "3", 10, MPFR_RNDN); mpfr_init(r24807); mpfr_init(r24808); mpfr_init(r24809); mpfr_init(r24810); mpfr_init(r24811); mpfr_init(r24812); } double f_fm(double x1, double x2, double x3, double x4) { mpfr_set_d(r24785, x1, MPFR_RNDN); mpfr_set_d(r24786, x4, MPFR_RNDN); mpfr_mul(r24787, r24785, r24786, MPFR_RNDN); mpfr_set_d(r24788, x2, MPFR_RNDN); mpfr_set_d(r24789, x3, MPFR_RNDN); mpfr_sub(r24790, r24789, r24786, MPFR_RNDN); mpfr_add(r24791, r24788, r24790, MPFR_RNDN); mpfr_neg(r24792, r24785, MPFR_RNDN); mpfr_add(r24793, r24791, r24792, MPFR_RNDN); mpfr_mul(r24794, r24787, r24793, MPFR_RNDN); mpfr_mul(r24795, r24788, r24786, MPFR_RNDN); mpfr_sub(r24796, r24785, r24795, MPFR_RNDN); mpfr_sub(r24797, r24790, r24788, MPFR_RNDN); mpfr_sub(r24798, r24796, r24797, MPFR_RNDN); mpfr_mul(r24799, r24798, r24789, MPFR_RNDN); mpfr_add(r24800, r24794, r24799, MPFR_RNDN); mpfr_add(r24801, r24789, r24786, MPFR_RNDN); mpfr_sub(r24802, r24785, r24788, MPFR_RNDN); mpfr_add(r24803, r24801, r24802, MPFR_RNDN); mpfr_mul(r24804, r24788, r24803, MPFR_RNDN); mpfr_add(r24805, r24800, r24804, MPFR_RNDN); ; mpfr_pow(r24807, r24805, r24806, MPFR_RNDN); mpfr_cbrt(r24808, r24807, MPFR_RNDN); mpfr_add(r24809, r24788, r24789, MPFR_RNDN); mpfr_mul(r24810, r24809, r24785, MPFR_RNDN); mpfr_add(r24811, r24810, r24786, MPFR_RNDN); mpfr_sub(r24812, r24808, r24811, MPFR_RNDN); return mpfr_get_d(r24812, MPFR_RNDN); } static mpfr_t r24813, r24814, r24815, r24816, r24817, r24818, r24819, r24820, r24821, r24822, r24823, r24824, r24825, r24826, r24827, r24828, r24829, r24830, r24831, r24832, r24833, r24834, r24835, r24836, r24837, r24838, r24839, r24840; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r24813); mpfr_init(r24814); mpfr_init(r24815); mpfr_init(r24816); mpfr_init(r24817); mpfr_init(r24818); mpfr_init(r24819); mpfr_init(r24820); mpfr_init(r24821); mpfr_init(r24822); mpfr_init(r24823); mpfr_init(r24824); mpfr_init(r24825); mpfr_init(r24826); mpfr_init(r24827); mpfr_init(r24828); mpfr_init(r24829); mpfr_init(r24830); mpfr_init(r24831); mpfr_init(r24832); mpfr_init(r24833); mpfr_init_set_str(r24834, "3", 10, MPFR_RNDN); mpfr_init(r24835); mpfr_init(r24836); mpfr_init(r24837); mpfr_init(r24838); mpfr_init(r24839); mpfr_init(r24840); } double f_dm(double x1, double x2, double x3, double x4) { mpfr_set_d(r24813, x1, MPFR_RNDN); mpfr_set_d(r24814, x4, MPFR_RNDN); mpfr_mul(r24815, r24813, r24814, MPFR_RNDN); mpfr_set_d(r24816, x2, MPFR_RNDN); mpfr_set_d(r24817, x3, MPFR_RNDN); mpfr_sub(r24818, r24817, r24814, MPFR_RNDN); mpfr_add(r24819, r24816, r24818, MPFR_RNDN); mpfr_neg(r24820, r24813, MPFR_RNDN); mpfr_add(r24821, r24819, r24820, MPFR_RNDN); mpfr_mul(r24822, r24815, r24821, MPFR_RNDN); mpfr_mul(r24823, r24816, r24814, MPFR_RNDN); mpfr_sub(r24824, r24813, r24823, MPFR_RNDN); mpfr_sub(r24825, r24818, r24816, MPFR_RNDN); mpfr_sub(r24826, r24824, r24825, MPFR_RNDN); mpfr_mul(r24827, r24826, r24817, MPFR_RNDN); mpfr_add(r24828, r24822, r24827, MPFR_RNDN); mpfr_add(r24829, r24817, r24814, MPFR_RNDN); mpfr_sub(r24830, r24813, r24816, MPFR_RNDN); mpfr_add(r24831, r24829, r24830, MPFR_RNDN); mpfr_mul(r24832, r24816, r24831, MPFR_RNDN); mpfr_add(r24833, r24828, r24832, MPFR_RNDN); ; mpfr_pow(r24835, r24833, r24834, MPFR_RNDN); mpfr_cbrt(r24836, r24835, MPFR_RNDN); mpfr_add(r24837, r24816, r24817, MPFR_RNDN); mpfr_mul(r24838, r24837, r24813, MPFR_RNDN); mpfr_add(r24839, r24838, r24814, MPFR_RNDN); mpfr_sub(r24840, r24836, r24839, MPFR_RNDN); return mpfr_get_d(r24840, MPFR_RNDN); }