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

char *name = "Quadratic roots, narrow range";

double f_if(float a, float b, float c) {
        float r21708 = b;
        float r21709 = -r21708;
        float r21710 = r21708 * r21708;
        float r21711 = 4;
        float r21712 = a;
        float r21713 = r21711 * r21712;
        float r21714 = c;
        float r21715 = r21713 * r21714;
        float r21716 = r21710 - r21715;
        float r21717 = sqrt(r21716);
        float r21718 = r21709 + r21717;
        float r21719 = 2;
        float r21720 = r21719 * r21712;
        float r21721 = r21718 / r21720;
        return r21721;
}

double f_id(double a, double b, double c) {
        double r21722 = b;
        double r21723 = -r21722;
        double r21724 = r21722 * r21722;
        double r21725 = 4;
        double r21726 = a;
        double r21727 = r21725 * r21726;
        double r21728 = c;
        double r21729 = r21727 * r21728;
        double r21730 = r21724 - r21729;
        double r21731 = sqrt(r21730);
        double r21732 = r21723 + r21731;
        double r21733 = 2;
        double r21734 = r21733 * r21726;
        double r21735 = r21732 / r21734;
        return r21735;
}


double f_of(float a, float b, float c) {
        float r21736 = 1;
        float r21737 = 2;
        float r21738 = a;
        float r21739 = r21737 * r21738;
        float r21740 = c;
        float r21741 = r21740 * r21738;
        float r21742 = 4;
        float r21743 = -r21742;
        float r21744 = r21741 * r21743;
        float r21745 = r21742 * r21738;
        float r21746 = -r21740;
        float r21747 = b;
        float r21748 = r21747 * r21747;
        float r21749 = fma(r21745, r21746, r21748);
        float r21750 = sqrt(r21749);
        float r21751 = r21750 + r21747;
        float r21752 = r21744 / r21751;
        float r21753 = r21739 / r21752;
        float r21754 = r21736 / r21753;
        return r21754;
}

double f_od(double a, double b, double c) {
        double r21755 = 1;
        double r21756 = 2;
        double r21757 = a;
        double r21758 = r21756 * r21757;
        double r21759 = c;
        double r21760 = r21759 * r21757;
        double r21761 = 4;
        double r21762 = -r21761;
        double r21763 = r21760 * r21762;
        double r21764 = r21761 * r21757;
        double r21765 = -r21759;
        double r21766 = b;
        double r21767 = r21766 * r21766;
        double r21768 = fma(r21764, r21765, r21767);
        double r21769 = sqrt(r21768);
        double r21770 = r21769 + r21766;
        double r21771 = r21763 / r21770;
        double r21772 = r21758 / r21771;
        double r21773 = r21755 / r21772;
        return r21773;
}

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 r21774, r21775, r21776, r21777, r21778, r21779, r21780, r21781, r21782, r21783, r21784, r21785, r21786, r21787;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(656);
        mpfr_init(r21774);
        mpfr_init(r21775);
        mpfr_init(r21776);
        mpfr_init_set_str(r21777, "4", 10, MPFR_RNDN);
        mpfr_init(r21778);
        mpfr_init(r21779);
        mpfr_init(r21780);
        mpfr_init(r21781);
        mpfr_init(r21782);
        mpfr_init(r21783);
        mpfr_init(r21784);
        mpfr_init_set_str(r21785, "2", 10, MPFR_RNDN);
        mpfr_init(r21786);
        mpfr_init(r21787);
}

double f_im(double a, double b, double c) {
        mpfr_set_d(r21774, b, MPFR_RNDN);
        mpfr_neg(r21775, r21774, MPFR_RNDN);
        mpfr_mul(r21776, r21774, r21774, MPFR_RNDN);
        ;
        mpfr_set_d(r21778, a, MPFR_RNDN);
        mpfr_mul(r21779, r21777, r21778, MPFR_RNDN);
        mpfr_set_d(r21780, c, MPFR_RNDN);
        mpfr_mul(r21781, r21779, r21780, MPFR_RNDN);
        mpfr_sub(r21782, r21776, r21781, MPFR_RNDN);
        mpfr_sqrt(r21783, r21782, MPFR_RNDN);
        mpfr_add(r21784, r21775, r21783, MPFR_RNDN);
        ;
        mpfr_mul(r21786, r21785, r21778, MPFR_RNDN);
        mpfr_div(r21787, r21784, r21786, MPFR_RNDN);
        return mpfr_get_d(r21787, MPFR_RNDN);
}

