#include #include #include #include #include char *name = "VandenBroeck and Keller, Equation (6)"; double f_if(float F, float l) { float r23894 = atan2(1.0, 0.0); float r23895 = l; float r23896 = r23894 * r23895; float r23897 = 1; float r23898 = F; float r23899 = r23898 * r23898; float r23900 = r23897 / r23899; float r23901 = tan(r23896); float r23902 = r23900 * r23901; float r23903 = r23896 - r23902; return r23903; } double f_id(double F, double l) { double r23904 = atan2(1.0, 0.0); double r23905 = l; double r23906 = r23904 * r23905; double r23907 = 1; double r23908 = F; double r23909 = r23908 * r23908; double r23910 = r23907 / r23909; double r23911 = tan(r23906); double r23912 = r23910 * r23911; double r23913 = r23906 - r23912; return r23913; } double f_of(float F, float l) { float r23914 = l; float r23915 = atan2(1.0, 0.0); float r23916 = r23914 * r23915; float r23917 = tan(r23916); float r23918 = cbrt(r23917); float r23919 = -r23918; float r23920 = F; float r23921 = r23919 / r23920; float r23922 = sin(r23916); float r23923 = cos(r23916); float r23924 = r23922 / r23923; float r23925 = cbrt(r23924); float r23926 = r23918 * r23918; float r23927 = r23926 * r23918; float r23928 = cbrt(r23927); float r23929 = r23920 / r23928; float r23930 = r23925 / r23929; float r23931 = fma(r23921, r23930, r23916); return r23931; } double f_od(double F, double l) { double r23932 = l; double r23933 = atan2(1.0, 0.0); double r23934 = r23932 * r23933; double r23935 = tan(r23934); double r23936 = cbrt(r23935); double r23937 = -r23936; double r23938 = F; double r23939 = r23937 / r23938; double r23940 = sin(r23934); double r23941 = cos(r23934); double r23942 = r23940 / r23941; double r23943 = cbrt(r23942); double r23944 = r23936 * r23936; double r23945 = r23944 * r23936; double r23946 = cbrt(r23945); double r23947 = r23938 / r23946; double r23948 = r23943 / r23947; double r23949 = fma(r23939, r23948, r23934); return r23949; } 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 r23950, r23951, r23952, r23953, r23954, r23955, r23956, r23957, r23958, r23959; void setup_mpfr_f_im() { mpfr_set_default_prec(3920); mpfr_init(r23950); mpfr_init(r23951); mpfr_init(r23952); mpfr_init_set_str(r23953, "1", 10, MPFR_RNDN); mpfr_init(r23954); mpfr_init(r23955); mpfr_init(r23956); mpfr_init(r23957); mpfr_init(r23958); mpfr_init(r23959); } double f_im(double F, double l) { mpfr_const_pi(r23950, MPFR_RNDN); mpfr_set_d(r23951, l, MPFR_RNDN); mpfr_mul(r23952, r23950, r23951, MPFR_RNDN); ; mpfr_set_d(r23954, F, MPFR_RNDN); mpfr_mul(r23955, r23954, r23954, MPFR_RNDN); mpfr_div(r23956, r23953, r23955, MPFR_RNDN); mpfr_tan(r23957, r23952, MPFR_RNDN); mpfr_mul(r23958, r23956, r23957, MPFR_RNDN); mpfr_sub(r23959, r23952, r23958, MPFR_RNDN); return mpfr_get_d(r23959, MPFR_RNDN); } static mpfr_t r23960, r23961, r23962, r23963, r23964, r23965, r23966, r23967, r23968, r23969, r23970, r23971, r23972, r23973, r23974, r23975, r23976, r23977; void setup_mpfr_f_fm() { mpfr_set_default_prec(3920); mpfr_init(r23960); mpfr_init(r23961); mpfr_init(r23962); mpfr_init(r23963); mpfr_init(r23964); mpfr_init(r23965); mpfr_init(r23966); mpfr_init(r23967); mpfr_init(r23968); mpfr_init(r23969); mpfr_init(r23970); mpfr_init(r23971); mpfr_init(r23972); mpfr_init(r23973); mpfr_init(r23974); mpfr_init(r23975); mpfr_init(r23976); mpfr_init(r23977); } double f_fm(double F, double l) { mpfr_set_d(r23960, l, MPFR_RNDN); mpfr_const_pi(r23961, MPFR_RNDN); mpfr_mul(r23962, r23960, r23961, MPFR_RNDN); mpfr_tan(r23963, r23962, MPFR_RNDN); mpfr_cbrt(r23964, r23963, MPFR_RNDN); mpfr_neg(r23965, r23964, MPFR_RNDN); mpfr_set_d(r23966, F, MPFR_RNDN); mpfr_div(r23967, r23965, r23966, MPFR_RNDN); mpfr_sin(r23968, r23962, MPFR_RNDN); mpfr_cos(r23969, r23962, MPFR_RNDN); mpfr_div(r23970, r23968, r23969, MPFR_RNDN); mpfr_cbrt(r23971, r23970, MPFR_RNDN); mpfr_mul(r23972, r23964, r23964, MPFR_RNDN); mpfr_mul(r23973, r23972, r23964, MPFR_RNDN); mpfr_cbrt(r23974, r23973, MPFR_RNDN); mpfr_div(r23975, r23966, r23974, MPFR_RNDN); mpfr_div(r23976, r23971, r23975, MPFR_RNDN); mpfr_fma(r23977, r23967, r23976, r23962, MPFR_RNDN); return mpfr_get_d(r23977, MPFR_RNDN); } static mpfr_t r23978, r23979, r23980, r23981, r23982, r23983, r23984, r23985, r23986, r23987, r23988, r23989, r23990, r23991, r23992, r23993, r23994, r23995; void setup_mpfr_f_dm() { mpfr_set_default_prec(3920); mpfr_init(r23978); mpfr_init(r23979); mpfr_init(r23980); mpfr_init(r23981); mpfr_init(r23982); mpfr_init(r23983); mpfr_init(r23984); mpfr_init(r23985); mpfr_init(r23986); mpfr_init(r23987); mpfr_init(r23988); mpfr_init(r23989); mpfr_init(r23990); mpfr_init(r23991); mpfr_init(r23992); mpfr_init(r23993); mpfr_init(r23994); mpfr_init(r23995); } double f_dm(double F, double l) { mpfr_set_d(r23978, l, MPFR_RNDN); mpfr_const_pi(r23979, MPFR_RNDN); mpfr_mul(r23980, r23978, r23979, MPFR_RNDN); mpfr_tan(r23981, r23980, MPFR_RNDN); mpfr_cbrt(r23982, r23981, MPFR_RNDN); mpfr_neg(r23983, r23982, MPFR_RNDN); mpfr_set_d(r23984, F, MPFR_RNDN); mpfr_div(r23985, r23983, r23984, MPFR_RNDN); mpfr_sin(r23986, r23980, MPFR_RNDN); mpfr_cos(r23987, r23980, MPFR_RNDN); mpfr_div(r23988, r23986, r23987, MPFR_RNDN); mpfr_cbrt(r23989, r23988, MPFR_RNDN); mpfr_mul(r23990, r23982, r23982, MPFR_RNDN); mpfr_mul(r23991, r23990, r23982, MPFR_RNDN); mpfr_cbrt(r23992, r23991, MPFR_RNDN); mpfr_div(r23993, r23984, r23992, MPFR_RNDN); mpfr_div(r23994, r23989, r23993, MPFR_RNDN); mpfr_fma(r23995, r23985, r23994, r23980, MPFR_RNDN); return mpfr_get_d(r23995, MPFR_RNDN); }