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

char *name = "arccos";

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

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


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

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

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

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

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

static mpfr_t r7500, r7501, r7502, r7503, r7504, r7505, r7506, r7507, r7508, r7509, r7510, r7511, r7512, r7513;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init(r7500);
        mpfr_init_set_str(r7501, "1", 10, MPFR_RNDN);
        mpfr_init(r7502);
        mpfr_init(r7503);
        mpfr_init_set_str(r7504, "3", 10, MPFR_RNDN);
        mpfr_init(r7505);
        mpfr_init(r7506);
        mpfr_init(r7507);
        mpfr_init(r7508);
        mpfr_init(r7509);
        mpfr_init(r7510);
        mpfr_init(r7511);
        mpfr_init_set_str(r7512, "2", 10, MPFR_RNDN);
        mpfr_init(r7513);
}

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

static mpfr_t r7514, r7515, r7516, r7517, r7518, r7519, r7520, r7521, r7522, r7523, r7524, r7525, r7526, r7527;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init(r7514);
        mpfr_init_set_str(r7515, "1", 10, MPFR_RNDN);
        mpfr_init(r7516);
        mpfr_init(r7517);
        mpfr_init_set_str(r7518, "3", 10, MPFR_RNDN);
        mpfr_init(r7519);
        mpfr_init(r7520);
        mpfr_init(r7521);
        mpfr_init(r7522);
        mpfr_init(r7523);
        mpfr_init(r7524);
        mpfr_init(r7525);
        mpfr_init_set_str(r7526, "2", 10, MPFR_RNDN);
        mpfr_init(r7527);
}

double f_dm(double x) {
        mpfr_set_d(r7514, x, MPFR_RNDN);
        ;
        mpfr_sub(r7516, r7515, r7514, MPFR_RNDN);
        mpfr_fma(r7517, r7514, r7514, r7516, MPFR_RNDN);
        ;
        mpfr_pow(r7519, r7515, r7518, MPFR_RNDN);
        mpfr_pow(r7520, r7514, r7518, MPFR_RNDN);
        mpfr_add(r7521, r7519, r7520, MPFR_RNDN);
        mpfr_div(r7522, r7516, r7521, MPFR_RNDN);
        mpfr_mul(r7523, r7517, r7522, MPFR_RNDN);
        mpfr_sqrt(r7524, r7523, MPFR_RNDN);
        mpfr_atan(r7525, r7524, MPFR_RNDN);
        ;
        mpfr_mul(r7527, r7525, r7526, MPFR_RNDN);
        return mpfr_get_d(r7527, MPFR_RNDN);
}

