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

char *name = "NMSE problem 3.2.1";

double f_if(float a, float b_2F2, float c) {
        float r22619 = b_2F2;
        float r22620 = -r22619;
        float r22621 = r22619 * r22619;
        float r22622 = a;
        float r22623 = c;
        float r22624 = r22622 * r22623;
        float r22625 = r22621 - r22624;
        float r22626 = sqrt(r22625);
        float r22627 = r22620 - r22626;
        float r22628 = r22627 / r22622;
        return r22628;
}

double f_id(double a, double b_2F2, double c) {
        double r22629 = b_2F2;
        double r22630 = -r22629;
        double r22631 = r22629 * r22629;
        double r22632 = a;
        double r22633 = c;
        double r22634 = r22632 * r22633;
        double r22635 = r22631 - r22634;
        double r22636 = sqrt(r22635);
        double r22637 = r22630 - r22636;
        double r22638 = r22637 / r22632;
        return r22638;
}


double f_of(float a, float b_2F2, float c) {
        float r22639 = b_2F2;
        float r22640 = -7.200245550033268e-104;
        bool r22641 = r22639 <= r22640;
        float r22642 = c;
        float r22643 = 1/2;
        float r22644 = a;
        float r22645 = r22643 * r22644;
        float r22646 = r22639 / r22642;
        float r22647 = r22645 / r22646;
        float r22648 = r22639 + r22639;
        float r22649 = r22647 - r22648;
        float r22650 = r22642 / r22649;
        float r22651 = 3.409351165504793e+65;
        bool r22652 = r22639 <= r22651;
        float r22653 = 1;
        float r22654 = -r22639;
        float r22655 = r22639 * r22639;
        float r22656 = r22644 * r22642;
        float r22657 = r22655 - r22656;
        float r22658 = sqrt(r22657);
        float r22659 = r22654 - r22658;
        float r22660 = r22644 / r22659;
        float r22661 = r22653 / r22660;
        float r22662 = -2;
        float r22663 = r22639 / r22644;
        float r22664 = r22662 * r22663;
        float r22665 = r22652 ? r22661 : r22664;
        float r22666 = r22641 ? r22650 : r22665;
        return r22666;
}

double f_od(double a, double b_2F2, double c) {
        double r22667 = b_2F2;
        double r22668 = -7.200245550033268e-104;
        bool r22669 = r22667 <= r22668;
        double r22670 = c;
        double r22671 = 1/2;
        double r22672 = a;
        double r22673 = r22671 * r22672;
        double r22674 = r22667 / r22670;
        double r22675 = r22673 / r22674;
        double r22676 = r22667 + r22667;
        double r22677 = r22675 - r22676;
        double r22678 = r22670 / r22677;
        double r22679 = 3.409351165504793e+65;
        bool r22680 = r22667 <= r22679;
        double r22681 = 1;
        double r22682 = -r22667;
        double r22683 = r22667 * r22667;
        double r22684 = r22672 * r22670;
        double r22685 = r22683 - r22684;
        double r22686 = sqrt(r22685);
        double r22687 = r22682 - r22686;
        double r22688 = r22672 / r22687;
        double r22689 = r22681 / r22688;
        double r22690 = -2;
        double r22691 = r22667 / r22672;
        double r22692 = r22690 * r22691;
        double r22693 = r22680 ? r22689 : r22692;
        double r22694 = r22669 ? r22678 : r22693;
        return r22694;
}

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 r22695, r22696, r22697, r22698, r22699, r22700, r22701, r22702, r22703, r22704;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(3472);
        mpfr_init(r22695);
        mpfr_init(r22696);
        mpfr_init(r22697);
        mpfr_init(r22698);
        mpfr_init(r22699);
        mpfr_init(r22700);
        mpfr_init(r22701);
        mpfr_init(r22702);
        mpfr_init(r22703);
        mpfr_init(r22704);
}

