#include #include #include #include #include char *name = "Toniolo and Linder, Equation (10+)"; double f_if(float t, float l, float k) { float r22782 = 2; float r22783 = t; float r22784 = 3; float r22785 = pow(r22783, r22784); float r22786 = l; float r22787 = r22786 * r22786; float r22788 = r22785 / r22787; float r22789 = k; float r22790 = sin(r22789); float r22791 = r22788 * r22790; float r22792 = tan(r22789); float r22793 = r22791 * r22792; float r22794 = 1; float r22795 = r22789 / r22783; float r22796 = pow(r22795, r22782); float r22797 = r22794 + r22796; float r22798 = r22797 + r22794; float r22799 = r22793 * r22798; float r22800 = r22782 / r22799; return r22800; } double f_id(double t, double l, double k) { double r22801 = 2; double r22802 = t; double r22803 = 3; double r22804 = pow(r22802, r22803); double r22805 = l; double r22806 = r22805 * r22805; double r22807 = r22804 / r22806; double r22808 = k; double r22809 = sin(r22808); double r22810 = r22807 * r22809; double r22811 = tan(r22808); double r22812 = r22810 * r22811; double r22813 = 1; double r22814 = r22808 / r22802; double r22815 = pow(r22814, r22801); double r22816 = r22813 + r22815; double r22817 = r22816 + r22813; double r22818 = r22812 * r22817; double r22819 = r22801 / r22818; return r22819; } double f_of(float t, float l, float k) { float r22820 = 2; float r22821 = k; float r22822 = sin(r22821); float r22823 = l; float r22824 = t; float r22825 = r22823 / r22824; float r22826 = r22822 / r22825; float r22827 = r22826 * r22824; float r22828 = 1; float r22829 = r22821 / r22824; float r22830 = pow(r22829, r22820); float r22831 = r22828 + r22830; float r22832 = r22831 + r22828; float r22833 = r22822 * r22832; float r22834 = r22827 * r22833; float r22835 = cos(r22821); float r22836 = r22825 * r22835; float r22837 = r22834 / r22836; float r22838 = r22820 / r22837; return r22838; } double f_od(double t, double l, double k) { double r22839 = 2; double r22840 = k; double r22841 = sin(r22840); double r22842 = l; double r22843 = t; double r22844 = r22842 / r22843; double r22845 = r22841 / r22844; double r22846 = r22845 * r22843; double r22847 = 1; double r22848 = r22840 / r22843; double r22849 = pow(r22848, r22839); double r22850 = r22847 + r22849; double r22851 = r22850 + r22847; double r22852 = r22841 * r22851; double r22853 = r22846 * r22852; double r22854 = cos(r22840); double r22855 = r22844 * r22854; double r22856 = r22853 / r22855; double r22857 = r22839 / r22856; return r22857; } 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 r22858, r22859, r22860, r22861, r22862, r22863, r22864, r22865, r22866, r22867, r22868, r22869, r22870, r22871, r22872, r22873, r22874, r22875, r22876; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r22858, "2", 10, MPFR_RNDN); mpfr_init(r22859); mpfr_init_set_str(r22860, "3", 10, MPFR_RNDN); mpfr_init(r22861); mpfr_init(r22862); mpfr_init(r22863); mpfr_init(r22864); mpfr_init(r22865); mpfr_init(r22866); mpfr_init(r22867); mpfr_init(r22868); mpfr_init(r22869); mpfr_init_set_str(r22870, "1", 10, MPFR_RNDN); mpfr_init(r22871); mpfr_init(r22872); mpfr_init(r22873); mpfr_init(r22874); mpfr_init(r22875); mpfr_init(r22876); } double f_im(double t, double l, double k) { ; mpfr_set_d(r22859, t, MPFR_RNDN); ; mpfr_pow(r22861, r22859, r22860, MPFR_RNDN); mpfr_set_d(r22862, l, MPFR_RNDN); mpfr_mul(r22863, r22862, r22862, MPFR_RNDN); mpfr_div(r22864, r22861, r22863, MPFR_RNDN); mpfr_set_d(r22865, k, MPFR_RNDN); mpfr_sin(r22866, r22865, MPFR_RNDN); mpfr_mul(r22867, r22864, r22866, MPFR_RNDN); mpfr_tan(r22868, r22865, MPFR_RNDN); mpfr_mul(r22869, r22867, r22868, MPFR_RNDN); ; mpfr_div(r22871, r22865, r22859, MPFR_RNDN); mpfr_pow(r22872, r22871, r22858, MPFR_RNDN); mpfr_add(r22873, r22870, r22872, MPFR_RNDN); mpfr_add(r22874, r22873, r22870, MPFR_RNDN); mpfr_mul(r22875, r22869, r22874, MPFR_RNDN); mpfr_div(r22876, r22858, r22875, MPFR_RNDN); return mpfr_get_d(r22876, MPFR_RNDN); } static mpfr_t r22877, r22878, r22879, r22880, r22881, r22882, r22883, r22884, r22885, r22886, r22887, r22888, r22889, r22890, r22891, r22892, r22893, r22894, r22895; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r22877, "2", 10, MPFR_RNDN); mpfr_init(r22878); mpfr_init(r22879); mpfr_init(r22880); mpfr_init(r22881); mpfr_init(r22882); mpfr_init(r22883); mpfr_init(r22884); mpfr_init_set_str(r22885, "1", 10, MPFR_RNDN); mpfr_init(r22886); mpfr_init(r22887); mpfr_init(r22888); mpfr_init(r22889); mpfr_init(r22890); mpfr_init(r22891); mpfr_init(r22892); mpfr_init(r22893); mpfr_init(r22894); mpfr_init(r22895); } double f_fm(double t, double l, double k) { ; mpfr_set_d(r22878, k, MPFR_RNDN); mpfr_sin(r22879, r22878, MPFR_RNDN); mpfr_set_d(r22880, l, MPFR_RNDN); mpfr_set_d(r22881, t, MPFR_RNDN); mpfr_div(r22882, r22880, r22881, MPFR_RNDN); mpfr_div(r22883, r22879, r22882, MPFR_RNDN); mpfr_mul(r22884, r22883, r22881, MPFR_RNDN); ; mpfr_div(r22886, r22878, r22881, MPFR_RNDN); mpfr_pow(r22887, r22886, r22877, MPFR_RNDN); mpfr_add(r22888, r22885, r22887, MPFR_RNDN); mpfr_add(r22889, r22888, r22885, MPFR_RNDN); mpfr_mul(r22890, r22879, r22889, MPFR_RNDN); mpfr_mul(r22891, r22884, r22890, MPFR_RNDN); mpfr_cos(r22892, r22878, MPFR_RNDN); mpfr_mul(r22893, r22882, r22892, MPFR_RNDN); mpfr_div(r22894, r22891, r22893, MPFR_RNDN); mpfr_div(r22895, r22877, r22894, MPFR_RNDN); return mpfr_get_d(r22895, MPFR_RNDN); } static mpfr_t r22896, r22897, r22898, r22899, r22900, r22901, r22902, r22903, r22904, r22905, r22906, r22907, r22908, r22909, r22910, r22911, r22912, r22913, r22914; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r22896, "2", 10, MPFR_RNDN); mpfr_init(r22897); mpfr_init(r22898); mpfr_init(r22899); mpfr_init(r22900); mpfr_init(r22901); mpfr_init(r22902); mpfr_init(r22903); mpfr_init_set_str(r22904, "1", 10, MPFR_RNDN); mpfr_init(r22905); mpfr_init(r22906); mpfr_init(r22907); mpfr_init(r22908); mpfr_init(r22909); mpfr_init(r22910); mpfr_init(r22911); mpfr_init(r22912); mpfr_init(r22913); mpfr_init(r22914); } double f_dm(double t, double l, double k) { ; mpfr_set_d(r22897, k, MPFR_RNDN); mpfr_sin(r22898, r22897, MPFR_RNDN); mpfr_set_d(r22899, l, MPFR_RNDN); mpfr_set_d(r22900, t, MPFR_RNDN); mpfr_div(r22901, r22899, r22900, MPFR_RNDN); mpfr_div(r22902, r22898, r22901, MPFR_RNDN); mpfr_mul(r22903, r22902, r22900, MPFR_RNDN); ; mpfr_div(r22905, r22897, r22900, MPFR_RNDN); mpfr_pow(r22906, r22905, r22896, MPFR_RNDN); mpfr_add(r22907, r22904, r22906, MPFR_RNDN); mpfr_add(r22908, r22907, r22904, MPFR_RNDN); mpfr_mul(r22909, r22898, r22908, MPFR_RNDN); mpfr_mul(r22910, r22903, r22909, MPFR_RNDN); mpfr_cos(r22911, r22897, MPFR_RNDN); mpfr_mul(r22912, r22901, r22911, MPFR_RNDN); mpfr_div(r22913, r22910, r22912, MPFR_RNDN); mpfr_div(r22914, r22896, r22913, MPFR_RNDN); return mpfr_get_d(r22914, MPFR_RNDN); }