#include #include #include #include #include char *name = "VandenBroeck and Keller, Equation (23)"; double f_if(float F, float B, float x) { float r19703 = x; float r19704 = 1; float r19705 = B; float r19706 = tan(r19705); float r19707 = r19704 / r19706; float r19708 = r19703 * r19707; float r19709 = -r19708; float r19710 = F; float r19711 = sin(r19705); float r19712 = r19710 / r19711; float r19713 = r19710 * r19710; float r19714 = 2; float r19715 = r19713 + r19714; float r19716 = r19714 * r19703; float r19717 = r19715 + r19716; float r19718 = r19704 / r19714; float r19719 = -r19718; float r19720 = pow(r19717, r19719); float r19721 = r19712 * r19720; float r19722 = r19709 + r19721; return r19722; } double f_id(double F, double B, double x) { double r19723 = x; double r19724 = 1; double r19725 = B; double r19726 = tan(r19725); double r19727 = r19724 / r19726; double r19728 = r19723 * r19727; double r19729 = -r19728; double r19730 = F; double r19731 = sin(r19725); double r19732 = r19730 / r19731; double r19733 = r19730 * r19730; double r19734 = 2; double r19735 = r19733 + r19734; double r19736 = r19734 * r19723; double r19737 = r19735 + r19736; double r19738 = r19724 / r19734; double r19739 = -r19738; double r19740 = pow(r19737, r19739); double r19741 = r19732 * r19740; double r19742 = r19729 + r19741; return r19742; } double f_of(float F, float B, float x) { float r19743 = x; float r19744 = -r19743; float r19745 = B; float r19746 = tan(r19745); float r19747 = r19744 / r19746; float r19748 = F; float r19749 = r19743 + r19743; float r19750 = 2; float r19751 = r19748 * r19748; float r19752 = r19750 + r19751; float r19753 = r19749 + r19752; float r19754 = 1; float r19755 = r19754 / r19750; float r19756 = pow(r19753, r19755); float r19757 = sin(r19745); float r19758 = r19756 * r19757; float r19759 = cbrt(r19758); float r19760 = r19759 * r19759; float r19761 = r19760 * r19759; float r19762 = r19748 / r19761; float r19763 = r19747 + r19762; return r19763; } double f_od(double F, double B, double x) { double r19764 = x; double r19765 = -r19764; double r19766 = B; double r19767 = tan(r19766); double r19768 = r19765 / r19767; double r19769 = F; double r19770 = r19764 + r19764; double r19771 = 2; double r19772 = r19769 * r19769; double r19773 = r19771 + r19772; double r19774 = r19770 + r19773; double r19775 = 1; double r19776 = r19775 / r19771; double r19777 = pow(r19774, r19776); double r19778 = sin(r19766); double r19779 = r19777 * r19778; double r19780 = cbrt(r19779); double r19781 = r19780 * r19780; double r19782 = r19781 * r19780; double r19783 = r19769 / r19782; double r19784 = r19768 + r19783; return r19784; } 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 r19785, r19786, r19787, r19788, r19789, r19790, r19791, r19792, r19793, r19794, r19795, r19796, r19797, r19798, r19799, r19800, r19801, r19802, r19803, r19804; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r19785); mpfr_init_set_str(r19786, "1", 10, MPFR_RNDN); mpfr_init(r19787); mpfr_init(r19788); mpfr_init(r19789); mpfr_init(r19790); mpfr_init(r19791); mpfr_init(r19792); mpfr_init(r19793); mpfr_init(r19794); mpfr_init(r19795); mpfr_init_set_str(r19796, "2", 10, MPFR_RNDN); mpfr_init(r19797); mpfr_init(r19798); mpfr_init(r19799); mpfr_init(r19800); mpfr_init(r19801); mpfr_init(r19802); mpfr_init(r19803); mpfr_init(r19804); } double f_im(double F, double B, double x) { mpfr_set_d(r19785, x, MPFR_RNDN); ; mpfr_set_d(r19787, B, MPFR_RNDN); mpfr_tan(r19788, r19787, MPFR_RNDN); mpfr_div(r19789, r19786, r19788, MPFR_RNDN); mpfr_mul(r19790, r19785, r19789, MPFR_RNDN); mpfr_neg(r19791, r19790, MPFR_RNDN); mpfr_set_d(r19792, F, MPFR_RNDN); mpfr_sin(r19793, r19787, MPFR_RNDN); mpfr_div(r19794, r19792, r19793, MPFR_RNDN); mpfr_mul(r19795, r19792, r19792, MPFR_RNDN); ; mpfr_add(r19797, r19795, r19796, MPFR_RNDN); mpfr_mul(r19798, r19796, r19785, MPFR_RNDN); mpfr_add(r19799, r19797, r19798, MPFR_RNDN); mpfr_div(r19800, r19786, r19796, MPFR_RNDN); mpfr_neg(r19801, r19800, MPFR_RNDN); mpfr_pow(r19802, r19799, r19801, MPFR_RNDN); mpfr_mul(r19803, r19794, r19802, MPFR_RNDN); mpfr_add(r19804, r19791, r19803, MPFR_RNDN); return mpfr_get_d(r19804, MPFR_RNDN); } static mpfr_t r19805, r19806, r19807, r19808, r19809, r19810, r19811, r19812, r19813, r19814, r19815, r19816, r19817, r19818, r19819, r19820, r19821, r19822, r19823, r19824, r19825; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r19805); mpfr_init(r19806); mpfr_init(r19807); mpfr_init(r19808); mpfr_init(r19809); mpfr_init(r19810); mpfr_init(r19811); mpfr_init_set_str(r19812, "2", 10, MPFR_RNDN); mpfr_init(r19813); mpfr_init(r19814); mpfr_init(r19815); mpfr_init_set_str(r19816, "1", 10, MPFR_RNDN); mpfr_init(r19817); mpfr_init(r19818); mpfr_init(r19819); mpfr_init(r19820); mpfr_init(r19821); mpfr_init(r19822); mpfr_init(r19823); mpfr_init(r19824); mpfr_init(r19825); } double f_fm(double F, double B, double x) { mpfr_set_d(r19805, x, MPFR_RNDN); mpfr_neg(r19806, r19805, MPFR_RNDN); mpfr_set_d(r19807, B, MPFR_RNDN); mpfr_tan(r19808, r19807, MPFR_RNDN); mpfr_div(r19809, r19806, r19808, MPFR_RNDN); mpfr_set_d(r19810, F, MPFR_RNDN); mpfr_add(r19811, r19805, r19805, MPFR_RNDN); ; mpfr_mul(r19813, r19810, r19810, MPFR_RNDN); mpfr_add(r19814, r19812, r19813, MPFR_RNDN); mpfr_add(r19815, r19811, r19814, MPFR_RNDN); ; mpfr_div(r19817, r19816, r19812, MPFR_RNDN); mpfr_pow(r19818, r19815, r19817, MPFR_RNDN); mpfr_sin(r19819, r19807, MPFR_RNDN); mpfr_mul(r19820, r19818, r19819, MPFR_RNDN); mpfr_cbrt(r19821, r19820, MPFR_RNDN); mpfr_mul(r19822, r19821, r19821, MPFR_RNDN); mpfr_mul(r19823, r19822, r19821, MPFR_RNDN); mpfr_div(r19824, r19810, r19823, MPFR_RNDN); mpfr_add(r19825, r19809, r19824, MPFR_RNDN); return mpfr_get_d(r19825, MPFR_RNDN); } static mpfr_t r19826, r19827, r19828, r19829, r19830, r19831, r19832, r19833, r19834, r19835, r19836, r19837, r19838, r19839, r19840, r19841, r19842, r19843, r19844, r19845, r19846; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r19826); mpfr_init(r19827); mpfr_init(r19828); mpfr_init(r19829); mpfr_init(r19830); mpfr_init(r19831); mpfr_init(r19832); mpfr_init_set_str(r19833, "2", 10, MPFR_RNDN); mpfr_init(r19834); mpfr_init(r19835); mpfr_init(r19836); mpfr_init_set_str(r19837, "1", 10, MPFR_RNDN); mpfr_init(r19838); mpfr_init(r19839); mpfr_init(r19840); mpfr_init(r19841); mpfr_init(r19842); mpfr_init(r19843); mpfr_init(r19844); mpfr_init(r19845); mpfr_init(r19846); } double f_dm(double F, double B, double x) { mpfr_set_d(r19826, x, MPFR_RNDN); mpfr_neg(r19827, r19826, MPFR_RNDN); mpfr_set_d(r19828, B, MPFR_RNDN); mpfr_tan(r19829, r19828, MPFR_RNDN); mpfr_div(r19830, r19827, r19829, MPFR_RNDN); mpfr_set_d(r19831, F, MPFR_RNDN); mpfr_add(r19832, r19826, r19826, MPFR_RNDN); ; mpfr_mul(r19834, r19831, r19831, MPFR_RNDN); mpfr_add(r19835, r19833, r19834, MPFR_RNDN); mpfr_add(r19836, r19832, r19835, MPFR_RNDN); ; mpfr_div(r19838, r19837, r19833, MPFR_RNDN); mpfr_pow(r19839, r19836, r19838, MPFR_RNDN); mpfr_sin(r19840, r19828, MPFR_RNDN); mpfr_mul(r19841, r19839, r19840, MPFR_RNDN); mpfr_cbrt(r19842, r19841, MPFR_RNDN); mpfr_mul(r19843, r19842, r19842, MPFR_RNDN); mpfr_mul(r19844, r19843, r19842, MPFR_RNDN); mpfr_div(r19845, r19831, r19844, MPFR_RNDN); mpfr_add(r19846, r19830, r19845, MPFR_RNDN); return mpfr_get_d(r19846, MPFR_RNDN); }