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