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

char *name = "Diagrams.ThreeD.Transform:aboutY from diagrams-lib-1.3.0.3";

double f_if(float x, float y, float z) {
        float r34343 = x;
        float r34344 = y;
        float r34345 = cos(r34344);
        float r34346 = r34343 * r34345;
        float r34347 = z;
        float r34348 = sin(r34344);
        float r34349 = r34347 * r34348;
        float r34350 = r34346 + r34349;
        return r34350;
}

double f_id(double x, double y, double z) {
        double r34351 = x;
        double r34352 = y;
        double r34353 = cos(r34352);
        double r34354 = r34351 * r34353;
        double r34355 = z;
        double r34356 = sin(r34352);
        double r34357 = r34355 * r34356;
        double r34358 = r34354 + r34357;
        return r34358;
}


double f_of(float x, float y, float z) {
        float r34359 = x;
        float r34360 = y;
        float r34361 = cos(r34360);
        float r34362 = r34361 * r34361;
        float r34363 = 1/3;
        float r34364 = pow(r34362, r34363);
        float r34365 = r34359 * r34364;
        float r34366 = cbrt(r34361);
        float r34367 = r34365 * r34366;
        float r34368 = z;
        float r34369 = sin(r34360);
        float r34370 = r34368 * r34369;
        float r34371 = r34367 + r34370;
        return r34371;
}

double f_od(double x, double y, double z) {
        double r34372 = x;
        double r34373 = y;
        double r34374 = cos(r34373);
        double r34375 = r34374 * r34374;
        double r34376 = 1/3;
        double r34377 = pow(r34375, r34376);
        double r34378 = r34372 * r34377;
        double r34379 = cbrt(r34374);
        double r34380 = r34378 * r34379;
        double r34381 = z;
        double r34382 = sin(r34373);
        double r34383 = r34381 * r34382;
        double r34384 = r34380 + r34383;
        return r34384;
}

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 r34385, r34386, r34387, r34388, r34389, r34390, r34391, r34392;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r34385);
        mpfr_init(r34386);
        mpfr_init(r34387);
        mpfr_init(r34388);
        mpfr_init(r34389);
        mpfr_init(r34390);
        mpfr_init(r34391);
        mpfr_init(r34392);
}

double f_im(double x, double y, double z) {
        mpfr_set_d(r34385, x, MPFR_RNDN);
        mpfr_set_d(r34386, y, MPFR_RNDN);
        mpfr_cos(r34387, r34386, MPFR_RNDN);
        mpfr_mul(r34388, r34385, r34387, MPFR_RNDN);
        mpfr_set_d(r34389, z, MPFR_RNDN);
        mpfr_sin(r34390, r34386, MPFR_RNDN);
        mpfr_mul(r34391, r34389, r34390, MPFR_RNDN);
        mpfr_add(r34392, r34388, r34391, MPFR_RNDN);
        return mpfr_get_d(r34392, MPFR_RNDN);
}

static mpfr_t r34393, r34394, r34395, r34396, r34397, r34398, r34399, r34400, r34401, r34402, r34403, r34404, r34405;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r34393);
        mpfr_init(r34394);
        mpfr_init(r34395);
        mpfr_init(r34396);
        mpfr_init_set_str(r34397, "1/3", 10, MPFR_RNDN);
        mpfr_init(r34398);
        mpfr_init(r34399);
        mpfr_init(r34400);
        mpfr_init(r34401);
        mpfr_init(r34402);
        mpfr_init(r34403);
        mpfr_init(r34404);
        mpfr_init(r34405);
}

double f_fm(double x, double y, double z) {
        mpfr_set_d(r34393, x, MPFR_RNDN);
        mpfr_set_d(r34394, y, MPFR_RNDN);
        mpfr_cos(r34395, r34394, MPFR_RNDN);
        mpfr_mul(r34396, r34395, r34395, MPFR_RNDN);
        ;
        mpfr_pow(r34398, r34396, r34397, MPFR_RNDN);
        mpfr_mul(r34399, r34393, r34398, MPFR_RNDN);
        mpfr_cbrt(r34400, r34395, MPFR_RNDN);
        mpfr_mul(r34401, r34399, r34400, MPFR_RNDN);
        mpfr_set_d(r34402, z, MPFR_RNDN);
        mpfr_sin(r34403, r34394, MPFR_RNDN);
        mpfr_mul(r34404, r34402, r34403, MPFR_RNDN);
        mpfr_add(r34405, r34401, r34404, MPFR_RNDN);
        return mpfr_get_d(r34405, MPFR_RNDN);
}

static mpfr_t r34406, r34407, r34408, r34409, r34410, r34411, r34412, r34413, r34414, r34415, r34416, r34417, r34418;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r34406);
        mpfr_init(r34407);
        mpfr_init(r34408);
        mpfr_init(r34409);
        mpfr_init_set_str(r34410, "1/3", 10, MPFR_RNDN);
        mpfr_init(r34411);
        mpfr_init(r34412);
        mpfr_init(r34413);
        mpfr_init(r34414);
        mpfr_init(r34415);
        mpfr_init(r34416);
        mpfr_init(r34417);
        mpfr_init(r34418);
}

double f_dm(double x, double y, double z) {
        mpfr_set_d(r34406, x, MPFR_RNDN);
        mpfr_set_d(r34407, y, MPFR_RNDN);
        mpfr_cos(r34408, r34407, MPFR_RNDN);
        mpfr_mul(r34409, r34408, r34408, MPFR_RNDN);
        ;
        mpfr_pow(r34411, r34409, r34410, MPFR_RNDN);
        mpfr_mul(r34412, r34406, r34411, MPFR_RNDN);
        mpfr_cbrt(r34413, r34408, MPFR_RNDN);
        mpfr_mul(r34414, r34412, r34413, MPFR_RNDN);
        mpfr_set_d(r34415, z, MPFR_RNDN);
        mpfr_sin(r34416, r34407, MPFR_RNDN);
        mpfr_mul(r34417, r34415, r34416, MPFR_RNDN);
        mpfr_add(r34418, r34414, r34417, MPFR_RNDN);
        return mpfr_get_d(r34418, MPFR_RNDN);
}