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

char *name = "Rectangular parallelepiped of dimension a×b×c";

double f_if() {
        float r41198 = 2;
        float r41199 = 1;
        float r41200 = 9;
        float r41201 = r41199 / r41200;
        float r41202 = r41199 * r41201;
        float r41203 = r41201 * r41201;
        float r41204 = r41202 + r41203;
        float r41205 = r41201 * r41199;
        float r41206 = r41204 + r41205;
        float r41207 = r41198 * r41206;
        return r41207;
}

double f_id() {
        double r41208 = 2;
        double r41209 = 1;
        double r41210 = 9;
        double r41211 = r41209 / r41210;
        double r41212 = r41209 * r41211;
        double r41213 = r41211 * r41211;
        double r41214 = r41212 + r41213;
        double r41215 = r41211 * r41209;
        double r41216 = r41214 + r41215;
        double r41217 = r41208 * r41216;
        return r41217;
}


double f_of() {
        float r41218 = 2;
        float r41219 = 1;
        float r41220 = 9;
        float r41221 = r41219 / r41220;
        float r41222 = r41219 * r41221;
        float r41223 = r41221 * r41221;
        float r41224 = r41222 + r41223;
        float r41225 = r41221 * r41219;
        float r41226 = r41224 + r41225;
        float r41227 = r41218 * r41226;
        return r41227;
}

double f_od() {
        double r41228 = 2;
        double r41229 = 1;
        double r41230 = 9;
        double r41231 = r41229 / r41230;
        double r41232 = r41229 * r41231;
        double r41233 = r41231 * r41231;
        double r41234 = r41232 + r41233;
        double r41235 = r41231 * r41229;
        double r41236 = r41234 + r41235;
        double r41237 = r41228 * r41236;
        return r41237;
}

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 r41238, r41239, r41240, r41241, r41242, r41243, r41244, r41245, r41246, r41247;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r41238, "2", 10, MPFR_RNDN);
        mpfr_init_set_str(r41239, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r41240, "9", 10, MPFR_RNDN);
        mpfr_init(r41241);
        mpfr_init(r41242);
        mpfr_init(r41243);
        mpfr_init(r41244);
        mpfr_init(r41245);
        mpfr_init(r41246);
        mpfr_init(r41247);
}

double f_im() {
        ;
        ;
        ;
        mpfr_div(r41241, r41239, r41240, MPFR_RNDN);
        mpfr_mul(r41242, r41239, r41241, MPFR_RNDN);
        mpfr_mul(r41243, r41241, r41241, MPFR_RNDN);
        mpfr_add(r41244, r41242, r41243, MPFR_RNDN);
        mpfr_mul(r41245, r41241, r41239, MPFR_RNDN);
        mpfr_add(r41246, r41244, r41245, MPFR_RNDN);
        mpfr_mul(r41247, r41238, r41246, MPFR_RNDN);
        return mpfr_get_d(r41247, MPFR_RNDN);
}

static mpfr_t r41248, r41249, r41250, r41251, r41252, r41253, r41254, r41255, r41256, r41257;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r41248, "2", 10, MPFR_RNDN);
        mpfr_init_set_str(r41249, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r41250, "9", 10, MPFR_RNDN);
        mpfr_init(r41251);
        mpfr_init(r41252);
        mpfr_init(r41253);
        mpfr_init(r41254);
        mpfr_init(r41255);
        mpfr_init(r41256);
        mpfr_init(r41257);
}

double f_fm() {
        ;
        ;
        ;
        mpfr_div(r41251, r41249, r41250, MPFR_RNDN);
        mpfr_mul(r41252, r41249, r41251, MPFR_RNDN);
        mpfr_mul(r41253, r41251, r41251, MPFR_RNDN);
        mpfr_add(r41254, r41252, r41253, MPFR_RNDN);
        mpfr_mul(r41255, r41251, r41249, MPFR_RNDN);
        mpfr_add(r41256, r41254, r41255, MPFR_RNDN);
        mpfr_mul(r41257, r41248, r41256, MPFR_RNDN);
        return mpfr_get_d(r41257, MPFR_RNDN);
}

static mpfr_t r41258, r41259, r41260, r41261, r41262, r41263, r41264, r41265, r41266, r41267;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r41258, "2", 10, MPFR_RNDN);
        mpfr_init_set_str(r41259, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r41260, "9", 10, MPFR_RNDN);
        mpfr_init(r41261);
        mpfr_init(r41262);
        mpfr_init(r41263);
        mpfr_init(r41264);
        mpfr_init(r41265);
        mpfr_init(r41266);
        mpfr_init(r41267);
}

double f_dm() {
        ;
        ;
        ;
        mpfr_div(r41261, r41259, r41260, MPFR_RNDN);
        mpfr_mul(r41262, r41259, r41261, MPFR_RNDN);
        mpfr_mul(r41263, r41261, r41261, MPFR_RNDN);
        mpfr_add(r41264, r41262, r41263, MPFR_RNDN);
        mpfr_mul(r41265, r41261, r41259, MPFR_RNDN);
        mpfr_add(r41266, r41264, r41265, MPFR_RNDN);
        mpfr_mul(r41267, r41258, r41266, MPFR_RNDN);
        return mpfr_get_d(r41267, MPFR_RNDN);
}

