#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 r10898 = x; float r10899 = 18.0; float r10900 = r10898 * r10899; float r10901 = y; float r10902 = r10900 * r10901; float r10903 = z; float r10904 = r10902 * r10903; float r10905 = t; float r10906 = r10904 * r10905; float r10907 = a; float r10908 = 4.0; float r10909 = r10907 * r10908; float r10910 = r10909 * r10905; float r10911 = r10906 - r10910; float r10912 = b; float r10913 = c; float r10914 = r10912 * r10913; float r10915 = r10911 + r10914; float r10916 = r10898 * r10908; float r10917 = i; float r10918 = r10916 * r10917; float r10919 = r10915 - r10918; float r10920 = j; float r10921 = 27.0; float r10922 = r10920 * r10921; float r10923 = k; float r10924 = r10922 * r10923; float r10925 = r10919 - r10924; return r10925; } double f_id(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) { double r10926 = x; double r10927 = 18.0; double r10928 = r10926 * r10927; double r10929 = y; double r10930 = r10928 * r10929; double r10931 = z; double r10932 = r10930 * r10931; double r10933 = t; double r10934 = r10932 * r10933; double r10935 = a; double r10936 = 4.0; double r10937 = r10935 * r10936; double r10938 = r10937 * r10933; double r10939 = r10934 - r10938; double r10940 = b; double r10941 = c; double r10942 = r10940 * r10941; double r10943 = r10939 + r10942; double r10944 = r10926 * r10936; double r10945 = i; double r10946 = r10944 * r10945; double r10947 = r10943 - r10946; double r10948 = j; double r10949 = 27.0; double r10950 = r10948 * r10949; double r10951 = k; double r10952 = r10950 * r10951; double r10953 = r10947 - r10952; return r10953; } double f_of(float x, float y, float z, float t, float a, float b, float c, float i, float j, float k) { float r10954 = c; float r10955 = b; float r10956 = r10954 * r10955; float r10957 = 27.0; float r10958 = j; float r10959 = k; float r10960 = r10958 * r10959; float r10961 = r10957 * r10960; float r10962 = r10956 - r10961; float r10963 = x; float r10964 = i; float r10965 = 4.0; float r10966 = r10964 * r10965; float r10967 = z; float r10968 = t; float r10969 = r10967 * r10968; float r10970 = y; float r10971 = 18.0; float r10972 = r10970 * r10971; float r10973 = r10969 * r10972; float r10974 = r10966 - r10973; float r10975 = r10963 * r10974; float r10976 = cbrt(r10975); float r10977 = r10976 * r10976; float r10978 = r10977 * r10976; float r10979 = r10962 - r10978; float r10980 = a; float r10981 = r10980 * r10965; float r10982 = -r10981; float r10983 = r10968 * r10982; float r10984 = r10979 + r10983; float r10985 = -4.398592172198208e+304; bool r10986 = r10984 <= r10985; float r10987 = r10963 * r10971; float r10988 = r10987 * r10970; float r10989 = r10988 * r10967; float r10990 = r10989 * r10968; float r10991 = r10981 * r10968; float r10992 = r10990 - r10991; float r10993 = r10955 * r10954; float r10994 = r10992 + r10993; float r10995 = r10963 * r10965; float r10996 = r10995 * r10964; float r10997 = r10994 - r10996; float r10998 = r10958 * r10957; float r10999 = r10998 * r10959; float r11000 = r10997 - r10999; float r11001 = 1.0563176570278051e+300; bool r11002 = r10984 <= r11001; float r11003 = r10969 * r10970; float r11004 = r11003 * r10971; float r11005 = r10966 - r11004; float r11006 = r10963 * r11005; float r11007 = r10962 - r11006; float r11008 = r11007 + r10983; float r11009 = r11002 ? r11008 : r11000; float r11010 = r10986 ? r11000 : r11009; return r11010; } double f_od(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) { double r11011 = c; double r11012 = b; double r11013 = r11011 * r11012; double r11014 = 27.0; double r11015 = j; double r11016 = k; double r11017 = r11015 * r11016; double r11018 = r11014 * r11017; double r11019 = r11013 - r11018; double r11020 = x; double r11021 = i; double r11022 = 4.0; double r11023 = r11021 * r11022; double r11024 = z; double r11025 = t; double r11026 = r11024 * r11025; double r11027 = y; double r11028 = 18.0; double r11029 = r11027 * r11028; double r11030 = r11026 * r11029; double r11031 = r11023 - r11030; double r11032 = r11020 * r11031; double r11033 = cbrt(r11032); double r11034 = r11033 * r11033; double r11035 = r11034 * r11033; double r11036 = r11019 - r11035; double r11037 = a; double r11038 = r11037 * r11022; double r11039 = -r11038; double r11040 = r11025 * r11039; double r11041 = r11036 + r11040; double r11042 = -4.398592172198208e+304; bool r11043 = r11041 <= r11042; double r11044 = r11020 * r11028; double r11045 = r11044 * r11027; double r11046 = r11045 * r11024; double r11047 = r11046 * r11025; double r11048 = r11038 * r11025; double r11049 = r11047 - r11048; double r11050 = r11012 * r11011; double r11051 = r11049 + r11050; double r11052 = r11020 * r11022; double r11053 = r11052 * r11021; double r11054 = r11051 - r11053; double r11055 = r11015 * r11014; double r11056 = r11055 * r11016; double r11057 = r11054 - r11056; double r11058 = 1.0563176570278051e+300; bool r11059 = r11041 <= r11058; double r11060 = r11026 * r11027; double r11061 = r11060 * r11028; double r11062 = r11023 - r11061; double r11063 = r11020 * r11062; double r11064 = r11019 - r11063; double r11065 = r11064 + r11040; double r11066 = r11059 ? r11065 : r11057; double r11067 = r11043 ? r11057 : r11066; return r11067; } 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 r11068, r11069, r11070, r11071, r11072, r11073, r11074, r11075, r11076, r11077, r11078, r11079, r11080, r11081, r11082, r11083, r11084, r11085, r11086, r11087, r11088, r11089, r11090, r11091, r11092, r11093, r11094, r11095; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r11068); mpfr_init_set_str(r11069, "18.0", 10, MPFR_RNDN); mpfr_init(r11070); mpfr_init(r11071); mpfr_init(r11072); mpfr_init(r11073); mpfr_init(r11074); mpfr_init(r11075); mpfr_init(r11076); mpfr_init(r11077); mpfr_init_set_str(r11078, "4.0", 10, MPFR_RNDN); mpfr_init(r11079); mpfr_init(r11080); mpfr_init(r11081); mpfr_init(r11082); mpfr_init(r11083); mpfr_init(r11084); mpfr_init(r11085); mpfr_init(r11086); mpfr_init(r11087); mpfr_init(r11088); mpfr_init(r11089); mpfr_init(r11090); mpfr_init_set_str(r11091, "27.0", 10, MPFR_RNDN); mpfr_init(r11092); mpfr_init(r11093); mpfr_init(r11094); mpfr_init(r11095); } 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(r11068, x, MPFR_RNDN); ; mpfr_mul(r11070, r11068, r11069, MPFR_RNDN); mpfr_set_d(r11071, y, MPFR_RNDN); mpfr_mul(r11072, r11070, r11071, MPFR_RNDN); mpfr_set_d(r11073, z, MPFR_RNDN); mpfr_mul(r11074, r11072, r11073, MPFR_RNDN); mpfr_set_d(r11075, t, MPFR_RNDN); mpfr_mul(r11076, r11074, r11075, MPFR_RNDN); mpfr_set_d(r11077, a, MPFR_RNDN); ; mpfr_mul(r11079, r11077, r11078, MPFR_RNDN); mpfr_mul(r11080, r11079, r11075, MPFR_RNDN); mpfr_sub(r11081, r11076, r11080, MPFR_RNDN); mpfr_set_d(r11082, b, MPFR_RNDN); mpfr_set_d(r11083, c, MPFR_RNDN); mpfr_mul(r11084, r11082, r11083, MPFR_RNDN); mpfr_add(r11085, r11081, r11084, MPFR_RNDN); mpfr_mul(r11086, r11068, r11078, MPFR_RNDN); mpfr_set_d(r11087, i, MPFR_RNDN); mpfr_mul(r11088, r11086, r11087, MPFR_RNDN); mpfr_sub(r11089, r11085, r11088, MPFR_RNDN); mpfr_set_d(r11090, j, MPFR_RNDN); ; mpfr_mul(r11092, r11090, r11091, MPFR_RNDN); mpfr_set_d(r11093, k, MPFR_RNDN); mpfr_mul(r11094, r11092, r11093, MPFR_RNDN); mpfr_sub(r11095, r11089, r11094, MPFR_RNDN); return mpfr_get_d(r11095, MPFR_RNDN); } static mpfr_t r11096, r11097, r11098, r11099, r11100, r11101, r11102, r11103, r11104, r11105, r11106, r11107, r11108, r11109, r11110, r11111, r11112, r11113, r11114, r11115, r11116, r11117, r11118, r11119, r11120, r11121, r11122, r11123, r11124, r11125, r11126, r11127, r11128, r11129, r11130, r11131, r11132, r11133, r11134, r11135, r11136, r11137, r11138, r11139, r11140, r11141, r11142, r11143, r11144, r11145, r11146, r11147, r11148, r11149, r11150, r11151, r11152; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r11096); mpfr_init(r11097); mpfr_init(r11098); mpfr_init_set_str(r11099, "27.0", 10, MPFR_RNDN); mpfr_init(r11100); mpfr_init(r11101); mpfr_init(r11102); mpfr_init(r11103); mpfr_init(r11104); mpfr_init(r11105); mpfr_init(r11106); mpfr_init_set_str(r11107, "4.0", 10, MPFR_RNDN); mpfr_init(r11108); mpfr_init(r11109); mpfr_init(r11110); mpfr_init(r11111); mpfr_init(r11112); mpfr_init_set_str(r11113, "18.0", 10, MPFR_RNDN); mpfr_init(r11114); mpfr_init(r11115); mpfr_init(r11116); mpfr_init(r11117); mpfr_init(r11118); mpfr_init(r11119); mpfr_init(r11120); mpfr_init(r11121); mpfr_init(r11122); mpfr_init(r11123); mpfr_init(r11124); mpfr_init(r11125); mpfr_init(r11126); mpfr_init_set_str(r11127, "-4.398592172198208e+304", 10, MPFR_RNDN); mpfr_init(r11128); mpfr_init(r11129); mpfr_init(r11130); mpfr_init(r11131); mpfr_init(r11132); mpfr_init(r11133); mpfr_init(r11134); mpfr_init(r11135); mpfr_init(r11136); mpfr_init(r11137); mpfr_init(r11138); mpfr_init(r11139); mpfr_init(r11140); mpfr_init(r11141); mpfr_init(r11142); mpfr_init_set_str(r11143, "1.0563176570278051e+300", 10, MPFR_RNDN); mpfr_init(r11144); mpfr_init(r11145); mpfr_init(r11146); mpfr_init(r11147); mpfr_init(r11148); mpfr_init(r11149); mpfr_init(r11150); mpfr_init(r11151); mpfr_init(r11152); } 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(r11096, c, MPFR_RNDN); mpfr_set_d(r11097, b, MPFR_RNDN); mpfr_mul(r11098, r11096, r11097, MPFR_RNDN); ; mpfr_set_d(r11100, j, MPFR_RNDN); mpfr_set_d(r11101, k, MPFR_RNDN); mpfr_mul(r11102, r11100, r11101, MPFR_RNDN); mpfr_mul(r11103, r11099, r11102, MPFR_RNDN); mpfr_sub(r11104, r11098, r11103, MPFR_RNDN); mpfr_set_d(r11105, x, MPFR_RNDN); mpfr_set_d(r11106, i, MPFR_RNDN); ; mpfr_mul(r11108, r11106, r11107, MPFR_RNDN); mpfr_set_d(r11109, z, MPFR_RNDN); mpfr_set_d(r11110, t, MPFR_RNDN); mpfr_mul(r11111, r11109, r11110, MPFR_RNDN); mpfr_set_d(r11112, y, MPFR_RNDN); ; mpfr_mul(r11114, r11112, r11113, MPFR_RNDN); mpfr_mul(r11115, r11111, r11114, MPFR_RNDN); mpfr_sub(r11116, r11108, r11115, MPFR_RNDN); mpfr_mul(r11117, r11105, r11116, MPFR_RNDN); mpfr_cbrt(r11118, r11117, MPFR_RNDN); mpfr_mul(r11119, r11118, r11118, MPFR_RNDN); mpfr_mul(r11120, r11119, r11118, MPFR_RNDN); mpfr_sub(r11121, r11104, r11120, MPFR_RNDN); mpfr_set_d(r11122, a, MPFR_RNDN); mpfr_mul(r11123, r11122, r11107, MPFR_RNDN); mpfr_neg(r11124, r11123, MPFR_RNDN); mpfr_mul(r11125, r11110, r11124, MPFR_RNDN); mpfr_add(r11126, r11121, r11125, MPFR_RNDN); ; mpfr_set_si(r11128, mpfr_cmp(r11126, r11127) <= 0, MPFR_RNDN); mpfr_mul(r11129, r11105, r11113, MPFR_RNDN); mpfr_mul(r11130, r11129, r11112, MPFR_RNDN); mpfr_mul(r11131, r11130, r11109, MPFR_RNDN); mpfr_mul(r11132, r11131, r11110, MPFR_RNDN); mpfr_mul(r11133, r11123, r11110, MPFR_RNDN); mpfr_sub(r11134, r11132, r11133, MPFR_RNDN); mpfr_mul(r11135, r11097, r11096, MPFR_RNDN); mpfr_add(r11136, r11134, r11135, MPFR_RNDN); mpfr_mul(r11137, r11105, r11107, MPFR_RNDN); mpfr_mul(r11138, r11137, r11106, MPFR_RNDN); mpfr_sub(r11139, r11136, r11138, MPFR_RNDN); mpfr_mul(r11140, r11100, r11099, MPFR_RNDN); mpfr_mul(r11141, r11140, r11101, MPFR_RNDN); mpfr_sub(r11142, r11139, r11141, MPFR_RNDN); ; mpfr_set_si(r11144, mpfr_cmp(r11126, r11143) <= 0, MPFR_RNDN); mpfr_mul(r11145, r11111, r11112, MPFR_RNDN); mpfr_mul(r11146, r11145, r11113, MPFR_RNDN); mpfr_sub(r11147, r11108, r11146, MPFR_RNDN); mpfr_mul(r11148, r11105, r11147, MPFR_RNDN); mpfr_sub(r11149, r11104, r11148, MPFR_RNDN); mpfr_add(r11150, r11149, r11125, MPFR_RNDN); if (mpfr_get_si(r11144, MPFR_RNDN)) { mpfr_set(r11151, r11150, MPFR_RNDN); } else { mpfr_set(r11151, r11142, MPFR_RNDN); }; if (mpfr_get_si(r11128, MPFR_RNDN)) { mpfr_set(r11152, r11142, MPFR_RNDN); } else { mpfr_set(r11152, r11151, MPFR_RNDN); }; return mpfr_get_d(r11152, MPFR_RNDN); } static mpfr_t r11153, r11154, r11155, r11156, r11157, r11158, r11159, r11160, r11161, r11162, r11163, r11164, r11165, r11166, r11167, r11168, r11169, r11170, r11171, r11172, r11173, r11174, r11175, r11176, r11177, r11178, r11179, r11180, r11181, r11182, r11183, r11184, r11185, r11186, r11187, r11188, r11189, r11190, r11191, r11192, r11193, r11194, r11195, r11196, r11197, r11198, r11199, r11200, r11201, r11202, r11203, r11204, r11205, r11206, r11207, r11208, r11209; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r11153); mpfr_init(r11154); mpfr_init(r11155); mpfr_init_set_str(r11156, "27.0", 10, MPFR_RNDN); mpfr_init(r11157); mpfr_init(r11158); mpfr_init(r11159); mpfr_init(r11160); mpfr_init(r11161); mpfr_init(r11162); mpfr_init(r11163); mpfr_init_set_str(r11164, "4.0", 10, MPFR_RNDN); mpfr_init(r11165); mpfr_init(r11166); mpfr_init(r11167); mpfr_init(r11168); mpfr_init(r11169); mpfr_init_set_str(r11170, "18.0", 10, MPFR_RNDN); mpfr_init(r11171); mpfr_init(r11172); mpfr_init(r11173); mpfr_init(r11174); mpfr_init(r11175); mpfr_init(r11176); mpfr_init(r11177); mpfr_init(r11178); mpfr_init(r11179); mpfr_init(r11180); mpfr_init(r11181); mpfr_init(r11182); mpfr_init(r11183); mpfr_init_set_str(r11184, "-4.398592172198208e+304", 10, MPFR_RNDN); mpfr_init(r11185); mpfr_init(r11186); mpfr_init(r11187); mpfr_init(r11188); mpfr_init(r11189); mpfr_init(r11190); mpfr_init(r11191); mpfr_init(r11192); mpfr_init(r11193); mpfr_init(r11194); mpfr_init(r11195); mpfr_init(r11196); mpfr_init(r11197); mpfr_init(r11198); mpfr_init(r11199); mpfr_init_set_str(r11200, "1.0563176570278051e+300", 10, MPFR_RNDN); mpfr_init(r11201); mpfr_init(r11202); mpfr_init(r11203); mpfr_init(r11204); mpfr_init(r11205); mpfr_init(r11206); mpfr_init(r11207); mpfr_init(r11208); mpfr_init(r11209); } 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(r11153, c, MPFR_RNDN); mpfr_set_d(r11154, b, MPFR_RNDN); mpfr_mul(r11155, r11153, r11154, MPFR_RNDN); ; mpfr_set_d(r11157, j, MPFR_RNDN); mpfr_set_d(r11158, k, MPFR_RNDN); mpfr_mul(r11159, r11157, r11158, MPFR_RNDN); mpfr_mul(r11160, r11156, r11159, MPFR_RNDN); mpfr_sub(r11161, r11155, r11160, MPFR_RNDN); mpfr_set_d(r11162, x, MPFR_RNDN); mpfr_set_d(r11163, i, MPFR_RNDN); ; mpfr_mul(r11165, r11163, r11164, MPFR_RNDN); mpfr_set_d(r11166, z, MPFR_RNDN); mpfr_set_d(r11167, t, MPFR_RNDN); mpfr_mul(r11168, r11166, r11167, MPFR_RNDN); mpfr_set_d(r11169, y, MPFR_RNDN); ; mpfr_mul(r11171, r11169, r11170, MPFR_RNDN); mpfr_mul(r11172, r11168, r11171, MPFR_RNDN); mpfr_sub(r11173, r11165, r11172, MPFR_RNDN); mpfr_mul(r11174, r11162, r11173, MPFR_RNDN); mpfr_cbrt(r11175, r11174, MPFR_RNDN); mpfr_mul(r11176, r11175, r11175, MPFR_RNDN); mpfr_mul(r11177, r11176, r11175, MPFR_RNDN); mpfr_sub(r11178, r11161, r11177, MPFR_RNDN); mpfr_set_d(r11179, a, MPFR_RNDN); mpfr_mul(r11180, r11179, r11164, MPFR_RNDN); mpfr_neg(r11181, r11180, MPFR_RNDN); mpfr_mul(r11182, r11167, r11181, MPFR_RNDN); mpfr_add(r11183, r11178, r11182, MPFR_RNDN); ; mpfr_set_si(r11185, mpfr_cmp(r11183, r11184) <= 0, MPFR_RNDN); mpfr_mul(r11186, r11162, r11170, MPFR_RNDN); mpfr_mul(r11187, r11186, r11169, MPFR_RNDN); mpfr_mul(r11188, r11187, r11166, MPFR_RNDN); mpfr_mul(r11189, r11188, r11167, MPFR_RNDN); mpfr_mul(r11190, r11180, r11167, MPFR_RNDN); mpfr_sub(r11191, r11189, r11190, MPFR_RNDN); mpfr_mul(r11192, r11154, r11153, MPFR_RNDN); mpfr_add(r11193, r11191, r11192, MPFR_RNDN); mpfr_mul(r11194, r11162, r11164, MPFR_RNDN); mpfr_mul(r11195, r11194, r11163, MPFR_RNDN); mpfr_sub(r11196, r11193, r11195, MPFR_RNDN); mpfr_mul(r11197, r11157, r11156, MPFR_RNDN); mpfr_mul(r11198, r11197, r11158, MPFR_RNDN); mpfr_sub(r11199, r11196, r11198, MPFR_RNDN); ; mpfr_set_si(r11201, mpfr_cmp(r11183, r11200) <= 0, MPFR_RNDN); mpfr_mul(r11202, r11168, r11169, MPFR_RNDN); mpfr_mul(r11203, r11202, r11170, MPFR_RNDN); mpfr_sub(r11204, r11165, r11203, MPFR_RNDN); mpfr_mul(r11205, r11162, r11204, MPFR_RNDN); mpfr_sub(r11206, r11161, r11205, MPFR_RNDN); mpfr_add(r11207, r11206, r11182, MPFR_RNDN); if (mpfr_get_si(r11201, MPFR_RNDN)) { mpfr_set(r11208, r11207, MPFR_RNDN); } else { mpfr_set(r11208, r11199, MPFR_RNDN); }; if (mpfr_get_si(r11185, MPFR_RNDN)) { mpfr_set(r11209, r11199, MPFR_RNDN); } else { mpfr_set(r11209, r11208, MPFR_RNDN); }; return mpfr_get_d(r11209, MPFR_RNDN); }