#include <tgmath.h>
#include <gmp.h>
#include <mpfr.h>
#include <stdio.h>
#include <stdbool.h>

char *name = "The quadratic formula (r2)";

double f_if(float a, float b, float c) {
        float r8793 = b;
        float r8794 = -r8793;
        float r8795 = r8793 * r8793;
        float r8796 = 4.0f;
        float r8797 = a;
        float r8798 = c;
        float r8799 = r8797 * r8798;
        float r8800 = r8796 * r8799;
        float r8801 = r8795 - r8800;
        float r8802 = sqrt(r8801);
        float r8803 = r8794 - r8802;
        float r8804 = 2.0f;
        float r8805 = r8804 * r8797;
        float r8806 = r8803 / r8805;
        return r8806;
}

double f_id(double a, double b, double c) {
        double r8807 = b;
        double r8808 = -r8807;
        double r8809 = r8807 * r8807;
        double r8810 = 4.0;
        double r8811 = a;
        double r8812 = c;
        double r8813 = r8811 * r8812;
        double r8814 = r8810 * r8813;
        double r8815 = r8809 - r8814;
        double r8816 = sqrt(r8815);
        double r8817 = r8808 - r8816;
        double r8818 = 2.0;
        double r8819 = r8818 * r8811;
        double r8820 = r8817 / r8819;
        return r8820;
}


double f_of(float a, float b, float c) {
        float r8821 = b;
        float r8822 = -r8821;
        float r8823 = -9.670079395517677e+115f;
        bool r8824 = r8822 <= r8823;
        float r8825 = a;
        float r8826 = r8821 / r8825;
        float r8827 = -r8826;
        float r8828 = 1.094348108854397e-160f;
        bool r8829 = r8822 <= r8828;
        float r8830 = 4.0f;
        float r8831 = r8825 * r8830;
        float r8832 = c;
        float r8833 = -r8832;
        float r8834 = r8821 * r8821;
        float r8835 = fma(r8831, r8833, r8834);
        float r8836 = sqrt(r8835);
        float r8837 = r8822 - r8836;
        float r8838 = 2.0f;
        float r8839 = r8825 * r8838;
        float r8840 = r8837 / r8839;
        float r8841 = 1.2760553527866898e+132f;
        bool r8842 = r8822 <= r8841;
        float r8843 = r8830 * r8832;
        float r8844 = r8843 / r8838;
        float r8845 = cbrt(r8844);
        float r8846 = -r8830;
        float r8847 = r8825 * r8832;
        float r8848 = fma(r8846, r8847, r8834);
        float r8849 = sqrt(r8848);
        float r8850 = r8849 - r8821;
        float r8851 = r8845 / r8850;
        float r8852 = r8838 / r8832;
        float r8853 = r8830 / r8852;
        float r8854 = cbrt(r8853);
        float r8855 = r8854 * r8854;
        float r8856 = r8851 * r8855;
        float r8857 = r8847 / r8821;
        float r8858 = -2.0f;
        float r8859 = r8857 * r8858;
        float r8860 = r8859 / r8839;
        float r8861 = r8842 ? r8856 : r8860;
        float r8862 = r8829 ? r8840 : r8861;
        float r8863 = r8824 ? r8827 : r8862;
        return r8863;
}

double f_od(double a, double b, double c) {
        double r8864 = b;
        double r8865 = -r8864;
        double r8866 = -9.670079395517677e+115;
        bool r8867 = r8865 <= r8866;
        double r8868 = a;
        double r8869 = r8864 / r8868;
        double r8870 = -r8869;
        double r8871 = 1.094348108854397e-160;
        bool r8872 = r8865 <= r8871;
        double r8873 = 4.0;
        double r8874 = r8868 * r8873;
        double r8875 = c;
        double r8876 = -r8875;
        double r8877 = r8864 * r8864;
        double r8878 = fma(r8874, r8876, r8877);
        double r8879 = sqrt(r8878);
        double r8880 = r8865 - r8879;
        double r8881 = 2.0;
        double r8882 = r8868 * r8881;
        double r8883 = r8880 / r8882;
        double r8884 = 1.2760553527866898e+132;
        bool r8885 = r8865 <= r8884;
        double r8886 = r8873 * r8875;
        double r8887 = r8886 / r8881;
        double r8888 = cbrt(r8887);
        double r8889 = -r8873;
        double r8890 = r8868 * r8875;
        double r8891 = fma(r8889, r8890, r8877);
        double r8892 = sqrt(r8891);
        double r8893 = r8892 - r8864;
        double r8894 = r8888 / r8893;
        double r8895 = r8881 / r8875;
        double r8896 = r8873 / r8895;
        double r8897 = cbrt(r8896);
        double r8898 = r8897 * r8897;
        double r8899 = r8894 * r8898;
        double r8900 = r8890 / r8864;
        double r8901 = -2.0;
        double r8902 = r8900 * r8901;
        double r8903 = r8902 / r8882;
        double r8904 = r8885 ? r8899 : r8903;
        double r8905 = r8872 ? r8883 : r8904;
        double r8906 = r8867 ? r8870 : r8905;
        return r8906;
}

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 r8907, r8908, r8909, r8910, r8911, r8912, r8913, r8914, r8915, r8916, r8917, r8918, r8919, r8920;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(3408);
        mpfr_init(r8907);
        mpfr_init(r8908);
        mpfr_init(r8909);
        mpfr_init_set_str(r8910, "4", 10, MPFR_RNDN);
        mpfr_init(r8911);
        mpfr_init(r8912);
        mpfr_init(r8913);
        mpfr_init(r8914);
        mpfr_init(r8915);
        mpfr_init(r8916);
        mpfr_init(r8917);
        mpfr_init_set_str(r8918, "2", 10, MPFR_RNDN);
        mpfr_init(r8919);
        mpfr_init(r8920);
}

