#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 r25528 = N;
        float r25529 = 1;
        float r25530 = r25528 + r25529;
        float r25531 = atan(r25530);
        float r25532 = atan(r25528);
        float r25533 = r25531 - r25532;
        return r25533;
}

double f_id(double N) {
        double r25534 = N;
        double r25535 = 1;
        double r25536 = r25534 + r25535;
        double r25537 = atan(r25536);
        double r25538 = atan(r25534);
        double r25539 = r25537 - r25538;
        return r25539;
}


double f_of(float N) {
        float r25540 = 1;
        float r25541 = N;
        float r25542 = r25541 + r25540;
        float r25543 = r25542 * r25541;
        float r25544 = r25540 + r25543;
        float r25545 = atan2(r25540, r25544);
        return r25545;
}

double f_od(double N) {
        double r25546 = 1;
        double r25547 = N;
        double r25548 = r25547 + r25546;
        double r25549 = r25548 * r25547;
        double r25550 = r25546 + r25549;
        double r25551 = atan2(r25546, r25550);
        return r25551;
}

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 r25552, r25553, r25554, r25555, r25556, r25557;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init(r25552);
        mpfr_init_set_str(r25553, "1", 10, MPFR_RNDN);
        mpfr_init(r25554);
        mpfr_init(r25555);
        mpfr_init(r25556);
        mpfr_init(r25557);
}

double f_im(double N) {
        mpfr_set_d(r25552, N, MPFR_RNDN);
        ;
        mpfr_add(r25554, r25552, r25553, MPFR_RNDN);
        mpfr_atan(r25555, r25554, MPFR_RNDN);
        mpfr_atan(r25556, r25552, MPFR_RNDN);
        mpfr_sub(r25557, r25555, r25556, MPFR_RNDN);
        return mpfr_get_d(r25557, MPFR_RNDN);
}

static mpfr_t r25558, r25559, r25560, r25561, r25562, r25563;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r25558, "1", 10, MPFR_RNDN);
        mpfr_init(r25559);
        mpfr_init(r25560);
        mpfr_init(r25561);
        mpfr_init(r25562);
        mpfr_init(r25563);
}

double f_fm(double N) {
        ;
        mpfr_set_d(r25559, N, MPFR_RNDN);
        mpfr_add(r25560, r25559, r25558, MPFR_RNDN);
        mpfr_mul(r25561, r25560, r25559, MPFR_RNDN);
        mpfr_add(r25562, r25558, r25561, MPFR_RNDN);
        mpfr_atan2(r25563, r25558, r25562, MPFR_RNDN);
        return mpfr_get_d(r25563, MPFR_RNDN);
}

static mpfr_t r25564, r25565, r25566, r25567, r25568, r25569;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r25564, "1", 10, MPFR_RNDN);
        mpfr_init(r25565);
        mpfr_init(r25566);
        mpfr_init(r25567);
        mpfr_init(r25568);
        mpfr_init(r25569);
}

double f_dm(double N) {
        ;
        mpfr_set_d(r25565, N, MPFR_RNDN);
        mpfr_add(r25566, r25565, r25564, MPFR_RNDN);
        mpfr_mul(r25567, r25566, r25565, MPFR_RNDN);
        mpfr_add(r25568, r25564, r25567, MPFR_RNDN);
        mpfr_atan2(r25569, r25564, r25568, MPFR_RNDN);
        return mpfr_get_d(r25569, MPFR_RNDN);
}