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

char *name = "Diagrams.ThreeD.Shapes:frustum from diagrams-lib-1.3.0.3, B";

double f_if(float x, float y, float z) {
        float r34359 = x;
        float r34360 = y;
        float r34361 = r34360 - r34359;
        float r34362 = z;
        float r34363 = r34361 * r34362;
        float r34364 = r34359 + r34363;
        return r34364;
}

double f_id(double x, double y, double z) {
        double r34365 = x;
        double r34366 = y;
        double r34367 = r34366 - r34365;
        double r34368 = z;
        double r34369 = r34367 * r34368;
        double r34370 = r34365 + r34369;
        return r34370;
}


double f_of(float x, float y, float z) {
        float r34371 = x;
        float r34372 = y;
        float r34373 = r34372 - r34371;
        float r34374 = z;
        float r34375 = r34373 * r34374;
        float r34376 = r34371 + r34375;
        return r34376;
}

double f_od(double x, double y, double z) {
        double r34377 = x;
        double r34378 = y;
        double r34379 = r34378 - r34377;
        double r34380 = z;
        double r34381 = r34379 * r34380;
        double r34382 = r34377 + r34381;
        return r34382;
}

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 r34383, r34384, r34385, r34386, r34387, r34388;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r34383);
        mpfr_init(r34384);
        mpfr_init(r34385);
        mpfr_init(r34386);
        mpfr_init(r34387);
        mpfr_init(r34388);
}

double f_im(double x, double y, double z) {
        mpfr_set_d(r34383, x, MPFR_RNDN);
        mpfr_set_d(r34384, y, MPFR_RNDN);
        mpfr_sub(r34385, r34384, r34383, MPFR_RNDN);
        mpfr_set_d(r34386, z, MPFR_RNDN);
        mpfr_mul(r34387, r34385, r34386, MPFR_RNDN);
        mpfr_add(r34388, r34383, r34387, MPFR_RNDN);
        return mpfr_get_d(r34388, MPFR_RNDN);
}

static mpfr_t r34389, r34390, r34391, r34392, r34393, r34394;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r34389);
        mpfr_init(r34390);
        mpfr_init(r34391);
        mpfr_init(r34392);
        mpfr_init(r34393);
        mpfr_init(r34394);
}

double f_fm(double x, double y, double z) {
        mpfr_set_d(r34389, x, MPFR_RNDN);
        mpfr_set_d(r34390, y, MPFR_RNDN);
        mpfr_sub(r34391, r34390, r34389, MPFR_RNDN);
        mpfr_set_d(r34392, z, MPFR_RNDN);
        mpfr_mul(r34393, r34391, r34392, MPFR_RNDN);
        mpfr_add(r34394, r34389, r34393, MPFR_RNDN);
        return mpfr_get_d(r34394, MPFR_RNDN);
}

static mpfr_t r34395, r34396, r34397, r34398, r34399, r34400;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r34395);
        mpfr_init(r34396);
        mpfr_init(r34397);
        mpfr_init(r34398);
        mpfr_init(r34399);
        mpfr_init(r34400);
}

double f_dm(double x, double y, double z) {
        mpfr_set_d(r34395, x, MPFR_RNDN);
        mpfr_set_d(r34396, y, MPFR_RNDN);
        mpfr_sub(r34397, r34396, r34395, MPFR_RNDN);
        mpfr_set_d(r34398, z, MPFR_RNDN);
        mpfr_mul(r34399, r34397, r34398, MPFR_RNDN);
        mpfr_add(r34400, r34395, r34399, MPFR_RNDN);
        return mpfr_get_d(r34400, MPFR_RNDN);
}