#include #include #include #include #include char *name = "Jmat.Real.dawson"; double f_if(float x) { float r25645 = 1; float r25646 = 0.1049934947; float r25647 = x; float r25648 = r25647 * r25647; float r25649 = r25646 * r25648; float r25650 = r25645 + r25649; float r25651 = 0.0424060604; float r25652 = r25648 * r25648; float r25653 = r25651 * r25652; float r25654 = r25650 + r25653; float r25655 = 0.0072644182; float r25656 = r25652 * r25648; float r25657 = r25655 * r25656; float r25658 = r25654 + r25657; float r25659 = 0.0005064034; float r25660 = r25656 * r25648; float r25661 = r25659 * r25660; float r25662 = r25658 + r25661; float r25663 = 0.0001789971; float r25664 = r25660 * r25648; float r25665 = r25663 * r25664; float r25666 = r25662 + r25665; float r25667 = 0.7715471019; float r25668 = r25667 * r25648; float r25669 = r25645 + r25668; float r25670 = 0.2909738639; float r25671 = r25670 * r25652; float r25672 = r25669 + r25671; float r25673 = 0.0694555761; float r25674 = r25673 * r25656; float r25675 = r25672 + r25674; float r25676 = 0.0140005442; float r25677 = r25676 * r25660; float r25678 = r25675 + r25677; float r25679 = 0.0008327945; float r25680 = r25679 * r25664; float r25681 = r25678 + r25680; float r25682 = 2; float r25683 = r25682 * r25663; float r25684 = r25664 * r25648; float r25685 = r25683 * r25684; float r25686 = r25681 + r25685; float r25687 = r25666 / r25686; float r25688 = r25687 * r25647; return r25688; } double f_id(double x) { double r25689 = 1; double r25690 = 0.1049934947; double r25691 = x; double r25692 = r25691 * r25691; double r25693 = r25690 * r25692; double r25694 = r25689 + r25693; double r25695 = 0.0424060604; double r25696 = r25692 * r25692; double r25697 = r25695 * r25696; double r25698 = r25694 + r25697; double r25699 = 0.0072644182; double r25700 = r25696 * r25692; double r25701 = r25699 * r25700; double r25702 = r25698 + r25701; double r25703 = 0.0005064034; double r25704 = r25700 * r25692; double r25705 = r25703 * r25704; double r25706 = r25702 + r25705; double r25707 = 0.0001789971; double r25708 = r25704 * r25692; double r25709 = r25707 * r25708; double r25710 = r25706 + r25709; double r25711 = 0.7715471019; double r25712 = r25711 * r25692; double r25713 = r25689 + r25712; double r25714 = 0.2909738639; double r25715 = r25714 * r25696; double r25716 = r25713 + r25715; double r25717 = 0.0694555761; double r25718 = r25717 * r25700; double r25719 = r25716 + r25718; double r25720 = 0.0140005442; double r25721 = r25720 * r25704; double r25722 = r25719 + r25721; double r25723 = 0.0008327945; double r25724 = r25723 * r25708; double r25725 = r25722 + r25724; double r25726 = 2; double r25727 = r25726 * r25707; double r25728 = r25708 * r25692; double r25729 = r25727 * r25728; double r25730 = r25725 + r25729; double r25731 = r25710 / r25730; double r25732 = r25731 * r25691; return r25732; } double f_of(float x) { float r25733 = x; float r25734 = -5.161343281950335e+18; bool r25735 = r25733 <= r25734; float r25736 = 0.5; float r25737 = r25736 / r25733; float r25738 = 0.2514179000665375; float r25739 = r25738 * r25733; float r25740 = 4; float r25741 = pow(r25733, r25740); float r25742 = r25739 / r25741; float r25743 = 0.15298196345929327; float r25744 = r25743 * r25733; float r25745 = 6; float r25746 = pow(r25733, r25745); float r25747 = r25744 / r25746; float r25748 = r25742 + r25747; float r25749 = r25737 + r25748; float r25750 = 510048.78143239743; bool r25751 = r25733 <= r25750; float r25752 = 0.1049934947; float r25753 = r25733 * r25752; float r25754 = r25753 * r25733; float r25755 = 1; float r25756 = 0.0424060604; float r25757 = r25756 * r25733; float r25758 = 3; float r25759 = pow(r25733, r25758); float r25760 = r25757 * r25759; float r25761 = r25755 + r25760; float r25762 = r25754 + r25761; float r25763 = r25733 * r25733; float r25764 = r25763 * r25763; float r25765 = 0.0072644182; float r25766 = r25763 * r25765; float r25767 = r25764 * r25766; float r25768 = 0.0005064034; float r25769 = r25768 * r25763; float r25770 = r25759 * r25759; float r25771 = r25769 * r25770; float r25772 = r25767 + r25771; float r25773 = r25764 * r25764; float r25774 = 0.0001789971; float r25775 = r25763 * r25774; float r25776 = r25773 * r25775; float r25777 = r25772 + r25776; float r25778 = r25762 + r25777; float r25779 = 0.0008327945; float r25780 = r25763 * r25779; float r25781 = r25773 * r25780; float r25782 = 2; float r25783 = r25774 * r25782; float r25784 = r25764 * r25783; float r25785 = r25773 * r25784; float r25786 = r25781 + r25785; float r25787 = 0.0140005442; float r25788 = r25787 * r25763; float r25789 = r25788 * r25770; float r25790 = 0.7715471019; float r25791 = r25790 * r25763; float r25792 = r25791 + r25755; float r25793 = r25789 + r25792; float r25794 = 0.2909738639; float r25795 = r25794 * r25763; float r25796 = 0.0694555761; float r25797 = r25796 * r25733; float r25798 = r25797 * r25759; float r25799 = r25795 + r25798; float r25800 = r25763 * r25799; float r25801 = r25793 + r25800; float r25802 = r25786 + r25801; float r25803 = r25778 / r25802; float r25804 = pow(r25803, r25758); float r25805 = cbrt(r25804); float r25806 = r25805 * r25733; float r25807 = r25751 ? r25806 : r25749; float r25808 = r25735 ? r25749 : r25807; return r25808; } double f_od(double x) { double r25809 = x; double r25810 = -5.161343281950335e+18; bool r25811 = r25809 <= r25810; double r25812 = 0.5; double r25813 = r25812 / r25809; double r25814 = 0.2514179000665375; double r25815 = r25814 * r25809; double r25816 = 4; double r25817 = pow(r25809, r25816); double r25818 = r25815 / r25817; double r25819 = 0.15298196345929327; double r25820 = r25819 * r25809; double r25821 = 6; double r25822 = pow(r25809, r25821); double r25823 = r25820 / r25822; double r25824 = r25818 + r25823; double r25825 = r25813 + r25824; double r25826 = 510048.78143239743; bool r25827 = r25809 <= r25826; double r25828 = 0.1049934947; double r25829 = r25809 * r25828; double r25830 = r25829 * r25809; double r25831 = 1; double r25832 = 0.0424060604; double r25833 = r25832 * r25809; double r25834 = 3; double r25835 = pow(r25809, r25834); double r25836 = r25833 * r25835; double r25837 = r25831 + r25836; double r25838 = r25830 + r25837; double r25839 = r25809 * r25809; double r25840 = r25839 * r25839; double r25841 = 0.0072644182; double r25842 = r25839 * r25841; double r25843 = r25840 * r25842; double r25844 = 0.0005064034; double r25845 = r25844 * r25839; double r25846 = r25835 * r25835; double r25847 = r25845 * r25846; double r25848 = r25843 + r25847; double r25849 = r25840 * r25840; double r25850 = 0.0001789971; double r25851 = r25839 * r25850; double r25852 = r25849 * r25851; double r25853 = r25848 + r25852; double r25854 = r25838 + r25853; double r25855 = 0.0008327945; double r25856 = r25839 * r25855; double r25857 = r25849 * r25856; double r25858 = 2; double r25859 = r25850 * r25858; double r25860 = r25840 * r25859; double r25861 = r25849 * r25860; double r25862 = r25857 + r25861; double r25863 = 0.0140005442; double r25864 = r25863 * r25839; double r25865 = r25864 * r25846; double r25866 = 0.7715471019; double r25867 = r25866 * r25839; double r25868 = r25867 + r25831; double r25869 = r25865 + r25868; double r25870 = 0.2909738639; double r25871 = r25870 * r25839; double r25872 = 0.0694555761; double r25873 = r25872 * r25809; double r25874 = r25873 * r25835; double r25875 = r25871 + r25874; double r25876 = r25839 * r25875; double r25877 = r25869 + r25876; double r25878 = r25862 + r25877; double r25879 = r25854 / r25878; double r25880 = pow(r25879, r25834); double r25881 = cbrt(r25880); double r25882 = r25881 * r25809; double r25883 = r25827 ? r25882 : r25825; double r25884 = r25811 ? r25825 : r25883; return r25884; } 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 r25885, r25886, r25887, r25888, r25889, r25890, r25891, r25892, r25893, r25894, r25895, r25896, r25897, r25898, r25899, r25900, r25901, r25902, r25903, r25904, r25905, r25906, r25907, r25908, r25909, r25910, r25911, r25912, r25913, r25914, r25915, r25916, r25917, r25918, r25919, r25920, r25921, r25922, r25923, r25924, r25925, r25926, r25927, r25928; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r25885, "1", 10, MPFR_RNDN); mpfr_init_set_str(r25886, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r25887); mpfr_init(r25888); mpfr_init(r25889); mpfr_init(r25890); mpfr_init_set_str(r25891, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r25892); mpfr_init(r25893); mpfr_init(r25894); mpfr_init_set_str(r25895, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r25896); mpfr_init(r25897); mpfr_init(r25898); mpfr_init_set_str(r25899, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r25900); mpfr_init(r25901); mpfr_init(r25902); mpfr_init_set_str(r25903, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r25904); mpfr_init(r25905); mpfr_init(r25906); mpfr_init_set_str(r25907, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r25908); mpfr_init(r25909); mpfr_init_set_str(r25910, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r25911); mpfr_init(r25912); mpfr_init_set_str(r25913, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r25914); mpfr_init(r25915); mpfr_init_set_str(r25916, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r25917); mpfr_init(r25918); mpfr_init_set_str(r25919, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r25920); mpfr_init(r25921); mpfr_init_set_str(r25922, "2", 10, MPFR_RNDN); mpfr_init(r25923); mpfr_init(r25924); mpfr_init(r25925); mpfr_init(r25926); mpfr_init(r25927); mpfr_init(r25928); } double f_im(double x) { ; ; mpfr_set_d(r25887, x, MPFR_RNDN); mpfr_mul(r25888, r25887, r25887, MPFR_RNDN); mpfr_mul(r25889, r25886, r25888, MPFR_RNDN); mpfr_add(r25890, r25885, r25889, MPFR_RNDN); ; mpfr_mul(r25892, r25888, r25888, MPFR_RNDN); mpfr_mul(r25893, r25891, r25892, MPFR_RNDN); mpfr_add(r25894, r25890, r25893, MPFR_RNDN); ; mpfr_mul(r25896, r25892, r25888, MPFR_RNDN); mpfr_mul(r25897, r25895, r25896, MPFR_RNDN); mpfr_add(r25898, r25894, r25897, MPFR_RNDN); ; mpfr_mul(r25900, r25896, r25888, MPFR_RNDN); mpfr_mul(r25901, r25899, r25900, MPFR_RNDN); mpfr_add(r25902, r25898, r25901, MPFR_RNDN); ; mpfr_mul(r25904, r25900, r25888, MPFR_RNDN); mpfr_mul(r25905, r25903, r25904, MPFR_RNDN); mpfr_add(r25906, r25902, r25905, MPFR_RNDN); ; mpfr_mul(r25908, r25907, r25888, MPFR_RNDN); mpfr_add(r25909, r25885, r25908, MPFR_RNDN); ; mpfr_mul(r25911, r25910, r25892, MPFR_RNDN); mpfr_add(r25912, r25909, r25911, MPFR_RNDN); ; mpfr_mul(r25914, r25913, r25896, MPFR_RNDN); mpfr_add(r25915, r25912, r25914, MPFR_RNDN); ; mpfr_mul(r25917, r25916, r25900, MPFR_RNDN); mpfr_add(r25918, r25915, r25917, MPFR_RNDN); ; mpfr_mul(r25920, r25919, r25904, MPFR_RNDN); mpfr_add(r25921, r25918, r25920, MPFR_RNDN); ; mpfr_mul(r25923, r25922, r25903, MPFR_RNDN); mpfr_mul(r25924, r25904, r25888, MPFR_RNDN); mpfr_mul(r25925, r25923, r25924, MPFR_RNDN); mpfr_add(r25926, r25921, r25925, MPFR_RNDN); mpfr_div(r25927, r25906, r25926, MPFR_RNDN); mpfr_mul(r25928, r25927, r25887, MPFR_RNDN); return mpfr_get_d(r25928, MPFR_RNDN); } static mpfr_t r25929, r25930, r25931, r25932, r25933, r25934, r25935, r25936, r25937, r25938, r25939, r25940, r25941, r25942, r25943, r25944, r25945, r25946, r25947, r25948, r25949, r25950, r25951, r25952, r25953, r25954, r25955, r25956, r25957, r25958, r25959, r25960, r25961, r25962, r25963, r25964, r25965, r25966, r25967, r25968, r25969, r25970, r25971, r25972, r25973, r25974, r25975, r25976, r25977, r25978, r25979, r25980, r25981, r25982, r25983, r25984, r25985, r25986, r25987, r25988, r25989, r25990, r25991, r25992, r25993, r25994, r25995, r25996, r25997, r25998, r25999, r26000, r26001, r26002, r26003, r26004; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r25929); mpfr_init_set_str(r25930, "-5.161343281950335e+18", 10, MPFR_RNDN); mpfr_init(r25931); mpfr_init_set_str(r25932, "0.5", 10, MPFR_RNDN); mpfr_init(r25933); mpfr_init_set_str(r25934, "0.2514179000665375", 10, MPFR_RNDN); mpfr_init(r25935); mpfr_init_set_str(r25936, "4", 10, MPFR_RNDN); mpfr_init(r25937); mpfr_init(r25938); mpfr_init_set_str(r25939, "0.15298196345929327", 10, MPFR_RNDN); mpfr_init(r25940); mpfr_init_set_str(r25941, "6", 10, MPFR_RNDN); mpfr_init(r25942); mpfr_init(r25943); mpfr_init(r25944); mpfr_init(r25945); mpfr_init_set_str(r25946, "510048.78143239743", 10, MPFR_RNDN); mpfr_init(r25947); mpfr_init_set_str(r25948, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r25949); mpfr_init(r25950); mpfr_init_set_str(r25951, "1", 10, MPFR_RNDN); mpfr_init_set_str(r25952, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r25953); mpfr_init_set_str(r25954, "3", 10, MPFR_RNDN); mpfr_init(r25955); mpfr_init(r25956); mpfr_init(r25957); mpfr_init(r25958); mpfr_init(r25959); mpfr_init(r25960); mpfr_init_set_str(r25961, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r25962); mpfr_init(r25963); mpfr_init_set_str(r25964, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r25965); mpfr_init(r25966); mpfr_init(r25967); mpfr_init(r25968); mpfr_init(r25969); mpfr_init_set_str(r25970, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r25971); mpfr_init(r25972); mpfr_init(r25973); mpfr_init(r25974); mpfr_init_set_str(r25975, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r25976); mpfr_init(r25977); mpfr_init_set_str(r25978, "2", 10, MPFR_RNDN); mpfr_init(r25979); mpfr_init(r25980); mpfr_init(r25981); mpfr_init(r25982); mpfr_init_set_str(r25983, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r25984); mpfr_init(r25985); mpfr_init_set_str(r25986, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r25987); mpfr_init(r25988); mpfr_init(r25989); mpfr_init_set_str(r25990, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r25991); mpfr_init_set_str(r25992, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r25993); mpfr_init(r25994); mpfr_init(r25995); mpfr_init(r25996); mpfr_init(r25997); mpfr_init(r25998); mpfr_init(r25999); mpfr_init(r26000); mpfr_init(r26001); mpfr_init(r26002); mpfr_init(r26003); mpfr_init(r26004); } double f_fm(double x) { mpfr_set_d(r25929, x, MPFR_RNDN); ; mpfr_set_si(r25931, mpfr_cmp(r25929, r25930) <= 0, MPFR_RNDN); ; mpfr_div(r25933, r25932, r25929, MPFR_RNDN); ; mpfr_mul(r25935, r25934, r25929, MPFR_RNDN); ; mpfr_pow(r25937, r25929, r25936, MPFR_RNDN); mpfr_div(r25938, r25935, r25937, MPFR_RNDN); ; mpfr_mul(r25940, r25939, r25929, MPFR_RNDN); ; mpfr_pow(r25942, r25929, r25941, MPFR_RNDN); mpfr_div(r25943, r25940, r25942, MPFR_RNDN); mpfr_add(r25944, r25938, r25943, MPFR_RNDN); mpfr_add(r25945, r25933, r25944, MPFR_RNDN); ; mpfr_set_si(r25947, mpfr_cmp(r25929, r25946) <= 0, MPFR_RNDN); ; mpfr_mul(r25949, r25929, r25948, MPFR_RNDN); mpfr_mul(r25950, r25949, r25929, MPFR_RNDN); ; ; mpfr_mul(r25953, r25952, r25929, MPFR_RNDN); ; mpfr_pow(r25955, r25929, r25954, MPFR_RNDN); mpfr_mul(r25956, r25953, r25955, MPFR_RNDN); mpfr_add(r25957, r25951, r25956, MPFR_RNDN); mpfr_add(r25958, r25950, r25957, MPFR_RNDN); mpfr_mul(r25959, r25929, r25929, MPFR_RNDN); mpfr_mul(r25960, r25959, r25959, MPFR_RNDN); ; mpfr_mul(r25962, r25959, r25961, MPFR_RNDN); mpfr_mul(r25963, r25960, r25962, MPFR_RNDN); ; mpfr_mul(r25965, r25964, r25959, MPFR_RNDN); mpfr_mul(r25966, r25955, r25955, MPFR_RNDN); mpfr_mul(r25967, r25965, r25966, MPFR_RNDN); mpfr_add(r25968, r25963, r25967, MPFR_RNDN); mpfr_mul(r25969, r25960, r25960, MPFR_RNDN); ; mpfr_mul(r25971, r25959, r25970, MPFR_RNDN); mpfr_mul(r25972, r25969, r25971, MPFR_RNDN); mpfr_add(r25973, r25968, r25972, MPFR_RNDN); mpfr_add(r25974, r25958, r25973, MPFR_RNDN); ; mpfr_mul(r25976, r25959, r25975, MPFR_RNDN); mpfr_mul(r25977, r25969, r25976, MPFR_RNDN); ; mpfr_mul(r25979, r25970, r25978, MPFR_RNDN); mpfr_mul(r25980, r25960, r25979, MPFR_RNDN); mpfr_mul(r25981, r25969, r25980, MPFR_RNDN); mpfr_add(r25982, r25977, r25981, MPFR_RNDN); ; mpfr_mul(r25984, r25983, r25959, MPFR_RNDN); mpfr_mul(r25985, r25984, r25966, MPFR_RNDN); ; mpfr_mul(r25987, r25986, r25959, MPFR_RNDN); mpfr_add(r25988, r25987, r25951, MPFR_RNDN); mpfr_add(r25989, r25985, r25988, MPFR_RNDN); ; mpfr_mul(r25991, r25990, r25959, MPFR_RNDN); ; mpfr_mul(r25993, r25992, r25929, MPFR_RNDN); mpfr_mul(r25994, r25993, r25955, MPFR_RNDN); mpfr_add(r25995, r25991, r25994, MPFR_RNDN); mpfr_mul(r25996, r25959, r25995, MPFR_RNDN); mpfr_add(r25997, r25989, r25996, MPFR_RNDN); mpfr_add(r25998, r25982, r25997, MPFR_RNDN); mpfr_div(r25999, r25974, r25998, MPFR_RNDN); mpfr_pow(r26000, r25999, r25954, MPFR_RNDN); mpfr_cbrt(r26001, r26000, MPFR_RNDN); mpfr_mul(r26002, r26001, r25929, MPFR_RNDN); if (mpfr_get_si(r25947, MPFR_RNDN)) { mpfr_set(r26003, r26002, MPFR_RNDN); } else { mpfr_set(r26003, r25945, MPFR_RNDN); }; if (mpfr_get_si(r25931, MPFR_RNDN)) { mpfr_set(r26004, r25945, MPFR_RNDN); } else { mpfr_set(r26004, r26003, MPFR_RNDN); }; return mpfr_get_d(r26004, MPFR_RNDN); } static mpfr_t r26005, r26006, r26007, r26008, r26009, r26010, r26011, r26012, r26013, r26014, r26015, r26016, r26017, r26018, r26019, r26020, r26021, r26022, r26023, r26024, r26025, r26026, r26027, r26028, r26029, r26030, r26031, r26032, r26033, r26034, r26035, r26036, r26037, r26038, r26039, r26040, r26041, r26042, r26043, r26044, r26045, r26046, r26047, r26048, r26049, r26050, r26051, r26052, r26053, r26054, r26055, r26056, r26057, r26058, r26059, r26060, r26061, r26062, r26063, r26064, r26065, r26066, r26067, r26068, r26069, r26070, r26071, r26072, r26073, r26074, r26075, r26076, r26077, r26078, r26079, r26080; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r26005); mpfr_init_set_str(r26006, "-5.161343281950335e+18", 10, MPFR_RNDN); mpfr_init(r26007); mpfr_init_set_str(r26008, "0.5", 10, MPFR_RNDN); mpfr_init(r26009); mpfr_init_set_str(r26010, "0.2514179000665375", 10, MPFR_RNDN); mpfr_init(r26011); mpfr_init_set_str(r26012, "4", 10, MPFR_RNDN); mpfr_init(r26013); mpfr_init(r26014); mpfr_init_set_str(r26015, "0.15298196345929327", 10, MPFR_RNDN); mpfr_init(r26016); mpfr_init_set_str(r26017, "6", 10, MPFR_RNDN); mpfr_init(r26018); mpfr_init(r26019); mpfr_init(r26020); mpfr_init(r26021); mpfr_init_set_str(r26022, "510048.78143239743", 10, MPFR_RNDN); mpfr_init(r26023); mpfr_init_set_str(r26024, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r26025); mpfr_init(r26026); mpfr_init_set_str(r26027, "1", 10, MPFR_RNDN); mpfr_init_set_str(r26028, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r26029); mpfr_init_set_str(r26030, "3", 10, MPFR_RNDN); mpfr_init(r26031); mpfr_init(r26032); mpfr_init(r26033); mpfr_init(r26034); mpfr_init(r26035); mpfr_init(r26036); mpfr_init_set_str(r26037, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r26038); mpfr_init(r26039); mpfr_init_set_str(r26040, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r26041); mpfr_init(r26042); mpfr_init(r26043); mpfr_init(r26044); mpfr_init(r26045); mpfr_init_set_str(r26046, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r26047); mpfr_init(r26048); mpfr_init(r26049); mpfr_init(r26050); mpfr_init_set_str(r26051, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r26052); mpfr_init(r26053); mpfr_init_set_str(r26054, "2", 10, MPFR_RNDN); mpfr_init(r26055); mpfr_init(r26056); mpfr_init(r26057); mpfr_init(r26058); mpfr_init_set_str(r26059, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r26060); mpfr_init(r26061); mpfr_init_set_str(r26062, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r26063); mpfr_init(r26064); mpfr_init(r26065); mpfr_init_set_str(r26066, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r26067); mpfr_init_set_str(r26068, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r26069); mpfr_init(r26070); mpfr_init(r26071); mpfr_init(r26072); mpfr_init(r26073); mpfr_init(r26074); mpfr_init(r26075); mpfr_init(r26076); mpfr_init(r26077); mpfr_init(r26078); mpfr_init(r26079); mpfr_init(r26080); } double f_dm(double x) { mpfr_set_d(r26005, x, MPFR_RNDN); ; mpfr_set_si(r26007, mpfr_cmp(r26005, r26006) <= 0, MPFR_RNDN); ; mpfr_div(r26009, r26008, r26005, MPFR_RNDN); ; mpfr_mul(r26011, r26010, r26005, MPFR_RNDN); ; mpfr_pow(r26013, r26005, r26012, MPFR_RNDN); mpfr_div(r26014, r26011, r26013, MPFR_RNDN); ; mpfr_mul(r26016, r26015, r26005, MPFR_RNDN); ; mpfr_pow(r26018, r26005, r26017, MPFR_RNDN); mpfr_div(r26019, r26016, r26018, MPFR_RNDN); mpfr_add(r26020, r26014, r26019, MPFR_RNDN); mpfr_add(r26021, r26009, r26020, MPFR_RNDN); ; mpfr_set_si(r26023, mpfr_cmp(r26005, r26022) <= 0, MPFR_RNDN); ; mpfr_mul(r26025, r26005, r26024, MPFR_RNDN); mpfr_mul(r26026, r26025, r26005, MPFR_RNDN); ; ; mpfr_mul(r26029, r26028, r26005, MPFR_RNDN); ; mpfr_pow(r26031, r26005, r26030, MPFR_RNDN); mpfr_mul(r26032, r26029, r26031, MPFR_RNDN); mpfr_add(r26033, r26027, r26032, MPFR_RNDN); mpfr_add(r26034, r26026, r26033, MPFR_RNDN); mpfr_mul(r26035, r26005, r26005, MPFR_RNDN); mpfr_mul(r26036, r26035, r26035, MPFR_RNDN); ; mpfr_mul(r26038, r26035, r26037, MPFR_RNDN); mpfr_mul(r26039, r26036, r26038, MPFR_RNDN); ; mpfr_mul(r26041, r26040, r26035, MPFR_RNDN); mpfr_mul(r26042, r26031, r26031, MPFR_RNDN); mpfr_mul(r26043, r26041, r26042, MPFR_RNDN); mpfr_add(r26044, r26039, r26043, MPFR_RNDN); mpfr_mul(r26045, r26036, r26036, MPFR_RNDN); ; mpfr_mul(r26047, r26035, r26046, MPFR_RNDN); mpfr_mul(r26048, r26045, r26047, MPFR_RNDN); mpfr_add(r26049, r26044, r26048, MPFR_RNDN); mpfr_add(r26050, r26034, r26049, MPFR_RNDN); ; mpfr_mul(r26052, r26035, r26051, MPFR_RNDN); mpfr_mul(r26053, r26045, r26052, MPFR_RNDN); ; mpfr_mul(r26055, r26046, r26054, MPFR_RNDN); mpfr_mul(r26056, r26036, r26055, MPFR_RNDN); mpfr_mul(r26057, r26045, r26056, MPFR_RNDN); mpfr_add(r26058, r26053, r26057, MPFR_RNDN); ; mpfr_mul(r26060, r26059, r26035, MPFR_RNDN); mpfr_mul(r26061, r26060, r26042, MPFR_RNDN); ; mpfr_mul(r26063, r26062, r26035, MPFR_RNDN); mpfr_add(r26064, r26063, r26027, MPFR_RNDN); mpfr_add(r26065, r26061, r26064, MPFR_RNDN); ; mpfr_mul(r26067, r26066, r26035, MPFR_RNDN); ; mpfr_mul(r26069, r26068, r26005, MPFR_RNDN); mpfr_mul(r26070, r26069, r26031, MPFR_RNDN); mpfr_add(r26071, r26067, r26070, MPFR_RNDN); mpfr_mul(r26072, r26035, r26071, MPFR_RNDN); mpfr_add(r26073, r26065, r26072, MPFR_RNDN); mpfr_add(r26074, r26058, r26073, MPFR_RNDN); mpfr_div(r26075, r26050, r26074, MPFR_RNDN); mpfr_pow(r26076, r26075, r26030, MPFR_RNDN); mpfr_cbrt(r26077, r26076, MPFR_RNDN); mpfr_mul(r26078, r26077, r26005, MPFR_RNDN); if (mpfr_get_si(r26023, MPFR_RNDN)) { mpfr_set(r26079, r26078, MPFR_RNDN); } else { mpfr_set(r26079, r26021, MPFR_RNDN); }; if (mpfr_get_si(r26007, MPFR_RNDN)) { mpfr_set(r26080, r26021, MPFR_RNDN); } else { mpfr_set(r26080, r26079, MPFR_RNDN); }; return mpfr_get_d(r26080, MPFR_RNDN); }