#include #include #include #include #include char *name = "Maksimov and Kolovsky, Equation (4)"; double f_if(float J, float l, float K, float U) { float r25911 = J; float r25912 = l; float r25913 = exp(r25912); float r25914 = -r25912; float r25915 = exp(r25914); float r25916 = r25913 - r25915; float r25917 = r25911 * r25916; float r25918 = K; float r25919 = 2; float r25920 = r25918 / r25919; float r25921 = cos(r25920); float r25922 = r25917 * r25921; float r25923 = U; float r25924 = r25922 + r25923; return r25924; } double f_id(double J, double l, double K, double U) { double r25925 = J; double r25926 = l; double r25927 = exp(r25926); double r25928 = -r25926; double r25929 = exp(r25928); double r25930 = r25927 - r25929; double r25931 = r25925 * r25930; double r25932 = K; double r25933 = 2; double r25934 = r25932 / r25933; double r25935 = cos(r25934); double r25936 = r25931 * r25935; double r25937 = U; double r25938 = r25936 + r25937; return r25938; } double f_of(float J, float l, float K, float U) { float r25939 = K; float r25940 = 1/2; float r25941 = r25939 * r25940; float r25942 = cos(r25941); float r25943 = l; float r25944 = r25943 * r25943; float r25945 = 1/3; float r25946 = r25944 * r25945; float r25947 = J; float r25948 = r25947 * r25943; float r25949 = 2; float r25950 = r25949 * r25947; float r25951 = r25950 * r25943; float r25952 = fma(r25946, r25948, r25951); float r25953 = U; float r25954 = fma(r25942, r25952, r25953); return r25954; } double f_od(double J, double l, double K, double U) { double r25955 = K; double r25956 = 1/2; double r25957 = r25955 * r25956; double r25958 = cos(r25957); double r25959 = l; double r25960 = r25959 * r25959; double r25961 = 1/3; double r25962 = r25960 * r25961; double r25963 = J; double r25964 = r25963 * r25959; double r25965 = 2; double r25966 = r25965 * r25963; double r25967 = r25966 * r25959; double r25968 = fma(r25962, r25964, r25967); double r25969 = U; double r25970 = fma(r25958, r25968, r25969); return r25970; } 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 r25971, r25972, r25973, r25974, r25975, r25976, r25977, r25978, r25979, r25980, r25981, r25982, r25983, r25984; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r25971); mpfr_init(r25972); mpfr_init(r25973); mpfr_init(r25974); mpfr_init(r25975); mpfr_init(r25976); mpfr_init(r25977); mpfr_init(r25978); mpfr_init_set_str(r25979, "2", 10, MPFR_RNDN); mpfr_init(r25980); mpfr_init(r25981); mpfr_init(r25982); mpfr_init(r25983); mpfr_init(r25984); } double f_im(double J, double l, double K, double U) { mpfr_set_d(r25971, J, MPFR_RNDN); mpfr_set_d(r25972, l, MPFR_RNDN); mpfr_exp(r25973, r25972, MPFR_RNDN); mpfr_neg(r25974, r25972, MPFR_RNDN); mpfr_exp(r25975, r25974, MPFR_RNDN); mpfr_sub(r25976, r25973, r25975, MPFR_RNDN); mpfr_mul(r25977, r25971, r25976, MPFR_RNDN); mpfr_set_d(r25978, K, MPFR_RNDN); ; mpfr_div(r25980, r25978, r25979, MPFR_RNDN); mpfr_cos(r25981, r25980, MPFR_RNDN); mpfr_mul(r25982, r25977, r25981, MPFR_RNDN); mpfr_set_d(r25983, U, MPFR_RNDN); mpfr_add(r25984, r25982, r25983, MPFR_RNDN); return mpfr_get_d(r25984, MPFR_RNDN); } static mpfr_t r25985, r25986, r25987, r25988, r25989, r25990, r25991, r25992, r25993, r25994, r25995, r25996, r25997, r25998, r25999, r26000; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r25985); mpfr_init_set_str(r25986, "1/2", 10, MPFR_RNDN); mpfr_init(r25987); mpfr_init(r25988); mpfr_init(r25989); mpfr_init(r25990); mpfr_init_set_str(r25991, "1/3", 10, MPFR_RNDN); mpfr_init(r25992); mpfr_init(r25993); mpfr_init(r25994); mpfr_init_set_str(r25995, "2", 10, MPFR_RNDN); mpfr_init(r25996); mpfr_init(r25997); mpfr_init(r25998); mpfr_init(r25999); mpfr_init(r26000); } double f_fm(double J, double l, double K, double U) { mpfr_set_d(r25985, K, MPFR_RNDN); ; mpfr_mul(r25987, r25985, r25986, MPFR_RNDN); mpfr_cos(r25988, r25987, MPFR_RNDN); mpfr_set_d(r25989, l, MPFR_RNDN); mpfr_mul(r25990, r25989, r25989, MPFR_RNDN); ; mpfr_mul(r25992, r25990, r25991, MPFR_RNDN); mpfr_set_d(r25993, J, MPFR_RNDN); mpfr_mul(r25994, r25993, r25989, MPFR_RNDN); ; mpfr_mul(r25996, r25995, r25993, MPFR_RNDN); mpfr_mul(r25997, r25996, r25989, MPFR_RNDN); mpfr_fma(r25998, r25992, r25994, r25997, MPFR_RNDN); mpfr_set_d(r25999, U, MPFR_RNDN); mpfr_fma(r26000, r25988, r25998, r25999, MPFR_RNDN); return mpfr_get_d(r26000, MPFR_RNDN); } static mpfr_t r26001, r26002, r26003, r26004, r26005, r26006, r26007, r26008, r26009, r26010, r26011, r26012, r26013, r26014, r26015, r26016; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r26001); mpfr_init_set_str(r26002, "1/2", 10, MPFR_RNDN); mpfr_init(r26003); mpfr_init(r26004); mpfr_init(r26005); mpfr_init(r26006); mpfr_init_set_str(r26007, "1/3", 10, MPFR_RNDN); mpfr_init(r26008); mpfr_init(r26009); mpfr_init(r26010); mpfr_init_set_str(r26011, "2", 10, MPFR_RNDN); mpfr_init(r26012); mpfr_init(r26013); mpfr_init(r26014); mpfr_init(r26015); mpfr_init(r26016); } double f_dm(double J, double l, double K, double U) { mpfr_set_d(r26001, K, MPFR_RNDN); ; mpfr_mul(r26003, r26001, r26002, MPFR_RNDN); mpfr_cos(r26004, r26003, MPFR_RNDN); mpfr_set_d(r26005, l, MPFR_RNDN); mpfr_mul(r26006, r26005, r26005, MPFR_RNDN); ; mpfr_mul(r26008, r26006, r26007, MPFR_RNDN); mpfr_set_d(r26009, J, MPFR_RNDN); mpfr_mul(r26010, r26009, r26005, MPFR_RNDN); ; mpfr_mul(r26012, r26011, r26009, MPFR_RNDN); mpfr_mul(r26013, r26012, r26005, MPFR_RNDN); mpfr_fma(r26014, r26008, r26010, r26013, MPFR_RNDN); mpfr_set_d(r26015, U, MPFR_RNDN); mpfr_fma(r26016, r26004, r26014, r26015, MPFR_RNDN); return mpfr_get_d(r26016, MPFR_RNDN); }