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

char *name = "2atan (example 3.5)";

double f_if(float N) {
        float r27437 = N;
        float r27438 = 1;
        float r27439 = r27437 + r27438;
        float r27440 = atan(r27439);
        float r27441 = atan(r27437);
        float r27442 = r27440 - r27441;
        return r27442;
}

double f_id(double N) {
        double r27443 = N;
        double r27444 = 1;
        double r27445 = r27443 + r27444;
        double r27446 = atan(r27445);
        double r27447 = atan(r27443);
        double r27448 = r27446 - r27447;
        return r27448;
}


double f_of(float N) {
        float r27449 = 1;
        float r27450 = N;
        float r27451 = r27450 + r27449;
        float r27452 = fma(r27451, r27450, r27449);
        float r27453 = atan2(r27449, r27452);
        return r27453;
}

double f_od(double N) {
        double r27454 = 1;
        double r27455 = N;
        double r27456 = r27455 + r27454;
        double r27457 = fma(r27456, r27455, r27454);
        double r27458 = atan2(r27454, r27457);
        return r27458;
}

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 r27459, r27460, r27461, r27462, r27463, r27464;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init(r27459);
        mpfr_init_set_str(r27460, "1", 10, MPFR_RNDN);
        mpfr_init(r27461);
        mpfr_init(r27462);
        mpfr_init(r27463);
        mpfr_init(r27464);
}

double f_im(double N) {
        mpfr_set_d(r27459, N, MPFR_RNDN);
        ;
        mpfr_add(r27461, r27459, r27460, MPFR_RNDN);
        mpfr_atan(r27462, r27461, MPFR_RNDN);
        mpfr_atan(r27463, r27459, MPFR_RNDN);
        mpfr_sub(r27464, r27462, r27463, MPFR_RNDN);
        return mpfr_get_d(r27464, MPFR_RNDN);
}

static mpfr_t r27465, r27466, r27467, r27468, r27469;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r27465, "1", 10, MPFR_RNDN);
        mpfr_init(r27466);
        mpfr_init(r27467);
        mpfr_init(r27468);
        mpfr_init(r27469);
}

double f_fm(double N) {
        ;
        mpfr_set_d(r27466, N, MPFR_RNDN);
        mpfr_add(r27467, r27466, r27465, MPFR_RNDN);
        mpfr_fma(r27468, r27467, r27466, r27465, MPFR_RNDN);
        mpfr_atan2(r27469, r27465, r27468, MPFR_RNDN);
        return mpfr_get_d(r27469, MPFR_RNDN);
}

static mpfr_t r27470, r27471, r27472, r27473, r27474;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r27470, "1", 10, MPFR_RNDN);
        mpfr_init(r27471);
        mpfr_init(r27472);
        mpfr_init(r27473);
        mpfr_init(r27474);
}

double f_dm(double N) {
        ;
        mpfr_set_d(r27471, N, MPFR_RNDN);
        mpfr_add(r27472, r27471, r27470, MPFR_RNDN);
        mpfr_fma(r27473, r27472, r27471, r27470, MPFR_RNDN);
        mpfr_atan2(r27474, r27470, r27473, MPFR_RNDN);
        return mpfr_get_d(r27474, MPFR_RNDN);
}