double f_im(double a, double b_2F2, double c) {
        mpfr_set_d(r22695, b_2F2, MPFR_RNDN);
        mpfr_neg(r22696, r22695, MPFR_RNDN);
        mpfr_mul(r22697, r22695, r22695, MPFR_RNDN);
        mpfr_set_d(r22698, a, MPFR_RNDN);
        mpfr_set_d(r22699, c, MPFR_RNDN);
        mpfr_mul(r22700, r22698, r22699, MPFR_RNDN);
        mpfr_sub(r22701, r22697, r22700, MPFR_RNDN);
        mpfr_sqrt(r22702, r22701, MPFR_RNDN);
        mpfr_sub(r22703, r22696, r22702, MPFR_RNDN);
        mpfr_div(r22704, r22703, r22698, MPFR_RNDN);
        return mpfr_get_d(r22704, MPFR_RNDN);
}

static mpfr_t r22705, r22706, r22707, r22708, r22709, r22710, r22711, r22712, r22713, r22714, r22715, r22716, r22717, r22718, r22719, r22720, r22721, r22722, r22723, r22724, r22725, r22726, r22727, r22728, r22729, r22730, r22731, r22732;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(3472);
        mpfr_init(r22705);
        mpfr_init_set_str(r22706, "-7.200245550033268e-104", 10, MPFR_RNDN);
        mpfr_init(r22707);
        mpfr_init(r22708);
        mpfr_init_set_str(r22709, "1/2", 10, MPFR_RNDN);
        mpfr_init(r22710);
        mpfr_init(r22711);
        mpfr_init(r22712);
        mpfr_init(r22713);
        mpfr_init(r22714);
        mpfr_init(r22715);
        mpfr_init(r22716);
        mpfr_init_set_str(r22717, "3.409351165504793e+65", 10, MPFR_RNDN);
        mpfr_init(r22718);
        mpfr_init_set_str(r22719, "1", 10, MPFR_RNDN);
        mpfr_init(r22720);
        mpfr_init(r22721);
        mpfr_init(r22722);
        mpfr_init(r22723);
        mpfr_init(r22724);
        mpfr_init(r22725);
        mpfr_init(r22726);
        mpfr_init(r22727);
        mpfr_init_set_str(r22728, "-2", 10, MPFR_RNDN);
        mpfr_init(r22729);
        mpfr_init(r22730);
        mpfr_init(r22731);
        mpfr_init(r22732);
}

double f_fm(double a, double b_2F2, double c) {
        mpfr_set_d(r22705, b_2F2, MPFR_RNDN);
        ;
        mpfr_set_si(r22707, mpfr_cmp(r22705, r22706) <= 0, MPFR_RNDN);
        mpfr_set_d(r22708, c, MPFR_RNDN);
        ;
        mpfr_set_d(r22710, a, MPFR_RNDN);
        mpfr_mul(r22711, r22709, r22710, MPFR_RNDN);
        mpfr_div(r22712, r22705, r22708, MPFR_RNDN);
        mpfr_div(r22713, r22711, r22712, MPFR_RNDN);
        mpfr_add(r22714, r22705, r22705, MPFR_RNDN);
        mpfr_sub(r22715, r22713, r22714, MPFR_RNDN);
        mpfr_div(r22716, r22708, r22715, MPFR_RNDN);
        ;
        mpfr_set_si(r22718, mpfr_cmp(r22705, r22717) <= 0, MPFR_RNDN);
        ;
        mpfr_neg(r22720, r22705, MPFR_RNDN);
        mpfr_mul(r22721, r22705, r22705, MPFR_RNDN);
        mpfr_mul(r22722, r22710, r22708, MPFR_RNDN);
        mpfr_sub(r22723, r22721, r22722, MPFR_RNDN);
        mpfr_sqrt(r22724, r22723, MPFR_RNDN);
        mpfr_sub(r22725, r22720, r22724, MPFR_RNDN);
        mpfr_div(r22726, r22710, r22725, MPFR_RNDN);
        mpfr_div(r22727, r22719, r22726, MPFR_RNDN);
        ;
        mpfr_div(r22729, r22705, r22710, MPFR_RNDN);
        mpfr_mul(r22730, r22728, r22729, MPFR_RNDN);
        if (mpfr_get_si(r22718, MPFR_RNDN)) { mpfr_set(r22731, r22727, MPFR_RNDN); } else { mpfr_set(r22731, r22730, MPFR_RNDN); };
        if (mpfr_get_si(r22707, MPFR_RNDN)) { mpfr_set(r22732, r22716, MPFR_RNDN); } else { mpfr_set(r22732, r22731, MPFR_RNDN); };
        return mpfr_get_d(r22732, MPFR_RNDN);
}

