#include #include #include #include #include char *name = "Diagrams.Solve.Polynomial:cubForm from diagrams-solve-0.1"; double f_if(float x, float y, float z, float t, float a, float b, float c, float i, float j, float k) { float r25232 = x; float r25233 = 18.0; float r25234 = r25232 * r25233; float r25235 = y; float r25236 = r25234 * r25235; float r25237 = z; float r25238 = r25236 * r25237; float r25239 = t; float r25240 = r25238 * r25239; float r25241 = a; float r25242 = 4.0; float r25243 = r25241 * r25242; float r25244 = r25243 * r25239; float r25245 = r25240 - r25244; float r25246 = b; float r25247 = c; float r25248 = r25246 * r25247; float r25249 = r25245 + r25248; float r25250 = r25232 * r25242; float r25251 = i; float r25252 = r25250 * r25251; float r25253 = r25249 - r25252; float r25254 = j; float r25255 = 27.0; float r25256 = r25254 * r25255; float r25257 = k; float r25258 = r25256 * r25257; float r25259 = r25253 - r25258; return r25259; } double f_id(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) { double r25260 = x; double r25261 = 18.0; double r25262 = r25260 * r25261; double r25263 = y; double r25264 = r25262 * r25263; double r25265 = z; double r25266 = r25264 * r25265; double r25267 = t; double r25268 = r25266 * r25267; double r25269 = a; double r25270 = 4.0; double r25271 = r25269 * r25270; double r25272 = r25271 * r25267; double r25273 = r25268 - r25272; double r25274 = b; double r25275 = c; double r25276 = r25274 * r25275; double r25277 = r25273 + r25276; double r25278 = r25260 * r25270; double r25279 = i; double r25280 = r25278 * r25279; double r25281 = r25277 - r25280; double r25282 = j; double r25283 = 27.0; double r25284 = r25282 * r25283; double r25285 = k; double r25286 = r25284 * r25285; double r25287 = r25281 - r25286; return r25287; } double f_of(float x, float y, float z, float t, float a, float b, float c, float i, float j, float k) { float r25288 = x; float r25289 = t; float r25290 = r25288 * r25289; float r25291 = y; float r25292 = 18.0; float r25293 = r25291 * r25292; float r25294 = r25290 * r25293; float r25295 = z; float r25296 = a; float r25297 = -r25296; float r25298 = 4.0; float r25299 = r25289 * r25298; float r25300 = c; float r25301 = b; float r25302 = r25300 * r25301; float r25303 = fma(r25297, r25299, r25302); float r25304 = fma(r25294, r25295, r25303); float r25305 = cbrt(r25304); float r25306 = r25305 * r25305; float r25307 = cbrt(r25306); float r25308 = cbrt(r25305); float r25309 = r25307 * r25308; float r25310 = r25309 * r25306; float r25311 = i; float r25312 = r25298 * r25288; float r25313 = j; float r25314 = k; float r25315 = 27.0; float r25316 = r25314 * r25315; float r25317 = r25313 * r25316; float r25318 = fma(r25311, r25312, r25317); float r25319 = r25310 - r25318; float r25320 = -inf.0; bool r25321 = r25319 <= r25320; float r25322 = r25288 * r25292; float r25323 = r25291 * r25295; float r25324 = r25322 * r25323; float r25325 = fma(r25289, r25324, r25303); float r25326 = r25313 * r25314; float r25327 = r25288 * r25311; float r25328 = r25327 * r25298; float r25329 = fma(r25326, r25315, r25328); float r25330 = r25325 - r25329; float r25331 = 1.290934848658723e+280; bool r25332 = r25319 <= r25331; float r25333 = r25290 * r25291; float r25334 = r25333 * r25292; float r25335 = fma(r25334, r25295, r25303); float r25336 = r25335 - r25318; float r25337 = r25289 * r25292; float r25338 = r25295 * r25288; float r25339 = r25337 * r25338; float r25340 = fma(r25299, r25297, r25302); float r25341 = fma(r25339, r25291, r25340); float r25342 = r25341 - r25329; float r25343 = r25332 ? r25336 : r25342; float r25344 = r25321 ? r25330 : r25343; return r25344; } double f_od(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) { double r25345 = x; double r25346 = t; double r25347 = r25345 * r25346; double r25348 = y; double r25349 = 18.0; double r25350 = r25348 * r25349; double r25351 = r25347 * r25350; double r25352 = z; double r25353 = a; double r25354 = -r25353; double r25355 = 4.0; double r25356 = r25346 * r25355; double r25357 = c; double r25358 = b; double r25359 = r25357 * r25358; double r25360 = fma(r25354, r25356, r25359); double r25361 = fma(r25351, r25352, r25360); double r25362 = cbrt(r25361); double r25363 = r25362 * r25362; double r25364 = cbrt(r25363); double r25365 = cbrt(r25362); double r25366 = r25364 * r25365; double r25367 = r25366 * r25363; double r25368 = i; double r25369 = r25355 * r25345; double r25370 = j; double r25371 = k; double r25372 = 27.0; double r25373 = r25371 * r25372; double r25374 = r25370 * r25373; double r25375 = fma(r25368, r25369, r25374); double r25376 = r25367 - r25375; double r25377 = -inf.0; bool r25378 = r25376 <= r25377; double r25379 = r25345 * r25349; double r25380 = r25348 * r25352; double r25381 = r25379 * r25380; double r25382 = fma(r25346, r25381, r25360); double r25383 = r25370 * r25371; double r25384 = r25345 * r25368; double r25385 = r25384 * r25355; double r25386 = fma(r25383, r25372, r25385); double r25387 = r25382 - r25386; double r25388 = 1.290934848658723e+280; bool r25389 = r25376 <= r25388; double r25390 = r25347 * r25348; double r25391 = r25390 * r25349; double r25392 = fma(r25391, r25352, r25360); double r25393 = r25392 - r25375; double r25394 = r25346 * r25349; double r25395 = r25352 * r25345; double r25396 = r25394 * r25395; double r25397 = fma(r25356, r25354, r25359); double r25398 = fma(r25396, r25348, r25397); double r25399 = r25398 - r25386; double r25400 = r25389 ? r25393 : r25399; double r25401 = r25378 ? r25387 : r25400; return r25401; } 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 r25402, r25403, r25404, r25405, r25406, r25407, r25408, r25409, r25410, r25411, r25412, r25413, r25414, r25415, r25416, r25417, r25418, r25419, r25420, r25421, r25422, r25423, r25424, r25425, r25426, r25427, r25428, r25429; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r25402); mpfr_init_set_str(r25403, "18.0", 10, MPFR_RNDN); mpfr_init(r25404); mpfr_init(r25405); mpfr_init(r25406); mpfr_init(r25407); mpfr_init(r25408); mpfr_init(r25409); mpfr_init(r25410); mpfr_init(r25411); mpfr_init_set_str(r25412, "4.0", 10, MPFR_RNDN); mpfr_init(r25413); mpfr_init(r25414); mpfr_init(r25415); mpfr_init(r25416); mpfr_init(r25417); mpfr_init(r25418); mpfr_init(r25419); mpfr_init(r25420); mpfr_init(r25421); mpfr_init(r25422); mpfr_init(r25423); mpfr_init(r25424); mpfr_init_set_str(r25425, "27.0", 10, MPFR_RNDN); mpfr_init(r25426); mpfr_init(r25427); mpfr_init(r25428); mpfr_init(r25429); } double f_im(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) { mpfr_set_d(r25402, x, MPFR_RNDN); ; mpfr_mul(r25404, r25402, r25403, MPFR_RNDN); mpfr_set_d(r25405, y, MPFR_RNDN); mpfr_mul(r25406, r25404, r25405, MPFR_RNDN); mpfr_set_d(r25407, z, MPFR_RNDN); mpfr_mul(r25408, r25406, r25407, MPFR_RNDN); mpfr_set_d(r25409, t, MPFR_RNDN); mpfr_mul(r25410, r25408, r25409, MPFR_RNDN); mpfr_set_d(r25411, a, MPFR_RNDN); ; mpfr_mul(r25413, r25411, r25412, MPFR_RNDN); mpfr_mul(r25414, r25413, r25409, MPFR_RNDN); mpfr_sub(r25415, r25410, r25414, MPFR_RNDN); mpfr_set_d(r25416, b, MPFR_RNDN); mpfr_set_d(r25417, c, MPFR_RNDN); mpfr_mul(r25418, r25416, r25417, MPFR_RNDN); mpfr_add(r25419, r25415, r25418, MPFR_RNDN); mpfr_mul(r25420, r25402, r25412, MPFR_RNDN); mpfr_set_d(r25421, i, MPFR_RNDN); mpfr_mul(r25422, r25420, r25421, MPFR_RNDN); mpfr_sub(r25423, r25419, r25422, MPFR_RNDN); mpfr_set_d(r25424, j, MPFR_RNDN); ; mpfr_mul(r25426, r25424, r25425, MPFR_RNDN); mpfr_set_d(r25427, k, MPFR_RNDN); mpfr_mul(r25428, r25426, r25427, MPFR_RNDN); mpfr_sub(r25429, r25423, r25428, MPFR_RNDN); return mpfr_get_d(r25429, MPFR_RNDN); } static mpfr_t r25430, r25431, r25432, r25433, r25434, r25435, r25436, r25437, r25438, r25439, r25440, r25441, r25442, r25443, r25444, r25445, r25446, r25447, r25448, r25449, r25450, r25451, r25452, r25453, r25454, r25455, r25456, r25457, r25458, r25459, r25460, r25461, r25462, r25463, r25464, r25465, r25466, r25467, r25468, r25469, r25470, r25471, r25472, r25473, r25474, r25475, r25476, r25477, r25478, r25479, r25480, r25481, r25482, r25483, r25484, r25485, r25486; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r25430); mpfr_init(r25431); mpfr_init(r25432); mpfr_init(r25433); mpfr_init_set_str(r25434, "18.0", 10, MPFR_RNDN); mpfr_init(r25435); mpfr_init(r25436); mpfr_init(r25437); mpfr_init(r25438); mpfr_init(r25439); mpfr_init_set_str(r25440, "4.0", 10, MPFR_RNDN); mpfr_init(r25441); mpfr_init(r25442); mpfr_init(r25443); mpfr_init(r25444); mpfr_init(r25445); mpfr_init(r25446); mpfr_init(r25447); mpfr_init(r25448); mpfr_init(r25449); mpfr_init(r25450); mpfr_init(r25451); mpfr_init(r25452); mpfr_init(r25453); mpfr_init(r25454); mpfr_init(r25455); mpfr_init(r25456); mpfr_init_set_str(r25457, "27.0", 10, MPFR_RNDN); mpfr_init(r25458); mpfr_init(r25459); mpfr_init(r25460); mpfr_init(r25461); mpfr_init_set_str(r25462, "-inf.0", 10, MPFR_RNDN); mpfr_init(r25463); mpfr_init(r25464); mpfr_init(r25465); mpfr_init(r25466); mpfr_init(r25467); mpfr_init(r25468); mpfr_init(r25469); mpfr_init(r25470); mpfr_init(r25471); mpfr_init(r25472); mpfr_init_set_str(r25473, "1.290934848658723e+280", 10, MPFR_RNDN); mpfr_init(r25474); mpfr_init(r25475); mpfr_init(r25476); mpfr_init(r25477); mpfr_init(r25478); mpfr_init(r25479); mpfr_init(r25480); mpfr_init(r25481); mpfr_init(r25482); mpfr_init(r25483); mpfr_init(r25484); mpfr_init(r25485); mpfr_init(r25486); } double f_fm(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) { mpfr_set_d(r25430, x, MPFR_RNDN); mpfr_set_d(r25431, t, MPFR_RNDN); mpfr_mul(r25432, r25430, r25431, MPFR_RNDN); mpfr_set_d(r25433, y, MPFR_RNDN); ; mpfr_mul(r25435, r25433, r25434, MPFR_RNDN); mpfr_mul(r25436, r25432, r25435, MPFR_RNDN); mpfr_set_d(r25437, z, MPFR_RNDN); mpfr_set_d(r25438, a, MPFR_RNDN); mpfr_neg(r25439, r25438, MPFR_RNDN); ; mpfr_mul(r25441, r25431, r25440, MPFR_RNDN); mpfr_set_d(r25442, c, MPFR_RNDN); mpfr_set_d(r25443, b, MPFR_RNDN); mpfr_mul(r25444, r25442, r25443, MPFR_RNDN); mpfr_fma(r25445, r25439, r25441, r25444, MPFR_RNDN); mpfr_fma(r25446, r25436, r25437, r25445, MPFR_RNDN); mpfr_cbrt(r25447, r25446, MPFR_RNDN); mpfr_mul(r25448, r25447, r25447, MPFR_RNDN); mpfr_cbrt(r25449, r25448, MPFR_RNDN); mpfr_cbrt(r25450, r25447, MPFR_RNDN); mpfr_mul(r25451, r25449, r25450, MPFR_RNDN); mpfr_mul(r25452, r25451, r25448, MPFR_RNDN); mpfr_set_d(r25453, i, MPFR_RNDN); mpfr_mul(r25454, r25440, r25430, MPFR_RNDN); mpfr_set_d(r25455, j, MPFR_RNDN); mpfr_set_d(r25456, k, MPFR_RNDN); ; mpfr_mul(r25458, r25456, r25457, MPFR_RNDN); mpfr_mul(r25459, r25455, r25458, MPFR_RNDN); mpfr_fma(r25460, r25453, r25454, r25459, MPFR_RNDN); mpfr_sub(r25461, r25452, r25460, MPFR_RNDN); ; mpfr_set_si(r25463, mpfr_cmp(r25461, r25462) <= 0, MPFR_RNDN); mpfr_mul(r25464, r25430, r25434, MPFR_RNDN); mpfr_mul(r25465, r25433, r25437, MPFR_RNDN); mpfr_mul(r25466, r25464, r25465, MPFR_RNDN); mpfr_fma(r25467, r25431, r25466, r25445, MPFR_RNDN); mpfr_mul(r25468, r25455, r25456, MPFR_RNDN); mpfr_mul(r25469, r25430, r25453, MPFR_RNDN); mpfr_mul(r25470, r25469, r25440, MPFR_RNDN); mpfr_fma(r25471, r25468, r25457, r25470, MPFR_RNDN); mpfr_sub(r25472, r25467, r25471, MPFR_RNDN); ; mpfr_set_si(r25474, mpfr_cmp(r25461, r25473) <= 0, MPFR_RNDN); mpfr_mul(r25475, r25432, r25433, MPFR_RNDN); mpfr_mul(r25476, r25475, r25434, MPFR_RNDN); mpfr_fma(r25477, r25476, r25437, r25445, MPFR_RNDN); mpfr_sub(r25478, r25477, r25460, MPFR_RNDN); mpfr_mul(r25479, r25431, r25434, MPFR_RNDN); mpfr_mul(r25480, r25437, r25430, MPFR_RNDN); mpfr_mul(r25481, r25479, r25480, MPFR_RNDN); mpfr_fma(r25482, r25441, r25439, r25444, MPFR_RNDN); mpfr_fma(r25483, r25481, r25433, r25482, MPFR_RNDN); mpfr_sub(r25484, r25483, r25471, MPFR_RNDN); if (mpfr_get_si(r25474, MPFR_RNDN)) { mpfr_set(r25485, r25478, MPFR_RNDN); } else { mpfr_set(r25485, r25484, MPFR_RNDN); }; if (mpfr_get_si(r25463, MPFR_RNDN)) { mpfr_set(r25486, r25472, MPFR_RNDN); } else { mpfr_set(r25486, r25485, MPFR_RNDN); }; return mpfr_get_d(r25486, MPFR_RNDN); } static mpfr_t r25487, r25488, r25489, r25490, r25491, r25492, r25493, r25494, r25495, r25496, r25497, r25498, r25499, r25500, r25501, r25502, r25503, r25504, r25505, r25506, r25507, r25508, r25509, r25510, r25511, r25512, r25513, r25514, r25515, r25516, r25517, r25518, r25519, r25520, r25521, r25522, r25523, r25524, r25525, r25526, r25527, r25528, r25529, r25530, r25531, r25532, r25533, r25534, r25535, r25536, r25537, r25538, r25539, r25540, r25541, r25542, r25543; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r25487); mpfr_init(r25488); mpfr_init(r25489); mpfr_init(r25490); mpfr_init_set_str(r25491, "18.0", 10, MPFR_RNDN); mpfr_init(r25492); mpfr_init(r25493); mpfr_init(r25494); mpfr_init(r25495); mpfr_init(r25496); mpfr_init_set_str(r25497, "4.0", 10, MPFR_RNDN); mpfr_init(r25498); mpfr_init(r25499); mpfr_init(r25500); mpfr_init(r25501); mpfr_init(r25502); mpfr_init(r25503); mpfr_init(r25504); mpfr_init(r25505); mpfr_init(r25506); mpfr_init(r25507); mpfr_init(r25508); mpfr_init(r25509); mpfr_init(r25510); mpfr_init(r25511); mpfr_init(r25512); mpfr_init(r25513); mpfr_init_set_str(r25514, "27.0", 10, MPFR_RNDN); mpfr_init(r25515); mpfr_init(r25516); mpfr_init(r25517); mpfr_init(r25518); mpfr_init_set_str(r25519, "-inf.0", 10, MPFR_RNDN); mpfr_init(r25520); mpfr_init(r25521); mpfr_init(r25522); mpfr_init(r25523); mpfr_init(r25524); mpfr_init(r25525); mpfr_init(r25526); mpfr_init(r25527); mpfr_init(r25528); mpfr_init(r25529); mpfr_init_set_str(r25530, "1.290934848658723e+280", 10, MPFR_RNDN); mpfr_init(r25531); mpfr_init(r25532); mpfr_init(r25533); mpfr_init(r25534); mpfr_init(r25535); mpfr_init(r25536); mpfr_init(r25537); mpfr_init(r25538); mpfr_init(r25539); mpfr_init(r25540); mpfr_init(r25541); mpfr_init(r25542); mpfr_init(r25543); } double f_dm(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) { mpfr_set_d(r25487, x, MPFR_RNDN); mpfr_set_d(r25488, t, MPFR_RNDN); mpfr_mul(r25489, r25487, r25488, MPFR_RNDN); mpfr_set_d(r25490, y, MPFR_RNDN); ; mpfr_mul(r25492, r25490, r25491, MPFR_RNDN); mpfr_mul(r25493, r25489, r25492, MPFR_RNDN); mpfr_set_d(r25494, z, MPFR_RNDN); mpfr_set_d(r25495, a, MPFR_RNDN); mpfr_neg(r25496, r25495, MPFR_RNDN); ; mpfr_mul(r25498, r25488, r25497, MPFR_RNDN); mpfr_set_d(r25499, c, MPFR_RNDN); mpfr_set_d(r25500, b, MPFR_RNDN); mpfr_mul(r25501, r25499, r25500, MPFR_RNDN); mpfr_fma(r25502, r25496, r25498, r25501, MPFR_RNDN); mpfr_fma(r25503, r25493, r25494, r25502, MPFR_RNDN); mpfr_cbrt(r25504, r25503, MPFR_RNDN); mpfr_mul(r25505, r25504, r25504, MPFR_RNDN); mpfr_cbrt(r25506, r25505, MPFR_RNDN); mpfr_cbrt(r25507, r25504, MPFR_RNDN); mpfr_mul(r25508, r25506, r25507, MPFR_RNDN); mpfr_mul(r25509, r25508, r25505, MPFR_RNDN); mpfr_set_d(r25510, i, MPFR_RNDN); mpfr_mul(r25511, r25497, r25487, MPFR_RNDN); mpfr_set_d(r25512, j, MPFR_RNDN); mpfr_set_d(r25513, k, MPFR_RNDN); ; mpfr_mul(r25515, r25513, r25514, MPFR_RNDN); mpfr_mul(r25516, r25512, r25515, MPFR_RNDN); mpfr_fma(r25517, r25510, r25511, r25516, MPFR_RNDN); mpfr_sub(r25518, r25509, r25517, MPFR_RNDN); ; mpfr_set_si(r25520, mpfr_cmp(r25518, r25519) <= 0, MPFR_RNDN); mpfr_mul(r25521, r25487, r25491, MPFR_RNDN); mpfr_mul(r25522, r25490, r25494, MPFR_RNDN); mpfr_mul(r25523, r25521, r25522, MPFR_RNDN); mpfr_fma(r25524, r25488, r25523, r25502, MPFR_RNDN); mpfr_mul(r25525, r25512, r25513, MPFR_RNDN); mpfr_mul(r25526, r25487, r25510, MPFR_RNDN); mpfr_mul(r25527, r25526, r25497, MPFR_RNDN); mpfr_fma(r25528, r25525, r25514, r25527, MPFR_RNDN); mpfr_sub(r25529, r25524, r25528, MPFR_RNDN); ; mpfr_set_si(r25531, mpfr_cmp(r25518, r25530) <= 0, MPFR_RNDN); mpfr_mul(r25532, r25489, r25490, MPFR_RNDN); mpfr_mul(r25533, r25532, r25491, MPFR_RNDN); mpfr_fma(r25534, r25533, r25494, r25502, MPFR_RNDN); mpfr_sub(r25535, r25534, r25517, MPFR_RNDN); mpfr_mul(r25536, r25488, r25491, MPFR_RNDN); mpfr_mul(r25537, r25494, r25487, MPFR_RNDN); mpfr_mul(r25538, r25536, r25537, MPFR_RNDN); mpfr_fma(r25539, r25498, r25496, r25501, MPFR_RNDN); mpfr_fma(r25540, r25538, r25490, r25539, MPFR_RNDN); mpfr_sub(r25541, r25540, r25528, MPFR_RNDN); if (mpfr_get_si(r25531, MPFR_RNDN)) { mpfr_set(r25542, r25535, MPFR_RNDN); } else { mpfr_set(r25542, r25541, MPFR_RNDN); }; if (mpfr_get_si(r25520, MPFR_RNDN)) { mpfr_set(r25543, r25529, MPFR_RNDN); } else { mpfr_set(r25543, r25542, MPFR_RNDN); }; return mpfr_get_d(r25543, MPFR_RNDN); }