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