#include #include #include #include #include char *name = "jeff quadratic root 1"; double f_if(float a, float b, float c) { float r20953 = b; float r20954 = 0; bool r20955 = r20953 >= r20954; float r20956 = -r20953; float r20957 = r20953 * r20953; float r20958 = 4; float r20959 = a; float r20960 = r20958 * r20959; float r20961 = c; float r20962 = r20960 * r20961; float r20963 = r20957 - r20962; float r20964 = sqrt(r20963); float r20965 = r20956 - r20964; float r20966 = 2; float r20967 = r20966 * r20959; float r20968 = r20965 / r20967; float r20969 = r20966 * r20961; float r20970 = r20956 + r20964; float r20971 = r20969 / r20970; float r20972 = r20955 ? r20968 : r20971; return r20972; } double f_id(double a, double b, double c) { double r20973 = b; double r20974 = 0; bool r20975 = r20973 >= r20974; double r20976 = -r20973; double r20977 = r20973 * r20973; double r20978 = 4; double r20979 = a; double r20980 = r20978 * r20979; double r20981 = c; double r20982 = r20980 * r20981; double r20983 = r20977 - r20982; double r20984 = sqrt(r20983); double r20985 = r20976 - r20984; double r20986 = 2; double r20987 = r20986 * r20979; double r20988 = r20985 / r20987; double r20989 = r20986 * r20981; double r20990 = r20976 + r20984; double r20991 = r20989 / r20990; double r20992 = r20975 ? r20988 : r20991; return r20992; } double f_of(float a, float b, float c) { float r20993 = b; float r20994 = -7.853185861915354e+117; bool r20995 = r20993 <= r20994; float r20996 = 0; bool r20997 = r20993 >= r20996; float r20998 = -r20993; float r20999 = r20993 * r20993; float r21000 = c; float r21001 = a; float r21002 = 4; float r21003 = r21001 * r21002; float r21004 = r21000 * r21003; float r21005 = r20999 - r21004; float r21006 = sqrt(r21005); float r21007 = r20998 - r21006; float r21008 = r21001 + r21001; float r21009 = r21007 / r21008; float r21010 = r21000 + r21000; float r21011 = 2; float r21012 = r21010 / r21011; float r21013 = r21000 / r20993; float r21014 = fma(r21001, r21013, r20998); float r21015 = r21012 / r21014; float r21016 = r20997 ? r21009 : r21015; float r21017 = -4.086812811986243e-293; bool r21018 = r20993 <= r21017; float r21019 = r21002 * r21001; float r21020 = r21019 * r21000; float r21021 = r20999 - r21020; float r21022 = sqrt(r21021); float r21023 = r20998 + r21022; float r21024 = r21004 / r21023; float r21025 = r21011 * r21001; float r21026 = r21024 / r21025; float r21027 = r21011 * r21000; float r21028 = r21027 / r21023; float r21029 = r20997 ? r21026 : r21028; float r21030 = 1.8918781781269702e+105; bool r21031 = r20993 <= r21030; float r21032 = r20998 - r21022; float r21033 = r21032 / r21025; float r21034 = r21001 * r21000; float r21035 = r21034 * r21002; float r21036 = r20999 - r21035; float r21037 = sqrt(r21036); float r21038 = r21037 * r21036; float r21039 = cbrt(r21038); float r21040 = r20998 + r21039; float r21041 = r21027 / r21040; float r21042 = r20997 ? r21033 : r21041; float r21043 = 1; float r21044 = r21013 / r21043; float r21045 = r20993 / r21001; float r21046 = r21044 - r21045; float r21047 = r21006 + r20998; float r21048 = r21010 / r21047; float r21049 = r20997 ? r21046 : r21048; float r21050 = r21031 ? r21042 : r21049; float r21051 = r21018 ? r21029 : r21050; float r21052 = r20995 ? r21016 : r21051; return r21052; } double f_od(double a, double b, double c) { double r21053 = b; double r21054 = -7.853185861915354e+117; bool r21055 = r21053 <= r21054; double r21056 = 0; bool r21057 = r21053 >= r21056; double r21058 = -r21053; double r21059 = r21053 * r21053; double r21060 = c; double r21061 = a; double r21062 = 4; double r21063 = r21061 * r21062; double r21064 = r21060 * r21063; double r21065 = r21059 - r21064; double r21066 = sqrt(r21065); double r21067 = r21058 - r21066; double r21068 = r21061 + r21061; double r21069 = r21067 / r21068; double r21070 = r21060 + r21060; double r21071 = 2; double r21072 = r21070 / r21071; double r21073 = r21060 / r21053; double r21074 = fma(r21061, r21073, r21058); double r21075 = r21072 / r21074; double r21076 = r21057 ? r21069 : r21075; double r21077 = -4.086812811986243e-293; bool r21078 = r21053 <= r21077; double r21079 = r21062 * r21061; double r21080 = r21079 * r21060; double r21081 = r21059 - r21080; double r21082 = sqrt(r21081); double r21083 = r21058 + r21082; double r21084 = r21064 / r21083; double r21085 = r21071 * r21061; double r21086 = r21084 / r21085; double r21087 = r21071 * r21060; double r21088 = r21087 / r21083; double r21089 = r21057 ? r21086 : r21088; double r21090 = 1.8918781781269702e+105; bool r21091 = r21053 <= r21090; double r21092 = r21058 - r21082; double r21093 = r21092 / r21085; double r21094 = r21061 * r21060; double r21095 = r21094 * r21062; double r21096 = r21059 - r21095; double r21097 = sqrt(r21096); double r21098 = r21097 * r21096; double r21099 = cbrt(r21098); double r21100 = r21058 + r21099; double r21101 = r21087 / r21100; double r21102 = r21057 ? r21093 : r21101; double r21103 = 1; double r21104 = r21073 / r21103; double r21105 = r21053 / r21061; double r21106 = r21104 - r21105; double r21107 = r21066 + r21058; double r21108 = r21070 / r21107; double r21109 = r21057 ? r21106 : r21108; double r21110 = r21091 ? r21102 : r21109; double r21111 = r21078 ? r21089 : r21110; double r21112 = r21055 ? r21076 : r21111; return r21112; } 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 r21113, r21114, r21115, r21116, r21117, r21118, r21119, r21120, r21121, r21122, r21123, r21124, r21125, r21126, r21127, r21128, r21129, r21130, r21131, r21132; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r21113); mpfr_init_set_str(r21114, "0", 10, MPFR_RNDN); mpfr_init(r21115); mpfr_init(r21116); mpfr_init(r21117); mpfr_init_set_str(r21118, "4", 10, MPFR_RNDN); mpfr_init(r21119); mpfr_init(r21120); mpfr_init(r21121); mpfr_init(r21122); mpfr_init(r21123); mpfr_init(r21124); mpfr_init(r21125); mpfr_init_set_str(r21126, "2", 10, MPFR_RNDN); mpfr_init(r21127); mpfr_init(r21128); mpfr_init(r21129); mpfr_init(r21130); mpfr_init(r21131); mpfr_init(r21132); } double f_im(double a, double b, double c) { mpfr_set_d(r21113, b, MPFR_RNDN); ; mpfr_set_si(r21115, mpfr_cmp(r21113, r21114) >= 0, MPFR_RNDN); mpfr_neg(r21116, r21113, MPFR_RNDN); mpfr_mul(r21117, r21113, r21113, MPFR_RNDN); ; mpfr_set_d(r21119, a, MPFR_RNDN); mpfr_mul(r21120, r21118, r21119, MPFR_RNDN); mpfr_set_d(r21121, c, MPFR_RNDN); mpfr_mul(r21122, r21120, r21121, MPFR_RNDN); mpfr_sub(r21123, r21117, r21122, MPFR_RNDN); mpfr_sqrt(r21124, r21123, MPFR_RNDN); mpfr_sub(r21125, r21116, r21124, MPFR_RNDN); ; mpfr_mul(r21127, r21126, r21119, MPFR_RNDN); mpfr_div(r21128, r21125, r21127, MPFR_RNDN); mpfr_mul(r21129, r21126, r21121, MPFR_RNDN); mpfr_add(r21130, r21116, r21124, MPFR_RNDN); mpfr_div(r21131, r21129, r21130, MPFR_RNDN); if (mpfr_get_si(r21115, MPFR_RNDN)) { mpfr_set(r21132, r21128, MPFR_RNDN); } else { mpfr_set(r21132, r21131, MPFR_RNDN); }; return mpfr_get_d(r21132, MPFR_RNDN); } static mpfr_t r21133, r21134, r21135, r21136, r21137, r21138, r21139, r21140, r21141, r21142, r21143, r21144, r21145, r21146, r21147, r21148, r21149, r21150, r21151, r21152, r21153, r21154, r21155, r21156, r21157, r21158, r21159, r21160, r21161, r21162, r21163, r21164, r21165, r21166, r21167, r21168, r21169, r21170, r21171, r21172, r21173, r21174, r21175, r21176, r21177, r21178, r21179, r21180, r21181, r21182, r21183, r21184, r21185, r21186, r21187, r21188, r21189, r21190, r21191, r21192; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r21133); mpfr_init_set_str(r21134, "-7.853185861915354e+117", 10, MPFR_RNDN); mpfr_init(r21135); mpfr_init_set_str(r21136, "0", 10, MPFR_RNDN); mpfr_init(r21137); mpfr_init(r21138); mpfr_init(r21139); mpfr_init(r21140); mpfr_init(r21141); mpfr_init_set_str(r21142, "4", 10, MPFR_RNDN); mpfr_init(r21143); mpfr_init(r21144); mpfr_init(r21145); mpfr_init(r21146); mpfr_init(r21147); mpfr_init(r21148); mpfr_init(r21149); mpfr_init(r21150); mpfr_init_set_str(r21151, "2", 10, MPFR_RNDN); mpfr_init(r21152); mpfr_init(r21153); mpfr_init(r21154); mpfr_init(r21155); mpfr_init(r21156); mpfr_init_set_str(r21157, "-4.086812811986243e-293", 10, MPFR_RNDN); mpfr_init(r21158); mpfr_init(r21159); mpfr_init(r21160); mpfr_init(r21161); mpfr_init(r21162); mpfr_init(r21163); mpfr_init(r21164); mpfr_init(r21165); mpfr_init(r21166); mpfr_init(r21167); mpfr_init(r21168); mpfr_init(r21169); mpfr_init_set_str(r21170, "1.8918781781269702e+105", 10, MPFR_RNDN); mpfr_init(r21171); mpfr_init(r21172); mpfr_init(r21173); mpfr_init(r21174); mpfr_init(r21175); mpfr_init(r21176); mpfr_init(r21177); mpfr_init(r21178); mpfr_init(r21179); mpfr_init(r21180); mpfr_init(r21181); mpfr_init(r21182); mpfr_init_set_str(r21183, "1", 10, MPFR_RNDN); mpfr_init(r21184); mpfr_init(r21185); mpfr_init(r21186); mpfr_init(r21187); mpfr_init(r21188); mpfr_init(r21189); mpfr_init(r21190); mpfr_init(r21191); mpfr_init(r21192); } double f_fm(double a, double b, double c) { mpfr_set_d(r21133, b, MPFR_RNDN); ; mpfr_set_si(r21135, mpfr_cmp(r21133, r21134) <= 0, MPFR_RNDN); ; mpfr_set_si(r21137, mpfr_cmp(r21133, r21136) >= 0, MPFR_RNDN); mpfr_neg(r21138, r21133, MPFR_RNDN); mpfr_mul(r21139, r21133, r21133, MPFR_RNDN); mpfr_set_d(r21140, c, MPFR_RNDN); mpfr_set_d(r21141, a, MPFR_RNDN); ; mpfr_mul(r21143, r21141, r21142, MPFR_RNDN); mpfr_mul(r21144, r21140, r21143, MPFR_RNDN); mpfr_sub(r21145, r21139, r21144, MPFR_RNDN); mpfr_sqrt(r21146, r21145, MPFR_RNDN); mpfr_sub(r21147, r21138, r21146, MPFR_RNDN); mpfr_add(r21148, r21141, r21141, MPFR_RNDN); mpfr_div(r21149, r21147, r21148, MPFR_RNDN); mpfr_add(r21150, r21140, r21140, MPFR_RNDN); ; mpfr_div(r21152, r21150, r21151, MPFR_RNDN); mpfr_div(r21153, r21140, r21133, MPFR_RNDN); mpfr_fma(r21154, r21141, r21153, r21138, MPFR_RNDN); mpfr_div(r21155, r21152, r21154, MPFR_RNDN); if (mpfr_get_si(r21137, MPFR_RNDN)) { mpfr_set(r21156, r21149, MPFR_RNDN); } else { mpfr_set(r21156, r21155, MPFR_RNDN); }; ; mpfr_set_si(r21158, mpfr_cmp(r21133, r21157) <= 0, MPFR_RNDN); mpfr_mul(r21159, r21142, r21141, MPFR_RNDN); mpfr_mul(r21160, r21159, r21140, MPFR_RNDN); mpfr_sub(r21161, r21139, r21160, MPFR_RNDN); mpfr_sqrt(r21162, r21161, MPFR_RNDN); mpfr_add(r21163, r21138, r21162, MPFR_RNDN); mpfr_div(r21164, r21144, r21163, MPFR_RNDN); mpfr_mul(r21165, r21151, r21141, MPFR_RNDN); mpfr_div(r21166, r21164, r21165, MPFR_RNDN); mpfr_mul(r21167, r21151, r21140, MPFR_RNDN); mpfr_div(r21168, r21167, r21163, MPFR_RNDN); if (mpfr_get_si(r21137, MPFR_RNDN)) { mpfr_set(r21169, r21166, MPFR_RNDN); } else { mpfr_set(r21169, r21168, MPFR_RNDN); }; ; mpfr_set_si(r21171, mpfr_cmp(r21133, r21170) <= 0, MPFR_RNDN); mpfr_sub(r21172, r21138, r21162, MPFR_RNDN); mpfr_div(r21173, r21172, r21165, MPFR_RNDN); mpfr_mul(r21174, r21141, r21140, MPFR_RNDN); mpfr_mul(r21175, r21174, r21142, MPFR_RNDN); mpfr_sub(r21176, r21139, r21175, MPFR_RNDN); mpfr_sqrt(r21177, r21176, MPFR_RNDN); mpfr_mul(r21178, r21177, r21176, MPFR_RNDN); mpfr_cbrt(r21179, r21178, MPFR_RNDN); mpfr_add(r21180, r21138, r21179, MPFR_RNDN); mpfr_div(r21181, r21167, r21180, MPFR_RNDN); if (mpfr_get_si(r21137, MPFR_RNDN)) { mpfr_set(r21182, r21173, MPFR_RNDN); } else { mpfr_set(r21182, r21181, MPFR_RNDN); }; ; mpfr_div(r21184, r21153, r21183, MPFR_RNDN); mpfr_div(r21185, r21133, r21141, MPFR_RNDN); mpfr_sub(r21186, r21184, r21185, MPFR_RNDN); mpfr_add(r21187, r21146, r21138, MPFR_RNDN); mpfr_div(r21188, r21150, r21187, MPFR_RNDN); if (mpfr_get_si(r21137, MPFR_RNDN)) { mpfr_set(r21189, r21186, MPFR_RNDN); } else { mpfr_set(r21189, r21188, MPFR_RNDN); }; if (mpfr_get_si(r21171, MPFR_RNDN)) { mpfr_set(r21190, r21182, MPFR_RNDN); } else { mpfr_set(r21190, r21189, MPFR_RNDN); }; if (mpfr_get_si(r21158, MPFR_RNDN)) { mpfr_set(r21191, r21169, MPFR_RNDN); } else { mpfr_set(r21191, r21190, MPFR_RNDN); }; if (mpfr_get_si(r21135, MPFR_RNDN)) { mpfr_set(r21192, r21156, MPFR_RNDN); } else { mpfr_set(r21192, r21191, MPFR_RNDN); }; return mpfr_get_d(r21192, MPFR_RNDN); } static mpfr_t r21193, r21194, r21195, r21196, r21197, r21198, r21199, r21200, r21201, r21202, r21203, r21204, r21205, r21206, r21207, r21208, r21209, r21210, r21211, r21212, r21213, r21214, r21215, r21216, r21217, r21218, r21219, r21220, r21221, r21222, r21223, r21224, r21225, r21226, r21227, r21228, r21229, r21230, r21231, r21232, r21233, r21234, r21235, r21236, r21237, r21238, r21239, r21240, r21241, r21242, r21243, r21244, r21245, r21246, r21247, r21248, r21249, r21250, r21251, r21252; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r21193); mpfr_init_set_str(r21194, "-7.853185861915354e+117", 10, MPFR_RNDN); mpfr_init(r21195); mpfr_init_set_str(r21196, "0", 10, MPFR_RNDN); mpfr_init(r21197); mpfr_init(r21198); mpfr_init(r21199); mpfr_init(r21200); mpfr_init(r21201); mpfr_init_set_str(r21202, "4", 10, MPFR_RNDN); mpfr_init(r21203); mpfr_init(r21204); mpfr_init(r21205); mpfr_init(r21206); mpfr_init(r21207); mpfr_init(r21208); mpfr_init(r21209); mpfr_init(r21210); mpfr_init_set_str(r21211, "2", 10, MPFR_RNDN); mpfr_init(r21212); mpfr_init(r21213); mpfr_init(r21214); mpfr_init(r21215); mpfr_init(r21216); mpfr_init_set_str(r21217, "-4.086812811986243e-293", 10, MPFR_RNDN); mpfr_init(r21218); mpfr_init(r21219); mpfr_init(r21220); mpfr_init(r21221); mpfr_init(r21222); mpfr_init(r21223); mpfr_init(r21224); mpfr_init(r21225); mpfr_init(r21226); mpfr_init(r21227); mpfr_init(r21228); mpfr_init(r21229); mpfr_init_set_str(r21230, "1.8918781781269702e+105", 10, MPFR_RNDN); mpfr_init(r21231); mpfr_init(r21232); mpfr_init(r21233); mpfr_init(r21234); mpfr_init(r21235); mpfr_init(r21236); mpfr_init(r21237); mpfr_init(r21238); mpfr_init(r21239); mpfr_init(r21240); mpfr_init(r21241); mpfr_init(r21242); mpfr_init_set_str(r21243, "1", 10, MPFR_RNDN); mpfr_init(r21244); mpfr_init(r21245); mpfr_init(r21246); mpfr_init(r21247); mpfr_init(r21248); mpfr_init(r21249); mpfr_init(r21250); mpfr_init(r21251); mpfr_init(r21252); } double f_dm(double a, double b, double c) { mpfr_set_d(r21193, b, MPFR_RNDN); ; mpfr_set_si(r21195, mpfr_cmp(r21193, r21194) <= 0, MPFR_RNDN); ; mpfr_set_si(r21197, mpfr_cmp(r21193, r21196) >= 0, MPFR_RNDN); mpfr_neg(r21198, r21193, MPFR_RNDN); mpfr_mul(r21199, r21193, r21193, MPFR_RNDN); mpfr_set_d(r21200, c, MPFR_RNDN); mpfr_set_d(r21201, a, MPFR_RNDN); ; mpfr_mul(r21203, r21201, r21202, MPFR_RNDN); mpfr_mul(r21204, r21200, r21203, MPFR_RNDN); mpfr_sub(r21205, r21199, r21204, MPFR_RNDN); mpfr_sqrt(r21206, r21205, MPFR_RNDN); mpfr_sub(r21207, r21198, r21206, MPFR_RNDN); mpfr_add(r21208, r21201, r21201, MPFR_RNDN); mpfr_div(r21209, r21207, r21208, MPFR_RNDN); mpfr_add(r21210, r21200, r21200, MPFR_RNDN); ; mpfr_div(r21212, r21210, r21211, MPFR_RNDN); mpfr_div(r21213, r21200, r21193, MPFR_RNDN); mpfr_fma(r21214, r21201, r21213, r21198, MPFR_RNDN); mpfr_div(r21215, r21212, r21214, MPFR_RNDN); if (mpfr_get_si(r21197, MPFR_RNDN)) { mpfr_set(r21216, r21209, MPFR_RNDN); } else { mpfr_set(r21216, r21215, MPFR_RNDN); }; ; mpfr_set_si(r21218, mpfr_cmp(r21193, r21217) <= 0, MPFR_RNDN); mpfr_mul(r21219, r21202, r21201, MPFR_RNDN); mpfr_mul(r21220, r21219, r21200, MPFR_RNDN); mpfr_sub(r21221, r21199, r21220, MPFR_RNDN); mpfr_sqrt(r21222, r21221, MPFR_RNDN); mpfr_add(r21223, r21198, r21222, MPFR_RNDN); mpfr_div(r21224, r21204, r21223, MPFR_RNDN); mpfr_mul(r21225, r21211, r21201, MPFR_RNDN); mpfr_div(r21226, r21224, r21225, MPFR_RNDN); mpfr_mul(r21227, r21211, r21200, MPFR_RNDN); mpfr_div(r21228, r21227, r21223, MPFR_RNDN); if (mpfr_get_si(r21197, MPFR_RNDN)) { mpfr_set(r21229, r21226, MPFR_RNDN); } else { mpfr_set(r21229, r21228, MPFR_RNDN); }; ; mpfr_set_si(r21231, mpfr_cmp(r21193, r21230) <= 0, MPFR_RNDN); mpfr_sub(r21232, r21198, r21222, MPFR_RNDN); mpfr_div(r21233, r21232, r21225, MPFR_RNDN); mpfr_mul(r21234, r21201, r21200, MPFR_RNDN); mpfr_mul(r21235, r21234, r21202, MPFR_RNDN); mpfr_sub(r21236, r21199, r21235, MPFR_RNDN); mpfr_sqrt(r21237, r21236, MPFR_RNDN); mpfr_mul(r21238, r21237, r21236, MPFR_RNDN); mpfr_cbrt(r21239, r21238, MPFR_RNDN); mpfr_add(r21240, r21198, r21239, MPFR_RNDN); mpfr_div(r21241, r21227, r21240, MPFR_RNDN); if (mpfr_get_si(r21197, MPFR_RNDN)) { mpfr_set(r21242, r21233, MPFR_RNDN); } else { mpfr_set(r21242, r21241, MPFR_RNDN); }; ; mpfr_div(r21244, r21213, r21243, MPFR_RNDN); mpfr_div(r21245, r21193, r21201, MPFR_RNDN); mpfr_sub(r21246, r21244, r21245, MPFR_RNDN); mpfr_add(r21247, r21206, r21198, MPFR_RNDN); mpfr_div(r21248, r21210, r21247, MPFR_RNDN); if (mpfr_get_si(r21197, MPFR_RNDN)) { mpfr_set(r21249, r21246, MPFR_RNDN); } else { mpfr_set(r21249, r21248, MPFR_RNDN); }; if (mpfr_get_si(r21231, MPFR_RNDN)) { mpfr_set(r21250, r21242, MPFR_RNDN); } else { mpfr_set(r21250, r21249, MPFR_RNDN); }; if (mpfr_get_si(r21218, MPFR_RNDN)) { mpfr_set(r21251, r21229, MPFR_RNDN); } else { mpfr_set(r21251, r21250, MPFR_RNDN); }; if (mpfr_get_si(r21195, MPFR_RNDN)) { mpfr_set(r21252, r21216, MPFR_RNDN); } else { mpfr_set(r21252, r21251, MPFR_RNDN); }; return mpfr_get_d(r21252, MPFR_RNDN); }