#include #include #include #include #include char *name = "Toniolo and Linder, Equation (10-)"; double f_if(float t, float l, float k) { float r22718 = 2; float r22719 = t; float r22720 = 3; float r22721 = pow(r22719, r22720); float r22722 = l; float r22723 = r22722 * r22722; float r22724 = r22721 / r22723; float r22725 = k; float r22726 = sin(r22725); float r22727 = r22724 * r22726; float r22728 = tan(r22725); float r22729 = r22727 * r22728; float r22730 = 1; float r22731 = r22725 / r22719; float r22732 = pow(r22731, r22718); float r22733 = r22730 + r22732; float r22734 = r22733 - r22730; float r22735 = r22729 * r22734; float r22736 = r22718 / r22735; return r22736; } double f_id(double t, double l, double k) { double r22737 = 2; double r22738 = t; double r22739 = 3; double r22740 = pow(r22738, r22739); double r22741 = l; double r22742 = r22741 * r22741; double r22743 = r22740 / r22742; double r22744 = k; double r22745 = sin(r22744); double r22746 = r22743 * r22745; double r22747 = tan(r22744); double r22748 = r22746 * r22747; double r22749 = 1; double r22750 = r22744 / r22738; double r22751 = pow(r22750, r22737); double r22752 = r22749 + r22751; double r22753 = r22752 - r22749; double r22754 = r22748 * r22753; double r22755 = r22737 / r22754; return r22755; } double f_of(float t, float l, float k) { float r22756 = 2; float r22757 = k; float r22758 = r22756 / r22757; float r22759 = 1; float r22760 = l; float r22761 = r22759 / r22760; float r22762 = r22758 / r22761; float r22763 = sin(r22757); float r22764 = -r22763; float r22765 = r22762 / r22764; float r22766 = cos(r22757); float r22767 = t; float r22768 = r22766 / r22767; float r22769 = -r22760; float r22770 = r22757 / r22769; float r22771 = r22768 / r22770; float r22772 = r22771 / r22763; float r22773 = r22765 * r22772; return r22773; } double f_od(double t, double l, double k) { double r22774 = 2; double r22775 = k; double r22776 = r22774 / r22775; double r22777 = 1; double r22778 = l; double r22779 = r22777 / r22778; double r22780 = r22776 / r22779; double r22781 = sin(r22775); double r22782 = -r22781; double r22783 = r22780 / r22782; double r22784 = cos(r22775); double r22785 = t; double r22786 = r22784 / r22785; double r22787 = -r22778; double r22788 = r22775 / r22787; double r22789 = r22786 / r22788; double r22790 = r22789 / r22781; double r22791 = r22783 * r22790; return r22791; } 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 r22792, r22793, r22794, r22795, r22796, r22797, r22798, r22799, r22800, r22801, r22802, r22803, r22804, r22805, r22806, r22807, r22808, r22809, r22810; void setup_mpfr_f_im() { mpfr_set_default_prec(4240); mpfr_init_set_str(r22792, "2", 10, MPFR_RNDN); mpfr_init(r22793); mpfr_init_set_str(r22794, "3", 10, MPFR_RNDN); mpfr_init(r22795); mpfr_init(r22796); mpfr_init(r22797); mpfr_init(r22798); mpfr_init(r22799); mpfr_init(r22800); mpfr_init(r22801); mpfr_init(r22802); mpfr_init(r22803); mpfr_init_set_str(r22804, "1", 10, MPFR_RNDN); mpfr_init(r22805); mpfr_init(r22806); mpfr_init(r22807); mpfr_init(r22808); mpfr_init(r22809); mpfr_init(r22810); } double f_im(double t, double l, double k) { ; mpfr_set_d(r22793, t, MPFR_RNDN); ; mpfr_pow(r22795, r22793, r22794, MPFR_RNDN); mpfr_set_d(r22796, l, MPFR_RNDN); mpfr_mul(r22797, r22796, r22796, MPFR_RNDN); mpfr_div(r22798, r22795, r22797, MPFR_RNDN); mpfr_set_d(r22799, k, MPFR_RNDN); mpfr_sin(r22800, r22799, MPFR_RNDN); mpfr_mul(r22801, r22798, r22800, MPFR_RNDN); mpfr_tan(r22802, r22799, MPFR_RNDN); mpfr_mul(r22803, r22801, r22802, MPFR_RNDN); ; mpfr_div(r22805, r22799, r22793, MPFR_RNDN); mpfr_pow(r22806, r22805, r22792, MPFR_RNDN); mpfr_add(r22807, r22804, r22806, MPFR_RNDN); mpfr_sub(r22808, r22807, r22804, MPFR_RNDN); mpfr_mul(r22809, r22803, r22808, MPFR_RNDN); mpfr_div(r22810, r22792, r22809, MPFR_RNDN); return mpfr_get_d(r22810, MPFR_RNDN); } static mpfr_t r22811, r22812, r22813, r22814, r22815, r22816, r22817, r22818, r22819, r22820, r22821, r22822, r22823, r22824, r22825, r22826, r22827, r22828; void setup_mpfr_f_fm() { mpfr_set_default_prec(4240); mpfr_init_set_str(r22811, "2", 10, MPFR_RNDN); mpfr_init(r22812); mpfr_init(r22813); mpfr_init_set_str(r22814, "1", 10, MPFR_RNDN); mpfr_init(r22815); mpfr_init(r22816); mpfr_init(r22817); mpfr_init(r22818); mpfr_init(r22819); mpfr_init(r22820); mpfr_init(r22821); mpfr_init(r22822); mpfr_init(r22823); mpfr_init(r22824); mpfr_init(r22825); mpfr_init(r22826); mpfr_init(r22827); mpfr_init(r22828); } double f_fm(double t, double l, double k) { ; mpfr_set_d(r22812, k, MPFR_RNDN); mpfr_div(r22813, r22811, r22812, MPFR_RNDN); ; mpfr_set_d(r22815, l, MPFR_RNDN); mpfr_div(r22816, r22814, r22815, MPFR_RNDN); mpfr_div(r22817, r22813, r22816, MPFR_RNDN); mpfr_sin(r22818, r22812, MPFR_RNDN); mpfr_neg(r22819, r22818, MPFR_RNDN); mpfr_div(r22820, r22817, r22819, MPFR_RNDN); mpfr_cos(r22821, r22812, MPFR_RNDN); mpfr_set_d(r22822, t, MPFR_RNDN); mpfr_div(r22823, r22821, r22822, MPFR_RNDN); mpfr_neg(r22824, r22815, MPFR_RNDN); mpfr_div(r22825, r22812, r22824, MPFR_RNDN); mpfr_div(r22826, r22823, r22825, MPFR_RNDN); mpfr_div(r22827, r22826, r22818, MPFR_RNDN); mpfr_mul(r22828, r22820, r22827, MPFR_RNDN); return mpfr_get_d(r22828, MPFR_RNDN); } static mpfr_t r22829, r22830, r22831, r22832, r22833, r22834, r22835, r22836, r22837, r22838, r22839, r22840, r22841, r22842, r22843, r22844, r22845, r22846; void setup_mpfr_f_dm() { mpfr_set_default_prec(4240); mpfr_init_set_str(r22829, "2", 10, MPFR_RNDN); mpfr_init(r22830); mpfr_init(r22831); mpfr_init_set_str(r22832, "1", 10, MPFR_RNDN); mpfr_init(r22833); mpfr_init(r22834); mpfr_init(r22835); mpfr_init(r22836); mpfr_init(r22837); mpfr_init(r22838); mpfr_init(r22839); mpfr_init(r22840); mpfr_init(r22841); mpfr_init(r22842); mpfr_init(r22843); mpfr_init(r22844); mpfr_init(r22845); mpfr_init(r22846); } double f_dm(double t, double l, double k) { ; mpfr_set_d(r22830, k, MPFR_RNDN); mpfr_div(r22831, r22829, r22830, MPFR_RNDN); ; mpfr_set_d(r22833, l, MPFR_RNDN); mpfr_div(r22834, r22832, r22833, MPFR_RNDN); mpfr_div(r22835, r22831, r22834, MPFR_RNDN); mpfr_sin(r22836, r22830, MPFR_RNDN); mpfr_neg(r22837, r22836, MPFR_RNDN); mpfr_div(r22838, r22835, r22837, MPFR_RNDN); mpfr_cos(r22839, r22830, MPFR_RNDN); mpfr_set_d(r22840, t, MPFR_RNDN); mpfr_div(r22841, r22839, r22840, MPFR_RNDN); mpfr_neg(r22842, r22833, MPFR_RNDN); mpfr_div(r22843, r22830, r22842, MPFR_RNDN); mpfr_div(r22844, r22841, r22843, MPFR_RNDN); mpfr_div(r22845, r22844, r22836, MPFR_RNDN); mpfr_mul(r22846, r22838, r22845, MPFR_RNDN); return mpfr_get_d(r22846, MPFR_RNDN); }