#include #include #include #include #include char *name = "Falkner and Boettcher, Equation (22+)"; double f_if(float v) { float r26873 = 4; float r26874 = 3; float r26875 = atan2(1.0, 0.0); float r26876 = r26874 * r26875; float r26877 = 1; float r26878 = v; float r26879 = r26878 * r26878; float r26880 = r26877 - r26879; float r26881 = r26876 * r26880; float r26882 = 2; float r26883 = 6; float r26884 = r26883 * r26879; float r26885 = r26882 - r26884; float r26886 = sqrt(r26885); float r26887 = r26881 * r26886; float r26888 = r26873 / r26887; return r26888; } double f_id(double v) { double r26889 = 4; double r26890 = 3; double r26891 = atan2(1.0, 0.0); double r26892 = r26890 * r26891; double r26893 = 1; double r26894 = v; double r26895 = r26894 * r26894; double r26896 = r26893 - r26895; double r26897 = r26892 * r26896; double r26898 = 2; double r26899 = 6; double r26900 = r26899 * r26895; double r26901 = r26898 - r26900; double r26902 = sqrt(r26901); double r26903 = r26897 * r26902; double r26904 = r26889 / r26903; return r26904; } double f_of(float v) { float r26905 = 4; float r26906 = 3; float r26907 = atan2(1.0, 0.0); float r26908 = r26906 * r26907; float r26909 = 1; float r26910 = v; float r26911 = r26910 * r26910; float r26912 = r26909 - r26911; float r26913 = r26908 * r26912; float r26914 = r26905 / r26913; float r26915 = 2; float r26916 = 6; float r26917 = r26916 * r26911; float r26918 = r26915 - r26917; float r26919 = sqrt(r26918); float r26920 = r26914 / r26919; return r26920; } double f_od(double v) { double r26921 = 4; double r26922 = 3; double r26923 = atan2(1.0, 0.0); double r26924 = r26922 * r26923; double r26925 = 1; double r26926 = v; double r26927 = r26926 * r26926; double r26928 = r26925 - r26927; double r26929 = r26924 * r26928; double r26930 = r26921 / r26929; double r26931 = 2; double r26932 = 6; double r26933 = r26932 * r26927; double r26934 = r26931 - r26933; double r26935 = sqrt(r26934); double r26936 = r26930 / r26935; return r26936; } 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 r26937, r26938, r26939, r26940, r26941, r26942, r26943, r26944, r26945, r26946, r26947, r26948, r26949, r26950, r26951, r26952; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init_set_str(r26937, "4", 10, MPFR_RNDN); mpfr_init_set_str(r26938, "3", 10, MPFR_RNDN); mpfr_init(r26939); mpfr_init(r26940); mpfr_init_set_str(r26941, "1", 10, MPFR_RNDN); mpfr_init(r26942); mpfr_init(r26943); mpfr_init(r26944); mpfr_init(r26945); mpfr_init_set_str(r26946, "2", 10, MPFR_RNDN); mpfr_init_set_str(r26947, "6", 10, MPFR_RNDN); mpfr_init(r26948); mpfr_init(r26949); mpfr_init(r26950); mpfr_init(r26951); mpfr_init(r26952); } double f_im(double v) { ; ; mpfr_const_pi(r26939, MPFR_RNDN); mpfr_mul(r26940, r26938, r26939, MPFR_RNDN); ; mpfr_set_d(r26942, v, MPFR_RNDN); mpfr_mul(r26943, r26942, r26942, MPFR_RNDN); mpfr_sub(r26944, r26941, r26943, MPFR_RNDN); mpfr_mul(r26945, r26940, r26944, MPFR_RNDN); ; ; mpfr_mul(r26948, r26947, r26943, MPFR_RNDN); mpfr_sub(r26949, r26946, r26948, MPFR_RNDN); mpfr_sqrt(r26950, r26949, MPFR_RNDN); mpfr_mul(r26951, r26945, r26950, MPFR_RNDN); mpfr_div(r26952, r26937, r26951, MPFR_RNDN); return mpfr_get_d(r26952, MPFR_RNDN); } static mpfr_t r26953, r26954, r26955, r26956, r26957, r26958, r26959, r26960, r26961, r26962, r26963, r26964, r26965, r26966, r26967, r26968; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init_set_str(r26953, "4", 10, MPFR_RNDN); mpfr_init_set_str(r26954, "3", 10, MPFR_RNDN); mpfr_init(r26955); mpfr_init(r26956); mpfr_init_set_str(r26957, "1", 10, MPFR_RNDN); mpfr_init(r26958); mpfr_init(r26959); mpfr_init(r26960); mpfr_init(r26961); mpfr_init(r26962); mpfr_init_set_str(r26963, "2", 10, MPFR_RNDN); mpfr_init_set_str(r26964, "6", 10, MPFR_RNDN); mpfr_init(r26965); mpfr_init(r26966); mpfr_init(r26967); mpfr_init(r26968); } double f_fm(double v) { ; ; mpfr_const_pi(r26955, MPFR_RNDN); mpfr_mul(r26956, r26954, r26955, MPFR_RNDN); ; mpfr_set_d(r26958, v, MPFR_RNDN); mpfr_mul(r26959, r26958, r26958, MPFR_RNDN); mpfr_sub(r26960, r26957, r26959, MPFR_RNDN); mpfr_mul(r26961, r26956, r26960, MPFR_RNDN); mpfr_div(r26962, r26953, r26961, MPFR_RNDN); ; ; mpfr_mul(r26965, r26964, r26959, MPFR_RNDN); mpfr_sub(r26966, r26963, r26965, MPFR_RNDN); mpfr_sqrt(r26967, r26966, MPFR_RNDN); mpfr_div(r26968, r26962, r26967, MPFR_RNDN); return mpfr_get_d(r26968, MPFR_RNDN); } static mpfr_t r26969, r26970, r26971, r26972, r26973, r26974, r26975, r26976, r26977, r26978, r26979, r26980, r26981, r26982, r26983, r26984; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init_set_str(r26969, "4", 10, MPFR_RNDN); mpfr_init_set_str(r26970, "3", 10, MPFR_RNDN); mpfr_init(r26971); mpfr_init(r26972); mpfr_init_set_str(r26973, "1", 10, MPFR_RNDN); mpfr_init(r26974); mpfr_init(r26975); mpfr_init(r26976); mpfr_init(r26977); mpfr_init(r26978); mpfr_init_set_str(r26979, "2", 10, MPFR_RNDN); mpfr_init_set_str(r26980, "6", 10, MPFR_RNDN); mpfr_init(r26981); mpfr_init(r26982); mpfr_init(r26983); mpfr_init(r26984); } double f_dm(double v) { ; ; mpfr_const_pi(r26971, MPFR_RNDN); mpfr_mul(r26972, r26970, r26971, MPFR_RNDN); ; mpfr_set_d(r26974, v, MPFR_RNDN); mpfr_mul(r26975, r26974, r26974, MPFR_RNDN); mpfr_sub(r26976, r26973, r26975, MPFR_RNDN); mpfr_mul(r26977, r26972, r26976, MPFR_RNDN); mpfr_div(r26978, r26969, r26977, MPFR_RNDN); ; ; mpfr_mul(r26981, r26980, r26975, MPFR_RNDN); mpfr_sub(r26982, r26979, r26981, MPFR_RNDN); mpfr_sqrt(r26983, r26982, MPFR_RNDN); mpfr_div(r26984, r26978, r26983, MPFR_RNDN); return mpfr_get_d(r26984, MPFR_RNDN); }