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

char *name = "(+ x (- (tan (+ y z)) (tan a)))";

double f_if(float x, float y, float z, float a) {
        float r27386 = x;
        float r27387 = y;
        float r27388 = z;
        float r27389 = r27387 + r27388;
        float r27390 = tan(r27389);
        float r27391 = a;
        float r27392 = tan(r27391);
        float r27393 = r27390 - r27392;
        float r27394 = r27386 + r27393;
        return r27394;
}

double f_id(double x, double y, double z, double a) {
        double r27395 = x;
        double r27396 = y;
        double r27397 = z;
        double r27398 = r27396 + r27397;
        double r27399 = tan(r27398);
        double r27400 = a;
        double r27401 = tan(r27400);
        double r27402 = r27399 - r27401;
        double r27403 = r27395 + r27402;
        return r27403;
}


double f_of(float x, float y, float z, float a) {
        float r27404 = x;
        float r27405 = y;
        float r27406 = tan(r27405);
        float r27407 = z;
        float r27408 = tan(r27407);
        float r27409 = r27406 + r27408;
        float r27410 = exp(r27409);
        float r27411 = log(r27410);
        float r27412 = 1;
        float r27413 = r27406 * r27408;
        float r27414 = r27412 - r27413;
        float r27415 = r27411 / r27414;
        float r27416 = a;
        float r27417 = tan(r27416);
        float r27418 = r27415 - r27417;
        float r27419 = r27404 + r27418;
        return r27419;
}

double f_od(double x, double y, double z, double a) {
        double r27420 = x;
        double r27421 = y;
        double r27422 = tan(r27421);
        double r27423 = z;
        double r27424 = tan(r27423);
        double r27425 = r27422 + r27424;
        double r27426 = exp(r27425);
        double r27427 = log(r27426);
        double r27428 = 1;
        double r27429 = r27422 * r27424;
        double r27430 = r27428 - r27429;
        double r27431 = r27427 / r27430;
        double r27432 = a;
        double r27433 = tan(r27432);
        double r27434 = r27431 - r27433;
        double r27435 = r27420 + r27434;
        return r27435;
}

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 r27436, r27437, r27438, r27439, r27440, r27441, r27442, r27443, r27444;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init(r27436);
        mpfr_init(r27437);
        mpfr_init(r27438);
        mpfr_init(r27439);
        mpfr_init(r27440);
        mpfr_init(r27441);
        mpfr_init(r27442);
        mpfr_init(r27443);
        mpfr_init(r27444);
}

double f_im(double x, double y, double z, double a) {
        mpfr_set_d(r27436, x, MPFR_RNDN);
        mpfr_set_d(r27437, y, MPFR_RNDN);
        mpfr_set_d(r27438, z, MPFR_RNDN);
        mpfr_add(r27439, r27437, r27438, MPFR_RNDN);
        mpfr_tan(r27440, r27439, MPFR_RNDN);
        mpfr_set_d(r27441, a, MPFR_RNDN);
        mpfr_tan(r27442, r27441, MPFR_RNDN);
        mpfr_sub(r27443, r27440, r27442, MPFR_RNDN);
        mpfr_add(r27444, r27436, r27443, MPFR_RNDN);
        return mpfr_get_d(r27444, MPFR_RNDN);
}

static mpfr_t r27445, r27446, r27447, r27448, r27449, r27450, r27451, r27452, r27453, r27454, r27455, r27456, r27457, r27458, r27459, r27460;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r27445);
        mpfr_init(r27446);
        mpfr_init(r27447);
        mpfr_init(r27448);
        mpfr_init(r27449);
        mpfr_init(r27450);
        mpfr_init(r27451);
        mpfr_init(r27452);
        mpfr_init_set_str(r27453, "1", 10, MPFR_RNDN);
        mpfr_init(r27454);
        mpfr_init(r27455);
        mpfr_init(r27456);
        mpfr_init(r27457);
        mpfr_init(r27458);
        mpfr_init(r27459);
        mpfr_init(r27460);
}

double f_fm(double x, double y, double z, double a) {
        mpfr_set_d(r27445, x, MPFR_RNDN);
        mpfr_set_d(r27446, y, MPFR_RNDN);
        mpfr_tan(r27447, r27446, MPFR_RNDN);
        mpfr_set_d(r27448, z, MPFR_RNDN);
        mpfr_tan(r27449, r27448, MPFR_RNDN);
        mpfr_add(r27450, r27447, r27449, MPFR_RNDN);
        mpfr_exp(r27451, r27450, MPFR_RNDN);
        mpfr_log(r27452, r27451, MPFR_RNDN);
        ;
        mpfr_mul(r27454, r27447, r27449, MPFR_RNDN);
        mpfr_sub(r27455, r27453, r27454, MPFR_RNDN);
        mpfr_div(r27456, r27452, r27455, MPFR_RNDN);
        mpfr_set_d(r27457, a, MPFR_RNDN);
        mpfr_tan(r27458, r27457, MPFR_RNDN);
        mpfr_sub(r27459, r27456, r27458, MPFR_RNDN);
        mpfr_add(r27460, r27445, r27459, MPFR_RNDN);
        return mpfr_get_d(r27460, MPFR_RNDN);
}

static mpfr_t r27461, r27462, r27463, r27464, r27465, r27466, r27467, r27468, r27469, r27470, r27471, r27472, r27473, r27474, r27475, r27476;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r27461);
        mpfr_init(r27462);
        mpfr_init(r27463);
        mpfr_init(r27464);
        mpfr_init(r27465);
        mpfr_init(r27466);
        mpfr_init(r27467);
        mpfr_init(r27468);
        mpfr_init_set_str(r27469, "1", 10, MPFR_RNDN);
        mpfr_init(r27470);
        mpfr_init(r27471);
        mpfr_init(r27472);
        mpfr_init(r27473);
        mpfr_init(r27474);
        mpfr_init(r27475);
        mpfr_init(r27476);
}

double f_dm(double x, double y, double z, double a) {
        mpfr_set_d(r27461, x, MPFR_RNDN);
        mpfr_set_d(r27462, y, MPFR_RNDN);
        mpfr_tan(r27463, r27462, MPFR_RNDN);
        mpfr_set_d(r27464, z, MPFR_RNDN);
        mpfr_tan(r27465, r27464, MPFR_RNDN);
        mpfr_add(r27466, r27463, r27465, MPFR_RNDN);
        mpfr_exp(r27467, r27466, MPFR_RNDN);
        mpfr_log(r27468, r27467, MPFR_RNDN);
        ;
        mpfr_mul(r27470, r27463, r27465, MPFR_RNDN);
        mpfr_sub(r27471, r27469, r27470, MPFR_RNDN);
        mpfr_div(r27472, r27468, r27471, MPFR_RNDN);
        mpfr_set_d(r27473, a, MPFR_RNDN);
        mpfr_tan(r27474, r27473, MPFR_RNDN);
        mpfr_sub(r27475, r27472, r27474, MPFR_RNDN);
        mpfr_add(r27476, r27461, r27475, MPFR_RNDN);
        return mpfr_get_d(r27476, MPFR_RNDN);
}