#include #include #include #include #include char *name = "kepler0"; double f_if(float x1, float x2, float x3, float x4, float x5, float x6) { float r24841 = x2; float r24842 = x5; float r24843 = r24841 * r24842; float r24844 = x3; float r24845 = x6; float r24846 = r24844 * r24845; float r24847 = r24843 + r24846; float r24848 = r24841 * r24844; float r24849 = r24847 - r24848; float r24850 = r24842 * r24845; float r24851 = r24849 - r24850; float r24852 = x1; float r24853 = -r24852; float r24854 = r24853 + r24841; float r24855 = r24854 + r24844; float r24856 = x4; float r24857 = r24855 - r24856; float r24858 = r24857 + r24842; float r24859 = r24858 + r24845; float r24860 = r24852 * r24859; float r24861 = r24851 + r24860; return r24861; } double f_id(double x1, double x2, double x3, double x4, double x5, double x6) { double r24862 = x2; double r24863 = x5; double r24864 = r24862 * r24863; double r24865 = x3; double r24866 = x6; double r24867 = r24865 * r24866; double r24868 = r24864 + r24867; double r24869 = r24862 * r24865; double r24870 = r24868 - r24869; double r24871 = r24863 * r24866; double r24872 = r24870 - r24871; double r24873 = x1; double r24874 = -r24873; double r24875 = r24874 + r24862; double r24876 = r24875 + r24865; double r24877 = x4; double r24878 = r24876 - r24877; double r24879 = r24878 + r24863; double r24880 = r24879 + r24866; double r24881 = r24873 * r24880; double r24882 = r24872 + r24881; return r24882; } double f_of(float x1, float x2, float x3, float x4, float x5, float x6) { float r24883 = x2; float r24884 = x4; float r24885 = x3; float r24886 = r24884 - r24885; float r24887 = r24883 - r24886; float r24888 = x1; float r24889 = -r24888; float r24890 = x5; float r24891 = x6; float r24892 = r24890 + r24891; float r24893 = r24889 + r24892; float r24894 = r24887 + r24893; float r24895 = r24894 * r24888; float r24896 = 3; float r24897 = pow(r24885, r24896); float r24898 = pow(r24890, r24896); float r24899 = r24897 - r24898; float r24900 = r24891 * r24891; float r24901 = r24883 * r24883; float r24902 = r24900 - r24901; float r24903 = r24899 * r24902; float r24904 = r24885 * r24885; float r24905 = r24885 + r24890; float r24906 = r24905 * r24890; float r24907 = r24904 + r24906; float r24908 = r24883 + r24891; float r24909 = r24907 * r24908; float r24910 = r24903 / r24909; float r24911 = r24895 + r24910; return r24911; } double f_od(double x1, double x2, double x3, double x4, double x5, double x6) { double r24912 = x2; double r24913 = x4; double r24914 = x3; double r24915 = r24913 - r24914; double r24916 = r24912 - r24915; double r24917 = x1; double r24918 = -r24917; double r24919 = x5; double r24920 = x6; double r24921 = r24919 + r24920; double r24922 = r24918 + r24921; double r24923 = r24916 + r24922; double r24924 = r24923 * r24917; double r24925 = 3; double r24926 = pow(r24914, r24925); double r24927 = pow(r24919, r24925); double r24928 = r24926 - r24927; double r24929 = r24920 * r24920; double r24930 = r24912 * r24912; double r24931 = r24929 - r24930; double r24932 = r24928 * r24931; double r24933 = r24914 * r24914; double r24934 = r24914 + r24919; double r24935 = r24934 * r24919; double r24936 = r24933 + r24935; double r24937 = r24912 + r24920; double r24938 = r24936 * r24937; double r24939 = r24932 / r24938; double r24940 = r24924 + r24939; return r24940; } 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 r24941, r24942, r24943, r24944, r24945, r24946, r24947, r24948, r24949, r24950, r24951, r24952, r24953, r24954, r24955, r24956, r24957, r24958, r24959, r24960, r24961; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r24941); mpfr_init(r24942); mpfr_init(r24943); mpfr_init(r24944); mpfr_init(r24945); mpfr_init(r24946); mpfr_init(r24947); mpfr_init(r24948); mpfr_init(r24949); mpfr_init(r24950); mpfr_init(r24951); mpfr_init(r24952); mpfr_init(r24953); mpfr_init(r24954); mpfr_init(r24955); mpfr_init(r24956); mpfr_init(r24957); mpfr_init(r24958); mpfr_init(r24959); mpfr_init(r24960); mpfr_init(r24961); } double f_im(double x1, double x2, double x3, double x4, double x5, double x6) { mpfr_set_d(r24941, x2, MPFR_RNDN); mpfr_set_d(r24942, x5, MPFR_RNDN); mpfr_mul(r24943, r24941, r24942, MPFR_RNDN); mpfr_set_d(r24944, x3, MPFR_RNDN); mpfr_set_d(r24945, x6, MPFR_RNDN); mpfr_mul(r24946, r24944, r24945, MPFR_RNDN); mpfr_add(r24947, r24943, r24946, MPFR_RNDN); mpfr_mul(r24948, r24941, r24944, MPFR_RNDN); mpfr_sub(r24949, r24947, r24948, MPFR_RNDN); mpfr_mul(r24950, r24942, r24945, MPFR_RNDN); mpfr_sub(r24951, r24949, r24950, MPFR_RNDN); mpfr_set_d(r24952, x1, MPFR_RNDN); mpfr_neg(r24953, r24952, MPFR_RNDN); mpfr_add(r24954, r24953, r24941, MPFR_RNDN); mpfr_add(r24955, r24954, r24944, MPFR_RNDN); mpfr_set_d(r24956, x4, MPFR_RNDN); mpfr_sub(r24957, r24955, r24956, MPFR_RNDN); mpfr_add(r24958, r24957, r24942, MPFR_RNDN); mpfr_add(r24959, r24958, r24945, MPFR_RNDN); mpfr_mul(r24960, r24952, r24959, MPFR_RNDN); mpfr_add(r24961, r24951, r24960, MPFR_RNDN); return mpfr_get_d(r24961, MPFR_RNDN); } static mpfr_t r24962, r24963, r24964, r24965, r24966, r24967, r24968, r24969, r24970, r24971, r24972, r24973, r24974, r24975, r24976, r24977, r24978, r24979, r24980, r24981, r24982, r24983, r24984, r24985, r24986, r24987, r24988, r24989, r24990; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r24962); mpfr_init(r24963); mpfr_init(r24964); mpfr_init(r24965); mpfr_init(r24966); mpfr_init(r24967); mpfr_init(r24968); mpfr_init(r24969); mpfr_init(r24970); mpfr_init(r24971); mpfr_init(r24972); mpfr_init(r24973); mpfr_init(r24974); mpfr_init_set_str(r24975, "3", 10, MPFR_RNDN); mpfr_init(r24976); mpfr_init(r24977); mpfr_init(r24978); mpfr_init(r24979); mpfr_init(r24980); mpfr_init(r24981); mpfr_init(r24982); mpfr_init(r24983); mpfr_init(r24984); mpfr_init(r24985); mpfr_init(r24986); mpfr_init(r24987); mpfr_init(r24988); mpfr_init(r24989); mpfr_init(r24990); } double f_fm(double x1, double x2, double x3, double x4, double x5, double x6) { mpfr_set_d(r24962, x2, MPFR_RNDN); mpfr_set_d(r24963, x4, MPFR_RNDN); mpfr_set_d(r24964, x3, MPFR_RNDN); mpfr_sub(r24965, r24963, r24964, MPFR_RNDN); mpfr_sub(r24966, r24962, r24965, MPFR_RNDN); mpfr_set_d(r24967, x1, MPFR_RNDN); mpfr_neg(r24968, r24967, MPFR_RNDN); mpfr_set_d(r24969, x5, MPFR_RNDN); mpfr_set_d(r24970, x6, MPFR_RNDN); mpfr_add(r24971, r24969, r24970, MPFR_RNDN); mpfr_add(r24972, r24968, r24971, MPFR_RNDN); mpfr_add(r24973, r24966, r24972, MPFR_RNDN); mpfr_mul(r24974, r24973, r24967, MPFR_RNDN); ; mpfr_pow(r24976, r24964, r24975, MPFR_RNDN); mpfr_pow(r24977, r24969, r24975, MPFR_RNDN); mpfr_sub(r24978, r24976, r24977, MPFR_RNDN); mpfr_mul(r24979, r24970, r24970, MPFR_RNDN); mpfr_mul(r24980, r24962, r24962, MPFR_RNDN); mpfr_sub(r24981, r24979, r24980, MPFR_RNDN); mpfr_mul(r24982, r24978, r24981, MPFR_RNDN); mpfr_mul(r24983, r24964, r24964, MPFR_RNDN); mpfr_add(r24984, r24964, r24969, MPFR_RNDN); mpfr_mul(r24985, r24984, r24969, MPFR_RNDN); mpfr_add(r24986, r24983, r24985, MPFR_RNDN); mpfr_add(r24987, r24962, r24970, MPFR_RNDN); mpfr_mul(r24988, r24986, r24987, MPFR_RNDN); mpfr_div(r24989, r24982, r24988, MPFR_RNDN); mpfr_add(r24990, r24974, r24989, MPFR_RNDN); return mpfr_get_d(r24990, MPFR_RNDN); } static mpfr_t r24991, r24992, r24993, r24994, r24995, r24996, r24997, r24998, r24999, r25000, r25001, r25002, r25003, r25004, r25005, r25006, r25007, r25008, r25009, r25010, r25011, r25012, r25013, r25014, r25015, r25016, r25017, r25018, r25019; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r24991); mpfr_init(r24992); mpfr_init(r24993); mpfr_init(r24994); mpfr_init(r24995); mpfr_init(r24996); mpfr_init(r24997); mpfr_init(r24998); mpfr_init(r24999); mpfr_init(r25000); mpfr_init(r25001); mpfr_init(r25002); mpfr_init(r25003); mpfr_init_set_str(r25004, "3", 10, MPFR_RNDN); mpfr_init(r25005); mpfr_init(r25006); mpfr_init(r25007); mpfr_init(r25008); mpfr_init(r25009); mpfr_init(r25010); mpfr_init(r25011); mpfr_init(r25012); mpfr_init(r25013); mpfr_init(r25014); mpfr_init(r25015); mpfr_init(r25016); mpfr_init(r25017); mpfr_init(r25018); mpfr_init(r25019); } double f_dm(double x1, double x2, double x3, double x4, double x5, double x6) { mpfr_set_d(r24991, x2, MPFR_RNDN); mpfr_set_d(r24992, x4, MPFR_RNDN); mpfr_set_d(r24993, x3, MPFR_RNDN); mpfr_sub(r24994, r24992, r24993, MPFR_RNDN); mpfr_sub(r24995, r24991, r24994, MPFR_RNDN); mpfr_set_d(r24996, x1, MPFR_RNDN); mpfr_neg(r24997, r24996, MPFR_RNDN); mpfr_set_d(r24998, x5, MPFR_RNDN); mpfr_set_d(r24999, x6, MPFR_RNDN); mpfr_add(r25000, r24998, r24999, MPFR_RNDN); mpfr_add(r25001, r24997, r25000, MPFR_RNDN); mpfr_add(r25002, r24995, r25001, MPFR_RNDN); mpfr_mul(r25003, r25002, r24996, MPFR_RNDN); ; mpfr_pow(r25005, r24993, r25004, MPFR_RNDN); mpfr_pow(r25006, r24998, r25004, MPFR_RNDN); mpfr_sub(r25007, r25005, r25006, MPFR_RNDN); mpfr_mul(r25008, r24999, r24999, MPFR_RNDN); mpfr_mul(r25009, r24991, r24991, MPFR_RNDN); mpfr_sub(r25010, r25008, r25009, MPFR_RNDN); mpfr_mul(r25011, r25007, r25010, MPFR_RNDN); mpfr_mul(r25012, r24993, r24993, MPFR_RNDN); mpfr_add(r25013, r24993, r24998, MPFR_RNDN); mpfr_mul(r25014, r25013, r24998, MPFR_RNDN); mpfr_add(r25015, r25012, r25014, MPFR_RNDN); mpfr_add(r25016, r24991, r24999, MPFR_RNDN); mpfr_mul(r25017, r25015, r25016, MPFR_RNDN); mpfr_div(r25018, r25011, r25017, MPFR_RNDN); mpfr_add(r25019, r25003, r25018, MPFR_RNDN); return mpfr_get_d(r25019, MPFR_RNDN); }