#include #include #include #include #include char *name = "Linear.Matrix:det44 from linear-1.19.1.3"; double f_if(float x, float y, float z, float t, float a, float b, float c, float i, float j, float k, float y0, float y1, float y2, float y3, float y4, float y5) { float r25840 = x; float r25841 = y; float r25842 = r25840 * r25841; float r25843 = z; float r25844 = t; float r25845 = r25843 * r25844; float r25846 = r25842 - r25845; float r25847 = a; float r25848 = b; float r25849 = r25847 * r25848; float r25850 = c; float r25851 = i; float r25852 = r25850 * r25851; float r25853 = r25849 - r25852; float r25854 = r25846 * r25853; float r25855 = j; float r25856 = r25840 * r25855; float r25857 = k; float r25858 = r25843 * r25857; float r25859 = r25856 - r25858; float r25860 = y0; float r25861 = r25860 * r25848; float r25862 = y1; float r25863 = r25862 * r25851; float r25864 = r25861 - r25863; float r25865 = r25859 * r25864; float r25866 = r25854 - r25865; float r25867 = y2; float r25868 = r25840 * r25867; float r25869 = y3; float r25870 = r25843 * r25869; float r25871 = r25868 - r25870; float r25872 = r25860 * r25850; float r25873 = r25862 * r25847; float r25874 = r25872 - r25873; float r25875 = r25871 * r25874; float r25876 = r25866 + r25875; float r25877 = r25844 * r25855; float r25878 = r25841 * r25857; float r25879 = r25877 - r25878; float r25880 = y4; float r25881 = r25880 * r25848; float r25882 = y5; float r25883 = r25882 * r25851; float r25884 = r25881 - r25883; float r25885 = r25879 * r25884; float r25886 = r25876 + r25885; float r25887 = r25844 * r25867; float r25888 = r25841 * r25869; float r25889 = r25887 - r25888; float r25890 = r25880 * r25850; float r25891 = r25882 * r25847; float r25892 = r25890 - r25891; float r25893 = r25889 * r25892; float r25894 = r25886 - r25893; float r25895 = r25857 * r25867; float r25896 = r25855 * r25869; float r25897 = r25895 - r25896; float r25898 = r25880 * r25862; float r25899 = r25882 * r25860; float r25900 = r25898 - r25899; float r25901 = r25897 * r25900; float r25902 = r25894 + r25901; return r25902; } double f_id(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k, double y0, double y1, double y2, double y3, double y4, double y5) { double r25903 = x; double r25904 = y; double r25905 = r25903 * r25904; double r25906 = z; double r25907 = t; double r25908 = r25906 * r25907; double r25909 = r25905 - r25908; double r25910 = a; double r25911 = b; double r25912 = r25910 * r25911; double r25913 = c; double r25914 = i; double r25915 = r25913 * r25914; double r25916 = r25912 - r25915; double r25917 = r25909 * r25916; double r25918 = j; double r25919 = r25903 * r25918; double r25920 = k; double r25921 = r25906 * r25920; double r25922 = r25919 - r25921; double r25923 = y0; double r25924 = r25923 * r25911; double r25925 = y1; double r25926 = r25925 * r25914; double r25927 = r25924 - r25926; double r25928 = r25922 * r25927; double r25929 = r25917 - r25928; double r25930 = y2; double r25931 = r25903 * r25930; double r25932 = y3; double r25933 = r25906 * r25932; double r25934 = r25931 - r25933; double r25935 = r25923 * r25913; double r25936 = r25925 * r25910; double r25937 = r25935 - r25936; double r25938 = r25934 * r25937; double r25939 = r25929 + r25938; double r25940 = r25907 * r25918; double r25941 = r25904 * r25920; double r25942 = r25940 - r25941; double r25943 = y4; double r25944 = r25943 * r25911; double r25945 = y5; double r25946 = r25945 * r25914; double r25947 = r25944 - r25946; double r25948 = r25942 * r25947; double r25949 = r25939 + r25948; double r25950 = r25907 * r25930; double r25951 = r25904 * r25932; double r25952 = r25950 - r25951; double r25953 = r25943 * r25913; double r25954 = r25945 * r25910; double r25955 = r25953 - r25954; double r25956 = r25952 * r25955; double r25957 = r25949 - r25956; double r25958 = r25920 * r25930; double r25959 = r25918 * r25932; double r25960 = r25958 - r25959; double r25961 = r25943 * r25925; double r25962 = r25945 * r25923; double r25963 = r25961 - r25962; double r25964 = r25960 * r25963; double r25965 = r25957 + r25964; return r25965; } double f_of(float x, float y, float z, float t, float a, float b, float c, float i, float j, float k, float y0, float y1, float y2, float y3, float y4, float y5) { float r25966 = x; float r25967 = y; float r25968 = r25966 * r25967; float r25969 = z; float r25970 = t; float r25971 = r25969 * r25970; float r25972 = r25968 - r25971; float r25973 = a; float r25974 = b; float r25975 = r25973 * r25974; float r25976 = c; float r25977 = i; float r25978 = r25976 * r25977; float r25979 = r25975 - r25978; float r25980 = r25972 * r25979; float r25981 = j; float r25982 = r25966 * r25981; float r25983 = k; float r25984 = r25969 * r25983; float r25985 = r25982 - r25984; float r25986 = y0; float r25987 = r25986 * r25974; float r25988 = y1; float r25989 = r25988 * r25977; float r25990 = r25987 - r25989; float r25991 = r25985 * r25990; float r25992 = r25980 - r25991; float r25993 = y2; float r25994 = r25966 * r25993; float r25995 = y3; float r25996 = r25969 * r25995; float r25997 = r25994 - r25996; float r25998 = r25986 * r25976; float r25999 = r25988 * r25973; float r26000 = r25998 - r25999; float r26001 = r25997 * r26000; float r26002 = r25992 + r26001; float r26003 = y5; float r26004 = r25983 * r25977; float r26005 = r26003 * r26004; float r26006 = r25967 * r26005; float r26007 = y4; float r26008 = r25974 * r25983; float r26009 = r25967 * r26008; float r26010 = r26007 * r26009; float r26011 = r25970 * r25977; float r26012 = r26003 * r26011; float r26013 = r25981 * r26012; float r26014 = r26010 + r26013; float r26015 = r26006 - r26014; float r26016 = r26002 + r26015; float r26017 = r25970 * r25993; float r26018 = r25967 * r25995; float r26019 = r26017 - r26018; float r26020 = r26007 * r25976; float r26021 = r26003 * r25973; float r26022 = r26020 - r26021; float r26023 = r26019 * r26022; float r26024 = r26016 - r26023; float r26025 = r25983 * r25993; float r26026 = r25981 * r25995; float r26027 = r26025 - r26026; float r26028 = r26007 * r25988; float r26029 = r26003 * r25986; float r26030 = r26028 - r26029; float r26031 = r26027 * r26030; float r26032 = r26024 + r26031; float r26033 = -2.982399978237135e+306; bool r26034 = r26032 <= r26033; float r26035 = r25988 * r26004; float r26036 = r25969 * r26035; float r26037 = r25986 * r25983; float r26038 = r25974 * r26037; float r26039 = r25969 * r26038; float r26040 = r25977 * r25966; float r26041 = r25981 * r26040; float r26042 = r25988 * r26041; float r26043 = r26039 + r26042; float r26044 = r26036 - r26043; float r26045 = r25980 - r26044; float r26046 = r26045 + r26001; float r26047 = r25970 * r25981; float r26048 = r25967 * r25983; float r26049 = r26047 - r26048; float r26050 = r26007 * r25974; float r26051 = r26003 * r25977; float r26052 = r26050 - r26051; float r26053 = r26049 * r26052; float r26054 = r26046 + r26053; float r26055 = r26054 - r26023; float r26056 = r26055 + r26031; float r26057 = -5.977989999399742e+190; bool r26058 = r26032 <= r26057; float r26059 = 5.20966244590318e+293; bool r26060 = r26032 <= r26059; float r26061 = r26002 + r26053; float r26062 = cbrt(r26023); float r26063 = r26062 * r26062; float r26064 = r26063 * r26062; float r26065 = r26061 - r26064; float r26066 = r26065 + r26031; float r26067 = r25993 * r25966; float r26068 = r26067 - r25996; float r26069 = r26068 * r26000; float r26070 = r25981 * r25970; float r26071 = r25983 * r25967; float r26072 = r26070 - r26071; float r26073 = r26072 * r26052; float r26074 = r25993 * r25970; float r26075 = r26074 - r26018; float r26076 = r26022 * r26075; float r26077 = r26073 - r26076; float r26078 = r25974 * r25973; float r26079 = r25977 * r25976; float r26080 = r26078 - r26079; float r26081 = r25967 * r25966; float r26082 = r26081 - r25971; float r26083 = r26080 * r26082; float r26084 = r25981 * r25966; float r26085 = r26084 - r25984; float r26086 = r25974 * r25986; float r26087 = r25977 * r25988; float r26088 = r26086 - r26087; float r26089 = r26085 * r26088; float r26090 = r26083 - r26089; float r26091 = r26077 + r26090; float r26092 = r26069 + r26091; float r26093 = r26060 ? r26066 : r26092; float r26094 = r26058 ? r26032 : r26093; float r26095 = r26034 ? r26056 : r26094; return r26095; } double f_od(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k, double y0, double y1, double y2, double y3, double y4, double y5) { double r26096 = x; double r26097 = y; double r26098 = r26096 * r26097; double r26099 = z; double r26100 = t; double r26101 = r26099 * r26100; double r26102 = r26098 - r26101; double r26103 = a; double r26104 = b; double r26105 = r26103 * r26104; double r26106 = c; double r26107 = i; double r26108 = r26106 * r26107; double r26109 = r26105 - r26108; double r26110 = r26102 * r26109; double r26111 = j; double r26112 = r26096 * r26111; double r26113 = k; double r26114 = r26099 * r26113; double r26115 = r26112 - r26114; double r26116 = y0; double r26117 = r26116 * r26104; double r26118 = y1; double r26119 = r26118 * r26107; double r26120 = r26117 - r26119; double r26121 = r26115 * r26120; double r26122 = r26110 - r26121; double r26123 = y2; double r26124 = r26096 * r26123; double r26125 = y3; double r26126 = r26099 * r26125; double r26127 = r26124 - r26126; double r26128 = r26116 * r26106; double r26129 = r26118 * r26103; double r26130 = r26128 - r26129; double r26131 = r26127 * r26130; double r26132 = r26122 + r26131; double r26133 = y5; double r26134 = r26113 * r26107; double r26135 = r26133 * r26134; double r26136 = r26097 * r26135; double r26137 = y4; double r26138 = r26104 * r26113; double r26139 = r26097 * r26138; double r26140 = r26137 * r26139; double r26141 = r26100 * r26107; double r26142 = r26133 * r26141; double r26143 = r26111 * r26142; double r26144 = r26140 + r26143; double r26145 = r26136 - r26144; double r26146 = r26132 + r26145; double r26147 = r26100 * r26123; double r26148 = r26097 * r26125; double r26149 = r26147 - r26148; double r26150 = r26137 * r26106; double r26151 = r26133 * r26103; double r26152 = r26150 - r26151; double r26153 = r26149 * r26152; double r26154 = r26146 - r26153; double r26155 = r26113 * r26123; double r26156 = r26111 * r26125; double r26157 = r26155 - r26156; double r26158 = r26137 * r26118; double r26159 = r26133 * r26116; double r26160 = r26158 - r26159; double r26161 = r26157 * r26160; double r26162 = r26154 + r26161; double r26163 = -2.982399978237135e+306; bool r26164 = r26162 <= r26163; double r26165 = r26118 * r26134; double r26166 = r26099 * r26165; double r26167 = r26116 * r26113; double r26168 = r26104 * r26167; double r26169 = r26099 * r26168; double r26170 = r26107 * r26096; double r26171 = r26111 * r26170; double r26172 = r26118 * r26171; double r26173 = r26169 + r26172; double r26174 = r26166 - r26173; double r26175 = r26110 - r26174; double r26176 = r26175 + r26131; double r26177 = r26100 * r26111; double r26178 = r26097 * r26113; double r26179 = r26177 - r26178; double r26180 = r26137 * r26104; double r26181 = r26133 * r26107; double r26182 = r26180 - r26181; double r26183 = r26179 * r26182; double r26184 = r26176 + r26183; double r26185 = r26184 - r26153; double r26186 = r26185 + r26161; double r26187 = -5.977989999399742e+190; bool r26188 = r26162 <= r26187; double r26189 = 5.20966244590318e+293; bool r26190 = r26162 <= r26189; double r26191 = r26132 + r26183; double r26192 = cbrt(r26153); double r26193 = r26192 * r26192; double r26194 = r26193 * r26192; double r26195 = r26191 - r26194; double r26196 = r26195 + r26161; double r26197 = r26123 * r26096; double r26198 = r26197 - r26126; double r26199 = r26198 * r26130; double r26200 = r26111 * r26100; double r26201 = r26113 * r26097; double r26202 = r26200 - r26201; double r26203 = r26202 * r26182; double r26204 = r26123 * r26100; double r26205 = r26204 - r26148; double r26206 = r26152 * r26205; double r26207 = r26203 - r26206; double r26208 = r26104 * r26103; double r26209 = r26107 * r26106; double r26210 = r26208 - r26209; double r26211 = r26097 * r26096; double r26212 = r26211 - r26101; double r26213 = r26210 * r26212; double r26214 = r26111 * r26096; double r26215 = r26214 - r26114; double r26216 = r26104 * r26116; double r26217 = r26107 * r26118; double r26218 = r26216 - r26217; double r26219 = r26215 * r26218; double r26220 = r26213 - r26219; double r26221 = r26207 + r26220; double r26222 = r26199 + r26221; double r26223 = r26190 ? r26196 : r26222; double r26224 = r26188 ? r26162 : r26223; double r26225 = r26164 ? r26186 : r26224; return r26225; } 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 r26226, r26227, r26228, r26229, r26230, r26231, r26232, r26233, r26234, r26235, r26236, r26237, r26238, r26239, r26240, r26241, r26242, r26243, r26244, r26245, r26246, r26247, r26248, r26249, r26250, r26251, r26252, r26253, r26254, r26255, r26256, r26257, r26258, r26259, r26260, r26261, r26262, r26263, r26264, r26265, r26266, r26267, r26268, r26269, r26270, r26271, r26272, r26273, r26274, r26275, r26276, r26277, r26278, r26279, r26280, r26281, r26282, r26283, r26284, r26285, r26286, r26287, r26288; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r26226); mpfr_init(r26227); mpfr_init(r26228); mpfr_init(r26229); mpfr_init(r26230); mpfr_init(r26231); mpfr_init(r26232); mpfr_init(r26233); mpfr_init(r26234); mpfr_init(r26235); mpfr_init(r26236); mpfr_init(r26237); mpfr_init(r26238); mpfr_init(r26239); mpfr_init(r26240); mpfr_init(r26241); mpfr_init(r26242); mpfr_init(r26243); mpfr_init(r26244); mpfr_init(r26245); mpfr_init(r26246); mpfr_init(r26247); mpfr_init(r26248); mpfr_init(r26249); mpfr_init(r26250); mpfr_init(r26251); mpfr_init(r26252); mpfr_init(r26253); mpfr_init(r26254); mpfr_init(r26255); mpfr_init(r26256); mpfr_init(r26257); mpfr_init(r26258); mpfr_init(r26259); mpfr_init(r26260); mpfr_init(r26261); mpfr_init(r26262); mpfr_init(r26263); mpfr_init(r26264); mpfr_init(r26265); mpfr_init(r26266); mpfr_init(r26267); mpfr_init(r26268); mpfr_init(r26269); mpfr_init(r26270); mpfr_init(r26271); mpfr_init(r26272); mpfr_init(r26273); mpfr_init(r26274); mpfr_init(r26275); mpfr_init(r26276); mpfr_init(r26277); mpfr_init(r26278); mpfr_init(r26279); mpfr_init(r26280); mpfr_init(r26281); mpfr_init(r26282); mpfr_init(r26283); mpfr_init(r26284); mpfr_init(r26285); mpfr_init(r26286); mpfr_init(r26287); mpfr_init(r26288); } double f_im(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k, double y0, double y1, double y2, double y3, double y4, double y5) { mpfr_set_d(r26226, x, MPFR_RNDN); mpfr_set_d(r26227, y, MPFR_RNDN); mpfr_mul(r26228, r26226, r26227, MPFR_RNDN); mpfr_set_d(r26229, z, MPFR_RNDN); mpfr_set_d(r26230, t, MPFR_RNDN); mpfr_mul(r26231, r26229, r26230, MPFR_RNDN); mpfr_sub(r26232, r26228, r26231, MPFR_RNDN); mpfr_set_d(r26233, a, MPFR_RNDN); mpfr_set_d(r26234, b, MPFR_RNDN); mpfr_mul(r26235, r26233, r26234, MPFR_RNDN); mpfr_set_d(r26236, c, MPFR_RNDN); mpfr_set_d(r26237, i, MPFR_RNDN); mpfr_mul(r26238, r26236, r26237, MPFR_RNDN); mpfr_sub(r26239, r26235, r26238, MPFR_RNDN); mpfr_mul(r26240, r26232, r26239, MPFR_RNDN); mpfr_set_d(r26241, j, MPFR_RNDN); mpfr_mul(r26242, r26226, r26241, MPFR_RNDN); mpfr_set_d(r26243, k, MPFR_RNDN); mpfr_mul(r26244, r26229, r26243, MPFR_RNDN); mpfr_sub(r26245, r26242, r26244, MPFR_RNDN); mpfr_set_d(r26246, y0, MPFR_RNDN); mpfr_mul(r26247, r26246, r26234, MPFR_RNDN); mpfr_set_d(r26248, y1, MPFR_RNDN); mpfr_mul(r26249, r26248, r26237, MPFR_RNDN); mpfr_sub(r26250, r26247, r26249, MPFR_RNDN); mpfr_mul(r26251, r26245, r26250, MPFR_RNDN); mpfr_sub(r26252, r26240, r26251, MPFR_RNDN); mpfr_set_d(r26253, y2, MPFR_RNDN); mpfr_mul(r26254, r26226, r26253, MPFR_RNDN); mpfr_set_d(r26255, y3, MPFR_RNDN); mpfr_mul(r26256, r26229, r26255, MPFR_RNDN); mpfr_sub(r26257, r26254, r26256, MPFR_RNDN); mpfr_mul(r26258, r26246, r26236, MPFR_RNDN); mpfr_mul(r26259, r26248, r26233, MPFR_RNDN); mpfr_sub(r26260, r26258, r26259, MPFR_RNDN); mpfr_mul(r26261, r26257, r26260, MPFR_RNDN); mpfr_add(r26262, r26252, r26261, MPFR_RNDN); mpfr_mul(r26263, r26230, r26241, MPFR_RNDN); mpfr_mul(r26264, r26227, r26243, MPFR_RNDN); mpfr_sub(r26265, r26263, r26264, MPFR_RNDN); mpfr_set_d(r26266, y4, MPFR_RNDN); mpfr_mul(r26267, r26266, r26234, MPFR_RNDN); mpfr_set_d(r26268, y5, MPFR_RNDN); mpfr_mul(r26269, r26268, r26237, MPFR_RNDN); mpfr_sub(r26270, r26267, r26269, MPFR_RNDN); mpfr_mul(r26271, r26265, r26270, MPFR_RNDN); mpfr_add(r26272, r26262, r26271, MPFR_RNDN); mpfr_mul(r26273, r26230, r26253, MPFR_RNDN); mpfr_mul(r26274, r26227, r26255, MPFR_RNDN); mpfr_sub(r26275, r26273, r26274, MPFR_RNDN); mpfr_mul(r26276, r26266, r26236, MPFR_RNDN); mpfr_mul(r26277, r26268, r26233, MPFR_RNDN); mpfr_sub(r26278, r26276, r26277, MPFR_RNDN); mpfr_mul(r26279, r26275, r26278, MPFR_RNDN); mpfr_sub(r26280, r26272, r26279, MPFR_RNDN); mpfr_mul(r26281, r26243, r26253, MPFR_RNDN); mpfr_mul(r26282, r26241, r26255, MPFR_RNDN); mpfr_sub(r26283, r26281, r26282, MPFR_RNDN); mpfr_mul(r26284, r26266, r26248, MPFR_RNDN); mpfr_mul(r26285, r26268, r26246, MPFR_RNDN); mpfr_sub(r26286, r26284, r26285, MPFR_RNDN); mpfr_mul(r26287, r26283, r26286, MPFR_RNDN); mpfr_add(r26288, r26280, r26287, MPFR_RNDN); return mpfr_get_d(r26288, MPFR_RNDN); } static mpfr_t r26289, r26290, r26291, r26292, r26293, r26294, r26295, r26296, r26297, r26298, r26299, r26300, r26301, r26302, r26303, r26304, r26305, r26306, r26307, r26308, r26309, r26310, r26311, r26312, r26313, r26314, r26315, r26316, r26317, r26318, r26319, r26320, r26321, r26322, r26323, r26324, r26325, r26326, r26327, r26328, r26329, r26330, r26331, r26332, r26333, r26334, r26335, r26336, r26337, r26338, r26339, r26340, r26341, r26342, r26343, r26344, r26345, r26346, r26347, r26348, r26349, r26350, r26351, r26352, r26353, r26354, r26355, r26356, r26357, r26358, r26359, r26360, r26361, r26362, r26363, r26364, r26365, r26366, r26367, r26368, r26369, r26370, r26371, r26372, r26373, r26374, r26375, r26376, r26377, r26378, r26379, r26380, r26381, r26382, r26383, r26384, r26385, r26386, r26387, r26388, r26389, r26390, r26391, r26392, r26393, r26394, r26395, r26396, r26397, r26398, r26399, r26400, r26401, r26402, r26403, r26404, r26405, r26406, r26407, r26408, r26409, r26410, r26411, r26412, r26413, r26414, r26415, r26416, r26417, r26418; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r26289); mpfr_init(r26290); mpfr_init(r26291); mpfr_init(r26292); mpfr_init(r26293); mpfr_init(r26294); mpfr_init(r26295); mpfr_init(r26296); mpfr_init(r26297); mpfr_init(r26298); mpfr_init(r26299); mpfr_init(r26300); mpfr_init(r26301); mpfr_init(r26302); mpfr_init(r26303); mpfr_init(r26304); mpfr_init(r26305); mpfr_init(r26306); mpfr_init(r26307); mpfr_init(r26308); mpfr_init(r26309); mpfr_init(r26310); mpfr_init(r26311); mpfr_init(r26312); mpfr_init(r26313); mpfr_init(r26314); mpfr_init(r26315); mpfr_init(r26316); mpfr_init(r26317); mpfr_init(r26318); mpfr_init(r26319); mpfr_init(r26320); mpfr_init(r26321); mpfr_init(r26322); mpfr_init(r26323); mpfr_init(r26324); mpfr_init(r26325); mpfr_init(r26326); mpfr_init(r26327); mpfr_init(r26328); mpfr_init(r26329); mpfr_init(r26330); mpfr_init(r26331); mpfr_init(r26332); mpfr_init(r26333); mpfr_init(r26334); mpfr_init(r26335); mpfr_init(r26336); mpfr_init(r26337); mpfr_init(r26338); mpfr_init(r26339); mpfr_init(r26340); mpfr_init(r26341); mpfr_init(r26342); mpfr_init(r26343); mpfr_init(r26344); mpfr_init(r26345); mpfr_init(r26346); mpfr_init(r26347); mpfr_init(r26348); mpfr_init(r26349); mpfr_init(r26350); mpfr_init(r26351); mpfr_init(r26352); mpfr_init(r26353); mpfr_init(r26354); mpfr_init(r26355); mpfr_init_set_str(r26356, "-2.982399978237135e+306", 10, MPFR_RNDN); mpfr_init(r26357); mpfr_init(r26358); mpfr_init(r26359); mpfr_init(r26360); mpfr_init(r26361); mpfr_init(r26362); mpfr_init(r26363); mpfr_init(r26364); mpfr_init(r26365); mpfr_init(r26366); mpfr_init(r26367); mpfr_init(r26368); mpfr_init(r26369); mpfr_init(r26370); mpfr_init(r26371); mpfr_init(r26372); mpfr_init(r26373); mpfr_init(r26374); mpfr_init(r26375); mpfr_init(r26376); mpfr_init(r26377); mpfr_init(r26378); mpfr_init(r26379); mpfr_init_set_str(r26380, "-5.977989999399742e+190", 10, MPFR_RNDN); mpfr_init(r26381); mpfr_init_set_str(r26382, "5.20966244590318e+293", 10, MPFR_RNDN); mpfr_init(r26383); mpfr_init(r26384); mpfr_init(r26385); mpfr_init(r26386); mpfr_init(r26387); mpfr_init(r26388); mpfr_init(r26389); mpfr_init(r26390); mpfr_init(r26391); mpfr_init(r26392); mpfr_init(r26393); mpfr_init(r26394); mpfr_init(r26395); mpfr_init(r26396); mpfr_init(r26397); mpfr_init(r26398); mpfr_init(r26399); mpfr_init(r26400); mpfr_init(r26401); mpfr_init(r26402); mpfr_init(r26403); mpfr_init(r26404); mpfr_init(r26405); mpfr_init(r26406); mpfr_init(r26407); mpfr_init(r26408); mpfr_init(r26409); mpfr_init(r26410); mpfr_init(r26411); mpfr_init(r26412); mpfr_init(r26413); mpfr_init(r26414); mpfr_init(r26415); mpfr_init(r26416); mpfr_init(r26417); mpfr_init(r26418); } double f_fm(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k, double y0, double y1, double y2, double y3, double y4, double y5) { mpfr_set_d(r26289, x, MPFR_RNDN); mpfr_set_d(r26290, y, MPFR_RNDN); mpfr_mul(r26291, r26289, r26290, MPFR_RNDN); mpfr_set_d(r26292, z, MPFR_RNDN); mpfr_set_d(r26293, t, MPFR_RNDN); mpfr_mul(r26294, r26292, r26293, MPFR_RNDN); mpfr_sub(r26295, r26291, r26294, MPFR_RNDN); mpfr_set_d(r26296, a, MPFR_RNDN); mpfr_set_d(r26297, b, MPFR_RNDN); mpfr_mul(r26298, r26296, r26297, MPFR_RNDN); mpfr_set_d(r26299, c, MPFR_RNDN); mpfr_set_d(r26300, i, MPFR_RNDN); mpfr_mul(r26301, r26299, r26300, MPFR_RNDN); mpfr_sub(r26302, r26298, r26301, MPFR_RNDN); mpfr_mul(r26303, r26295, r26302, MPFR_RNDN); mpfr_set_d(r26304, j, MPFR_RNDN); mpfr_mul(r26305, r26289, r26304, MPFR_RNDN); mpfr_set_d(r26306, k, MPFR_RNDN); mpfr_mul(r26307, r26292, r26306, MPFR_RNDN); mpfr_sub(r26308, r26305, r26307, MPFR_RNDN); mpfr_set_d(r26309, y0, MPFR_RNDN); mpfr_mul(r26310, r26309, r26297, MPFR_RNDN); mpfr_set_d(r26311, y1, MPFR_RNDN); mpfr_mul(r26312, r26311, r26300, MPFR_RNDN); mpfr_sub(r26313, r26310, r26312, MPFR_RNDN); mpfr_mul(r26314, r26308, r26313, MPFR_RNDN); mpfr_sub(r26315, r26303, r26314, MPFR_RNDN); mpfr_set_d(r26316, y2, MPFR_RNDN); mpfr_mul(r26317, r26289, r26316, MPFR_RNDN); mpfr_set_d(r26318, y3, MPFR_RNDN); mpfr_mul(r26319, r26292, r26318, MPFR_RNDN); mpfr_sub(r26320, r26317, r26319, MPFR_RNDN); mpfr_mul(r26321, r26309, r26299, MPFR_RNDN); mpfr_mul(r26322, r26311, r26296, MPFR_RNDN); mpfr_sub(r26323, r26321, r26322, MPFR_RNDN); mpfr_mul(r26324, r26320, r26323, MPFR_RNDN); mpfr_add(r26325, r26315, r26324, MPFR_RNDN); mpfr_set_d(r26326, y5, MPFR_RNDN); mpfr_mul(r26327, r26306, r26300, MPFR_RNDN); mpfr_mul(r26328, r26326, r26327, MPFR_RNDN); mpfr_mul(r26329, r26290, r26328, MPFR_RNDN); mpfr_set_d(r26330, y4, MPFR_RNDN); mpfr_mul(r26331, r26297, r26306, MPFR_RNDN); mpfr_mul(r26332, r26290, r26331, MPFR_RNDN); mpfr_mul(r26333, r26330, r26332, MPFR_RNDN); mpfr_mul(r26334, r26293, r26300, MPFR_RNDN); mpfr_mul(r26335, r26326, r26334, MPFR_RNDN); mpfr_mul(r26336, r26304, r26335, MPFR_RNDN); mpfr_add(r26337, r26333, r26336, MPFR_RNDN); mpfr_sub(r26338, r26329, r26337, MPFR_RNDN); mpfr_add(r26339, r26325, r26338, MPFR_RNDN); mpfr_mul(r26340, r26293, r26316, MPFR_RNDN); mpfr_mul(r26341, r26290, r26318, MPFR_RNDN); mpfr_sub(r26342, r26340, r26341, MPFR_RNDN); mpfr_mul(r26343, r26330, r26299, MPFR_RNDN); mpfr_mul(r26344, r26326, r26296, MPFR_RNDN); mpfr_sub(r26345, r26343, r26344, MPFR_RNDN); mpfr_mul(r26346, r26342, r26345, MPFR_RNDN); mpfr_sub(r26347, r26339, r26346, MPFR_RNDN); mpfr_mul(r26348, r26306, r26316, MPFR_RNDN); mpfr_mul(r26349, r26304, r26318, MPFR_RNDN); mpfr_sub(r26350, r26348, r26349, MPFR_RNDN); mpfr_mul(r26351, r26330, r26311, MPFR_RNDN); mpfr_mul(r26352, r26326, r26309, MPFR_RNDN); mpfr_sub(r26353, r26351, r26352, MPFR_RNDN); mpfr_mul(r26354, r26350, r26353, MPFR_RNDN); mpfr_add(r26355, r26347, r26354, MPFR_RNDN); ; mpfr_set_si(r26357, mpfr_cmp(r26355, r26356) <= 0, MPFR_RNDN); mpfr_mul(r26358, r26311, r26327, MPFR_RNDN); mpfr_mul(r26359, r26292, r26358, MPFR_RNDN); mpfr_mul(r26360, r26309, r26306, MPFR_RNDN); mpfr_mul(r26361, r26297, r26360, MPFR_RNDN); mpfr_mul(r26362, r26292, r26361, MPFR_RNDN); mpfr_mul(r26363, r26300, r26289, MPFR_RNDN); mpfr_mul(r26364, r26304, r26363, MPFR_RNDN); mpfr_mul(r26365, r26311, r26364, MPFR_RNDN); mpfr_add(r26366, r26362, r26365, MPFR_RNDN); mpfr_sub(r26367, r26359, r26366, MPFR_RNDN); mpfr_sub(r26368, r26303, r26367, MPFR_RNDN); mpfr_add(r26369, r26368, r26324, MPFR_RNDN); mpfr_mul(r26370, r26293, r26304, MPFR_RNDN); mpfr_mul(r26371, r26290, r26306, MPFR_RNDN); mpfr_sub(r26372, r26370, r26371, MPFR_RNDN); mpfr_mul(r26373, r26330, r26297, MPFR_RNDN); mpfr_mul(r26374, r26326, r26300, MPFR_RNDN); mpfr_sub(r26375, r26373, r26374, MPFR_RNDN); mpfr_mul(r26376, r26372, r26375, MPFR_RNDN); mpfr_add(r26377, r26369, r26376, MPFR_RNDN); mpfr_sub(r26378, r26377, r26346, MPFR_RNDN); mpfr_add(r26379, r26378, r26354, MPFR_RNDN); ; mpfr_set_si(r26381, mpfr_cmp(r26355, r26380) <= 0, MPFR_RNDN); ; mpfr_set_si(r26383, mpfr_cmp(r26355, r26382) <= 0, MPFR_RNDN); mpfr_add(r26384, r26325, r26376, MPFR_RNDN); mpfr_cbrt(r26385, r26346, MPFR_RNDN); mpfr_mul(r26386, r26385, r26385, MPFR_RNDN); mpfr_mul(r26387, r26386, r26385, MPFR_RNDN); mpfr_sub(r26388, r26384, r26387, MPFR_RNDN); mpfr_add(r26389, r26388, r26354, MPFR_RNDN); mpfr_mul(r26390, r26316, r26289, MPFR_RNDN); mpfr_sub(r26391, r26390, r26319, MPFR_RNDN); mpfr_mul(r26392, r26391, r26323, MPFR_RNDN); mpfr_mul(r26393, r26304, r26293, MPFR_RNDN); mpfr_mul(r26394, r26306, r26290, MPFR_RNDN); mpfr_sub(r26395, r26393, r26394, MPFR_RNDN); mpfr_mul(r26396, r26395, r26375, MPFR_RNDN); mpfr_mul(r26397, r26316, r26293, MPFR_RNDN); mpfr_sub(r26398, r26397, r26341, MPFR_RNDN); mpfr_mul(r26399, r26345, r26398, MPFR_RNDN); mpfr_sub(r26400, r26396, r26399, MPFR_RNDN); mpfr_mul(r26401, r26297, r26296, MPFR_RNDN); mpfr_mul(r26402, r26300, r26299, MPFR_RNDN); mpfr_sub(r26403, r26401, r26402, MPFR_RNDN); mpfr_mul(r26404, r26290, r26289, MPFR_RNDN); mpfr_sub(r26405, r26404, r26294, MPFR_RNDN); mpfr_mul(r26406, r26403, r26405, MPFR_RNDN); mpfr_mul(r26407, r26304, r26289, MPFR_RNDN); mpfr_sub(r26408, r26407, r26307, MPFR_RNDN); mpfr_mul(r26409, r26297, r26309, MPFR_RNDN); mpfr_mul(r26410, r26300, r26311, MPFR_RNDN); mpfr_sub(r26411, r26409, r26410, MPFR_RNDN); mpfr_mul(r26412, r26408, r26411, MPFR_RNDN); mpfr_sub(r26413, r26406, r26412, MPFR_RNDN); mpfr_add(r26414, r26400, r26413, MPFR_RNDN); mpfr_add(r26415, r26392, r26414, MPFR_RNDN); if (mpfr_get_si(r26383, MPFR_RNDN)) { mpfr_set(r26416, r26389, MPFR_RNDN); } else { mpfr_set(r26416, r26415, MPFR_RNDN); }; if (mpfr_get_si(r26381, MPFR_RNDN)) { mpfr_set(r26417, r26355, MPFR_RNDN); } else { mpfr_set(r26417, r26416, MPFR_RNDN); }; if (mpfr_get_si(r26357, MPFR_RNDN)) { mpfr_set(r26418, r26379, MPFR_RNDN); } else { mpfr_set(r26418, r26417, MPFR_RNDN); }; return mpfr_get_d(r26418, MPFR_RNDN); } static mpfr_t r26419, r26420, r26421, r26422, r26423, r26424, r26425, r26426, r26427, r26428, r26429, r26430, r26431, r26432, r26433, r26434, r26435, r26436, r26437, r26438, r26439, r26440, r26441, r26442, r26443, r26444, r26445, r26446, r26447, r26448, r26449, r26450, r26451, r26452, r26453, r26454, r26455, r26456, r26457, r26458, r26459, r26460, r26461, r26462, r26463, r26464, r26465, r26466, r26467, r26468, r26469, r26470, r26471, r26472, r26473, r26474, r26475, r26476, r26477, r26478, r26479, r26480, r26481, r26482, r26483, r26484, r26485, r26486, r26487, r26488, r26489, r26490, r26491, r26492, r26493, r26494, r26495, r26496, r26497, r26498, r26499, r26500, r26501, r26502, r26503, r26504, r26505, r26506, r26507, r26508, r26509, r26510, r26511, r26512, r26513, r26514, r26515, r26516, r26517, r26518, r26519, r26520, r26521, r26522, r26523, r26524, r26525, r26526, r26527, r26528, r26529, r26530, r26531, r26532, r26533, r26534, r26535, r26536, r26537, r26538, r26539, r26540, r26541, r26542, r26543, r26544, r26545, r26546, r26547, r26548; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r26419); mpfr_init(r26420); mpfr_init(r26421); mpfr_init(r26422); mpfr_init(r26423); mpfr_init(r26424); mpfr_init(r26425); mpfr_init(r26426); mpfr_init(r26427); mpfr_init(r26428); mpfr_init(r26429); mpfr_init(r26430); mpfr_init(r26431); mpfr_init(r26432); mpfr_init(r26433); mpfr_init(r26434); mpfr_init(r26435); mpfr_init(r26436); mpfr_init(r26437); mpfr_init(r26438); mpfr_init(r26439); mpfr_init(r26440); mpfr_init(r26441); mpfr_init(r26442); mpfr_init(r26443); mpfr_init(r26444); mpfr_init(r26445); mpfr_init(r26446); mpfr_init(r26447); mpfr_init(r26448); mpfr_init(r26449); mpfr_init(r26450); mpfr_init(r26451); mpfr_init(r26452); mpfr_init(r26453); mpfr_init(r26454); mpfr_init(r26455); mpfr_init(r26456); mpfr_init(r26457); mpfr_init(r26458); mpfr_init(r26459); mpfr_init(r26460); mpfr_init(r26461); mpfr_init(r26462); mpfr_init(r26463); mpfr_init(r26464); mpfr_init(r26465); mpfr_init(r26466); mpfr_init(r26467); mpfr_init(r26468); mpfr_init(r26469); mpfr_init(r26470); mpfr_init(r26471); mpfr_init(r26472); mpfr_init(r26473); mpfr_init(r26474); mpfr_init(r26475); mpfr_init(r26476); mpfr_init(r26477); mpfr_init(r26478); mpfr_init(r26479); mpfr_init(r26480); mpfr_init(r26481); mpfr_init(r26482); mpfr_init(r26483); mpfr_init(r26484); mpfr_init(r26485); mpfr_init_set_str(r26486, "-2.982399978237135e+306", 10, MPFR_RNDN); mpfr_init(r26487); mpfr_init(r26488); mpfr_init(r26489); mpfr_init(r26490); mpfr_init(r26491); mpfr_init(r26492); mpfr_init(r26493); mpfr_init(r26494); mpfr_init(r26495); mpfr_init(r26496); mpfr_init(r26497); mpfr_init(r26498); mpfr_init(r26499); mpfr_init(r26500); mpfr_init(r26501); mpfr_init(r26502); mpfr_init(r26503); mpfr_init(r26504); mpfr_init(r26505); mpfr_init(r26506); mpfr_init(r26507); mpfr_init(r26508); mpfr_init(r26509); mpfr_init_set_str(r26510, "-5.977989999399742e+190", 10, MPFR_RNDN); mpfr_init(r26511); mpfr_init_set_str(r26512, "5.20966244590318e+293", 10, MPFR_RNDN); mpfr_init(r26513); mpfr_init(r26514); mpfr_init(r26515); mpfr_init(r26516); mpfr_init(r26517); mpfr_init(r26518); mpfr_init(r26519); mpfr_init(r26520); mpfr_init(r26521); mpfr_init(r26522); mpfr_init(r26523); mpfr_init(r26524); mpfr_init(r26525); mpfr_init(r26526); mpfr_init(r26527); mpfr_init(r26528); mpfr_init(r26529); mpfr_init(r26530); mpfr_init(r26531); mpfr_init(r26532); mpfr_init(r26533); mpfr_init(r26534); mpfr_init(r26535); mpfr_init(r26536); mpfr_init(r26537); mpfr_init(r26538); mpfr_init(r26539); mpfr_init(r26540); mpfr_init(r26541); mpfr_init(r26542); mpfr_init(r26543); mpfr_init(r26544); mpfr_init(r26545); mpfr_init(r26546); mpfr_init(r26547); mpfr_init(r26548); } double f_dm(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k, double y0, double y1, double y2, double y3, double y4, double y5) { mpfr_set_d(r26419, x, MPFR_RNDN); mpfr_set_d(r26420, y, MPFR_RNDN); mpfr_mul(r26421, r26419, r26420, MPFR_RNDN); mpfr_set_d(r26422, z, MPFR_RNDN); mpfr_set_d(r26423, t, MPFR_RNDN); mpfr_mul(r26424, r26422, r26423, MPFR_RNDN); mpfr_sub(r26425, r26421, r26424, MPFR_RNDN); mpfr_set_d(r26426, a, MPFR_RNDN); mpfr_set_d(r26427, b, MPFR_RNDN); mpfr_mul(r26428, r26426, r26427, MPFR_RNDN); mpfr_set_d(r26429, c, MPFR_RNDN); mpfr_set_d(r26430, i, MPFR_RNDN); mpfr_mul(r26431, r26429, r26430, MPFR_RNDN); mpfr_sub(r26432, r26428, r26431, MPFR_RNDN); mpfr_mul(r26433, r26425, r26432, MPFR_RNDN); mpfr_set_d(r26434, j, MPFR_RNDN); mpfr_mul(r26435, r26419, r26434, MPFR_RNDN); mpfr_set_d(r26436, k, MPFR_RNDN); mpfr_mul(r26437, r26422, r26436, MPFR_RNDN); mpfr_sub(r26438, r26435, r26437, MPFR_RNDN); mpfr_set_d(r26439, y0, MPFR_RNDN); mpfr_mul(r26440, r26439, r26427, MPFR_RNDN); mpfr_set_d(r26441, y1, MPFR_RNDN); mpfr_mul(r26442, r26441, r26430, MPFR_RNDN); mpfr_sub(r26443, r26440, r26442, MPFR_RNDN); mpfr_mul(r26444, r26438, r26443, MPFR_RNDN); mpfr_sub(r26445, r26433, r26444, MPFR_RNDN); mpfr_set_d(r26446, y2, MPFR_RNDN); mpfr_mul(r26447, r26419, r26446, MPFR_RNDN); mpfr_set_d(r26448, y3, MPFR_RNDN); mpfr_mul(r26449, r26422, r26448, MPFR_RNDN); mpfr_sub(r26450, r26447, r26449, MPFR_RNDN); mpfr_mul(r26451, r26439, r26429, MPFR_RNDN); mpfr_mul(r26452, r26441, r26426, MPFR_RNDN); mpfr_sub(r26453, r26451, r26452, MPFR_RNDN); mpfr_mul(r26454, r26450, r26453, MPFR_RNDN); mpfr_add(r26455, r26445, r26454, MPFR_RNDN); mpfr_set_d(r26456, y5, MPFR_RNDN); mpfr_mul(r26457, r26436, r26430, MPFR_RNDN); mpfr_mul(r26458, r26456, r26457, MPFR_RNDN); mpfr_mul(r26459, r26420, r26458, MPFR_RNDN); mpfr_set_d(r26460, y4, MPFR_RNDN); mpfr_mul(r26461, r26427, r26436, MPFR_RNDN); mpfr_mul(r26462, r26420, r26461, MPFR_RNDN); mpfr_mul(r26463, r26460, r26462, MPFR_RNDN); mpfr_mul(r26464, r26423, r26430, MPFR_RNDN); mpfr_mul(r26465, r26456, r26464, MPFR_RNDN); mpfr_mul(r26466, r26434, r26465, MPFR_RNDN); mpfr_add(r26467, r26463, r26466, MPFR_RNDN); mpfr_sub(r26468, r26459, r26467, MPFR_RNDN); mpfr_add(r26469, r26455, r26468, MPFR_RNDN); mpfr_mul(r26470, r26423, r26446, MPFR_RNDN); mpfr_mul(r26471, r26420, r26448, MPFR_RNDN); mpfr_sub(r26472, r26470, r26471, MPFR_RNDN); mpfr_mul(r26473, r26460, r26429, MPFR_RNDN); mpfr_mul(r26474, r26456, r26426, MPFR_RNDN); mpfr_sub(r26475, r26473, r26474, MPFR_RNDN); mpfr_mul(r26476, r26472, r26475, MPFR_RNDN); mpfr_sub(r26477, r26469, r26476, MPFR_RNDN); mpfr_mul(r26478, r26436, r26446, MPFR_RNDN); mpfr_mul(r26479, r26434, r26448, MPFR_RNDN); mpfr_sub(r26480, r26478, r26479, MPFR_RNDN); mpfr_mul(r26481, r26460, r26441, MPFR_RNDN); mpfr_mul(r26482, r26456, r26439, MPFR_RNDN); mpfr_sub(r26483, r26481, r26482, MPFR_RNDN); mpfr_mul(r26484, r26480, r26483, MPFR_RNDN); mpfr_add(r26485, r26477, r26484, MPFR_RNDN); ; mpfr_set_si(r26487, mpfr_cmp(r26485, r26486) <= 0, MPFR_RNDN); mpfr_mul(r26488, r26441, r26457, MPFR_RNDN); mpfr_mul(r26489, r26422, r26488, MPFR_RNDN); mpfr_mul(r26490, r26439, r26436, MPFR_RNDN); mpfr_mul(r26491, r26427, r26490, MPFR_RNDN); mpfr_mul(r26492, r26422, r26491, MPFR_RNDN); mpfr_mul(r26493, r26430, r26419, MPFR_RNDN); mpfr_mul(r26494, r26434, r26493, MPFR_RNDN); mpfr_mul(r26495, r26441, r26494, MPFR_RNDN); mpfr_add(r26496, r26492, r26495, MPFR_RNDN); mpfr_sub(r26497, r26489, r26496, MPFR_RNDN); mpfr_sub(r26498, r26433, r26497, MPFR_RNDN); mpfr_add(r26499, r26498, r26454, MPFR_RNDN); mpfr_mul(r26500, r26423, r26434, MPFR_RNDN); mpfr_mul(r26501, r26420, r26436, MPFR_RNDN); mpfr_sub(r26502, r26500, r26501, MPFR_RNDN); mpfr_mul(r26503, r26460, r26427, MPFR_RNDN); mpfr_mul(r26504, r26456, r26430, MPFR_RNDN); mpfr_sub(r26505, r26503, r26504, MPFR_RNDN); mpfr_mul(r26506, r26502, r26505, MPFR_RNDN); mpfr_add(r26507, r26499, r26506, MPFR_RNDN); mpfr_sub(r26508, r26507, r26476, MPFR_RNDN); mpfr_add(r26509, r26508, r26484, MPFR_RNDN); ; mpfr_set_si(r26511, mpfr_cmp(r26485, r26510) <= 0, MPFR_RNDN); ; mpfr_set_si(r26513, mpfr_cmp(r26485, r26512) <= 0, MPFR_RNDN); mpfr_add(r26514, r26455, r26506, MPFR_RNDN); mpfr_cbrt(r26515, r26476, MPFR_RNDN); mpfr_mul(r26516, r26515, r26515, MPFR_RNDN); mpfr_mul(r26517, r26516, r26515, MPFR_RNDN); mpfr_sub(r26518, r26514, r26517, MPFR_RNDN); mpfr_add(r26519, r26518, r26484, MPFR_RNDN); mpfr_mul(r26520, r26446, r26419, MPFR_RNDN); mpfr_sub(r26521, r26520, r26449, MPFR_RNDN); mpfr_mul(r26522, r26521, r26453, MPFR_RNDN); mpfr_mul(r26523, r26434, r26423, MPFR_RNDN); mpfr_mul(r26524, r26436, r26420, MPFR_RNDN); mpfr_sub(r26525, r26523, r26524, MPFR_RNDN); mpfr_mul(r26526, r26525, r26505, MPFR_RNDN); mpfr_mul(r26527, r26446, r26423, MPFR_RNDN); mpfr_sub(r26528, r26527, r26471, MPFR_RNDN); mpfr_mul(r26529, r26475, r26528, MPFR_RNDN); mpfr_sub(r26530, r26526, r26529, MPFR_RNDN); mpfr_mul(r26531, r26427, r26426, MPFR_RNDN); mpfr_mul(r26532, r26430, r26429, MPFR_RNDN); mpfr_sub(r26533, r26531, r26532, MPFR_RNDN); mpfr_mul(r26534, r26420, r26419, MPFR_RNDN); mpfr_sub(r26535, r26534, r26424, MPFR_RNDN); mpfr_mul(r26536, r26533, r26535, MPFR_RNDN); mpfr_mul(r26537, r26434, r26419, MPFR_RNDN); mpfr_sub(r26538, r26537, r26437, MPFR_RNDN); mpfr_mul(r26539, r26427, r26439, MPFR_RNDN); mpfr_mul(r26540, r26430, r26441, MPFR_RNDN); mpfr_sub(r26541, r26539, r26540, MPFR_RNDN); mpfr_mul(r26542, r26538, r26541, MPFR_RNDN); mpfr_sub(r26543, r26536, r26542, MPFR_RNDN); mpfr_add(r26544, r26530, r26543, MPFR_RNDN); mpfr_add(r26545, r26522, r26544, MPFR_RNDN); if (mpfr_get_si(r26513, MPFR_RNDN)) { mpfr_set(r26546, r26519, MPFR_RNDN); } else { mpfr_set(r26546, r26545, MPFR_RNDN); }; if (mpfr_get_si(r26511, MPFR_RNDN)) { mpfr_set(r26547, r26485, MPFR_RNDN); } else { mpfr_set(r26547, r26546, MPFR_RNDN); }; if (mpfr_get_si(r26487, MPFR_RNDN)) { mpfr_set(r26548, r26509, MPFR_RNDN); } else { mpfr_set(r26548, r26547, MPFR_RNDN); }; return mpfr_get_d(r26548, MPFR_RNDN); }