#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 r26929 = x; float r26930 = 18.0; float r26931 = r26929 * r26930; float r26932 = y; float r26933 = r26931 * r26932; float r26934 = z; float r26935 = r26933 * r26934; float r26936 = t; float r26937 = r26935 * r26936; float r26938 = a; float r26939 = 4.0; float r26940 = r26938 * r26939; float r26941 = r26940 * r26936; float r26942 = r26937 - r26941; float r26943 = b; float r26944 = c; float r26945 = r26943 * r26944; float r26946 = r26942 + r26945; float r26947 = r26929 * r26939; float r26948 = i; float r26949 = r26947 * r26948; float r26950 = r26946 - r26949; float r26951 = j; float r26952 = 27.0; float r26953 = r26951 * r26952; float r26954 = k; float r26955 = r26953 * r26954; float r26956 = r26950 - r26955; return r26956; } double f_id(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) { double r26957 = x; double r26958 = 18.0; double r26959 = r26957 * r26958; double r26960 = y; double r26961 = r26959 * r26960; double r26962 = z; double r26963 = r26961 * r26962; double r26964 = t; double r26965 = r26963 * r26964; double r26966 = a; double r26967 = 4.0; double r26968 = r26966 * r26967; double r26969 = r26968 * r26964; double r26970 = r26965 - r26969; double r26971 = b; double r26972 = c; double r26973 = r26971 * r26972; double r26974 = r26970 + r26973; double r26975 = r26957 * r26967; double r26976 = i; double r26977 = r26975 * r26976; double r26978 = r26974 - r26977; double r26979 = j; double r26980 = 27.0; double r26981 = r26979 * r26980; double r26982 = k; double r26983 = r26981 * r26982; double r26984 = r26978 - r26983; return r26984; } double f_of(float x, float y, float z, float t, float a, float b, float c, float i, float j, float k) { float r26985 = c; float r26986 = b; float r26987 = r26985 * r26986; float r26988 = i; float r26989 = x; float r26990 = r26988 * r26989; float r26991 = 4.0; float r26992 = r26990 * r26991; float r26993 = j; float r26994 = 27.0; float r26995 = k; float r26996 = r26994 * r26995; float r26997 = r26993 * r26996; float r26998 = r26992 + r26997; float r26999 = z; float r27000 = r26999 * r26989; float r27001 = 18.0; float r27002 = t; float r27003 = y; float r27004 = r27002 * r27003; float r27005 = r27001 * r27004; float r27006 = r27000 * r27005; float r27007 = r26998 - r27006; float r27008 = r26987 - r27007; float r27009 = a; float r27010 = r26991 * r27009; float r27011 = -r27010; float r27012 = r27002 * r27011; float r27013 = r27008 + r27012; float r27014 = -inf.0; bool r27015 = r27013 <= r27014; float r27016 = r27002 * r26989; float r27017 = r27003 * r27016; float r27018 = r26999 * r27017; float r27019 = r27001 * r27018; float r27020 = r27009 * r27002; float r27021 = r26991 * r27020; float r27022 = r26993 * r26995; float r27023 = r26994 * r27022; float r27024 = r27021 + r27023; float r27025 = r27019 - r27024; float r27026 = 1.6747781086077075e+293; bool r27027 = r27013 <= r27026; float r27028 = r26995 * r26993; float r27029 = r27028 * r26994; float r27030 = r26987 - r27029; float r27031 = r26988 * r26991; float r27032 = r26989 * r27031; float r27033 = r27002 * r27001; float r27034 = r26999 * r27003; float r27035 = r27033 * r27034; float r27036 = -r27035; float r27037 = r26989 * r27036; float r27038 = r27032 + r27037; float r27039 = r27030 - r27038; float r27040 = r27039 + r27012; float r27041 = r27027 ? r27013 : r27040; float r27042 = r27015 ? r27025 : r27041; return r27042; } double f_od(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) { double r27043 = c; double r27044 = b; double r27045 = r27043 * r27044; double r27046 = i; double r27047 = x; double r27048 = r27046 * r27047; double r27049 = 4.0; double r27050 = r27048 * r27049; double r27051 = j; double r27052 = 27.0; double r27053 = k; double r27054 = r27052 * r27053; double r27055 = r27051 * r27054; double r27056 = r27050 + r27055; double r27057 = z; double r27058 = r27057 * r27047; double r27059 = 18.0; double r27060 = t; double r27061 = y; double r27062 = r27060 * r27061; double r27063 = r27059 * r27062; double r27064 = r27058 * r27063; double r27065 = r27056 - r27064; double r27066 = r27045 - r27065; double r27067 = a; double r27068 = r27049 * r27067; double r27069 = -r27068; double r27070 = r27060 * r27069; double r27071 = r27066 + r27070; double r27072 = -inf.0; bool r27073 = r27071 <= r27072; double r27074 = r27060 * r27047; double r27075 = r27061 * r27074; double r27076 = r27057 * r27075; double r27077 = r27059 * r27076; double r27078 = r27067 * r27060; double r27079 = r27049 * r27078; double r27080 = r27051 * r27053; double r27081 = r27052 * r27080; double r27082 = r27079 + r27081; double r27083 = r27077 - r27082; double r27084 = 1.6747781086077075e+293; bool r27085 = r27071 <= r27084; double r27086 = r27053 * r27051; double r27087 = r27086 * r27052; double r27088 = r27045 - r27087; double r27089 = r27046 * r27049; double r27090 = r27047 * r27089; double r27091 = r27060 * r27059; double r27092 = r27057 * r27061; double r27093 = r27091 * r27092; double r27094 = -r27093; double r27095 = r27047 * r27094; double r27096 = r27090 + r27095; double r27097 = r27088 - r27096; double r27098 = r27097 + r27070; double r27099 = r27085 ? r27071 : r27098; double r27100 = r27073 ? r27083 : r27099; return r27100; } 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 r27101, r27102, r27103, r27104, r27105, r27106, r27107, r27108, r27109, r27110, r27111, r27112, r27113, r27114, r27115, r27116, r27117, r27118, r27119, r27120, r27121, r27122, r27123, r27124, r27125, r27126, r27127, r27128; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r27101); mpfr_init_set_str(r27102, "18.0", 10, MPFR_RNDN); mpfr_init(r27103); mpfr_init(r27104); mpfr_init(r27105); mpfr_init(r27106); mpfr_init(r27107); mpfr_init(r27108); mpfr_init(r27109); mpfr_init(r27110); mpfr_init_set_str(r27111, "4.0", 10, MPFR_RNDN); mpfr_init(r27112); mpfr_init(r27113); mpfr_init(r27114); mpfr_init(r27115); mpfr_init(r27116); mpfr_init(r27117); mpfr_init(r27118); mpfr_init(r27119); mpfr_init(r27120); mpfr_init(r27121); mpfr_init(r27122); mpfr_init(r27123); mpfr_init_set_str(r27124, "27.0", 10, MPFR_RNDN); mpfr_init(r27125); mpfr_init(r27126); mpfr_init(r27127); mpfr_init(r27128); } 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(r27101, x, MPFR_RNDN); ; mpfr_mul(r27103, r27101, r27102, MPFR_RNDN); mpfr_set_d(r27104, y, MPFR_RNDN); mpfr_mul(r27105, r27103, r27104, MPFR_RNDN); mpfr_set_d(r27106, z, MPFR_RNDN); mpfr_mul(r27107, r27105, r27106, MPFR_RNDN); mpfr_set_d(r27108, t, MPFR_RNDN); mpfr_mul(r27109, r27107, r27108, MPFR_RNDN); mpfr_set_d(r27110, a, MPFR_RNDN); ; mpfr_mul(r27112, r27110, r27111, MPFR_RNDN); mpfr_mul(r27113, r27112, r27108, MPFR_RNDN); mpfr_sub(r27114, r27109, r27113, MPFR_RNDN); mpfr_set_d(r27115, b, MPFR_RNDN); mpfr_set_d(r27116, c, MPFR_RNDN); mpfr_mul(r27117, r27115, r27116, MPFR_RNDN); mpfr_add(r27118, r27114, r27117, MPFR_RNDN); mpfr_mul(r27119, r27101, r27111, MPFR_RNDN); mpfr_set_d(r27120, i, MPFR_RNDN); mpfr_mul(r27121, r27119, r27120, MPFR_RNDN); mpfr_sub(r27122, r27118, r27121, MPFR_RNDN); mpfr_set_d(r27123, j, MPFR_RNDN); ; mpfr_mul(r27125, r27123, r27124, MPFR_RNDN); mpfr_set_d(r27126, k, MPFR_RNDN); mpfr_mul(r27127, r27125, r27126, MPFR_RNDN); mpfr_sub(r27128, r27122, r27127, MPFR_RNDN); return mpfr_get_d(r27128, MPFR_RNDN); } static mpfr_t r27129, r27130, r27131, r27132, r27133, r27134, r27135, r27136, r27137, r27138, r27139, r27140, r27141, r27142, r27143, r27144, r27145, r27146, r27147, r27148, r27149, r27150, r27151, r27152, r27153, r27154, r27155, r27156, r27157, r27158, r27159, r27160, r27161, r27162, r27163, r27164, r27165, r27166, r27167, r27168, r27169, r27170, r27171, r27172, r27173, r27174, r27175, r27176, r27177, r27178, r27179, r27180, r27181, r27182, r27183, r27184, r27185, r27186; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r27129); mpfr_init(r27130); mpfr_init(r27131); mpfr_init(r27132); mpfr_init(r27133); mpfr_init(r27134); mpfr_init_set_str(r27135, "4.0", 10, MPFR_RNDN); mpfr_init(r27136); mpfr_init(r27137); mpfr_init_set_str(r27138, "27.0", 10, MPFR_RNDN); mpfr_init(r27139); mpfr_init(r27140); mpfr_init(r27141); mpfr_init(r27142); mpfr_init(r27143); mpfr_init(r27144); mpfr_init_set_str(r27145, "18.0", 10, MPFR_RNDN); mpfr_init(r27146); mpfr_init(r27147); mpfr_init(r27148); mpfr_init(r27149); mpfr_init(r27150); mpfr_init(r27151); mpfr_init(r27152); mpfr_init(r27153); mpfr_init(r27154); mpfr_init(r27155); mpfr_init(r27156); mpfr_init(r27157); mpfr_init_set_str(r27158, "-inf.0", 10, MPFR_RNDN); mpfr_init(r27159); mpfr_init(r27160); mpfr_init(r27161); mpfr_init(r27162); mpfr_init(r27163); mpfr_init(r27164); mpfr_init(r27165); mpfr_init(r27166); mpfr_init(r27167); mpfr_init(r27168); mpfr_init(r27169); mpfr_init_set_str(r27170, "1.6747781086077075e+293", 10, MPFR_RNDN); mpfr_init(r27171); mpfr_init(r27172); mpfr_init(r27173); mpfr_init(r27174); mpfr_init(r27175); mpfr_init(r27176); mpfr_init(r27177); mpfr_init(r27178); mpfr_init(r27179); mpfr_init(r27180); mpfr_init(r27181); mpfr_init(r27182); mpfr_init(r27183); mpfr_init(r27184); mpfr_init(r27185); mpfr_init(r27186); } 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(r27129, c, MPFR_RNDN); mpfr_set_d(r27130, b, MPFR_RNDN); mpfr_mul(r27131, r27129, r27130, MPFR_RNDN); mpfr_set_d(r27132, i, MPFR_RNDN); mpfr_set_d(r27133, x, MPFR_RNDN); mpfr_mul(r27134, r27132, r27133, MPFR_RNDN); ; mpfr_mul(r27136, r27134, r27135, MPFR_RNDN); mpfr_set_d(r27137, j, MPFR_RNDN); ; mpfr_set_d(r27139, k, MPFR_RNDN); mpfr_mul(r27140, r27138, r27139, MPFR_RNDN); mpfr_mul(r27141, r27137, r27140, MPFR_RNDN); mpfr_add(r27142, r27136, r27141, MPFR_RNDN); mpfr_set_d(r27143, z, MPFR_RNDN); mpfr_mul(r27144, r27143, r27133, MPFR_RNDN); ; mpfr_set_d(r27146, t, MPFR_RNDN); mpfr_set_d(r27147, y, MPFR_RNDN); mpfr_mul(r27148, r27146, r27147, MPFR_RNDN); mpfr_mul(r27149, r27145, r27148, MPFR_RNDN); mpfr_mul(r27150, r27144, r27149, MPFR_RNDN); mpfr_sub(r27151, r27142, r27150, MPFR_RNDN); mpfr_sub(r27152, r27131, r27151, MPFR_RNDN); mpfr_set_d(r27153, a, MPFR_RNDN); mpfr_mul(r27154, r27135, r27153, MPFR_RNDN); mpfr_neg(r27155, r27154, MPFR_RNDN); mpfr_mul(r27156, r27146, r27155, MPFR_RNDN); mpfr_add(r27157, r27152, r27156, MPFR_RNDN); ; mpfr_set_si(r27159, mpfr_cmp(r27157, r27158) <= 0, MPFR_RNDN); mpfr_mul(r27160, r27146, r27133, MPFR_RNDN); mpfr_mul(r27161, r27147, r27160, MPFR_RNDN); mpfr_mul(r27162, r27143, r27161, MPFR_RNDN); mpfr_mul(r27163, r27145, r27162, MPFR_RNDN); mpfr_mul(r27164, r27153, r27146, MPFR_RNDN); mpfr_mul(r27165, r27135, r27164, MPFR_RNDN); mpfr_mul(r27166, r27137, r27139, MPFR_RNDN); mpfr_mul(r27167, r27138, r27166, MPFR_RNDN); mpfr_add(r27168, r27165, r27167, MPFR_RNDN); mpfr_sub(r27169, r27163, r27168, MPFR_RNDN); ; mpfr_set_si(r27171, mpfr_cmp(r27157, r27170) <= 0, MPFR_RNDN); mpfr_mul(r27172, r27139, r27137, MPFR_RNDN); mpfr_mul(r27173, r27172, r27138, MPFR_RNDN); mpfr_sub(r27174, r27131, r27173, MPFR_RNDN); mpfr_mul(r27175, r27132, r27135, MPFR_RNDN); mpfr_mul(r27176, r27133, r27175, MPFR_RNDN); mpfr_mul(r27177, r27146, r27145, MPFR_RNDN); mpfr_mul(r27178, r27143, r27147, MPFR_RNDN); mpfr_mul(r27179, r27177, r27178, MPFR_RNDN); mpfr_neg(r27180, r27179, MPFR_RNDN); mpfr_mul(r27181, r27133, r27180, MPFR_RNDN); mpfr_add(r27182, r27176, r27181, MPFR_RNDN); mpfr_sub(r27183, r27174, r27182, MPFR_RNDN); mpfr_add(r27184, r27183, r27156, MPFR_RNDN); if (mpfr_get_si(r27171, MPFR_RNDN)) { mpfr_set(r27185, r27157, MPFR_RNDN); } else { mpfr_set(r27185, r27184, MPFR_RNDN); }; if (mpfr_get_si(r27159, MPFR_RNDN)) { mpfr_set(r27186, r27169, MPFR_RNDN); } else { mpfr_set(r27186, r27185, MPFR_RNDN); }; return mpfr_get_d(r27186, MPFR_RNDN); } static mpfr_t r27187, r27188, r27189, r27190, r27191, r27192, r27193, r27194, r27195, r27196, r27197, r27198, r27199, r27200, r27201, r27202, r27203, r27204, r27205, r27206, r27207, r27208, r27209, r27210, r27211, r27212, r27213, r27214, r27215, r27216, r27217, r27218, r27219, r27220, r27221, r27222, r27223, r27224, r27225, r27226, r27227, r27228, r27229, r27230, r27231, r27232, r27233, r27234, r27235, r27236, r27237, r27238, r27239, r27240, r27241, r27242, r27243, r27244; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r27187); mpfr_init(r27188); mpfr_init(r27189); mpfr_init(r27190); mpfr_init(r27191); mpfr_init(r27192); mpfr_init_set_str(r27193, "4.0", 10, MPFR_RNDN); mpfr_init(r27194); mpfr_init(r27195); mpfr_init_set_str(r27196, "27.0", 10, MPFR_RNDN); mpfr_init(r27197); mpfr_init(r27198); mpfr_init(r27199); mpfr_init(r27200); mpfr_init(r27201); mpfr_init(r27202); mpfr_init_set_str(r27203, "18.0", 10, MPFR_RNDN); mpfr_init(r27204); mpfr_init(r27205); mpfr_init(r27206); mpfr_init(r27207); mpfr_init(r27208); mpfr_init(r27209); mpfr_init(r27210); mpfr_init(r27211); mpfr_init(r27212); mpfr_init(r27213); mpfr_init(r27214); mpfr_init(r27215); mpfr_init_set_str(r27216, "-inf.0", 10, MPFR_RNDN); mpfr_init(r27217); mpfr_init(r27218); mpfr_init(r27219); mpfr_init(r27220); mpfr_init(r27221); mpfr_init(r27222); mpfr_init(r27223); mpfr_init(r27224); mpfr_init(r27225); mpfr_init(r27226); mpfr_init(r27227); mpfr_init_set_str(r27228, "1.6747781086077075e+293", 10, MPFR_RNDN); mpfr_init(r27229); mpfr_init(r27230); mpfr_init(r27231); mpfr_init(r27232); mpfr_init(r27233); mpfr_init(r27234); mpfr_init(r27235); mpfr_init(r27236); mpfr_init(r27237); mpfr_init(r27238); mpfr_init(r27239); mpfr_init(r27240); mpfr_init(r27241); mpfr_init(r27242); mpfr_init(r27243); mpfr_init(r27244); } 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(r27187, c, MPFR_RNDN); mpfr_set_d(r27188, b, MPFR_RNDN); mpfr_mul(r27189, r27187, r27188, MPFR_RNDN); mpfr_set_d(r27190, i, MPFR_RNDN); mpfr_set_d(r27191, x, MPFR_RNDN); mpfr_mul(r27192, r27190, r27191, MPFR_RNDN); ; mpfr_mul(r27194, r27192, r27193, MPFR_RNDN); mpfr_set_d(r27195, j, MPFR_RNDN); ; mpfr_set_d(r27197, k, MPFR_RNDN); mpfr_mul(r27198, r27196, r27197, MPFR_RNDN); mpfr_mul(r27199, r27195, r27198, MPFR_RNDN); mpfr_add(r27200, r27194, r27199, MPFR_RNDN); mpfr_set_d(r27201, z, MPFR_RNDN); mpfr_mul(r27202, r27201, r27191, MPFR_RNDN); ; mpfr_set_d(r27204, t, MPFR_RNDN); mpfr_set_d(r27205, y, MPFR_RNDN); mpfr_mul(r27206, r27204, r27205, MPFR_RNDN); mpfr_mul(r27207, r27203, r27206, MPFR_RNDN); mpfr_mul(r27208, r27202, r27207, MPFR_RNDN); mpfr_sub(r27209, r27200, r27208, MPFR_RNDN); mpfr_sub(r27210, r27189, r27209, MPFR_RNDN); mpfr_set_d(r27211, a, MPFR_RNDN); mpfr_mul(r27212, r27193, r27211, MPFR_RNDN); mpfr_neg(r27213, r27212, MPFR_RNDN); mpfr_mul(r27214, r27204, r27213, MPFR_RNDN); mpfr_add(r27215, r27210, r27214, MPFR_RNDN); ; mpfr_set_si(r27217, mpfr_cmp(r27215, r27216) <= 0, MPFR_RNDN); mpfr_mul(r27218, r27204, r27191, MPFR_RNDN); mpfr_mul(r27219, r27205, r27218, MPFR_RNDN); mpfr_mul(r27220, r27201, r27219, MPFR_RNDN); mpfr_mul(r27221, r27203, r27220, MPFR_RNDN); mpfr_mul(r27222, r27211, r27204, MPFR_RNDN); mpfr_mul(r27223, r27193, r27222, MPFR_RNDN); mpfr_mul(r27224, r27195, r27197, MPFR_RNDN); mpfr_mul(r27225, r27196, r27224, MPFR_RNDN); mpfr_add(r27226, r27223, r27225, MPFR_RNDN); mpfr_sub(r27227, r27221, r27226, MPFR_RNDN); ; mpfr_set_si(r27229, mpfr_cmp(r27215, r27228) <= 0, MPFR_RNDN); mpfr_mul(r27230, r27197, r27195, MPFR_RNDN); mpfr_mul(r27231, r27230, r27196, MPFR_RNDN); mpfr_sub(r27232, r27189, r27231, MPFR_RNDN); mpfr_mul(r27233, r27190, r27193, MPFR_RNDN); mpfr_mul(r27234, r27191, r27233, MPFR_RNDN); mpfr_mul(r27235, r27204, r27203, MPFR_RNDN); mpfr_mul(r27236, r27201, r27205, MPFR_RNDN); mpfr_mul(r27237, r27235, r27236, MPFR_RNDN); mpfr_neg(r27238, r27237, MPFR_RNDN); mpfr_mul(r27239, r27191, r27238, MPFR_RNDN); mpfr_add(r27240, r27234, r27239, MPFR_RNDN); mpfr_sub(r27241, r27232, r27240, MPFR_RNDN); mpfr_add(r27242, r27241, r27214, MPFR_RNDN); if (mpfr_get_si(r27229, MPFR_RNDN)) { mpfr_set(r27243, r27215, MPFR_RNDN); } else { mpfr_set(r27243, r27242, MPFR_RNDN); }; if (mpfr_get_si(r27217, MPFR_RNDN)) { mpfr_set(r27244, r27227, MPFR_RNDN); } else { mpfr_set(r27244, r27243, MPFR_RNDN); }; return mpfr_get_d(r27244, MPFR_RNDN); }