#include #include #include #include #include char *name = "Toniolo and Linder, Equation (3a)"; double f_if(float l, float Om, float kx, float ky) { float r20828 = 1.0f; float r20829 = 2.0f; float r20830 = r20828 / r20829; float r20831 = l; float r20832 = r20829 * r20831; float r20833 = Om; float r20834 = r20832 / r20833; float r20835 = r20834 * r20834; float r20836 = kx; float r20837 = sin(r20836); float r20838 = r20837 * r20837; float r20839 = ky; float r20840 = sin(r20839); float r20841 = r20840 * r20840; float r20842 = r20838 + r20841; float r20843 = r20835 * r20842; float r20844 = r20828 + r20843; float r20845 = sqrt(r20844); float r20846 = 1.0/r20845; float r20847 = r20828 + r20846; float r20848 = r20830 * r20847; float r20849 = sqrt(r20848); return r20849; } double f_id(double l, double Om, double kx, double ky) { double r20850 = 1.0; double r20851 = 2.0; double r20852 = r20850 / r20851; double r20853 = l; double r20854 = r20851 * r20853; double r20855 = Om; double r20856 = r20854 / r20855; double r20857 = r20856 * r20856; double r20858 = kx; double r20859 = sin(r20858); double r20860 = r20859 * r20859; double r20861 = ky; double r20862 = sin(r20861); double r20863 = r20862 * r20862; double r20864 = r20860 + r20863; double r20865 = r20857 * r20864; double r20866 = r20850 + r20865; double r20867 = sqrt(r20866); double r20868 = 1.0/r20867; double r20869 = r20850 + r20868; double r20870 = r20852 * r20869; double r20871 = sqrt(r20870); return r20871; } double f_of(float l, float Om, float kx, float ky) { float r20872 = 1.0f; float r20873 = kx; float r20874 = sin(r20873); float r20875 = r20874 * r20874; float r20876 = ky; float r20877 = sin(r20876); float r20878 = r20877 * r20877; float r20879 = r20875 + r20878; float r20880 = l; float r20881 = r20880 + r20880; float r20882 = Om; float r20883 = r20881 / r20882; float r20884 = r20883 * r20883; float r20885 = r20879 * r20884; float r20886 = r20885 + r20872; float r20887 = sqrt(r20886); float r20888 = 1.0/r20887; float r20889 = r20872 + r20888; float r20890 = 2.0f; float r20891 = r20889 / r20890; float r20892 = sqrt(r20891); return r20892; } double f_od(double l, double Om, double kx, double ky) { double r20893 = 1.0; double r20894 = kx; double r20895 = sin(r20894); double r20896 = r20895 * r20895; double r20897 = ky; double r20898 = sin(r20897); double r20899 = r20898 * r20898; double r20900 = r20896 + r20899; double r20901 = l; double r20902 = r20901 + r20901; double r20903 = Om; double r20904 = r20902 / r20903; double r20905 = r20904 * r20904; double r20906 = r20900 * r20905; double r20907 = r20906 + r20893; double r20908 = sqrt(r20907); double r20909 = 1.0/r20908; double r20910 = r20893 + r20909; double r20911 = 2.0; double r20912 = r20910 / r20911; double r20913 = sqrt(r20912); return r20913; } 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 r20914, r20915, r20916, r20917, r20918, r20919, r20920, r20921, r20922, r20923, r20924, r20925, r20926, r20927, r20928, r20929, r20930, r20931, r20932, r20933, r20934, r20935; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r20914, "1", 10, MPFR_RNDN); mpfr_init_set_str(r20915, "2", 10, MPFR_RNDN); mpfr_init(r20916); mpfr_init(r20917); mpfr_init(r20918); mpfr_init(r20919); mpfr_init(r20920); mpfr_init(r20921); mpfr_init(r20922); mpfr_init(r20923); mpfr_init(r20924); mpfr_init(r20925); mpfr_init(r20926); mpfr_init(r20927); mpfr_init(r20928); mpfr_init(r20929); mpfr_init(r20930); mpfr_init(r20931); mpfr_init(r20932); mpfr_init(r20933); mpfr_init(r20934); mpfr_init(r20935); } double f_im(double l, double Om, double kx, double ky) { ; ; mpfr_div(r20916, r20914, r20915, MPFR_RNDN); mpfr_set_d(r20917, l, MPFR_RNDN); mpfr_mul(r20918, r20915, r20917, MPFR_RNDN); mpfr_set_d(r20919, Om, MPFR_RNDN); mpfr_div(r20920, r20918, r20919, MPFR_RNDN); mpfr_sqr(r20921, r20920, MPFR_RNDN); mpfr_set_d(r20922, kx, MPFR_RNDN); mpfr_sin(r20923, r20922, MPFR_RNDN); mpfr_sqr(r20924, r20923, MPFR_RNDN); mpfr_set_d(r20925, ky, MPFR_RNDN); mpfr_sin(r20926, r20925, MPFR_RNDN); mpfr_sqr(r20927, r20926, MPFR_RNDN); mpfr_add(r20928, r20924, r20927, MPFR_RNDN); mpfr_mul(r20929, r20921, r20928, MPFR_RNDN); mpfr_add(r20930, r20914, r20929, MPFR_RNDN); mpfr_sqrt(r20931, r20930, MPFR_RNDN); mpfr_ui_div(r20932, 1, r20931, MPFR_RNDN); mpfr_add(r20933, r20914, r20932, MPFR_RNDN); mpfr_mul(r20934, r20916, r20933, MPFR_RNDN); mpfr_sqrt(r20935, r20934, MPFR_RNDN); return mpfr_get_d(r20935, MPFR_RNDN); } static mpfr_t r20936, r20937, r20938, r20939, r20940, r20941, r20942, r20943, r20944, r20945, r20946, r20947, r20948, r20949, r20950, r20951, r20952, r20953, r20954, r20955, r20956; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r20936, "1", 10, MPFR_RNDN); mpfr_init(r20937); mpfr_init(r20938); mpfr_init(r20939); mpfr_init(r20940); mpfr_init(r20941); mpfr_init(r20942); mpfr_init(r20943); mpfr_init(r20944); mpfr_init(r20945); mpfr_init(r20946); mpfr_init(r20947); mpfr_init(r20948); mpfr_init(r20949); mpfr_init(r20950); mpfr_init(r20951); mpfr_init(r20952); mpfr_init(r20953); mpfr_init_set_str(r20954, "2", 10, MPFR_RNDN); mpfr_init(r20955); mpfr_init(r20956); } double f_fm(double l, double Om, double kx, double ky) { ; mpfr_set_d(r20937, kx, MPFR_RNDN); mpfr_sin(r20938, r20937, MPFR_RNDN); mpfr_sqr(r20939, r20938, MPFR_RNDN); mpfr_set_d(r20940, ky, MPFR_RNDN); mpfr_sin(r20941, r20940, MPFR_RNDN); mpfr_sqr(r20942, r20941, MPFR_RNDN); mpfr_add(r20943, r20939, r20942, MPFR_RNDN); mpfr_set_d(r20944, l, MPFR_RNDN); mpfr_add(r20945, r20944, r20944, MPFR_RNDN); mpfr_set_d(r20946, Om, MPFR_RNDN); mpfr_div(r20947, r20945, r20946, MPFR_RNDN); mpfr_sqr(r20948, r20947, MPFR_RNDN); mpfr_mul(r20949, r20943, r20948, MPFR_RNDN); mpfr_add(r20950, r20949, r20936, MPFR_RNDN); mpfr_sqrt(r20951, r20950, MPFR_RNDN); mpfr_ui_div(r20952, 1, r20951, MPFR_RNDN); mpfr_add(r20953, r20936, r20952, MPFR_RNDN); ; mpfr_div(r20955, r20953, r20954, MPFR_RNDN); mpfr_sqrt(r20956, r20955, MPFR_RNDN); return mpfr_get_d(r20956, MPFR_RNDN); } static mpfr_t r20957, r20958, r20959, r20960, r20961, r20962, r20963, r20964, r20965, r20966, r20967, r20968, r20969, r20970, r20971, r20972, r20973, r20974, r20975, r20976, r20977; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r20957, "1", 10, MPFR_RNDN); mpfr_init(r20958); mpfr_init(r20959); mpfr_init(r20960); mpfr_init(r20961); mpfr_init(r20962); mpfr_init(r20963); mpfr_init(r20964); mpfr_init(r20965); mpfr_init(r20966); mpfr_init(r20967); mpfr_init(r20968); mpfr_init(r20969); mpfr_init(r20970); mpfr_init(r20971); mpfr_init(r20972); mpfr_init(r20973); mpfr_init(r20974); mpfr_init_set_str(r20975, "2", 10, MPFR_RNDN); mpfr_init(r20976); mpfr_init(r20977); } double f_dm(double l, double Om, double kx, double ky) { ; mpfr_set_d(r20958, kx, MPFR_RNDN); mpfr_sin(r20959, r20958, MPFR_RNDN); mpfr_sqr(r20960, r20959, MPFR_RNDN); mpfr_set_d(r20961, ky, MPFR_RNDN); mpfr_sin(r20962, r20961, MPFR_RNDN); mpfr_sqr(r20963, r20962, MPFR_RNDN); mpfr_add(r20964, r20960, r20963, MPFR_RNDN); mpfr_set_d(r20965, l, MPFR_RNDN); mpfr_add(r20966, r20965, r20965, MPFR_RNDN); mpfr_set_d(r20967, Om, MPFR_RNDN); mpfr_div(r20968, r20966, r20967, MPFR_RNDN); mpfr_sqr(r20969, r20968, MPFR_RNDN); mpfr_mul(r20970, r20964, r20969, MPFR_RNDN); mpfr_add(r20971, r20970, r20957, MPFR_RNDN); mpfr_sqrt(r20972, r20971, MPFR_RNDN); mpfr_ui_div(r20973, 1, r20972, MPFR_RNDN); mpfr_add(r20974, r20957, r20973, MPFR_RNDN); ; mpfr_div(r20976, r20974, r20975, MPFR_RNDN); mpfr_sqrt(r20977, r20976, MPFR_RNDN); return mpfr_get_d(r20977, MPFR_RNDN); }