#include #include #include #include #include char *name = "Jmat.Real.dawson"; double f_if(float x) { float r10628 = 1.0f; float r10629 = 0.1049934947f; float r10630 = x; float r10631 = r10630 * r10630; float r10632 = r10629 * r10631; float r10633 = r10628 + r10632; float r10634 = 0.0424060604f; float r10635 = r10631 * r10631; float r10636 = r10634 * r10635; float r10637 = r10633 + r10636; float r10638 = 0.0072644182f; float r10639 = r10635 * r10631; float r10640 = r10638 * r10639; float r10641 = r10637 + r10640; float r10642 = 0.0005064034f; float r10643 = r10639 * r10631; float r10644 = r10642 * r10643; float r10645 = r10641 + r10644; float r10646 = 0.0001789971f; float r10647 = r10643 * r10631; float r10648 = r10646 * r10647; float r10649 = r10645 + r10648; float r10650 = 0.7715471019f; float r10651 = r10650 * r10631; float r10652 = r10628 + r10651; float r10653 = 0.2909738639f; float r10654 = r10653 * r10635; float r10655 = r10652 + r10654; float r10656 = 0.0694555761f; float r10657 = r10656 * r10639; float r10658 = r10655 + r10657; float r10659 = 0.0140005442f; float r10660 = r10659 * r10643; float r10661 = r10658 + r10660; float r10662 = 0.0008327945f; float r10663 = r10662 * r10647; float r10664 = r10661 + r10663; float r10665 = 2.0f; float r10666 = r10665 * r10646; float r10667 = r10647 * r10631; float r10668 = r10666 * r10667; float r10669 = r10664 + r10668; float r10670 = r10649 / r10669; float r10671 = r10670 * r10630; return r10671; } double f_id(double x) { double r10672 = 1.0; double r10673 = 0.1049934947; double r10674 = x; double r10675 = r10674 * r10674; double r10676 = r10673 * r10675; double r10677 = r10672 + r10676; double r10678 = 0.0424060604; double r10679 = r10675 * r10675; double r10680 = r10678 * r10679; double r10681 = r10677 + r10680; double r10682 = 0.0072644182; double r10683 = r10679 * r10675; double r10684 = r10682 * r10683; double r10685 = r10681 + r10684; double r10686 = 0.0005064034; double r10687 = r10683 * r10675; double r10688 = r10686 * r10687; double r10689 = r10685 + r10688; double r10690 = 0.0001789971; double r10691 = r10687 * r10675; double r10692 = r10690 * r10691; double r10693 = r10689 + r10692; double r10694 = 0.7715471019; double r10695 = r10694 * r10675; double r10696 = r10672 + r10695; double r10697 = 0.2909738639; double r10698 = r10697 * r10679; double r10699 = r10696 + r10698; double r10700 = 0.0694555761; double r10701 = r10700 * r10683; double r10702 = r10699 + r10701; double r10703 = 0.0140005442; double r10704 = r10703 * r10687; double r10705 = r10702 + r10704; double r10706 = 0.0008327945; double r10707 = r10706 * r10691; double r10708 = r10705 + r10707; double r10709 = 2.0; double r10710 = r10709 * r10690; double r10711 = r10691 * r10675; double r10712 = r10710 * r10711; double r10713 = r10708 + r10712; double r10714 = r10693 / r10713; double r10715 = r10714 * r10674; return r10715; } double f_of(float x) { float r10716 = x; float r10717 = -4456723949048.314f; bool r10718 = r10716 <= r10717; float r10719 = 896.8339123666102f; bool r10720 = r10716 <= r10719; float r10721 = !r10720; bool r10722 = r10718 || r10721; float r10723 = 0.2514179000665375f; float r10724 = r10723 / r10716; float r10725 = r10716 * r10716; float r10726 = r10724 / r10725; float r10727 = 0.5f; float r10728 = r10727 / r10716; float r10729 = 0.15298196345928972f; float r10730 = 5.0f; float r10731 = pow(r10716, r10730); float r10732 = r10729 / r10731; float r10733 = r10728 + r10732; float r10734 = r10726 + r10733; float r10735 = 1.0f; float r10736 = r10725 * r10725; float r10737 = r10736 * r10725; float r10738 = r10737 * r10725; float r10739 = r10725 * r10738; float r10740 = r10739 * r10725; float r10741 = 2.0f; float r10742 = 0.0001789971f; float r10743 = r10741 * r10742; float r10744 = r10740 * r10743; float r10745 = 0.0008327945f; float r10746 = r10739 * r10745; float r10747 = 0.7715471019f; float r10748 = r10725 * r10747; float r10749 = r10748 + r10735; float r10750 = 0.2909738639f; float r10751 = r10750 * r10736; float r10752 = r10749 + r10751; float r10753 = 0.0694555761f; float r10754 = r10737 * r10753; float r10755 = r10752 + r10754; float r10756 = 0.0140005442f; float r10757 = r10738 * r10756; float r10758 = r10755 + r10757; float r10759 = r10746 + r10758; float r10760 = r10744 + r10759; float r10761 = r10735 / r10760; float r10762 = r10739 * r10742; float r10763 = 0.0072644182f; float r10764 = r10763 * r10737; float r10765 = 0.1049934947f; float r10766 = r10725 * r10765; float r10767 = r10735 + r10766; float r10768 = 0.0424060604f; float r10769 = r10768 * r10736; float r10770 = r10767 + r10769; float r10771 = r10764 + r10770; float r10772 = 0.0005064034f; float r10773 = r10738 * r10772; float r10774 = r10771 + r10773; float r10775 = r10762 + r10774; float r10776 = r10761 * r10775; float r10777 = r10716 * r10776; float r10778 = r10722 ? r10734 : r10777; return r10778; } double f_od(double x) { double r10779 = x; double r10780 = -4456723949048.314; bool r10781 = r10779 <= r10780; double r10782 = 896.8339123666102; bool r10783 = r10779 <= r10782; double r10784 = !r10783; bool r10785 = r10781 || r10784; double r10786 = 0.2514179000665375; double r10787 = r10786 / r10779; double r10788 = r10779 * r10779; double r10789 = r10787 / r10788; double r10790 = 0.5; double r10791 = r10790 / r10779; double r10792 = 0.15298196345928972; double r10793 = 5.0; double r10794 = pow(r10779, r10793); double r10795 = r10792 / r10794; double r10796 = r10791 + r10795; double r10797 = r10789 + r10796; double r10798 = 1.0; double r10799 = r10788 * r10788; double r10800 = r10799 * r10788; double r10801 = r10800 * r10788; double r10802 = r10788 * r10801; double r10803 = r10802 * r10788; double r10804 = 2.0; double r10805 = 0.0001789971; double r10806 = r10804 * r10805; double r10807 = r10803 * r10806; double r10808 = 0.0008327945; double r10809 = r10802 * r10808; double r10810 = 0.7715471019; double r10811 = r10788 * r10810; double r10812 = r10811 + r10798; double r10813 = 0.2909738639; double r10814 = r10813 * r10799; double r10815 = r10812 + r10814; double r10816 = 0.0694555761; double r10817 = r10800 * r10816; double r10818 = r10815 + r10817; double r10819 = 0.0140005442; double r10820 = r10801 * r10819; double r10821 = r10818 + r10820; double r10822 = r10809 + r10821; double r10823 = r10807 + r10822; double r10824 = r10798 / r10823; double r10825 = r10802 * r10805; double r10826 = 0.0072644182; double r10827 = r10826 * r10800; double r10828 = 0.1049934947; double r10829 = r10788 * r10828; double r10830 = r10798 + r10829; double r10831 = 0.0424060604; double r10832 = r10831 * r10799; double r10833 = r10830 + r10832; double r10834 = r10827 + r10833; double r10835 = 0.0005064034; double r10836 = r10801 * r10835; double r10837 = r10834 + r10836; double r10838 = r10825 + r10837; double r10839 = r10824 * r10838; double r10840 = r10779 * r10839; double r10841 = r10785 ? r10797 : r10840; return r10841; } 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 r10842, r10843, r10844, r10845, r10846, r10847, r10848, r10849, r10850, r10851, r10852, r10853, r10854, r10855, r10856, r10857, r10858, r10859, r10860, r10861, r10862, r10863, r10864, r10865, r10866, r10867, r10868, r10869, r10870, r10871, r10872, r10873, r10874, r10875, r10876, r10877, r10878, r10879, r10880, r10881, r10882, r10883, r10884, r10885; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init_set_str(r10842, "1", 10, MPFR_RNDN); mpfr_init_set_str(r10843, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r10844); mpfr_init(r10845); mpfr_init(r10846); mpfr_init(r10847); mpfr_init_set_str(r10848, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r10849); mpfr_init(r10850); mpfr_init(r10851); mpfr_init_set_str(r10852, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r10853); mpfr_init(r10854); mpfr_init(r10855); mpfr_init_set_str(r10856, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r10857); mpfr_init(r10858); mpfr_init(r10859); mpfr_init_set_str(r10860, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r10861); mpfr_init(r10862); mpfr_init(r10863); mpfr_init_set_str(r10864, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r10865); mpfr_init(r10866); mpfr_init_set_str(r10867, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r10868); mpfr_init(r10869); mpfr_init_set_str(r10870, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r10871); mpfr_init(r10872); mpfr_init_set_str(r10873, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r10874); mpfr_init(r10875); mpfr_init_set_str(r10876, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r10877); mpfr_init(r10878); mpfr_init_set_str(r10879, "2", 10, MPFR_RNDN); mpfr_init(r10880); mpfr_init(r10881); mpfr_init(r10882); mpfr_init(r10883); mpfr_init(r10884); mpfr_init(r10885); } double f_im(double x) { ; ; mpfr_set_d(r10844, x, MPFR_RNDN); mpfr_mul(r10845, r10844, r10844, MPFR_RNDN); mpfr_mul(r10846, r10843, r10845, MPFR_RNDN); mpfr_add(r10847, r10842, r10846, MPFR_RNDN); ; mpfr_mul(r10849, r10845, r10845, MPFR_RNDN); mpfr_mul(r10850, r10848, r10849, MPFR_RNDN); mpfr_add(r10851, r10847, r10850, MPFR_RNDN); ; mpfr_mul(r10853, r10849, r10845, MPFR_RNDN); mpfr_mul(r10854, r10852, r10853, MPFR_RNDN); mpfr_add(r10855, r10851, r10854, MPFR_RNDN); ; mpfr_mul(r10857, r10853, r10845, MPFR_RNDN); mpfr_mul(r10858, r10856, r10857, MPFR_RNDN); mpfr_add(r10859, r10855, r10858, MPFR_RNDN); ; mpfr_mul(r10861, r10857, r10845, MPFR_RNDN); mpfr_mul(r10862, r10860, r10861, MPFR_RNDN); mpfr_add(r10863, r10859, r10862, MPFR_RNDN); ; mpfr_mul(r10865, r10864, r10845, MPFR_RNDN); mpfr_add(r10866, r10842, r10865, MPFR_RNDN); ; mpfr_mul(r10868, r10867, r10849, MPFR_RNDN); mpfr_add(r10869, r10866, r10868, MPFR_RNDN); ; mpfr_mul(r10871, r10870, r10853, MPFR_RNDN); mpfr_add(r10872, r10869, r10871, MPFR_RNDN); ; mpfr_mul(r10874, r10873, r10857, MPFR_RNDN); mpfr_add(r10875, r10872, r10874, MPFR_RNDN); ; mpfr_mul(r10877, r10876, r10861, MPFR_RNDN); mpfr_add(r10878, r10875, r10877, MPFR_RNDN); ; mpfr_mul(r10880, r10879, r10860, MPFR_RNDN); mpfr_mul(r10881, r10861, r10845, MPFR_RNDN); mpfr_mul(r10882, r10880, r10881, MPFR_RNDN); mpfr_add(r10883, r10878, r10882, MPFR_RNDN); mpfr_div(r10884, r10863, r10883, MPFR_RNDN); mpfr_mul(r10885, r10884, r10844, MPFR_RNDN); return mpfr_get_d(r10885, MPFR_RNDN); } static mpfr_t r10886, r10887, r10888, r10889, r10890, r10891, r10892, r10893, r10894, r10895, r10896, r10897, r10898, r10899, r10900, r10901, r10902, r10903, r10904, r10905, r10906, r10907, r10908, r10909, r10910, r10911, r10912, r10913, r10914, r10915, r10916, r10917, r10918, r10919, r10920, r10921, r10922, r10923, r10924, r10925, r10926, r10927, r10928, r10929, r10930, r10931, r10932, r10933, r10934, r10935, r10936, r10937, r10938, r10939, r10940, r10941, r10942, r10943, r10944, r10945, r10946, r10947, r10948; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init(r10886); mpfr_init_set_str(r10887, "-4456723949048.314", 10, MPFR_RNDN); mpfr_init(r10888); mpfr_init_set_str(r10889, "896.8339123666102", 10, MPFR_RNDN); mpfr_init(r10890); mpfr_init(r10891); mpfr_init(r10892); mpfr_init_set_str(r10893, "0.2514179000665375", 10, MPFR_RNDN); mpfr_init(r10894); mpfr_init(r10895); mpfr_init(r10896); mpfr_init_set_str(r10897, "0.5", 10, MPFR_RNDN); mpfr_init(r10898); mpfr_init_set_str(r10899, "0.15298196345928972", 10, MPFR_RNDN); mpfr_init_set_str(r10900, "5", 10, MPFR_RNDN); mpfr_init(r10901); mpfr_init(r10902); mpfr_init(r10903); mpfr_init(r10904); mpfr_init_set_str(r10905, "1", 10, MPFR_RNDN); mpfr_init(r10906); mpfr_init(r10907); mpfr_init(r10908); mpfr_init(r10909); mpfr_init(r10910); mpfr_init_set_str(r10911, "2", 10, MPFR_RNDN); mpfr_init_set_str(r10912, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r10913); mpfr_init(r10914); mpfr_init_set_str(r10915, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r10916); mpfr_init_set_str(r10917, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r10918); mpfr_init(r10919); mpfr_init_set_str(r10920, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r10921); mpfr_init(r10922); mpfr_init_set_str(r10923, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r10924); mpfr_init(r10925); mpfr_init_set_str(r10926, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r10927); mpfr_init(r10928); mpfr_init(r10929); mpfr_init(r10930); mpfr_init(r10931); mpfr_init(r10932); mpfr_init_set_str(r10933, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r10934); mpfr_init_set_str(r10935, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r10936); mpfr_init(r10937); mpfr_init_set_str(r10938, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r10939); mpfr_init(r10940); mpfr_init(r10941); mpfr_init_set_str(r10942, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r10943); mpfr_init(r10944); mpfr_init(r10945); mpfr_init(r10946); mpfr_init(r10947); mpfr_init(r10948); } double f_fm(double x) { mpfr_set_d(r10886, x, MPFR_RNDN); ; mpfr_set_si(r10888, mpfr_cmp(r10886, r10887) <= 0, MPFR_RNDN); ; mpfr_set_si(r10890, mpfr_cmp(r10886, r10889) <= 0, MPFR_RNDN); mpfr_set_si(r10891, !mpfr_get_si(r10890, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r10892, mpfr_get_si(r10888, MPFR_RNDN) || mpfr_get_si(r10891, MPFR_RNDN), MPFR_RNDN); ; mpfr_div(r10894, r10893, r10886, MPFR_RNDN); mpfr_mul(r10895, r10886, r10886, MPFR_RNDN); mpfr_div(r10896, r10894, r10895, MPFR_RNDN); ; mpfr_div(r10898, r10897, r10886, MPFR_RNDN); ; ; mpfr_pow(r10901, r10886, r10900, MPFR_RNDN); mpfr_div(r10902, r10899, r10901, MPFR_RNDN); mpfr_add(r10903, r10898, r10902, MPFR_RNDN); mpfr_add(r10904, r10896, r10903, MPFR_RNDN); ; mpfr_mul(r10906, r10895, r10895, MPFR_RNDN); mpfr_mul(r10907, r10906, r10895, MPFR_RNDN); mpfr_mul(r10908, r10907, r10895, MPFR_RNDN); mpfr_mul(r10909, r10895, r10908, MPFR_RNDN); mpfr_mul(r10910, r10909, r10895, MPFR_RNDN); ; ; mpfr_mul(r10913, r10911, r10912, MPFR_RNDN); mpfr_mul(r10914, r10910, r10913, MPFR_RNDN); ; mpfr_mul(r10916, r10909, r10915, MPFR_RNDN); ; mpfr_mul(r10918, r10895, r10917, MPFR_RNDN); mpfr_add(r10919, r10918, r10905, MPFR_RNDN); ; mpfr_mul(r10921, r10920, r10906, MPFR_RNDN); mpfr_add(r10922, r10919, r10921, MPFR_RNDN); ; mpfr_mul(r10924, r10907, r10923, MPFR_RNDN); mpfr_add(r10925, r10922, r10924, MPFR_RNDN); ; mpfr_mul(r10927, r10908, r10926, MPFR_RNDN); mpfr_add(r10928, r10925, r10927, MPFR_RNDN); mpfr_add(r10929, r10916, r10928, MPFR_RNDN); mpfr_add(r10930, r10914, r10929, MPFR_RNDN); mpfr_div(r10931, r10905, r10930, MPFR_RNDN); mpfr_mul(r10932, r10909, r10912, MPFR_RNDN); ; mpfr_mul(r10934, r10933, r10907, MPFR_RNDN); ; mpfr_mul(r10936, r10895, r10935, MPFR_RNDN); mpfr_add(r10937, r10905, r10936, MPFR_RNDN); ; mpfr_mul(r10939, r10938, r10906, MPFR_RNDN); mpfr_add(r10940, r10937, r10939, MPFR_RNDN); mpfr_add(r10941, r10934, r10940, MPFR_RNDN); ; mpfr_mul(r10943, r10908, r10942, MPFR_RNDN); mpfr_add(r10944, r10941, r10943, MPFR_RNDN); mpfr_add(r10945, r10932, r10944, MPFR_RNDN); mpfr_mul(r10946, r10931, r10945, MPFR_RNDN); mpfr_mul(r10947, r10886, r10946, MPFR_RNDN); if (mpfr_get_si(r10892, MPFR_RNDN)) { mpfr_set(r10948, r10904, MPFR_RNDN); } else { mpfr_set(r10948, r10947, MPFR_RNDN); }; return mpfr_get_d(r10948, MPFR_RNDN); } static mpfr_t r10949, r10950, r10951, r10952, r10953, r10954, r10955, r10956, r10957, r10958, r10959, r10960, r10961, r10962, r10963, r10964, r10965, r10966, r10967, r10968, r10969, r10970, r10971, r10972, r10973, r10974, r10975, r10976, r10977, r10978, r10979, r10980, r10981, r10982, r10983, r10984, r10985, r10986, r10987, r10988, r10989, r10990, r10991, r10992, r10993, r10994, r10995, r10996, r10997, r10998, r10999, r11000, r11001, r11002, r11003, r11004, r11005, r11006, r11007, r11008, r11009, r11010, r11011; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init(r10949); mpfr_init_set_str(r10950, "-4456723949048.314", 10, MPFR_RNDN); mpfr_init(r10951); mpfr_init_set_str(r10952, "896.8339123666102", 10, MPFR_RNDN); mpfr_init(r10953); mpfr_init(r10954); mpfr_init(r10955); mpfr_init_set_str(r10956, "0.2514179000665375", 10, MPFR_RNDN); mpfr_init(r10957); mpfr_init(r10958); mpfr_init(r10959); mpfr_init_set_str(r10960, "0.5", 10, MPFR_RNDN); mpfr_init(r10961); mpfr_init_set_str(r10962, "0.15298196345928972", 10, MPFR_RNDN); mpfr_init_set_str(r10963, "5", 10, MPFR_RNDN); mpfr_init(r10964); mpfr_init(r10965); mpfr_init(r10966); mpfr_init(r10967); mpfr_init_set_str(r10968, "1", 10, MPFR_RNDN); mpfr_init(r10969); mpfr_init(r10970); mpfr_init(r10971); mpfr_init(r10972); mpfr_init(r10973); mpfr_init_set_str(r10974, "2", 10, MPFR_RNDN); mpfr_init_set_str(r10975, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r10976); mpfr_init(r10977); mpfr_init_set_str(r10978, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r10979); mpfr_init_set_str(r10980, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r10981); mpfr_init(r10982); mpfr_init_set_str(r10983, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r10984); mpfr_init(r10985); mpfr_init_set_str(r10986, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r10987); mpfr_init(r10988); mpfr_init_set_str(r10989, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r10990); mpfr_init(r10991); mpfr_init(r10992); mpfr_init(r10993); mpfr_init(r10994); mpfr_init(r10995); mpfr_init_set_str(r10996, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r10997); mpfr_init_set_str(r10998, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r10999); mpfr_init(r11000); mpfr_init_set_str(r11001, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r11002); mpfr_init(r11003); mpfr_init(r11004); mpfr_init_set_str(r11005, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r11006); mpfr_init(r11007); mpfr_init(r11008); mpfr_init(r11009); mpfr_init(r11010); mpfr_init(r11011); } double f_dm(double x) { mpfr_set_d(r10949, x, MPFR_RNDN); ; mpfr_set_si(r10951, mpfr_cmp(r10949, r10950) <= 0, MPFR_RNDN); ; mpfr_set_si(r10953, mpfr_cmp(r10949, r10952) <= 0, MPFR_RNDN); mpfr_set_si(r10954, !mpfr_get_si(r10953, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r10955, mpfr_get_si(r10951, MPFR_RNDN) || mpfr_get_si(r10954, MPFR_RNDN), MPFR_RNDN); ; mpfr_div(r10957, r10956, r10949, MPFR_RNDN); mpfr_mul(r10958, r10949, r10949, MPFR_RNDN); mpfr_div(r10959, r10957, r10958, MPFR_RNDN); ; mpfr_div(r10961, r10960, r10949, MPFR_RNDN); ; ; mpfr_pow(r10964, r10949, r10963, MPFR_RNDN); mpfr_div(r10965, r10962, r10964, MPFR_RNDN); mpfr_add(r10966, r10961, r10965, MPFR_RNDN); mpfr_add(r10967, r10959, r10966, MPFR_RNDN); ; mpfr_mul(r10969, r10958, r10958, MPFR_RNDN); mpfr_mul(r10970, r10969, r10958, MPFR_RNDN); mpfr_mul(r10971, r10970, r10958, MPFR_RNDN); mpfr_mul(r10972, r10958, r10971, MPFR_RNDN); mpfr_mul(r10973, r10972, r10958, MPFR_RNDN); ; ; mpfr_mul(r10976, r10974, r10975, MPFR_RNDN); mpfr_mul(r10977, r10973, r10976, MPFR_RNDN); ; mpfr_mul(r10979, r10972, r10978, MPFR_RNDN); ; mpfr_mul(r10981, r10958, r10980, MPFR_RNDN); mpfr_add(r10982, r10981, r10968, MPFR_RNDN); ; mpfr_mul(r10984, r10983, r10969, MPFR_RNDN); mpfr_add(r10985, r10982, r10984, MPFR_RNDN); ; mpfr_mul(r10987, r10970, r10986, MPFR_RNDN); mpfr_add(r10988, r10985, r10987, MPFR_RNDN); ; mpfr_mul(r10990, r10971, r10989, MPFR_RNDN); mpfr_add(r10991, r10988, r10990, MPFR_RNDN); mpfr_add(r10992, r10979, r10991, MPFR_RNDN); mpfr_add(r10993, r10977, r10992, MPFR_RNDN); mpfr_div(r10994, r10968, r10993, MPFR_RNDN); mpfr_mul(r10995, r10972, r10975, MPFR_RNDN); ; mpfr_mul(r10997, r10996, r10970, MPFR_RNDN); ; mpfr_mul(r10999, r10958, r10998, MPFR_RNDN); mpfr_add(r11000, r10968, r10999, MPFR_RNDN); ; mpfr_mul(r11002, r11001, r10969, MPFR_RNDN); mpfr_add(r11003, r11000, r11002, MPFR_RNDN); mpfr_add(r11004, r10997, r11003, MPFR_RNDN); ; mpfr_mul(r11006, r10971, r11005, MPFR_RNDN); mpfr_add(r11007, r11004, r11006, MPFR_RNDN); mpfr_add(r11008, r10995, r11007, MPFR_RNDN); mpfr_mul(r11009, r10994, r11008, MPFR_RNDN); mpfr_mul(r11010, r10949, r11009, MPFR_RNDN); if (mpfr_get_si(r10955, MPFR_RNDN)) { mpfr_set(r11011, r10967, MPFR_RNDN); } else { mpfr_set(r11011, r11010, MPFR_RNDN); }; return mpfr_get_d(r11011, MPFR_RNDN); }