#include #include #include #include #include char *name = "Maksimov and Kolovsky, Equation (3)"; double f_if(float J, float K, float U) { float r25847 = -2; float r25848 = J; float r25849 = r25847 * r25848; float r25850 = K; float r25851 = 2; float r25852 = r25850 / r25851; float r25853 = cos(r25852); float r25854 = r25849 * r25853; float r25855 = 1; float r25856 = U; float r25857 = r25851 * r25848; float r25858 = r25857 * r25853; float r25859 = r25856 / r25858; float r25860 = pow(r25859, r25851); float r25861 = r25855 + r25860; float r25862 = sqrt(r25861); float r25863 = r25854 * r25862; return r25863; } double f_id(double J, double K, double U) { double r25864 = -2; double r25865 = J; double r25866 = r25864 * r25865; double r25867 = K; double r25868 = 2; double r25869 = r25867 / r25868; double r25870 = cos(r25869); double r25871 = r25866 * r25870; double r25872 = 1; double r25873 = U; double r25874 = r25868 * r25865; double r25875 = r25874 * r25870; double r25876 = r25873 / r25875; double r25877 = pow(r25876, r25868); double r25878 = r25872 + r25877; double r25879 = sqrt(r25878); double r25880 = r25871 * r25879; return r25880; } double f_of(float J, float K, float U) { float r25881 = J; float r25882 = -2; float r25883 = r25881 * r25882; float r25884 = K; float r25885 = 2; float r25886 = r25884 / r25885; float r25887 = cos(r25886); float r25888 = 1; float r25889 = U; float r25890 = r25889 / r25885; float r25891 = r25890 / r25881; float r25892 = r25891 / r25887; float r25893 = hypot(r25888, r25892); float r25894 = r25887 * r25893; float r25895 = r25883 * r25894; float r25896 = -5.2693998927071626e+306; bool r25897 = r25895 <= r25896; float r25898 = 1/2; float r25899 = r25898 * r25889; float r25900 = r25882 * r25899; float r25901 = 3.668814568454269e+307; bool r25902 = r25895 <= r25901; float r25903 = r25902 ? r25895 : r25900; float r25904 = r25897 ? r25900 : r25903; return r25904; } double f_od(double J, double K, double U) { double r25905 = J; double r25906 = -2; double r25907 = r25905 * r25906; double r25908 = K; double r25909 = 2; double r25910 = r25908 / r25909; double r25911 = cos(r25910); double r25912 = 1; double r25913 = U; double r25914 = r25913 / r25909; double r25915 = r25914 / r25905; double r25916 = r25915 / r25911; double r25917 = hypot(r25912, r25916); double r25918 = r25911 * r25917; double r25919 = r25907 * r25918; double r25920 = -5.2693998927071626e+306; bool r25921 = r25919 <= r25920; double r25922 = 1/2; double r25923 = r25922 * r25913; double r25924 = r25906 * r25923; double r25925 = 3.668814568454269e+307; bool r25926 = r25919 <= r25925; double r25927 = r25926 ? r25919 : r25924; double r25928 = r25921 ? r25924 : r25927; return r25928; } 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 r25929, r25930, r25931, r25932, r25933, r25934, r25935, r25936, r25937, r25938, r25939, r25940, r25941, r25942, r25943, r25944, r25945; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init_set_str(r25929, "-2", 10, MPFR_RNDN); mpfr_init(r25930); mpfr_init(r25931); mpfr_init(r25932); mpfr_init_set_str(r25933, "2", 10, MPFR_RNDN); mpfr_init(r25934); mpfr_init(r25935); mpfr_init(r25936); mpfr_init_set_str(r25937, "1", 10, MPFR_RNDN); mpfr_init(r25938); mpfr_init(r25939); mpfr_init(r25940); mpfr_init(r25941); mpfr_init(r25942); mpfr_init(r25943); mpfr_init(r25944); mpfr_init(r25945); } double f_im(double J, double K, double U) { ; mpfr_set_d(r25930, J, MPFR_RNDN); mpfr_mul(r25931, r25929, r25930, MPFR_RNDN); mpfr_set_d(r25932, K, MPFR_RNDN); ; mpfr_div(r25934, r25932, r25933, MPFR_RNDN); mpfr_cos(r25935, r25934, MPFR_RNDN); mpfr_mul(r25936, r25931, r25935, MPFR_RNDN); ; mpfr_set_d(r25938, U, MPFR_RNDN); mpfr_mul(r25939, r25933, r25930, MPFR_RNDN); mpfr_mul(r25940, r25939, r25935, MPFR_RNDN); mpfr_div(r25941, r25938, r25940, MPFR_RNDN); mpfr_pow(r25942, r25941, r25933, MPFR_RNDN); mpfr_add(r25943, r25937, r25942, MPFR_RNDN); mpfr_sqrt(r25944, r25943, MPFR_RNDN); mpfr_mul(r25945, r25936, r25944, MPFR_RNDN); return mpfr_get_d(r25945, MPFR_RNDN); } static mpfr_t r25946, r25947, r25948, r25949, r25950, r25951, r25952, r25953, r25954, r25955, r25956, r25957, r25958, r25959, r25960, r25961, r25962, r25963, r25964, r25965, r25966, r25967, r25968, r25969; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init(r25946); mpfr_init_set_str(r25947, "-2", 10, MPFR_RNDN); mpfr_init(r25948); mpfr_init(r25949); mpfr_init_set_str(r25950, "2", 10, MPFR_RNDN); mpfr_init(r25951); mpfr_init(r25952); mpfr_init_set_str(r25953, "1", 10, MPFR_RNDN); mpfr_init(r25954); mpfr_init(r25955); mpfr_init(r25956); mpfr_init(r25957); mpfr_init(r25958); mpfr_init(r25959); mpfr_init(r25960); mpfr_init_set_str(r25961, "-5.2693998927071626e+306", 10, MPFR_RNDN); mpfr_init(r25962); mpfr_init_set_str(r25963, "1/2", 10, MPFR_RNDN); mpfr_init(r25964); mpfr_init(r25965); mpfr_init_set_str(r25966, "3.668814568454269e+307", 10, MPFR_RNDN); mpfr_init(r25967); mpfr_init(r25968); mpfr_init(r25969); } double f_fm(double J, double K, double U) { mpfr_set_d(r25946, J, MPFR_RNDN); ; mpfr_mul(r25948, r25946, r25947, MPFR_RNDN); mpfr_set_d(r25949, K, MPFR_RNDN); ; mpfr_div(r25951, r25949, r25950, MPFR_RNDN); mpfr_cos(r25952, r25951, MPFR_RNDN); ; mpfr_set_d(r25954, U, MPFR_RNDN); mpfr_div(r25955, r25954, r25950, MPFR_RNDN); mpfr_div(r25956, r25955, r25946, MPFR_RNDN); mpfr_div(r25957, r25956, r25952, MPFR_RNDN); mpfr_hypot(r25958, r25953, r25957, MPFR_RNDN); mpfr_mul(r25959, r25952, r25958, MPFR_RNDN); mpfr_mul(r25960, r25948, r25959, MPFR_RNDN); ; mpfr_set_si(r25962, mpfr_cmp(r25960, r25961) <= 0, MPFR_RNDN); ; mpfr_mul(r25964, r25963, r25954, MPFR_RNDN); mpfr_mul(r25965, r25947, r25964, MPFR_RNDN); ; mpfr_set_si(r25967, mpfr_cmp(r25960, r25966) <= 0, MPFR_RNDN); if (mpfr_get_si(r25967, MPFR_RNDN)) { mpfr_set(r25968, r25960, MPFR_RNDN); } else { mpfr_set(r25968, r25965, MPFR_RNDN); }; if (mpfr_get_si(r25962, MPFR_RNDN)) { mpfr_set(r25969, r25965, MPFR_RNDN); } else { mpfr_set(r25969, r25968, MPFR_RNDN); }; return mpfr_get_d(r25969, MPFR_RNDN); } static mpfr_t r25970, r25971, r25972, r25973, r25974, r25975, r25976, r25977, r25978, r25979, r25980, r25981, r25982, r25983, r25984, r25985, r25986, r25987, r25988, r25989, r25990, r25991, r25992, r25993; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init(r25970); mpfr_init_set_str(r25971, "-2", 10, MPFR_RNDN); mpfr_init(r25972); mpfr_init(r25973); mpfr_init_set_str(r25974, "2", 10, MPFR_RNDN); mpfr_init(r25975); mpfr_init(r25976); mpfr_init_set_str(r25977, "1", 10, MPFR_RNDN); mpfr_init(r25978); mpfr_init(r25979); mpfr_init(r25980); mpfr_init(r25981); mpfr_init(r25982); mpfr_init(r25983); mpfr_init(r25984); mpfr_init_set_str(r25985, "-5.2693998927071626e+306", 10, MPFR_RNDN); mpfr_init(r25986); mpfr_init_set_str(r25987, "1/2", 10, MPFR_RNDN); mpfr_init(r25988); mpfr_init(r25989); mpfr_init_set_str(r25990, "3.668814568454269e+307", 10, MPFR_RNDN); mpfr_init(r25991); mpfr_init(r25992); mpfr_init(r25993); } double f_dm(double J, double K, double U) { mpfr_set_d(r25970, J, MPFR_RNDN); ; mpfr_mul(r25972, r25970, r25971, MPFR_RNDN); mpfr_set_d(r25973, K, MPFR_RNDN); ; mpfr_div(r25975, r25973, r25974, MPFR_RNDN); mpfr_cos(r25976, r25975, MPFR_RNDN); ; mpfr_set_d(r25978, U, MPFR_RNDN); mpfr_div(r25979, r25978, r25974, MPFR_RNDN); mpfr_div(r25980, r25979, r25970, MPFR_RNDN); mpfr_div(r25981, r25980, r25976, MPFR_RNDN); mpfr_hypot(r25982, r25977, r25981, MPFR_RNDN); mpfr_mul(r25983, r25976, r25982, MPFR_RNDN); mpfr_mul(r25984, r25972, r25983, MPFR_RNDN); ; mpfr_set_si(r25986, mpfr_cmp(r25984, r25985) <= 0, MPFR_RNDN); ; mpfr_mul(r25988, r25987, r25978, MPFR_RNDN); mpfr_mul(r25989, r25971, r25988, MPFR_RNDN); ; mpfr_set_si(r25991, mpfr_cmp(r25984, r25990) <= 0, MPFR_RNDN); if (mpfr_get_si(r25991, MPFR_RNDN)) { mpfr_set(r25992, r25984, MPFR_RNDN); } else { mpfr_set(r25992, r25989, MPFR_RNDN); }; if (mpfr_get_si(r25986, MPFR_RNDN)) { mpfr_set(r25993, r25989, MPFR_RNDN); } else { mpfr_set(r25993, r25992, MPFR_RNDN); }; return mpfr_get_d(r25993, MPFR_RNDN); }