#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 r25920 = x; float r25921 = y; float r25922 = r25920 * r25921; float r25923 = z; float r25924 = t; float r25925 = r25923 * r25924; float r25926 = r25922 - r25925; float r25927 = a; float r25928 = b; float r25929 = r25927 * r25928; float r25930 = c; float r25931 = i; float r25932 = r25930 * r25931; float r25933 = r25929 - r25932; float r25934 = r25926 * r25933; float r25935 = j; float r25936 = r25920 * r25935; float r25937 = k; float r25938 = r25923 * r25937; float r25939 = r25936 - r25938; float r25940 = y0; float r25941 = r25940 * r25928; float r25942 = y1; float r25943 = r25942 * r25931; float r25944 = r25941 - r25943; float r25945 = r25939 * r25944; float r25946 = r25934 - r25945; float r25947 = y2; float r25948 = r25920 * r25947; float r25949 = y3; float r25950 = r25923 * r25949; float r25951 = r25948 - r25950; float r25952 = r25940 * r25930; float r25953 = r25942 * r25927; float r25954 = r25952 - r25953; float r25955 = r25951 * r25954; float r25956 = r25946 + r25955; float r25957 = r25924 * r25935; float r25958 = r25921 * r25937; float r25959 = r25957 - r25958; float r25960 = y4; float r25961 = r25960 * r25928; float r25962 = y5; float r25963 = r25962 * r25931; float r25964 = r25961 - r25963; float r25965 = r25959 * r25964; float r25966 = r25956 + r25965; float r25967 = r25924 * r25947; float r25968 = r25921 * r25949; float r25969 = r25967 - r25968; float r25970 = r25960 * r25930; float r25971 = r25962 * r25927; float r25972 = r25970 - r25971; float r25973 = r25969 * r25972; float r25974 = r25966 - r25973; float r25975 = r25937 * r25947; float r25976 = r25935 * r25949; float r25977 = r25975 - r25976; float r25978 = r25960 * r25942; float r25979 = r25962 * r25940; float r25980 = r25978 - r25979; float r25981 = r25977 * r25980; float r25982 = r25974 + r25981; return r25982; } 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 r25983 = x; double r25984 = y; double r25985 = r25983 * r25984; double r25986 = z; double r25987 = t; double r25988 = r25986 * r25987; double r25989 = r25985 - r25988; double r25990 = a; double r25991 = b; double r25992 = r25990 * r25991; double r25993 = c; double r25994 = i; double r25995 = r25993 * r25994; double r25996 = r25992 - r25995; double r25997 = r25989 * r25996; double r25998 = j; double r25999 = r25983 * r25998; double r26000 = k; double r26001 = r25986 * r26000; double r26002 = r25999 - r26001; double r26003 = y0; double r26004 = r26003 * r25991; double r26005 = y1; double r26006 = r26005 * r25994; double r26007 = r26004 - r26006; double r26008 = r26002 * r26007; double r26009 = r25997 - r26008; double r26010 = y2; double r26011 = r25983 * r26010; double r26012 = y3; double r26013 = r25986 * r26012; double r26014 = r26011 - r26013; double r26015 = r26003 * r25993; double r26016 = r26005 * r25990; double r26017 = r26015 - r26016; double r26018 = r26014 * r26017; double r26019 = r26009 + r26018; double r26020 = r25987 * r25998; double r26021 = r25984 * r26000; double r26022 = r26020 - r26021; double r26023 = y4; double r26024 = r26023 * r25991; double r26025 = y5; double r26026 = r26025 * r25994; double r26027 = r26024 - r26026; double r26028 = r26022 * r26027; double r26029 = r26019 + r26028; double r26030 = r25987 * r26010; double r26031 = r25984 * r26012; double r26032 = r26030 - r26031; double r26033 = r26023 * r25993; double r26034 = r26025 * r25990; double r26035 = r26033 - r26034; double r26036 = r26032 * r26035; double r26037 = r26029 - r26036; double r26038 = r26000 * r26010; double r26039 = r25998 * r26012; double r26040 = r26038 - r26039; double r26041 = r26023 * r26005; double r26042 = r26025 * r26003; double r26043 = r26041 - r26042; double r26044 = r26040 * r26043; double r26045 = r26037 + r26044; return r26045; } 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 r26046 = y0; float r26047 = c; float r26048 = r26046 * r26047; float r26049 = y1; float r26050 = a; float r26051 = r26049 * r26050; float r26052 = r26048 - r26051; float r26053 = y2; float r26054 = x; float r26055 = r26053 * r26054; float r26056 = z; float r26057 = y3; float r26058 = r26056 * r26057; float r26059 = r26055 - r26058; float r26060 = b; float r26061 = r26060 * r26050; float r26062 = i; float r26063 = r26062 * r26047; float r26064 = r26061 - r26063; float r26065 = y; float r26066 = r26065 * r26054; float r26067 = t; float r26068 = r26056 * r26067; float r26069 = r26066 - r26068; float r26070 = r26064 * r26069; float r26071 = fma(r26052, r26059, r26070); float r26072 = j; float r26073 = r26072 * r26054; float r26074 = k; float r26075 = r26056 * r26074; float r26076 = r26073 - r26075; float r26077 = r26060 * r26046; float r26078 = r26062 * r26049; float r26079 = r26077 - r26078; float r26080 = r26076 * r26079; float r26081 = r26071 - r26080; float r26082 = r26062 * r26074; float r26083 = y5; float r26084 = r26083 * r26065; float r26085 = r26082 * r26084; float r26086 = r26083 * r26067; float r26087 = r26062 * r26086; float r26088 = r26074 * r26060; float r26089 = y4; float r26090 = r26089 * r26065; float r26091 = r26088 * r26090; float r26092 = fma(r26072, r26087, r26091); float r26093 = r26085 - r26092; float r26094 = r26047 * r26089; float r26095 = r26050 * r26083; float r26096 = r26094 - r26095; float r26097 = r26053 * r26067; float r26098 = r26065 * r26057; float r26099 = r26097 - r26098; float r26100 = r26096 * r26099; float r26101 = r26049 * r26089; float r26102 = r26046 * r26083; float r26103 = r26101 - r26102; float r26104 = r26074 * r26053; float r26105 = r26057 * r26072; float r26106 = r26104 - r26105; float r26107 = r26103 * r26106; float r26108 = r26100 - r26107; float r26109 = r26093 - r26108; float r26110 = r26081 + r26109; float r26111 = -1.0735482180163745e+306; bool r26112 = r26110 <= r26111; float r26113 = r26053 * r26074; float r26114 = r26072 * r26057; float r26115 = r26113 - r26114; float r26116 = r26103 * r26115; float r26117 = r26060 * r26089; float r26118 = r26062 * r26083; float r26119 = r26117 - r26118; float r26120 = r26067 * r26072; float r26121 = r26074 * r26065; float r26122 = r26120 - r26121; float r26123 = r26119 * r26122; float r26124 = r26083 * r26050; float r26125 = r26094 - r26124; float r26126 = r26067 * r26053; float r26127 = r26057 * r26065; float r26128 = r26126 - r26127; float r26129 = r26125 * r26128; float r26130 = r26123 - r26129; float r26131 = r26116 + r26130; float r26132 = r26050 * r26049; float r26133 = r26048 - r26132; float r26134 = r26054 * r26053; float r26135 = r26057 * r26056; float r26136 = r26134 - r26135; float r26137 = fma(r26133, r26136, r26070); float r26138 = r26078 - r26077; float r26139 = r26074 * r26138; float r26140 = r26054 * r26062; float r26141 = r26140 * r26049; float r26142 = -r26072; float r26143 = r26141 * r26142; float r26144 = fma(r26056, r26139, r26143); float r26145 = r26137 - r26144; float r26146 = r26131 + r26145; float r26147 = -4.723058607936372e+200; bool r26148 = r26110 <= r26147; float r26149 = 1.4214353239211479e+87; bool r26150 = r26110 <= r26149; float r26151 = r26054 * r26065; float r26152 = r26151 - r26068; float r26153 = r26050 * r26060; float r26154 = r26047 * r26062; float r26155 = r26153 - r26154; float r26156 = r26152 * r26155; float r26157 = r26054 * r26072; float r26158 = r26157 - r26075; float r26159 = r26046 * r26060; float r26160 = r26049 * r26062; float r26161 = r26159 - r26160; float r26162 = r26158 * r26161; float r26163 = r26156 - r26162; float r26164 = r26134 - r26058; float r26165 = r26164 * r26052; float r26166 = r26163 + r26165; float r26167 = r26065 * r26074; float r26168 = r26120 - r26167; float r26169 = r26089 * r26060; float r26170 = r26083 * r26062; float r26171 = r26169 - r26170; float r26172 = r26168 * r26171; float r26173 = r26166 + r26172; float r26174 = r26126 - r26098; float r26175 = r26089 * r26047; float r26176 = r26175 - r26124; float r26177 = r26174 * r26176; float r26178 = r26173 - r26177; float r26179 = r26104 - r26114; float r26180 = cbrt(r26179); float r26181 = r26180 * r26180; float r26182 = r26089 * r26049; float r26183 = r26083 * r26046; float r26184 = r26182 - r26183; float r26185 = r26180 * r26184; float r26186 = r26181 * r26185; float r26187 = r26178 + r26186; float r26188 = 1.9943670921598213e+287; bool r26189 = r26110 <= r26188; float r26190 = r26072 * r26067; float r26191 = r26190 - r26121; float r26192 = r26176 * r26099; float r26193 = r26080 + r26192; float r26194 = r26071 - r26193; float r26195 = fma(r26191, r26171, r26194); float r26196 = r26189 ? r26110 : r26195; float r26197 = r26150 ? r26187 : r26196; float r26198 = r26148 ? r26110 : r26197; float r26199 = r26112 ? r26146 : r26198; return r26199; } 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 r26200 = y0; double r26201 = c; double r26202 = r26200 * r26201; double r26203 = y1; double r26204 = a; double r26205 = r26203 * r26204; double r26206 = r26202 - r26205; double r26207 = y2; double r26208 = x; double r26209 = r26207 * r26208; double r26210 = z; double r26211 = y3; double r26212 = r26210 * r26211; double r26213 = r26209 - r26212; double r26214 = b; double r26215 = r26214 * r26204; double r26216 = i; double r26217 = r26216 * r26201; double r26218 = r26215 - r26217; double r26219 = y; double r26220 = r26219 * r26208; double r26221 = t; double r26222 = r26210 * r26221; double r26223 = r26220 - r26222; double r26224 = r26218 * r26223; double r26225 = fma(r26206, r26213, r26224); double r26226 = j; double r26227 = r26226 * r26208; double r26228 = k; double r26229 = r26210 * r26228; double r26230 = r26227 - r26229; double r26231 = r26214 * r26200; double r26232 = r26216 * r26203; double r26233 = r26231 - r26232; double r26234 = r26230 * r26233; double r26235 = r26225 - r26234; double r26236 = r26216 * r26228; double r26237 = y5; double r26238 = r26237 * r26219; double r26239 = r26236 * r26238; double r26240 = r26237 * r26221; double r26241 = r26216 * r26240; double r26242 = r26228 * r26214; double r26243 = y4; double r26244 = r26243 * r26219; double r26245 = r26242 * r26244; double r26246 = fma(r26226, r26241, r26245); double r26247 = r26239 - r26246; double r26248 = r26201 * r26243; double r26249 = r26204 * r26237; double r26250 = r26248 - r26249; double r26251 = r26207 * r26221; double r26252 = r26219 * r26211; double r26253 = r26251 - r26252; double r26254 = r26250 * r26253; double r26255 = r26203 * r26243; double r26256 = r26200 * r26237; double r26257 = r26255 - r26256; double r26258 = r26228 * r26207; double r26259 = r26211 * r26226; double r26260 = r26258 - r26259; double r26261 = r26257 * r26260; double r26262 = r26254 - r26261; double r26263 = r26247 - r26262; double r26264 = r26235 + r26263; double r26265 = -1.0735482180163745e+306; bool r26266 = r26264 <= r26265; double r26267 = r26207 * r26228; double r26268 = r26226 * r26211; double r26269 = r26267 - r26268; double r26270 = r26257 * r26269; double r26271 = r26214 * r26243; double r26272 = r26216 * r26237; double r26273 = r26271 - r26272; double r26274 = r26221 * r26226; double r26275 = r26228 * r26219; double r26276 = r26274 - r26275; double r26277 = r26273 * r26276; double r26278 = r26237 * r26204; double r26279 = r26248 - r26278; double r26280 = r26221 * r26207; double r26281 = r26211 * r26219; double r26282 = r26280 - r26281; double r26283 = r26279 * r26282; double r26284 = r26277 - r26283; double r26285 = r26270 + r26284; double r26286 = r26204 * r26203; double r26287 = r26202 - r26286; double r26288 = r26208 * r26207; double r26289 = r26211 * r26210; double r26290 = r26288 - r26289; double r26291 = fma(r26287, r26290, r26224); double r26292 = r26232 - r26231; double r26293 = r26228 * r26292; double r26294 = r26208 * r26216; double r26295 = r26294 * r26203; double r26296 = -r26226; double r26297 = r26295 * r26296; double r26298 = fma(r26210, r26293, r26297); double r26299 = r26291 - r26298; double r26300 = r26285 + r26299; double r26301 = -4.723058607936372e+200; bool r26302 = r26264 <= r26301; double r26303 = 1.4214353239211479e+87; bool r26304 = r26264 <= r26303; double r26305 = r26208 * r26219; double r26306 = r26305 - r26222; double r26307 = r26204 * r26214; double r26308 = r26201 * r26216; double r26309 = r26307 - r26308; double r26310 = r26306 * r26309; double r26311 = r26208 * r26226; double r26312 = r26311 - r26229; double r26313 = r26200 * r26214; double r26314 = r26203 * r26216; double r26315 = r26313 - r26314; double r26316 = r26312 * r26315; double r26317 = r26310 - r26316; double r26318 = r26288 - r26212; double r26319 = r26318 * r26206; double r26320 = r26317 + r26319; double r26321 = r26219 * r26228; double r26322 = r26274 - r26321; double r26323 = r26243 * r26214; double r26324 = r26237 * r26216; double r26325 = r26323 - r26324; double r26326 = r26322 * r26325; double r26327 = r26320 + r26326; double r26328 = r26280 - r26252; double r26329 = r26243 * r26201; double r26330 = r26329 - r26278; double r26331 = r26328 * r26330; double r26332 = r26327 - r26331; double r26333 = r26258 - r26268; double r26334 = cbrt(r26333); double r26335 = r26334 * r26334; double r26336 = r26243 * r26203; double r26337 = r26237 * r26200; double r26338 = r26336 - r26337; double r26339 = r26334 * r26338; double r26340 = r26335 * r26339; double r26341 = r26332 + r26340; double r26342 = 1.9943670921598213e+287; bool r26343 = r26264 <= r26342; double r26344 = r26226 * r26221; double r26345 = r26344 - r26275; double r26346 = r26330 * r26253; double r26347 = r26234 + r26346; double r26348 = r26225 - r26347; double r26349 = fma(r26345, r26325, r26348); double r26350 = r26343 ? r26264 : r26349; double r26351 = r26304 ? r26341 : r26350; double r26352 = r26302 ? r26264 : r26351; double r26353 = r26266 ? r26300 : r26352; return r26353; } 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 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; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r26354); mpfr_init(r26355); mpfr_init(r26356); 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(r26380); mpfr_init(r26381); mpfr_init(r26382); 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); } 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(r26354, x, MPFR_RNDN); mpfr_set_d(r26355, y, MPFR_RNDN); mpfr_mul(r26356, r26354, r26355, MPFR_RNDN); mpfr_set_d(r26357, z, MPFR_RNDN); mpfr_set_d(r26358, t, MPFR_RNDN); mpfr_mul(r26359, r26357, r26358, MPFR_RNDN); mpfr_sub(r26360, r26356, r26359, MPFR_RNDN); mpfr_set_d(r26361, a, MPFR_RNDN); mpfr_set_d(r26362, b, MPFR_RNDN); mpfr_mul(r26363, r26361, r26362, MPFR_RNDN); mpfr_set_d(r26364, c, MPFR_RNDN); mpfr_set_d(r26365, i, MPFR_RNDN); mpfr_mul(r26366, r26364, r26365, MPFR_RNDN); mpfr_sub(r26367, r26363, r26366, MPFR_RNDN); mpfr_mul(r26368, r26360, r26367, MPFR_RNDN); mpfr_set_d(r26369, j, MPFR_RNDN); mpfr_mul(r26370, r26354, r26369, MPFR_RNDN); mpfr_set_d(r26371, k, MPFR_RNDN); mpfr_mul(r26372, r26357, r26371, MPFR_RNDN); mpfr_sub(r26373, r26370, r26372, MPFR_RNDN); mpfr_set_d(r26374, y0, MPFR_RNDN); mpfr_mul(r26375, r26374, r26362, MPFR_RNDN); mpfr_set_d(r26376, y1, MPFR_RNDN); mpfr_mul(r26377, r26376, r26365, MPFR_RNDN); mpfr_sub(r26378, r26375, r26377, MPFR_RNDN); mpfr_mul(r26379, r26373, r26378, MPFR_RNDN); mpfr_sub(r26380, r26368, r26379, MPFR_RNDN); mpfr_set_d(r26381, y2, MPFR_RNDN); mpfr_mul(r26382, r26354, r26381, MPFR_RNDN); mpfr_set_d(r26383, y3, MPFR_RNDN); mpfr_mul(r26384, r26357, r26383, MPFR_RNDN); mpfr_sub(r26385, r26382, r26384, MPFR_RNDN); mpfr_mul(r26386, r26374, r26364, MPFR_RNDN); mpfr_mul(r26387, r26376, r26361, MPFR_RNDN); mpfr_sub(r26388, r26386, r26387, MPFR_RNDN); mpfr_mul(r26389, r26385, r26388, MPFR_RNDN); mpfr_add(r26390, r26380, r26389, MPFR_RNDN); mpfr_mul(r26391, r26358, r26369, MPFR_RNDN); mpfr_mul(r26392, r26355, r26371, MPFR_RNDN); mpfr_sub(r26393, r26391, r26392, MPFR_RNDN); mpfr_set_d(r26394, y4, MPFR_RNDN); mpfr_mul(r26395, r26394, r26362, MPFR_RNDN); mpfr_set_d(r26396, y5, MPFR_RNDN); mpfr_mul(r26397, r26396, r26365, MPFR_RNDN); mpfr_sub(r26398, r26395, r26397, MPFR_RNDN); mpfr_mul(r26399, r26393, r26398, MPFR_RNDN); mpfr_add(r26400, r26390, r26399, MPFR_RNDN); mpfr_mul(r26401, r26358, r26381, MPFR_RNDN); mpfr_mul(r26402, r26355, r26383, MPFR_RNDN); mpfr_sub(r26403, r26401, r26402, MPFR_RNDN); mpfr_mul(r26404, r26394, r26364, MPFR_RNDN); mpfr_mul(r26405, r26396, r26361, MPFR_RNDN); mpfr_sub(r26406, r26404, r26405, MPFR_RNDN); mpfr_mul(r26407, r26403, r26406, MPFR_RNDN); mpfr_sub(r26408, r26400, r26407, MPFR_RNDN); mpfr_mul(r26409, r26371, r26381, MPFR_RNDN); mpfr_mul(r26410, r26369, r26383, MPFR_RNDN); mpfr_sub(r26411, r26409, r26410, MPFR_RNDN); mpfr_mul(r26412, r26394, r26376, MPFR_RNDN); mpfr_mul(r26413, r26396, r26374, MPFR_RNDN); mpfr_sub(r26414, r26412, r26413, MPFR_RNDN); mpfr_mul(r26415, r26411, r26414, MPFR_RNDN); mpfr_add(r26416, r26408, r26415, MPFR_RNDN); return mpfr_get_d(r26416, MPFR_RNDN); } static mpfr_t r26417, r26418, 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, r26549, r26550, r26551, r26552, r26553, r26554, r26555, r26556, r26557, r26558, r26559, r26560, r26561, r26562, r26563, r26564, r26565, r26566, r26567, r26568, r26569, r26570; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r26417); mpfr_init(r26418); 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_set_str(r26482, "-1.0735482180163745e+306", 10, MPFR_RNDN); mpfr_init(r26483); mpfr_init(r26484); mpfr_init(r26485); mpfr_init(r26486); 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(r26510); mpfr_init(r26511); mpfr_init(r26512); mpfr_init(r26513); mpfr_init(r26514); mpfr_init(r26515); mpfr_init(r26516); mpfr_init(r26517); mpfr_init_set_str(r26518, "-4.723058607936372e+200", 10, MPFR_RNDN); mpfr_init(r26519); mpfr_init_set_str(r26520, "1.4214353239211479e+87", 10, MPFR_RNDN); 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); mpfr_init(r26549); mpfr_init(r26550); mpfr_init(r26551); mpfr_init(r26552); mpfr_init(r26553); mpfr_init(r26554); mpfr_init(r26555); mpfr_init(r26556); mpfr_init(r26557); mpfr_init(r26558); mpfr_init_set_str(r26559, "1.9943670921598213e+287", 10, MPFR_RNDN); mpfr_init(r26560); mpfr_init(r26561); mpfr_init(r26562); mpfr_init(r26563); mpfr_init(r26564); mpfr_init(r26565); mpfr_init(r26566); mpfr_init(r26567); mpfr_init(r26568); mpfr_init(r26569); mpfr_init(r26570); } 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(r26417, y0, MPFR_RNDN); mpfr_set_d(r26418, c, MPFR_RNDN); mpfr_mul(r26419, r26417, r26418, MPFR_RNDN); mpfr_set_d(r26420, y1, MPFR_RNDN); mpfr_set_d(r26421, a, MPFR_RNDN); mpfr_mul(r26422, r26420, r26421, MPFR_RNDN); mpfr_sub(r26423, r26419, r26422, MPFR_RNDN); mpfr_set_d(r26424, y2, MPFR_RNDN); mpfr_set_d(r26425, x, MPFR_RNDN); mpfr_mul(r26426, r26424, r26425, MPFR_RNDN); mpfr_set_d(r26427, z, MPFR_RNDN); mpfr_set_d(r26428, y3, MPFR_RNDN); mpfr_mul(r26429, r26427, r26428, MPFR_RNDN); mpfr_sub(r26430, r26426, r26429, MPFR_RNDN); mpfr_set_d(r26431, b, MPFR_RNDN); mpfr_mul(r26432, r26431, r26421, MPFR_RNDN); mpfr_set_d(r26433, i, MPFR_RNDN); mpfr_mul(r26434, r26433, r26418, MPFR_RNDN); mpfr_sub(r26435, r26432, r26434, MPFR_RNDN); mpfr_set_d(r26436, y, MPFR_RNDN); mpfr_mul(r26437, r26436, r26425, MPFR_RNDN); mpfr_set_d(r26438, t, MPFR_RNDN); mpfr_mul(r26439, r26427, r26438, MPFR_RNDN); mpfr_sub(r26440, r26437, r26439, MPFR_RNDN); mpfr_mul(r26441, r26435, r26440, MPFR_RNDN); mpfr_fma(r26442, r26423, r26430, r26441, MPFR_RNDN); mpfr_set_d(r26443, j, MPFR_RNDN); mpfr_mul(r26444, r26443, r26425, MPFR_RNDN); mpfr_set_d(r26445, k, MPFR_RNDN); mpfr_mul(r26446, r26427, r26445, MPFR_RNDN); mpfr_sub(r26447, r26444, r26446, MPFR_RNDN); mpfr_mul(r26448, r26431, r26417, MPFR_RNDN); mpfr_mul(r26449, r26433, r26420, MPFR_RNDN); mpfr_sub(r26450, r26448, r26449, MPFR_RNDN); mpfr_mul(r26451, r26447, r26450, MPFR_RNDN); mpfr_sub(r26452, r26442, r26451, MPFR_RNDN); mpfr_mul(r26453, r26433, r26445, MPFR_RNDN); mpfr_set_d(r26454, y5, MPFR_RNDN); mpfr_mul(r26455, r26454, r26436, MPFR_RNDN); mpfr_mul(r26456, r26453, r26455, MPFR_RNDN); mpfr_mul(r26457, r26454, r26438, MPFR_RNDN); mpfr_mul(r26458, r26433, r26457, MPFR_RNDN); mpfr_mul(r26459, r26445, r26431, MPFR_RNDN); mpfr_set_d(r26460, y4, MPFR_RNDN); mpfr_mul(r26461, r26460, r26436, MPFR_RNDN); mpfr_mul(r26462, r26459, r26461, MPFR_RNDN); mpfr_fma(r26463, r26443, r26458, r26462, MPFR_RNDN); mpfr_sub(r26464, r26456, r26463, MPFR_RNDN); mpfr_mul(r26465, r26418, r26460, MPFR_RNDN); mpfr_mul(r26466, r26421, r26454, MPFR_RNDN); mpfr_sub(r26467, r26465, r26466, MPFR_RNDN); mpfr_mul(r26468, r26424, r26438, MPFR_RNDN); mpfr_mul(r26469, r26436, r26428, MPFR_RNDN); mpfr_sub(r26470, r26468, r26469, MPFR_RNDN); mpfr_mul(r26471, r26467, r26470, MPFR_RNDN); mpfr_mul(r26472, r26420, r26460, MPFR_RNDN); mpfr_mul(r26473, r26417, r26454, MPFR_RNDN); mpfr_sub(r26474, r26472, r26473, MPFR_RNDN); mpfr_mul(r26475, r26445, r26424, MPFR_RNDN); mpfr_mul(r26476, r26428, r26443, MPFR_RNDN); mpfr_sub(r26477, r26475, r26476, MPFR_RNDN); mpfr_mul(r26478, r26474, r26477, MPFR_RNDN); mpfr_sub(r26479, r26471, r26478, MPFR_RNDN); mpfr_sub(r26480, r26464, r26479, MPFR_RNDN); mpfr_add(r26481, r26452, r26480, MPFR_RNDN); ; mpfr_set_si(r26483, mpfr_cmp(r26481, r26482) <= 0, MPFR_RNDN); mpfr_mul(r26484, r26424, r26445, MPFR_RNDN); mpfr_mul(r26485, r26443, r26428, MPFR_RNDN); mpfr_sub(r26486, r26484, r26485, MPFR_RNDN); mpfr_mul(r26487, r26474, r26486, MPFR_RNDN); mpfr_mul(r26488, r26431, r26460, MPFR_RNDN); mpfr_mul(r26489, r26433, r26454, MPFR_RNDN); mpfr_sub(r26490, r26488, r26489, MPFR_RNDN); mpfr_mul(r26491, r26438, r26443, MPFR_RNDN); mpfr_mul(r26492, r26445, r26436, MPFR_RNDN); mpfr_sub(r26493, r26491, r26492, MPFR_RNDN); mpfr_mul(r26494, r26490, r26493, MPFR_RNDN); mpfr_mul(r26495, r26454, r26421, MPFR_RNDN); mpfr_sub(r26496, r26465, r26495, MPFR_RNDN); mpfr_mul(r26497, r26438, r26424, MPFR_RNDN); mpfr_mul(r26498, r26428, r26436, MPFR_RNDN); mpfr_sub(r26499, r26497, r26498, MPFR_RNDN); mpfr_mul(r26500, r26496, r26499, MPFR_RNDN); mpfr_sub(r26501, r26494, r26500, MPFR_RNDN); mpfr_add(r26502, r26487, r26501, MPFR_RNDN); mpfr_mul(r26503, r26421, r26420, MPFR_RNDN); mpfr_sub(r26504, r26419, r26503, MPFR_RNDN); mpfr_mul(r26505, r26425, r26424, MPFR_RNDN); mpfr_mul(r26506, r26428, r26427, MPFR_RNDN); mpfr_sub(r26507, r26505, r26506, MPFR_RNDN); mpfr_fma(r26508, r26504, r26507, r26441, MPFR_RNDN); mpfr_sub(r26509, r26449, r26448, MPFR_RNDN); mpfr_mul(r26510, r26445, r26509, MPFR_RNDN); mpfr_mul(r26511, r26425, r26433, MPFR_RNDN); mpfr_mul(r26512, r26511, r26420, MPFR_RNDN); mpfr_neg(r26513, r26443, MPFR_RNDN); mpfr_mul(r26514, r26512, r26513, MPFR_RNDN); mpfr_fma(r26515, r26427, r26510, r26514, MPFR_RNDN); mpfr_sub(r26516, r26508, r26515, MPFR_RNDN); mpfr_add(r26517, r26502, r26516, MPFR_RNDN); ; mpfr_set_si(r26519, mpfr_cmp(r26481, r26518) <= 0, MPFR_RNDN); ; mpfr_set_si(r26521, mpfr_cmp(r26481, r26520) <= 0, MPFR_RNDN); mpfr_mul(r26522, r26425, r26436, MPFR_RNDN); mpfr_sub(r26523, r26522, r26439, MPFR_RNDN); mpfr_mul(r26524, r26421, r26431, MPFR_RNDN); mpfr_mul(r26525, r26418, r26433, MPFR_RNDN); mpfr_sub(r26526, r26524, r26525, MPFR_RNDN); mpfr_mul(r26527, r26523, r26526, MPFR_RNDN); mpfr_mul(r26528, r26425, r26443, MPFR_RNDN); mpfr_sub(r26529, r26528, r26446, MPFR_RNDN); mpfr_mul(r26530, r26417, r26431, MPFR_RNDN); mpfr_mul(r26531, r26420, r26433, MPFR_RNDN); mpfr_sub(r26532, r26530, r26531, MPFR_RNDN); mpfr_mul(r26533, r26529, r26532, MPFR_RNDN); mpfr_sub(r26534, r26527, r26533, MPFR_RNDN); mpfr_sub(r26535, r26505, r26429, MPFR_RNDN); mpfr_mul(r26536, r26535, r26423, MPFR_RNDN); mpfr_add(r26537, r26534, r26536, MPFR_RNDN); mpfr_mul(r26538, r26436, r26445, MPFR_RNDN); mpfr_sub(r26539, r26491, r26538, MPFR_RNDN); mpfr_mul(r26540, r26460, r26431, MPFR_RNDN); mpfr_mul(r26541, r26454, r26433, MPFR_RNDN); mpfr_sub(r26542, r26540, r26541, MPFR_RNDN); mpfr_mul(r26543, r26539, r26542, MPFR_RNDN); mpfr_add(r26544, r26537, r26543, MPFR_RNDN); mpfr_sub(r26545, r26497, r26469, MPFR_RNDN); mpfr_mul(r26546, r26460, r26418, MPFR_RNDN); mpfr_sub(r26547, r26546, r26495, MPFR_RNDN); mpfr_mul(r26548, r26545, r26547, MPFR_RNDN); mpfr_sub(r26549, r26544, r26548, MPFR_RNDN); mpfr_sub(r26550, r26475, r26485, MPFR_RNDN); mpfr_cbrt(r26551, r26550, MPFR_RNDN); mpfr_mul(r26552, r26551, r26551, MPFR_RNDN); mpfr_mul(r26553, r26460, r26420, MPFR_RNDN); mpfr_mul(r26554, r26454, r26417, MPFR_RNDN); mpfr_sub(r26555, r26553, r26554, MPFR_RNDN); mpfr_mul(r26556, r26551, r26555, MPFR_RNDN); mpfr_mul(r26557, r26552, r26556, MPFR_RNDN); mpfr_add(r26558, r26549, r26557, MPFR_RNDN); ; mpfr_set_si(r26560, mpfr_cmp(r26481, r26559) <= 0, MPFR_RNDN); mpfr_mul(r26561, r26443, r26438, MPFR_RNDN); mpfr_sub(r26562, r26561, r26492, MPFR_RNDN); mpfr_mul(r26563, r26547, r26470, MPFR_RNDN); mpfr_add(r26564, r26451, r26563, MPFR_RNDN); mpfr_sub(r26565, r26442, r26564, MPFR_RNDN); mpfr_fma(r26566, r26562, r26542, r26565, MPFR_RNDN); if (mpfr_get_si(r26560, MPFR_RNDN)) { mpfr_set(r26567, r26481, MPFR_RNDN); } else { mpfr_set(r26567, r26566, MPFR_RNDN); }; if (mpfr_get_si(r26521, MPFR_RNDN)) { mpfr_set(r26568, r26558, MPFR_RNDN); } else { mpfr_set(r26568, r26567, MPFR_RNDN); }; if (mpfr_get_si(r26519, MPFR_RNDN)) { mpfr_set(r26569, r26481, MPFR_RNDN); } else { mpfr_set(r26569, r26568, MPFR_RNDN); }; if (mpfr_get_si(r26483, MPFR_RNDN)) { mpfr_set(r26570, r26517, MPFR_RNDN); } else { mpfr_set(r26570, r26569, MPFR_RNDN); }; return mpfr_get_d(r26570, MPFR_RNDN); } static mpfr_t r26571, r26572, r26573, r26574, r26575, r26576, r26577, r26578, r26579, r26580, r26581, r26582, r26583, r26584, r26585, r26586, r26587, r26588, r26589, r26590, r26591, r26592, r26593, r26594, r26595, r26596, r26597, r26598, r26599, r26600, r26601, r26602, r26603, r26604, r26605, r26606, r26607, r26608, r26609, r26610, r26611, r26612, r26613, r26614, r26615, r26616, r26617, r26618, r26619, r26620, r26621, r26622, r26623, r26624, r26625, r26626, r26627, r26628, r26629, r26630, r26631, r26632, r26633, r26634, r26635, r26636, r26637, r26638, r26639, r26640, r26641, r26642, r26643, r26644, r26645, r26646, r26647, r26648, r26649, r26650, r26651, r26652, r26653, r26654, r26655, r26656, r26657, r26658, r26659, r26660, r26661, r26662, r26663, r26664, r26665, r26666, r26667, r26668, r26669, r26670, r26671, r26672, r26673, r26674, r26675, r26676, r26677, r26678, r26679, r26680, r26681, r26682, r26683, r26684, r26685, r26686, r26687, r26688, r26689, r26690, r26691, r26692, r26693, r26694, r26695, r26696, r26697, r26698, r26699, r26700, r26701, r26702, r26703, r26704, r26705, r26706, r26707, r26708, r26709, r26710, r26711, r26712, r26713, r26714, r26715, r26716, r26717, r26718, r26719, r26720, r26721, r26722, r26723, r26724; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r26571); mpfr_init(r26572); mpfr_init(r26573); mpfr_init(r26574); mpfr_init(r26575); mpfr_init(r26576); mpfr_init(r26577); mpfr_init(r26578); mpfr_init(r26579); mpfr_init(r26580); mpfr_init(r26581); mpfr_init(r26582); mpfr_init(r26583); mpfr_init(r26584); mpfr_init(r26585); mpfr_init(r26586); mpfr_init(r26587); mpfr_init(r26588); mpfr_init(r26589); mpfr_init(r26590); mpfr_init(r26591); mpfr_init(r26592); mpfr_init(r26593); mpfr_init(r26594); mpfr_init(r26595); mpfr_init(r26596); mpfr_init(r26597); mpfr_init(r26598); mpfr_init(r26599); mpfr_init(r26600); mpfr_init(r26601); mpfr_init(r26602); mpfr_init(r26603); mpfr_init(r26604); mpfr_init(r26605); mpfr_init(r26606); mpfr_init(r26607); mpfr_init(r26608); mpfr_init(r26609); mpfr_init(r26610); mpfr_init(r26611); mpfr_init(r26612); mpfr_init(r26613); mpfr_init(r26614); mpfr_init(r26615); mpfr_init(r26616); mpfr_init(r26617); mpfr_init(r26618); mpfr_init(r26619); mpfr_init(r26620); mpfr_init(r26621); mpfr_init(r26622); mpfr_init(r26623); mpfr_init(r26624); mpfr_init(r26625); mpfr_init(r26626); mpfr_init(r26627); mpfr_init(r26628); mpfr_init(r26629); mpfr_init(r26630); mpfr_init(r26631); mpfr_init(r26632); mpfr_init(r26633); mpfr_init(r26634); mpfr_init(r26635); mpfr_init_set_str(r26636, "-1.0735482180163745e+306", 10, MPFR_RNDN); mpfr_init(r26637); mpfr_init(r26638); mpfr_init(r26639); mpfr_init(r26640); mpfr_init(r26641); mpfr_init(r26642); mpfr_init(r26643); mpfr_init(r26644); mpfr_init(r26645); mpfr_init(r26646); mpfr_init(r26647); mpfr_init(r26648); mpfr_init(r26649); mpfr_init(r26650); mpfr_init(r26651); mpfr_init(r26652); mpfr_init(r26653); mpfr_init(r26654); mpfr_init(r26655); mpfr_init(r26656); mpfr_init(r26657); mpfr_init(r26658); mpfr_init(r26659); mpfr_init(r26660); mpfr_init(r26661); mpfr_init(r26662); mpfr_init(r26663); mpfr_init(r26664); mpfr_init(r26665); mpfr_init(r26666); mpfr_init(r26667); mpfr_init(r26668); mpfr_init(r26669); mpfr_init(r26670); mpfr_init(r26671); mpfr_init_set_str(r26672, "-4.723058607936372e+200", 10, MPFR_RNDN); mpfr_init(r26673); mpfr_init_set_str(r26674, "1.4214353239211479e+87", 10, MPFR_RNDN); mpfr_init(r26675); mpfr_init(r26676); mpfr_init(r26677); mpfr_init(r26678); mpfr_init(r26679); mpfr_init(r26680); mpfr_init(r26681); mpfr_init(r26682); mpfr_init(r26683); mpfr_init(r26684); mpfr_init(r26685); mpfr_init(r26686); mpfr_init(r26687); mpfr_init(r26688); mpfr_init(r26689); mpfr_init(r26690); mpfr_init(r26691); mpfr_init(r26692); mpfr_init(r26693); mpfr_init(r26694); mpfr_init(r26695); mpfr_init(r26696); mpfr_init(r26697); mpfr_init(r26698); mpfr_init(r26699); mpfr_init(r26700); mpfr_init(r26701); mpfr_init(r26702); mpfr_init(r26703); mpfr_init(r26704); mpfr_init(r26705); mpfr_init(r26706); mpfr_init(r26707); mpfr_init(r26708); mpfr_init(r26709); mpfr_init(r26710); mpfr_init(r26711); mpfr_init(r26712); mpfr_init_set_str(r26713, "1.9943670921598213e+287", 10, MPFR_RNDN); mpfr_init(r26714); mpfr_init(r26715); mpfr_init(r26716); mpfr_init(r26717); mpfr_init(r26718); mpfr_init(r26719); mpfr_init(r26720); mpfr_init(r26721); mpfr_init(r26722); mpfr_init(r26723); mpfr_init(r26724); } 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(r26571, y0, MPFR_RNDN); mpfr_set_d(r26572, c, MPFR_RNDN); mpfr_mul(r26573, r26571, r26572, MPFR_RNDN); mpfr_set_d(r26574, y1, MPFR_RNDN); mpfr_set_d(r26575, a, MPFR_RNDN); mpfr_mul(r26576, r26574, r26575, MPFR_RNDN); mpfr_sub(r26577, r26573, r26576, MPFR_RNDN); mpfr_set_d(r26578, y2, MPFR_RNDN); mpfr_set_d(r26579, x, MPFR_RNDN); mpfr_mul(r26580, r26578, r26579, MPFR_RNDN); mpfr_set_d(r26581, z, MPFR_RNDN); mpfr_set_d(r26582, y3, MPFR_RNDN); mpfr_mul(r26583, r26581, r26582, MPFR_RNDN); mpfr_sub(r26584, r26580, r26583, MPFR_RNDN); mpfr_set_d(r26585, b, MPFR_RNDN); mpfr_mul(r26586, r26585, r26575, MPFR_RNDN); mpfr_set_d(r26587, i, MPFR_RNDN); mpfr_mul(r26588, r26587, r26572, MPFR_RNDN); mpfr_sub(r26589, r26586, r26588, MPFR_RNDN); mpfr_set_d(r26590, y, MPFR_RNDN); mpfr_mul(r26591, r26590, r26579, MPFR_RNDN); mpfr_set_d(r26592, t, MPFR_RNDN); mpfr_mul(r26593, r26581, r26592, MPFR_RNDN); mpfr_sub(r26594, r26591, r26593, MPFR_RNDN); mpfr_mul(r26595, r26589, r26594, MPFR_RNDN); mpfr_fma(r26596, r26577, r26584, r26595, MPFR_RNDN); mpfr_set_d(r26597, j, MPFR_RNDN); mpfr_mul(r26598, r26597, r26579, MPFR_RNDN); mpfr_set_d(r26599, k, MPFR_RNDN); mpfr_mul(r26600, r26581, r26599, MPFR_RNDN); mpfr_sub(r26601, r26598, r26600, MPFR_RNDN); mpfr_mul(r26602, r26585, r26571, MPFR_RNDN); mpfr_mul(r26603, r26587, r26574, MPFR_RNDN); mpfr_sub(r26604, r26602, r26603, MPFR_RNDN); mpfr_mul(r26605, r26601, r26604, MPFR_RNDN); mpfr_sub(r26606, r26596, r26605, MPFR_RNDN); mpfr_mul(r26607, r26587, r26599, MPFR_RNDN); mpfr_set_d(r26608, y5, MPFR_RNDN); mpfr_mul(r26609, r26608, r26590, MPFR_RNDN); mpfr_mul(r26610, r26607, r26609, MPFR_RNDN); mpfr_mul(r26611, r26608, r26592, MPFR_RNDN); mpfr_mul(r26612, r26587, r26611, MPFR_RNDN); mpfr_mul(r26613, r26599, r26585, MPFR_RNDN); mpfr_set_d(r26614, y4, MPFR_RNDN); mpfr_mul(r26615, r26614, r26590, MPFR_RNDN); mpfr_mul(r26616, r26613, r26615, MPFR_RNDN); mpfr_fma(r26617, r26597, r26612, r26616, MPFR_RNDN); mpfr_sub(r26618, r26610, r26617, MPFR_RNDN); mpfr_mul(r26619, r26572, r26614, MPFR_RNDN); mpfr_mul(r26620, r26575, r26608, MPFR_RNDN); mpfr_sub(r26621, r26619, r26620, MPFR_RNDN); mpfr_mul(r26622, r26578, r26592, MPFR_RNDN); mpfr_mul(r26623, r26590, r26582, MPFR_RNDN); mpfr_sub(r26624, r26622, r26623, MPFR_RNDN); mpfr_mul(r26625, r26621, r26624, MPFR_RNDN); mpfr_mul(r26626, r26574, r26614, MPFR_RNDN); mpfr_mul(r26627, r26571, r26608, MPFR_RNDN); mpfr_sub(r26628, r26626, r26627, MPFR_RNDN); mpfr_mul(r26629, r26599, r26578, MPFR_RNDN); mpfr_mul(r26630, r26582, r26597, MPFR_RNDN); mpfr_sub(r26631, r26629, r26630, MPFR_RNDN); mpfr_mul(r26632, r26628, r26631, MPFR_RNDN); mpfr_sub(r26633, r26625, r26632, MPFR_RNDN); mpfr_sub(r26634, r26618, r26633, MPFR_RNDN); mpfr_add(r26635, r26606, r26634, MPFR_RNDN); ; mpfr_set_si(r26637, mpfr_cmp(r26635, r26636) <= 0, MPFR_RNDN); mpfr_mul(r26638, r26578, r26599, MPFR_RNDN); mpfr_mul(r26639, r26597, r26582, MPFR_RNDN); mpfr_sub(r26640, r26638, r26639, MPFR_RNDN); mpfr_mul(r26641, r26628, r26640, MPFR_RNDN); mpfr_mul(r26642, r26585, r26614, MPFR_RNDN); mpfr_mul(r26643, r26587, r26608, MPFR_RNDN); mpfr_sub(r26644, r26642, r26643, MPFR_RNDN); mpfr_mul(r26645, r26592, r26597, MPFR_RNDN); mpfr_mul(r26646, r26599, r26590, MPFR_RNDN); mpfr_sub(r26647, r26645, r26646, MPFR_RNDN); mpfr_mul(r26648, r26644, r26647, MPFR_RNDN); mpfr_mul(r26649, r26608, r26575, MPFR_RNDN); mpfr_sub(r26650, r26619, r26649, MPFR_RNDN); mpfr_mul(r26651, r26592, r26578, MPFR_RNDN); mpfr_mul(r26652, r26582, r26590, MPFR_RNDN); mpfr_sub(r26653, r26651, r26652, MPFR_RNDN); mpfr_mul(r26654, r26650, r26653, MPFR_RNDN); mpfr_sub(r26655, r26648, r26654, MPFR_RNDN); mpfr_add(r26656, r26641, r26655, MPFR_RNDN); mpfr_mul(r26657, r26575, r26574, MPFR_RNDN); mpfr_sub(r26658, r26573, r26657, MPFR_RNDN); mpfr_mul(r26659, r26579, r26578, MPFR_RNDN); mpfr_mul(r26660, r26582, r26581, MPFR_RNDN); mpfr_sub(r26661, r26659, r26660, MPFR_RNDN); mpfr_fma(r26662, r26658, r26661, r26595, MPFR_RNDN); mpfr_sub(r26663, r26603, r26602, MPFR_RNDN); mpfr_mul(r26664, r26599, r26663, MPFR_RNDN); mpfr_mul(r26665, r26579, r26587, MPFR_RNDN); mpfr_mul(r26666, r26665, r26574, MPFR_RNDN); mpfr_neg(r26667, r26597, MPFR_RNDN); mpfr_mul(r26668, r26666, r26667, MPFR_RNDN); mpfr_fma(r26669, r26581, r26664, r26668, MPFR_RNDN); mpfr_sub(r26670, r26662, r26669, MPFR_RNDN); mpfr_add(r26671, r26656, r26670, MPFR_RNDN); ; mpfr_set_si(r26673, mpfr_cmp(r26635, r26672) <= 0, MPFR_RNDN); ; mpfr_set_si(r26675, mpfr_cmp(r26635, r26674) <= 0, MPFR_RNDN); mpfr_mul(r26676, r26579, r26590, MPFR_RNDN); mpfr_sub(r26677, r26676, r26593, MPFR_RNDN); mpfr_mul(r26678, r26575, r26585, MPFR_RNDN); mpfr_mul(r26679, r26572, r26587, MPFR_RNDN); mpfr_sub(r26680, r26678, r26679, MPFR_RNDN); mpfr_mul(r26681, r26677, r26680, MPFR_RNDN); mpfr_mul(r26682, r26579, r26597, MPFR_RNDN); mpfr_sub(r26683, r26682, r26600, MPFR_RNDN); mpfr_mul(r26684, r26571, r26585, MPFR_RNDN); mpfr_mul(r26685, r26574, r26587, MPFR_RNDN); mpfr_sub(r26686, r26684, r26685, MPFR_RNDN); mpfr_mul(r26687, r26683, r26686, MPFR_RNDN); mpfr_sub(r26688, r26681, r26687, MPFR_RNDN); mpfr_sub(r26689, r26659, r26583, MPFR_RNDN); mpfr_mul(r26690, r26689, r26577, MPFR_RNDN); mpfr_add(r26691, r26688, r26690, MPFR_RNDN); mpfr_mul(r26692, r26590, r26599, MPFR_RNDN); mpfr_sub(r26693, r26645, r26692, MPFR_RNDN); mpfr_mul(r26694, r26614, r26585, MPFR_RNDN); mpfr_mul(r26695, r26608, r26587, MPFR_RNDN); mpfr_sub(r26696, r26694, r26695, MPFR_RNDN); mpfr_mul(r26697, r26693, r26696, MPFR_RNDN); mpfr_add(r26698, r26691, r26697, MPFR_RNDN); mpfr_sub(r26699, r26651, r26623, MPFR_RNDN); mpfr_mul(r26700, r26614, r26572, MPFR_RNDN); mpfr_sub(r26701, r26700, r26649, MPFR_RNDN); mpfr_mul(r26702, r26699, r26701, MPFR_RNDN); mpfr_sub(r26703, r26698, r26702, MPFR_RNDN); mpfr_sub(r26704, r26629, r26639, MPFR_RNDN); mpfr_cbrt(r26705, r26704, MPFR_RNDN); mpfr_mul(r26706, r26705, r26705, MPFR_RNDN); mpfr_mul(r26707, r26614, r26574, MPFR_RNDN); mpfr_mul(r26708, r26608, r26571, MPFR_RNDN); mpfr_sub(r26709, r26707, r26708, MPFR_RNDN); mpfr_mul(r26710, r26705, r26709, MPFR_RNDN); mpfr_mul(r26711, r26706, r26710, MPFR_RNDN); mpfr_add(r26712, r26703, r26711, MPFR_RNDN); ; mpfr_set_si(r26714, mpfr_cmp(r26635, r26713) <= 0, MPFR_RNDN); mpfr_mul(r26715, r26597, r26592, MPFR_RNDN); mpfr_sub(r26716, r26715, r26646, MPFR_RNDN); mpfr_mul(r26717, r26701, r26624, MPFR_RNDN); mpfr_add(r26718, r26605, r26717, MPFR_RNDN); mpfr_sub(r26719, r26596, r26718, MPFR_RNDN); mpfr_fma(r26720, r26716, r26696, r26719, MPFR_RNDN); if (mpfr_get_si(r26714, MPFR_RNDN)) { mpfr_set(r26721, r26635, MPFR_RNDN); } else { mpfr_set(r26721, r26720, MPFR_RNDN); }; if (mpfr_get_si(r26675, MPFR_RNDN)) { mpfr_set(r26722, r26712, MPFR_RNDN); } else { mpfr_set(r26722, r26721, MPFR_RNDN); }; if (mpfr_get_si(r26673, MPFR_RNDN)) { mpfr_set(r26723, r26635, MPFR_RNDN); } else { mpfr_set(r26723, r26722, MPFR_RNDN); }; if (mpfr_get_si(r26637, MPFR_RNDN)) { mpfr_set(r26724, r26671, MPFR_RNDN); } else { mpfr_set(r26724, r26723, MPFR_RNDN); }; return mpfr_get_d(r26724, MPFR_RNDN); }