double f_im(double a, double b, double c) {
        mpfr_set_d(r8907, b, MPFR_RNDN);
        mpfr_neg(r8908, r8907, MPFR_RNDN);
        mpfr_mul(r8909, r8907, r8907, MPFR_RNDN);
        ;
        mpfr_set_d(r8911, a, MPFR_RNDN);
        mpfr_set_d(r8912, c, MPFR_RNDN);
        mpfr_mul(r8913, r8911, r8912, MPFR_RNDN);
        mpfr_mul(r8914, r8910, r8913, MPFR_RNDN);
        mpfr_sub(r8915, r8909, r8914, MPFR_RNDN);
        mpfr_sqrt(r8916, r8915, MPFR_RNDN);
        mpfr_sub(r8917, r8908, r8916, MPFR_RNDN);
        ;
        mpfr_mul(r8919, r8918, r8911, MPFR_RNDN);
        mpfr_div(r8920, r8917, r8919, MPFR_RNDN);
        return mpfr_get_d(r8920, MPFR_RNDN);
}

static mpfr_t r8921, r8922, r8923, r8924, r8925, r8926, r8927, r8928, r8929, r8930, r8931, r8932, r8933, r8934, r8935, r8936, r8937, r8938, r8939, r8940, r8941, r8942, r8943, r8944, r8945, r8946, r8947, r8948, r8949, r8950, r8951, r8952, r8953, r8954, r8955, r8956, r8957, r8958, r8959, r8960, r8961, r8962, r8963;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(3408);
        mpfr_init(r8921);
        mpfr_init(r8922);
        mpfr_init_set_str(r8923, "-9.670079395517677e+115", 10, MPFR_RNDN);
        mpfr_init(r8924);
        mpfr_init(r8925);
        mpfr_init(r8926);
        mpfr_init(r8927);
        mpfr_init_set_str(r8928, "1.094348108854397e-160", 10, MPFR_RNDN);
        mpfr_init(r8929);
        mpfr_init_set_str(r8930, "4", 10, MPFR_RNDN);
        mpfr_init(r8931);
        mpfr_init(r8932);
        mpfr_init(r8933);
        mpfr_init(r8934);
        mpfr_init(r8935);
        mpfr_init(r8936);
        mpfr_init(r8937);
        mpfr_init_set_str(r8938, "2", 10, MPFR_RNDN);
        mpfr_init(r8939);
        mpfr_init(r8940);
        mpfr_init_set_str(r8941, "1.2760553527866898e+132", 10, MPFR_RNDN);
        mpfr_init(r8942);
        mpfr_init(r8943);
        mpfr_init(r8944);
        mpfr_init(r8945);
        mpfr_init(r8946);
        mpfr_init(r8947);
        mpfr_init(r8948);
        mpfr_init(r8949);
        mpfr_init(r8950);
        mpfr_init(r8951);
        mpfr_init(r8952);
        mpfr_init(r8953);
        mpfr_init(r8954);
        mpfr_init(r8955);
        mpfr_init(r8956);
        mpfr_init(r8957);
        mpfr_init_set_str(r8958, "-2", 10, MPFR_RNDN);
        mpfr_init(r8959);
        mpfr_init(r8960);
        mpfr_init(r8961);
        mpfr_init(r8962);
        mpfr_init(r8963);
}

