#include #include #include #include #include char *name = "Toniolo and Linder, Equation (3a)"; double f_if(float l, float Om, float kx, float ky) { float r21750 = 1; float r21751 = 2; float r21752 = r21750 / r21751; float r21753 = l; float r21754 = r21751 * r21753; float r21755 = Om; float r21756 = r21754 / r21755; float r21757 = pow(r21756, r21751); float r21758 = kx; float r21759 = sin(r21758); float r21760 = pow(r21759, r21751); float r21761 = ky; float r21762 = sin(r21761); float r21763 = pow(r21762, r21751); float r21764 = r21760 + r21763; float r21765 = r21757 * r21764; float r21766 = r21750 + r21765; float r21767 = sqrt(r21766); float r21768 = r21750 / r21767; float r21769 = r21750 + r21768; float r21770 = r21752 * r21769; float r21771 = sqrt(r21770); return r21771; } double f_id(double l, double Om, double kx, double ky) { double r21772 = 1; double r21773 = 2; double r21774 = r21772 / r21773; double r21775 = l; double r21776 = r21773 * r21775; double r21777 = Om; double r21778 = r21776 / r21777; double r21779 = pow(r21778, r21773); double r21780 = kx; double r21781 = sin(r21780); double r21782 = pow(r21781, r21773); double r21783 = ky; double r21784 = sin(r21783); double r21785 = pow(r21784, r21773); double r21786 = r21782 + r21785; double r21787 = r21779 * r21786; double r21788 = r21772 + r21787; double r21789 = sqrt(r21788); double r21790 = r21772 / r21789; double r21791 = r21772 + r21790; double r21792 = r21774 * r21791; double r21793 = sqrt(r21792); return r21793; } double f_of(float l, float Om, float kx, float ky) { float r21794 = 1; float r21795 = 2; float r21796 = r21794 / r21795; float r21797 = l; float r21798 = r21795 * r21797; float r21799 = Om; float r21800 = r21798 / r21799; float r21801 = pow(r21800, r21795); float r21802 = kx; float r21803 = sin(r21802); float r21804 = pow(r21803, r21795); float r21805 = ky; float r21806 = sin(r21805); float r21807 = pow(r21806, r21795); float r21808 = r21804 + r21807; float r21809 = r21801 * r21808; float r21810 = r21794 + r21809; float r21811 = sqrt(r21810); float r21812 = r21794 / r21811; float r21813 = exp(r21812); float r21814 = log(r21813); float r21815 = r21794 + r21814; float r21816 = r21796 * r21815; float r21817 = sqrt(r21816); return r21817; } double f_od(double l, double Om, double kx, double ky) { double r21818 = 1; double r21819 = 2; double r21820 = r21818 / r21819; double r21821 = l; double r21822 = r21819 * r21821; double r21823 = Om; double r21824 = r21822 / r21823; double r21825 = pow(r21824, r21819); double r21826 = kx; double r21827 = sin(r21826); double r21828 = pow(r21827, r21819); double r21829 = ky; double r21830 = sin(r21829); double r21831 = pow(r21830, r21819); double r21832 = r21828 + r21831; double r21833 = r21825 * r21832; double r21834 = r21818 + r21833; double r21835 = sqrt(r21834); double r21836 = r21818 / r21835; double r21837 = exp(r21836); double r21838 = log(r21837); double r21839 = r21818 + r21838; double r21840 = r21820 * r21839; double r21841 = sqrt(r21840); return r21841; } 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 r21842, r21843, r21844, r21845, r21846, r21847, r21848, r21849, r21850, r21851, r21852, r21853, r21854, r21855, r21856, r21857, r21858, r21859, r21860, r21861, r21862, r21863; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r21842, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21843, "2", 10, MPFR_RNDN); mpfr_init(r21844); mpfr_init(r21845); mpfr_init(r21846); mpfr_init(r21847); mpfr_init(r21848); mpfr_init(r21849); mpfr_init(r21850); mpfr_init(r21851); mpfr_init(r21852); mpfr_init(r21853); mpfr_init(r21854); mpfr_init(r21855); mpfr_init(r21856); mpfr_init(r21857); mpfr_init(r21858); mpfr_init(r21859); mpfr_init(r21860); mpfr_init(r21861); mpfr_init(r21862); mpfr_init(r21863); } double f_im(double l, double Om, double kx, double ky) { ; ; mpfr_div(r21844, r21842, r21843, MPFR_RNDN); mpfr_set_d(r21845, l, MPFR_RNDN); mpfr_mul(r21846, r21843, r21845, MPFR_RNDN); mpfr_set_d(r21847, Om, MPFR_RNDN); mpfr_div(r21848, r21846, r21847, MPFR_RNDN); mpfr_pow(r21849, r21848, r21843, MPFR_RNDN); mpfr_set_d(r21850, kx, MPFR_RNDN); mpfr_sin(r21851, r21850, MPFR_RNDN); mpfr_pow(r21852, r21851, r21843, MPFR_RNDN); mpfr_set_d(r21853, ky, MPFR_RNDN); mpfr_sin(r21854, r21853, MPFR_RNDN); mpfr_pow(r21855, r21854, r21843, MPFR_RNDN); mpfr_add(r21856, r21852, r21855, MPFR_RNDN); mpfr_mul(r21857, r21849, r21856, MPFR_RNDN); mpfr_add(r21858, r21842, r21857, MPFR_RNDN); mpfr_sqrt(r21859, r21858, MPFR_RNDN); mpfr_div(r21860, r21842, r21859, MPFR_RNDN); mpfr_add(r21861, r21842, r21860, MPFR_RNDN); mpfr_mul(r21862, r21844, r21861, MPFR_RNDN); mpfr_sqrt(r21863, r21862, MPFR_RNDN); return mpfr_get_d(r21863, MPFR_RNDN); } static mpfr_t r21864, r21865, r21866, r21867, r21868, r21869, r21870, r21871, r21872, r21873, r21874, r21875, r21876, r21877, r21878, r21879, r21880, r21881, r21882, r21883, r21884, r21885, r21886, r21887; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r21864, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21865, "2", 10, MPFR_RNDN); mpfr_init(r21866); mpfr_init(r21867); mpfr_init(r21868); mpfr_init(r21869); mpfr_init(r21870); mpfr_init(r21871); mpfr_init(r21872); mpfr_init(r21873); mpfr_init(r21874); mpfr_init(r21875); mpfr_init(r21876); mpfr_init(r21877); mpfr_init(r21878); mpfr_init(r21879); mpfr_init(r21880); mpfr_init(r21881); mpfr_init(r21882); mpfr_init(r21883); mpfr_init(r21884); mpfr_init(r21885); mpfr_init(r21886); mpfr_init(r21887); } double f_fm(double l, double Om, double kx, double ky) { ; ; mpfr_div(r21866, r21864, r21865, MPFR_RNDN); mpfr_set_d(r21867, l, MPFR_RNDN); mpfr_mul(r21868, r21865, r21867, MPFR_RNDN); mpfr_set_d(r21869, Om, MPFR_RNDN); mpfr_div(r21870, r21868, r21869, MPFR_RNDN); mpfr_pow(r21871, r21870, r21865, MPFR_RNDN); mpfr_set_d(r21872, kx, MPFR_RNDN); mpfr_sin(r21873, r21872, MPFR_RNDN); mpfr_pow(r21874, r21873, r21865, MPFR_RNDN); mpfr_set_d(r21875, ky, MPFR_RNDN); mpfr_sin(r21876, r21875, MPFR_RNDN); mpfr_pow(r21877, r21876, r21865, MPFR_RNDN); mpfr_add(r21878, r21874, r21877, MPFR_RNDN); mpfr_mul(r21879, r21871, r21878, MPFR_RNDN); mpfr_add(r21880, r21864, r21879, MPFR_RNDN); mpfr_sqrt(r21881, r21880, MPFR_RNDN); mpfr_div(r21882, r21864, r21881, MPFR_RNDN); mpfr_exp(r21883, r21882, MPFR_RNDN); mpfr_log(r21884, r21883, MPFR_RNDN); mpfr_add(r21885, r21864, r21884, MPFR_RNDN); mpfr_mul(r21886, r21866, r21885, MPFR_RNDN); mpfr_sqrt(r21887, r21886, MPFR_RNDN); return mpfr_get_d(r21887, MPFR_RNDN); } static mpfr_t r21888, r21889, r21890, r21891, r21892, r21893, r21894, r21895, r21896, r21897, r21898, r21899, r21900, r21901, r21902, r21903, r21904, r21905, r21906, r21907, r21908, r21909, r21910, r21911; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r21888, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21889, "2", 10, MPFR_RNDN); mpfr_init(r21890); mpfr_init(r21891); mpfr_init(r21892); 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); } double f_dm(double l, double Om, double kx, double ky) { ; ; mpfr_div(r21890, r21888, r21889, MPFR_RNDN); mpfr_set_d(r21891, l, MPFR_RNDN); mpfr_mul(r21892, r21889, r21891, MPFR_RNDN); mpfr_set_d(r21893, Om, MPFR_RNDN); mpfr_div(r21894, r21892, r21893, MPFR_RNDN); mpfr_pow(r21895, r21894, r21889, MPFR_RNDN); mpfr_set_d(r21896, kx, MPFR_RNDN); mpfr_sin(r21897, r21896, MPFR_RNDN); mpfr_pow(r21898, r21897, r21889, MPFR_RNDN); mpfr_set_d(r21899, ky, MPFR_RNDN); mpfr_sin(r21900, r21899, MPFR_RNDN); mpfr_pow(r21901, r21900, r21889, MPFR_RNDN); mpfr_add(r21902, r21898, r21901, MPFR_RNDN); mpfr_mul(r21903, r21895, r21902, MPFR_RNDN); mpfr_add(r21904, r21888, r21903, MPFR_RNDN); mpfr_sqrt(r21905, r21904, MPFR_RNDN); mpfr_div(r21906, r21888, r21905, MPFR_RNDN); mpfr_exp(r21907, r21906, MPFR_RNDN); mpfr_log(r21908, r21907, MPFR_RNDN); mpfr_add(r21909, r21888, r21908, MPFR_RNDN); mpfr_mul(r21910, r21890, r21909, MPFR_RNDN); mpfr_sqrt(r21911, r21910, MPFR_RNDN); return mpfr_get_d(r21911, MPFR_RNDN); }