static mpfr_t r22733, r22734, r22735, r22736, r22737, r22738, r22739, r22740, r22741, r22742, r22743, r22744, r22745, r22746, r22747, r22748, r22749, r22750, r22751, r22752, r22753, r22754, r22755, r22756, r22757, r22758, r22759, r22760;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(3472);
        mpfr_init(r22733);
        mpfr_init_set_str(r22734, "-7.200245550033268e-104", 10, MPFR_RNDN);
        mpfr_init(r22735);
        mpfr_init(r22736);
        mpfr_init_set_str(r22737, "1/2", 10, MPFR_RNDN);
        mpfr_init(r22738);
        mpfr_init(r22739);
        mpfr_init(r22740);
        mpfr_init(r22741);
        mpfr_init(r22742);
        mpfr_init(r22743);
        mpfr_init(r22744);
        mpfr_init_set_str(r22745, "3.409351165504793e+65", 10, MPFR_RNDN);
        mpfr_init(r22746);
        mpfr_init_set_str(r22747, "1", 10, MPFR_RNDN);
        mpfr_init(r22748);
        mpfr_init(r22749);
        mpfr_init(r22750);
        mpfr_init(r22751);
        mpfr_init(r22752);
        mpfr_init(r22753);
        mpfr_init(r22754);
        mpfr_init(r22755);
        mpfr_init_set_str(r22756, "-2", 10, MPFR_RNDN);
        mpfr_init(r22757);
        mpfr_init(r22758);
        mpfr_init(r22759);
        mpfr_init(r22760);
}

double f_dm(double a, double b_2F2, double c) {
        mpfr_set_d(r22733, b_2F2, MPFR_RNDN);
        ;
        mpfr_set_si(r22735, mpfr_cmp(r22733, r22734) <= 0, MPFR_RNDN);
        mpfr_set_d(r22736, c, MPFR_RNDN);
        ;
        mpfr_set_d(r22738, a, MPFR_RNDN);
        mpfr_mul(r22739, r22737, r22738, MPFR_RNDN);
        mpfr_div(r22740, r22733, r22736, MPFR_RNDN);
        mpfr_div(r22741, r22739, r22740, MPFR_RNDN);
        mpfr_add(r22742, r22733, r22733, MPFR_RNDN);
        mpfr_sub(r22743, r22741, r22742, MPFR_RNDN);
        mpfr_div(r22744, r22736, r22743, MPFR_RNDN);
        ;
        mpfr_set_si(r22746, mpfr_cmp(r22733, r22745) <= 0, MPFR_RNDN);
        ;
        mpfr_neg(r22748, r22733, MPFR_RNDN);
        mpfr_mul(r22749, r22733, r22733, MPFR_RNDN);
        mpfr_mul(r22750, r22738, r22736, MPFR_RNDN);
        mpfr_sub(r22751, r22749, r22750, MPFR_RNDN);
        mpfr_sqrt(r22752, r22751, MPFR_RNDN);
        mpfr_sub(r22753, r22748, r22752, MPFR_RNDN);
        mpfr_div(r22754, r22738, r22753, MPFR_RNDN);
        mpfr_div(r22755, r22747, r22754, MPFR_RNDN);
        ;
        mpfr_div(r22757, r22733, r22738, MPFR_RNDN);
        mpfr_mul(r22758, r22756, r22757, MPFR_RNDN);
        if (mpfr_get_si(r22746, MPFR_RNDN)) { mpfr_set(r22759, r22755, MPFR_RNDN); } else { mpfr_set(r22759, r22758, MPFR_RNDN); };
        if (mpfr_get_si(r22735, MPFR_RNDN)) { mpfr_set(r22760, r22744, MPFR_RNDN); } else { mpfr_set(r22760, r22759, MPFR_RNDN); };
        return mpfr_get_d(r22760, MPFR_RNDN);
}