#include #include #include #include #include char *name = "Hyperbolic arcsine"; double f_if(float x) { float r19742 = x; float r19743 = r19742 * r19742; float r19744 = 1.0f; float r19745 = r19743 + r19744; float r19746 = sqrt(r19745); float r19747 = r19742 + r19746; float r19748 = log(r19747); return r19748; } double f_id(double x) { double r19749 = x; double r19750 = r19749 * r19749; double r19751 = 1.0; double r19752 = r19750 + r19751; double r19753 = sqrt(r19752); double r19754 = r19749 + r19753; double r19755 = log(r19754); return r19755; } double f_of(float x) { float r19756 = x; float r19757 = -0.11215751184115397f; bool r19758 = r19756 <= r19757; float r19759 = 0.125f; float r19760 = r19756 * (r19756 * r19756); float r19761 = r19759 / r19760; float r19762 = 0.5f; float r19763 = r19762 / r19756; float r19764 = r19761 - r19763; float r19765 = log(r19764); float r19766 = 7.988163388129748f; bool r19767 = r19756 <= r19766; float r19768 = 0.075f; float r19769 = 5.0f; float r19770 = pow(r19756, r19769); float r19771 = r19768 * r19770; float r19772 = r19771 + r19756; float r19773 = 0.16666666666666666f; float r19774 = 3.0f; float r19775 = pow(r19756, r19774); float r19776 = r19773 * r19775; float r19777 = r19772 - r19776; float r19778 = r19756 + r19756; float r19779 = r19778 + r19763; float r19780 = r19779 - r19761; float r19781 = log(r19780); float r19782 = r19767 ? r19777 : r19781; float r19783 = r19758 ? r19765 : r19782; return r19783; } double f_od(double x) { double r19784 = x; double r19785 = -0.11215751184115397; bool r19786 = r19784 <= r19785; double r19787 = 0.125; double r19788 = r19784 * (r19784 * r19784); double r19789 = r19787 / r19788; double r19790 = 0.5; double r19791 = r19790 / r19784; double r19792 = r19789 - r19791; double r19793 = log(r19792); double r19794 = 7.988163388129748; bool r19795 = r19784 <= r19794; double r19796 = 0.075; double r19797 = 5.0; double r19798 = pow(r19784, r19797); double r19799 = r19796 * r19798; double r19800 = r19799 + r19784; double r19801 = 0.16666666666666666; double r19802 = 3.0; double r19803 = pow(r19784, r19802); double r19804 = r19801 * r19803; double r19805 = r19800 - r19804; double r19806 = r19784 + r19784; double r19807 = r19806 + r19791; double r19808 = r19807 - r19789; double r19809 = log(r19808); double r19810 = r19795 ? r19805 : r19809; double r19811 = r19786 ? r19793 : r19810; return r19811; } 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 r19812, r19813, r19814, r19815, r19816, r19817, r19818; void setup_mpfr_f_im() { mpfr_set_default_prec(2448); mpfr_init(r19812); mpfr_init(r19813); mpfr_init_set_str(r19814, "1", 10, MPFR_RNDN); mpfr_init(r19815); mpfr_init(r19816); mpfr_init(r19817); mpfr_init(r19818); } double f_im(double x) { mpfr_set_d(r19812, x, MPFR_RNDN); mpfr_sqr(r19813, r19812, MPFR_RNDN); ; mpfr_add(r19815, r19813, r19814, MPFR_RNDN); mpfr_sqrt(r19816, r19815, MPFR_RNDN); mpfr_add(r19817, r19812, r19816, MPFR_RNDN); mpfr_log(r19818, r19817, MPFR_RNDN); return mpfr_get_d(r19818, MPFR_RNDN); } static mpfr_t r19819, r19820, r19821, r19822, r19823, r19824, r19825, 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_fm() { mpfr_set_default_prec(2448); mpfr_init(r19819); mpfr_init_set_str(r19820, "-0.11215751184115397", 10, MPFR_RNDN); mpfr_init(r19821); mpfr_init_set_str(r19822, "1/8", 10, MPFR_RNDN); mpfr_init(r19823); mpfr_init(r19824); mpfr_init_set_str(r19825, "1/2", 10, MPFR_RNDN); mpfr_init(r19826); mpfr_init(r19827); mpfr_init(r19828); mpfr_init_set_str(r19829, "7.988163388129748", 10, MPFR_RNDN); mpfr_init(r19830); mpfr_init_set_str(r19831, "3/40", 10, MPFR_RNDN); mpfr_init_set_str(r19832, "5", 10, MPFR_RNDN); mpfr_init(r19833); mpfr_init(r19834); mpfr_init(r19835); mpfr_init_set_str(r19836, "1/6", 10, MPFR_RNDN); mpfr_init_set_str(r19837, "3", 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_fm(double x) { mpfr_set_d(r19819, x, MPFR_RNDN); ; mpfr_set_si(r19821, mpfr_cmp(r19819, r19820) <= 0, MPFR_RNDN); ; mpfr_mul(r19823, r19819, r19819, MPFR_RNDN); mpfr_mul(r19823, r19823, r19819, MPFR_RNDN); mpfr_div(r19824, r19822, r19823, MPFR_RNDN); ; mpfr_div(r19826, r19825, r19819, MPFR_RNDN); mpfr_sub(r19827, r19824, r19826, MPFR_RNDN); mpfr_log(r19828, r19827, MPFR_RNDN); ; mpfr_set_si(r19830, mpfr_cmp(r19819, r19829) <= 0, MPFR_RNDN); ; ; mpfr_pow(r19833, r19819, r19832, MPFR_RNDN); mpfr_mul(r19834, r19831, r19833, MPFR_RNDN); mpfr_add(r19835, r19834, r19819, MPFR_RNDN); ; ; mpfr_pow(r19838, r19819, r19837, MPFR_RNDN); mpfr_mul(r19839, r19836, r19838, MPFR_RNDN); mpfr_sub(r19840, r19835, r19839, MPFR_RNDN); mpfr_add(r19841, r19819, r19819, MPFR_RNDN); mpfr_add(r19842, r19841, r19826, MPFR_RNDN); mpfr_sub(r19843, r19842, r19824, MPFR_RNDN); mpfr_log(r19844, r19843, MPFR_RNDN); if (mpfr_get_si(r19830, MPFR_RNDN)) { mpfr_set(r19845, r19840, MPFR_RNDN); } else { mpfr_set(r19845, r19844, MPFR_RNDN); }; if (mpfr_get_si(r19821, MPFR_RNDN)) { mpfr_set(r19846, r19828, MPFR_RNDN); } else { mpfr_set(r19846, r19845, MPFR_RNDN); }; return mpfr_get_d(r19846, MPFR_RNDN); } static mpfr_t r19847, r19848, r19849, r19850, r19851, r19852, r19853, r19854, r19855, r19856, r19857, r19858, r19859, r19860, r19861, r19862, r19863, r19864, r19865, r19866, r19867, r19868, r19869, r19870, r19871, r19872, r19873, r19874; void setup_mpfr_f_dm() { mpfr_set_default_prec(2448); mpfr_init(r19847); mpfr_init_set_str(r19848, "-0.11215751184115397", 10, MPFR_RNDN); mpfr_init(r19849); mpfr_init_set_str(r19850, "1/8", 10, MPFR_RNDN); mpfr_init(r19851); mpfr_init(r19852); mpfr_init_set_str(r19853, "1/2", 10, MPFR_RNDN); mpfr_init(r19854); mpfr_init(r19855); mpfr_init(r19856); mpfr_init_set_str(r19857, "7.988163388129748", 10, MPFR_RNDN); mpfr_init(r19858); mpfr_init_set_str(r19859, "3/40", 10, MPFR_RNDN); mpfr_init_set_str(r19860, "5", 10, MPFR_RNDN); mpfr_init(r19861); mpfr_init(r19862); mpfr_init(r19863); mpfr_init_set_str(r19864, "1/6", 10, MPFR_RNDN); mpfr_init_set_str(r19865, "3", 10, MPFR_RNDN); mpfr_init(r19866); mpfr_init(r19867); mpfr_init(r19868); mpfr_init(r19869); mpfr_init(r19870); mpfr_init(r19871); mpfr_init(r19872); mpfr_init(r19873); mpfr_init(r19874); } double f_dm(double x) { mpfr_set_d(r19847, x, MPFR_RNDN); ; mpfr_set_si(r19849, mpfr_cmp(r19847, r19848) <= 0, MPFR_RNDN); ; mpfr_mul(r19851, r19847, r19847, MPFR_RNDN); mpfr_mul(r19851, r19851, r19847, MPFR_RNDN); mpfr_div(r19852, r19850, r19851, MPFR_RNDN); ; mpfr_div(r19854, r19853, r19847, MPFR_RNDN); mpfr_sub(r19855, r19852, r19854, MPFR_RNDN); mpfr_log(r19856, r19855, MPFR_RNDN); ; mpfr_set_si(r19858, mpfr_cmp(r19847, r19857) <= 0, MPFR_RNDN); ; ; mpfr_pow(r19861, r19847, r19860, MPFR_RNDN); mpfr_mul(r19862, r19859, r19861, MPFR_RNDN); mpfr_add(r19863, r19862, r19847, MPFR_RNDN); ; ; mpfr_pow(r19866, r19847, r19865, MPFR_RNDN); mpfr_mul(r19867, r19864, r19866, MPFR_RNDN); mpfr_sub(r19868, r19863, r19867, MPFR_RNDN); mpfr_add(r19869, r19847, r19847, MPFR_RNDN); mpfr_add(r19870, r19869, r19854, MPFR_RNDN); mpfr_sub(r19871, r19870, r19852, MPFR_RNDN); mpfr_log(r19872, r19871, MPFR_RNDN); if (mpfr_get_si(r19858, MPFR_RNDN)) { mpfr_set(r19873, r19868, MPFR_RNDN); } else { mpfr_set(r19873, r19872, MPFR_RNDN); }; if (mpfr_get_si(r19849, MPFR_RNDN)) { mpfr_set(r19874, r19856, MPFR_RNDN); } else { mpfr_set(r19874, r19873, MPFR_RNDN); }; return mpfr_get_d(r19874, MPFR_RNDN); }