static mpfr_t r21788, r21789, r21790, r21791, r21792, r21793, r21794, r21795, r21796, r21797, r21798, r21799, r21800, r21801, r21802, r21803, r21804, r21805, r21806;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(656);
        mpfr_init_set_str(r21788, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r21789, "2", 10, MPFR_RNDN);
        mpfr_init(r21790);
        mpfr_init(r21791);
        mpfr_init(r21792);
        mpfr_init(r21793);
        mpfr_init_set_str(r21794, "4", 10, MPFR_RNDN);
        mpfr_init(r21795);
        mpfr_init(r21796);
        mpfr_init(r21797);
        mpfr_init(r21798);
        mpfr_init(r21799);
        mpfr_init(r21800);
        mpfr_init(r21801);
        mpfr_init(r21802);
        mpfr_init(r21803);
        mpfr_init(r21804);
        mpfr_init(r21805);
        mpfr_init(r21806);
}

double f_fm(double a, double b, double c) {
        ;
        ;
        mpfr_set_d(r21790, a, MPFR_RNDN);
        mpfr_mul(r21791, r21789, r21790, MPFR_RNDN);
        mpfr_set_d(r21792, c, MPFR_RNDN);
        mpfr_mul(r21793, r21792, r21790, MPFR_RNDN);
        ;
        mpfr_neg(r21795, r21794, MPFR_RNDN);
        mpfr_mul(r21796, r21793, r21795, MPFR_RNDN);
        mpfr_mul(r21797, r21794, r21790, MPFR_RNDN);
        mpfr_neg(r21798, r21792, MPFR_RNDN);
        mpfr_set_d(r21799, b, MPFR_RNDN);
        mpfr_mul(r21800, r21799, r21799, MPFR_RNDN);
        mpfr_fma(r21801, r21797, r21798, r21800, MPFR_RNDN);
        mpfr_sqrt(r21802, r21801, MPFR_RNDN);
        mpfr_add(r21803, r21802, r21799, MPFR_RNDN);
        mpfr_div(r21804, r21796, r21803, MPFR_RNDN);
        mpfr_div(r21805, r21791, r21804, MPFR_RNDN);
        mpfr_div(r21806, r21788, r21805, MPFR_RNDN);
        return mpfr_get_d(r21806, MPFR_RNDN);
}

static mpfr_t r21807, r21808, r21809, r21810, r21811, r21812, r21813, r21814, r21815, r21816, r21817, r21818, r21819, r21820, r21821, r21822, r21823, r21824, r21825;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(656);
        mpfr_init_set_str(r21807, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r21808, "2", 10, MPFR_RNDN);
        mpfr_init(r21809);
        mpfr_init(r21810);
        mpfr_init(r21811);
        mpfr_init(r21812);
        mpfr_init_set_str(r21813, "4", 10, MPFR_RNDN);
        mpfr_init(r21814);
        mpfr_init(r21815);
        mpfr_init(r21816);
        mpfr_init(r21817);
        mpfr_init(r21818);
        mpfr_init(r21819);
        mpfr_init(r21820);
        mpfr_init(r21821);
        mpfr_init(r21822);
        mpfr_init(r21823);
        mpfr_init(r21824);
        mpfr_init(r21825);
}

double f_dm(double a, double b, double c) {
        ;
        ;
        mpfr_set_d(r21809, a, MPFR_RNDN);
        mpfr_mul(r21810, r21808, r21809, MPFR_RNDN);
        mpfr_set_d(r21811, c, MPFR_RNDN);
        mpfr_mul(r21812, r21811, r21809, MPFR_RNDN);
        ;
        mpfr_neg(r21814, r21813, MPFR_RNDN);
        mpfr_mul(r21815, r21812, r21814, MPFR_RNDN);
        mpfr_mul(r21816, r21813, r21809, MPFR_RNDN);
        mpfr_neg(r21817, r21811, MPFR_RNDN);
        mpfr_set_d(r21818, b, MPFR_RNDN);
        mpfr_mul(r21819, r21818, r21818, MPFR_RNDN);
        mpfr_fma(r21820, r21816, r21817, r21819, MPFR_RNDN);
        mpfr_sqrt(r21821, r21820, MPFR_RNDN);
        mpfr_add(r21822, r21821, r21818, MPFR_RNDN);
        mpfr_div(r21823, r21815, r21822, MPFR_RNDN);
        mpfr_div(r21824, r21810, r21823, MPFR_RNDN);
        mpfr_div(r21825, r21807, r21824, MPFR_RNDN);
        return mpfr_get_d(r21825, MPFR_RNDN);
}