double f_fm(double a, double b, double c) {
        mpfr_set_d(r8921, b, MPFR_RNDN);
        mpfr_neg(r8922, r8921, MPFR_RNDN);
        ;
        mpfr_set_si(r8924, mpfr_cmp(r8922, r8923) <= 0, MPFR_RNDN);
        mpfr_set_d(r8925, a, MPFR_RNDN);
        mpfr_div(r8926, r8921, r8925, MPFR_RNDN);
        mpfr_neg(r8927, r8926, MPFR_RNDN);
        ;
        mpfr_set_si(r8929, mpfr_cmp(r8922, r8928) <= 0, MPFR_RNDN);
        ;
        mpfr_mul(r8931, r8925, r8930, MPFR_RNDN);
        mpfr_set_d(r8932, c, MPFR_RNDN);
        mpfr_neg(r8933, r8932, MPFR_RNDN);
        mpfr_mul(r8934, r8921, r8921, MPFR_RNDN);
        mpfr_fma(r8935, r8931, r8933, r8934, MPFR_RNDN);
        mpfr_sqrt(r8936, r8935, MPFR_RNDN);
        mpfr_sub(r8937, r8922, r8936, MPFR_RNDN);
        ;
        mpfr_mul(r8939, r8925, r8938, MPFR_RNDN);
        mpfr_div(r8940, r8937, r8939, MPFR_RNDN);
        ;
        mpfr_set_si(r8942, mpfr_cmp(r8922, r8941) <= 0, MPFR_RNDN);
        mpfr_mul(r8943, r8930, r8932, MPFR_RNDN);
        mpfr_div(r8944, r8943, r8938, MPFR_RNDN);
        mpfr_cbrt(r8945, r8944, MPFR_RNDN);
        mpfr_neg(r8946, r8930, MPFR_RNDN);
        mpfr_mul(r8947, r8925, r8932, MPFR_RNDN);
        mpfr_fma(r8948, r8946, r8947, r8934, MPFR_RNDN);
        mpfr_sqrt(r8949, r8948, MPFR_RNDN);
        mpfr_sub(r8950, r8949, r8921, MPFR_RNDN);
        mpfr_div(r8951, r8945, r8950, MPFR_RNDN);
        mpfr_div(r8952, r8938, r8932, MPFR_RNDN);
        mpfr_div(r8953, r8930, r8952, MPFR_RNDN);
        mpfr_cbrt(r8954, r8953, MPFR_RNDN);
        mpfr_mul(r8955, r8954, r8954, MPFR_RNDN);
        mpfr_mul(r8956, r8951, r8955, MPFR_RNDN);
        mpfr_div(r8957, r8947, r8921, MPFR_RNDN);
        ;
        mpfr_mul(r8959, r8957, r8958, MPFR_RNDN);
        mpfr_div(r8960, r8959, r8939, MPFR_RNDN);
        if (mpfr_get_si(r8942, MPFR_RNDN)) { mpfr_set(r8961, r8956, MPFR_RNDN); } else { mpfr_set(r8961, r8960, MPFR_RNDN); };
        if (mpfr_get_si(r8929, MPFR_RNDN)) { mpfr_set(r8962, r8940, MPFR_RNDN); } else { mpfr_set(r8962, r8961, MPFR_RNDN); };
        if (mpfr_get_si(r8924, MPFR_RNDN)) { mpfr_set(r8963, r8927, MPFR_RNDN); } else { mpfr_set(r8963, r8962, MPFR_RNDN); };
        return mpfr_get_d(r8963, MPFR_RNDN);
}

static mpfr_t r8964, r8965, r8966, r8967, r8968, r8969, r8970, r8971, r8972, r8973, r8974, r8975, r8976, r8977, r8978, r8979, r8980, r8981, r8982, r8983, r8984, r8985, r8986, r8987, r8988, r8989, r8990, r8991, r8992, r8993, r8994, r8995, r8996, r8997, r8998, r8999, r9000, r9001, r9002, r9003, r9004, r9005, r9006;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(3408);
        mpfr_init(r8964);
        mpfr_init(r8965);
        mpfr_init_set_str(r8966, "-9.670079395517677e+115", 10, MPFR_RNDN);
        mpfr_init(r8967);
        mpfr_init(r8968);
        mpfr_init(r8969);
        mpfr_init(r8970);
        mpfr_init_set_str(r8971, "1.094348108854397e-160", 10, MPFR_RNDN);
        mpfr_init(r8972);
        mpfr_init_set_str(r8973, "4", 10, MPFR_RNDN);
        mpfr_init(r8974);
        mpfr_init(r8975);
        mpfr_init(r8976);
        mpfr_init(r8977);
        mpfr_init(r8978);
        mpfr_init(r8979);
        mpfr_init(r8980);
        mpfr_init_set_str(r8981, "2", 10, MPFR_RNDN);
        mpfr_init(r8982);
        mpfr_init(r8983);
        mpfr_init_set_str(r8984, "1.2760553527866898e+132", 10, MPFR_RNDN);
        mpfr_init(r8985);
        mpfr_init(r8986);
        mpfr_init(r8987);
        mpfr_init(r8988);
        mpfr_init(r8989);
        mpfr_init(r8990);
        mpfr_init(r8991);
        mpfr_init(r8992);
        mpfr_init(r8993);
        mpfr_init(r8994);
        mpfr_init(r8995);
        mpfr_init(r8996);
        mpfr_init(r8997);
        mpfr_init(r8998);
        mpfr_init(r8999);
        mpfr_init(r9000);
        mpfr_init_set_str(r9001, "-2", 10, MPFR_RNDN);
        mpfr_init(r9002);
        mpfr_init(r9003);
        mpfr_init(r9004);
        mpfr_init(r9005);
        mpfr_init(r9006);
}

