#include #include #include #include #include char *name = "Toniolo and Linder, Equation (3a)"; double f_if(float l, float Om, float kx, float ky) { float r21803 = 1; float r21804 = 2; float r21805 = r21803 / r21804; float r21806 = l; float r21807 = r21804 * r21806; float r21808 = Om; float r21809 = r21807 / r21808; float r21810 = pow(r21809, r21804); float r21811 = kx; float r21812 = sin(r21811); float r21813 = pow(r21812, r21804); float r21814 = ky; float r21815 = sin(r21814); float r21816 = pow(r21815, r21804); float r21817 = r21813 + r21816; float r21818 = r21810 * r21817; float r21819 = r21803 + r21818; float r21820 = sqrt(r21819); float r21821 = r21803 / r21820; float r21822 = r21803 + r21821; float r21823 = r21805 * r21822; float r21824 = sqrt(r21823); return r21824; } double f_id(double l, double Om, double kx, double ky) { double r21825 = 1; double r21826 = 2; double r21827 = r21825 / r21826; double r21828 = l; double r21829 = r21826 * r21828; double r21830 = Om; double r21831 = r21829 / r21830; double r21832 = pow(r21831, r21826); double r21833 = kx; double r21834 = sin(r21833); double r21835 = pow(r21834, r21826); double r21836 = ky; double r21837 = sin(r21836); double r21838 = pow(r21837, r21826); double r21839 = r21835 + r21838; double r21840 = r21832 * r21839; double r21841 = r21825 + r21840; double r21842 = sqrt(r21841); double r21843 = r21825 / r21842; double r21844 = r21825 + r21843; double r21845 = r21827 * r21844; double r21846 = sqrt(r21845); return r21846; } double f_of(float l, float Om, float kx, float ky) { float r21847 = 1; float r21848 = 2; float r21849 = r21847 / r21848; float r21850 = l; float r21851 = r21848 * r21850; float r21852 = Om; float r21853 = r21851 / r21852; float r21854 = kx; float r21855 = sin(r21854); float r21856 = pow(r21855, r21848); float r21857 = ky; float r21858 = sin(r21857); float r21859 = pow(r21858, r21848); float r21860 = r21856 + r21859; float r21861 = r21853 * r21860; float r21862 = r21853 * r21861; float r21863 = r21847 + r21862; float r21864 = sqrt(r21863); float r21865 = r21847 / r21864; float r21866 = r21847 + r21865; float r21867 = r21849 * r21866; float r21868 = sqrt(r21867); return r21868; } double f_od(double l, double Om, double kx, double ky) { double r21869 = 1; double r21870 = 2; double r21871 = r21869 / r21870; double r21872 = l; double r21873 = r21870 * r21872; double r21874 = Om; double r21875 = r21873 / r21874; double r21876 = kx; double r21877 = sin(r21876); double r21878 = pow(r21877, r21870); double r21879 = ky; double r21880 = sin(r21879); double r21881 = pow(r21880, r21870); double r21882 = r21878 + r21881; double r21883 = r21875 * r21882; double r21884 = r21875 * r21883; double r21885 = r21869 + r21884; double r21886 = sqrt(r21885); double r21887 = r21869 / r21886; double r21888 = r21869 + r21887; double r21889 = r21871 * r21888; double r21890 = sqrt(r21889); return r21890; } 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 r21891, r21892, r21893, r21894, r21895, r21896, r21897, r21898, r21899, r21900, r21901, r21902, r21903, r21904, r21905, r21906, r21907, r21908, r21909, r21910, r21911, r21912; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r21891, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21892, "2", 10, MPFR_RNDN); mpfr_init(r21893); mpfr_init(r21894); mpfr_init(r21895); mpfr_init(r21896); mpfr_init(r21897); mpfr_init(r21898); mpfr_init(r21899); mpfr_init(r21900); mpfr_init(r21901); mpfr_init(r21902); mpfr_init(r21903); mpfr_init(r21904); mpfr_init(r21905); mpfr_init(r21906); mpfr_init(r21907); mpfr_init(r21908); mpfr_init(r21909); mpfr_init(r21910); mpfr_init(r21911); mpfr_init(r21912); } double f_im(double l, double Om, double kx, double ky) { ; ; mpfr_div(r21893, r21891, r21892, MPFR_RNDN); mpfr_set_d(r21894, l, MPFR_RNDN); mpfr_mul(r21895, r21892, r21894, MPFR_RNDN); mpfr_set_d(r21896, Om, MPFR_RNDN); mpfr_div(r21897, r21895, r21896, MPFR_RNDN); mpfr_pow(r21898, r21897, r21892, MPFR_RNDN); mpfr_set_d(r21899, kx, MPFR_RNDN); mpfr_sin(r21900, r21899, MPFR_RNDN); mpfr_pow(r21901, r21900, r21892, MPFR_RNDN); mpfr_set_d(r21902, ky, MPFR_RNDN); mpfr_sin(r21903, r21902, MPFR_RNDN); mpfr_pow(r21904, r21903, r21892, MPFR_RNDN); mpfr_add(r21905, r21901, r21904, MPFR_RNDN); mpfr_mul(r21906, r21898, r21905, MPFR_RNDN); mpfr_add(r21907, r21891, r21906, MPFR_RNDN); mpfr_sqrt(r21908, r21907, MPFR_RNDN); mpfr_div(r21909, r21891, r21908, MPFR_RNDN); mpfr_add(r21910, r21891, r21909, MPFR_RNDN); mpfr_mul(r21911, r21893, r21910, MPFR_RNDN); mpfr_sqrt(r21912, r21911, MPFR_RNDN); return mpfr_get_d(r21912, MPFR_RNDN); } static mpfr_t r21913, r21914, r21915, r21916, r21917, r21918, r21919, r21920, r21921, r21922, r21923, r21924, r21925, r21926, r21927, r21928, r21929, r21930, r21931, r21932, r21933, r21934; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r21913, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21914, "2", 10, MPFR_RNDN); mpfr_init(r21915); mpfr_init(r21916); mpfr_init(r21917); mpfr_init(r21918); mpfr_init(r21919); mpfr_init(r21920); mpfr_init(r21921); mpfr_init(r21922); mpfr_init(r21923); mpfr_init(r21924); mpfr_init(r21925); mpfr_init(r21926); mpfr_init(r21927); mpfr_init(r21928); mpfr_init(r21929); mpfr_init(r21930); mpfr_init(r21931); mpfr_init(r21932); mpfr_init(r21933); mpfr_init(r21934); } double f_fm(double l, double Om, double kx, double ky) { ; ; mpfr_div(r21915, r21913, r21914, MPFR_RNDN); mpfr_set_d(r21916, l, MPFR_RNDN); mpfr_mul(r21917, r21914, r21916, MPFR_RNDN); mpfr_set_d(r21918, Om, MPFR_RNDN); mpfr_div(r21919, r21917, r21918, MPFR_RNDN); mpfr_set_d(r21920, kx, MPFR_RNDN); mpfr_sin(r21921, r21920, MPFR_RNDN); mpfr_pow(r21922, r21921, r21914, MPFR_RNDN); mpfr_set_d(r21923, ky, MPFR_RNDN); mpfr_sin(r21924, r21923, MPFR_RNDN); mpfr_pow(r21925, r21924, r21914, MPFR_RNDN); mpfr_add(r21926, r21922, r21925, MPFR_RNDN); mpfr_mul(r21927, r21919, r21926, MPFR_RNDN); mpfr_mul(r21928, r21919, r21927, MPFR_RNDN); mpfr_add(r21929, r21913, r21928, MPFR_RNDN); mpfr_sqrt(r21930, r21929, MPFR_RNDN); mpfr_div(r21931, r21913, r21930, MPFR_RNDN); mpfr_add(r21932, r21913, r21931, MPFR_RNDN); mpfr_mul(r21933, r21915, r21932, MPFR_RNDN); mpfr_sqrt(r21934, r21933, MPFR_RNDN); return mpfr_get_d(r21934, MPFR_RNDN); } static mpfr_t r21935, r21936, r21937, r21938, r21939, r21940, r21941, r21942, r21943, r21944, r21945, r21946, r21947, r21948, r21949, r21950, r21951, r21952, r21953, r21954, r21955, r21956; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r21935, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21936, "2", 10, MPFR_RNDN); mpfr_init(r21937); mpfr_init(r21938); mpfr_init(r21939); mpfr_init(r21940); mpfr_init(r21941); mpfr_init(r21942); mpfr_init(r21943); mpfr_init(r21944); mpfr_init(r21945); mpfr_init(r21946); mpfr_init(r21947); mpfr_init(r21948); mpfr_init(r21949); mpfr_init(r21950); mpfr_init(r21951); mpfr_init(r21952); mpfr_init(r21953); mpfr_init(r21954); mpfr_init(r21955); mpfr_init(r21956); } double f_dm(double l, double Om, double kx, double ky) { ; ; mpfr_div(r21937, r21935, r21936, MPFR_RNDN); mpfr_set_d(r21938, l, MPFR_RNDN); mpfr_mul(r21939, r21936, r21938, MPFR_RNDN); mpfr_set_d(r21940, Om, MPFR_RNDN); mpfr_div(r21941, r21939, r21940, MPFR_RNDN); mpfr_set_d(r21942, kx, MPFR_RNDN); mpfr_sin(r21943, r21942, MPFR_RNDN); mpfr_pow(r21944, r21943, r21936, MPFR_RNDN); mpfr_set_d(r21945, ky, MPFR_RNDN); mpfr_sin(r21946, r21945, MPFR_RNDN); mpfr_pow(r21947, r21946, r21936, MPFR_RNDN); mpfr_add(r21948, r21944, r21947, MPFR_RNDN); mpfr_mul(r21949, r21941, r21948, MPFR_RNDN); mpfr_mul(r21950, r21941, r21949, MPFR_RNDN); mpfr_add(r21951, r21935, r21950, MPFR_RNDN); mpfr_sqrt(r21952, r21951, MPFR_RNDN); mpfr_div(r21953, r21935, r21952, MPFR_RNDN); mpfr_add(r21954, r21935, r21953, MPFR_RNDN); mpfr_mul(r21955, r21937, r21954, MPFR_RNDN); mpfr_sqrt(r21956, r21955, MPFR_RNDN); return mpfr_get_d(r21956, MPFR_RNDN); }