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

char *name = "sqrtexp (problem 3.4.4)";

double f_if(float x) {
        float r7445 = 2.0f;
        float r7446 = x;
        float r7447 = r7445 * r7446;
        float r7448 = exp(r7447);
        float r7449 = 1.0f;
        float r7450 = r7448 - r7449;
        float r7451 = exp(r7446);
        float r7452 = r7451 - r7449;
        float r7453 = r7450 / r7452;
        float r7454 = sqrt(r7453);
        return r7454;
}

double f_id(double x) {
        double r7455 = 2.0;
        double r7456 = x;
        double r7457 = r7455 * r7456;
        double r7458 = exp(r7457);
        double r7459 = 1.0;
        double r7460 = r7458 - r7459;
        double r7461 = exp(r7456);
        double r7462 = r7461 - r7459;
        double r7463 = r7460 / r7462;
        double r7464 = sqrt(r7463);
        return r7464;
}


double f_of(float x) {
        float r7465 = x;
        float r7466 = exp(r7465);
        float r7467 = 3.0f;
        float r7468 = pow(r7466, r7467);
        float r7469 = 1.0f;
        float r7470 = pow(r7469, r7467);
        float r7471 = r7468 + r7470;
        float r7472 = r7466 * r7466;
        float r7473 = r7469 - r7466;
        float r7474 = r7472 + r7473;
        float r7475 = r7471 / r7474;
        float r7476 = sqrt(r7475);
        return r7476;
}

double f_od(double x) {
        double r7477 = x;
        double r7478 = exp(r7477);
        double r7479 = 3.0;
        double r7480 = pow(r7478, r7479);
        double r7481 = 1.0;
        double r7482 = pow(r7481, r7479);
        double r7483 = r7480 + r7482;
        double r7484 = r7478 * r7478;
        double r7485 = r7481 - r7478;
        double r7486 = r7484 + r7485;
        double r7487 = r7483 / r7486;
        double r7488 = sqrt(r7487);
        return r7488;
}

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 r7489, r7490, r7491, r7492, r7493, r7494, r7495, r7496, r7497, r7498;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r7489, "2", 10, MPFR_RNDN);
        mpfr_init(r7490);
        mpfr_init(r7491);
        mpfr_init(r7492);
        mpfr_init_set_str(r7493, "1", 10, MPFR_RNDN);
        mpfr_init(r7494);
        mpfr_init(r7495);
        mpfr_init(r7496);
        mpfr_init(r7497);
        mpfr_init(r7498);
}

double f_im(double x) {
        ;
        mpfr_set_d(r7490, x, MPFR_RNDN);
        mpfr_mul(r7491, r7489, r7490, MPFR_RNDN);
        mpfr_exp(r7492, r7491, MPFR_RNDN);
        ;
        mpfr_sub(r7494, r7492, r7493, MPFR_RNDN);
        mpfr_exp(r7495, r7490, MPFR_RNDN);
        mpfr_sub(r7496, r7495, r7493, MPFR_RNDN);
        mpfr_div(r7497, r7494, r7496, MPFR_RNDN);
        mpfr_sqrt(r7498, r7497, MPFR_RNDN);
        return mpfr_get_d(r7498, MPFR_RNDN);
}

static mpfr_t r7499, r7500, r7501, r7502, r7503, r7504, r7505, r7506, r7507, r7508, r7509, r7510;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r7499);
        mpfr_init(r7500);
        mpfr_init_set_str(r7501, "3", 10, MPFR_RNDN);
        mpfr_init(r7502);
        mpfr_init_set_str(r7503, "1", 10, MPFR_RNDN);
        mpfr_init(r7504);
        mpfr_init(r7505);
        mpfr_init(r7506);
        mpfr_init(r7507);
        mpfr_init(r7508);
        mpfr_init(r7509);
        mpfr_init(r7510);
}

double f_fm(double x) {
        mpfr_set_d(r7499, x, MPFR_RNDN);
        mpfr_exp(r7500, r7499, MPFR_RNDN);
        ;
        mpfr_pow(r7502, r7500, r7501, MPFR_RNDN);
        ;
        mpfr_pow(r7504, r7503, r7501, MPFR_RNDN);
        mpfr_add(r7505, r7502, r7504, MPFR_RNDN);
        mpfr_mul(r7506, r7500, r7500, MPFR_RNDN);
        mpfr_sub(r7507, r7503, r7500, MPFR_RNDN);
        mpfr_add(r7508, r7506, r7507, MPFR_RNDN);
        mpfr_div(r7509, r7505, r7508, MPFR_RNDN);
        mpfr_sqrt(r7510, r7509, MPFR_RNDN);
        return mpfr_get_d(r7510, MPFR_RNDN);
}

static mpfr_t r7511, r7512, r7513, r7514, r7515, r7516, r7517, r7518, r7519, r7520, r7521, r7522;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r7511);
        mpfr_init(r7512);
        mpfr_init_set_str(r7513, "3", 10, MPFR_RNDN);
        mpfr_init(r7514);
        mpfr_init_set_str(r7515, "1", 10, MPFR_RNDN);
        mpfr_init(r7516);
        mpfr_init(r7517);
        mpfr_init(r7518);
        mpfr_init(r7519);
        mpfr_init(r7520);
        mpfr_init(r7521);
        mpfr_init(r7522);
}

double f_dm(double x) {
        mpfr_set_d(r7511, x, MPFR_RNDN);
        mpfr_exp(r7512, r7511, MPFR_RNDN);
        ;
        mpfr_pow(r7514, r7512, r7513, MPFR_RNDN);
        ;
        mpfr_pow(r7516, r7515, r7513, MPFR_RNDN);
        mpfr_add(r7517, r7514, r7516, MPFR_RNDN);
        mpfr_mul(r7518, r7512, r7512, MPFR_RNDN);
        mpfr_sub(r7519, r7515, r7512, MPFR_RNDN);
        mpfr_add(r7520, r7518, r7519, MPFR_RNDN);
        mpfr_div(r7521, r7517, r7520, MPFR_RNDN);
        mpfr_sqrt(r7522, r7521, MPFR_RNDN);
        return mpfr_get_d(r7522, MPFR_RNDN);
}