double f_dm(double a, double b, double c) {
        mpfr_set_d(r8964, b, MPFR_RNDN);
        mpfr_neg(r8965, r8964, MPFR_RNDN);
        ;
        mpfr_set_si(r8967, mpfr_cmp(r8965, r8966) <= 0, MPFR_RNDN);
        mpfr_set_d(r8968, a, MPFR_RNDN);
        mpfr_div(r8969, r8964, r8968, MPFR_RNDN);
        mpfr_neg(r8970, r8969, MPFR_RNDN);
        ;
        mpfr_set_si(r8972, mpfr_cmp(r8965, r8971) <= 0, MPFR_RNDN);
        ;
        mpfr_mul(r8974, r8968, r8973, MPFR_RNDN);
        mpfr_set_d(r8975, c, MPFR_RNDN);
        mpfr_neg(r8976, r8975, MPFR_RNDN);
        mpfr_mul(r8977, r8964, r8964, MPFR_RNDN);
        mpfr_fma(r8978, r8974, r8976, r8977, MPFR_RNDN);
        mpfr_sqrt(r8979, r8978, MPFR_RNDN);
        mpfr_sub(r8980, r8965, r8979, MPFR_RNDN);
        ;
        mpfr_mul(r8982, r8968, r8981, MPFR_RNDN);
        mpfr_div(r8983, r8980, r8982, MPFR_RNDN);
        ;
        mpfr_set_si(r8985, mpfr_cmp(r8965, r8984) <= 0, MPFR_RNDN);
        mpfr_mul(r8986, r8973, r8975, MPFR_RNDN);
        mpfr_div(r8987, r8986, r8981, MPFR_RNDN);
        mpfr_cbrt(r8988, r8987, MPFR_RNDN);
        mpfr_neg(r8989, r8973, MPFR_RNDN);
        mpfr_mul(r8990, r8968, r8975, MPFR_RNDN);
        mpfr_fma(r8991, r8989, r8990, r8977, MPFR_RNDN);
        mpfr_sqrt(r8992, r8991, MPFR_RNDN);
        mpfr_sub(r8993, r8992, r8964, MPFR_RNDN);
        mpfr_div(r8994, r8988, r8993, MPFR_RNDN);
        mpfr_div(r8995, r8981, r8975, MPFR_RNDN);
        mpfr_div(r8996, r8973, r8995, MPFR_RNDN);
        mpfr_cbrt(r8997, r8996, MPFR_RNDN);
        mpfr_mul(r8998, r8997, r8997, MPFR_RNDN);
        mpfr_mul(r8999, r8994, r8998, MPFR_RNDN);
        mpfr_div(r9000, r8990, r8964, MPFR_RNDN);
        ;
        mpfr_mul(r9002, r9000, r9001, MPFR_RNDN);
        mpfr_div(r9003, r9002, r8982, MPFR_RNDN);
        if (mpfr_get_si(r8985, MPFR_RNDN)) { mpfr_set(r9004, r8999, MPFR_RNDN); } else { mpfr_set(r9004, r9003, MPFR_RNDN); };
        if (mpfr_get_si(r8972, MPFR_RNDN)) { mpfr_set(r9005, r8983, MPFR_RNDN); } else { mpfr_set(r9005, r9004, MPFR_RNDN); };
        if (mpfr_get_si(r8967, MPFR_RNDN)) { mpfr_set(r9006, r8970, MPFR_RNDN); } else { mpfr_set(r9006, r9005, MPFR_RNDN); };
        return mpfr_get_d(r9006, MPFR_RNDN);
}