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

char *name = "Quotient of products";

double f_if(float a1, float a2, float b1, float b2) {
        float r24373 = a1;
        float r24374 = a2;
        float r24375 = r24373 * r24374;
        float r24376 = b1;
        float r24377 = b2;
        float r24378 = r24376 * r24377;
        float r24379 = r24375 / r24378;
        return r24379;
}

double f_id(double a1, double a2, double b1, double b2) {
        double r24380 = a1;
        double r24381 = a2;
        double r24382 = r24380 * r24381;
        double r24383 = b1;
        double r24384 = b2;
        double r24385 = r24383 * r24384;
        double r24386 = r24382 / r24385;
        return r24386;
}


double f_of(float a1, float a2, float b1, float b2) {
        float r24387 = a1;
        float r24388 = b1;
        float r24389 = r24387 / r24388;
        float r24390 = a2;
        float r24391 = b2;
        float r24392 = r24390 / r24391;
        float r24393 = r24389 * r24392;
        float r24394 = -inf.0;
        bool r24395 = r24393 <= r24394;
        float r24396 = r24388 / r24390;
        float r24397 = r24387 / r24396;
        float r24398 = r24397 / r24391;
        float r24399 = -6.1104555447917e-315;
        bool r24400 = r24393 <= r24399;
        float r24401 = 1.2496181882040065e-305;
        bool r24402 = r24393 <= r24401;
        float r24403 = 1;
        float r24404 = r24388 * r24391;
        float r24405 = r24387 * r24390;
        float r24406 = r24404 / r24405;
        float r24407 = r24403 / r24406;
        float r24408 = 5.4538264489329024e+302;
        bool r24409 = r24393 <= r24408;
        float r24410 = r24409 ? r24393 : r24407;
        float r24411 = r24402 ? r24407 : r24410;
        float r24412 = r24400 ? r24393 : r24411;
        float r24413 = r24395 ? r24398 : r24412;
        return r24413;
}

double f_od(double a1, double a2, double b1, double b2) {
        double r24414 = a1;
        double r24415 = b1;
        double r24416 = r24414 / r24415;
        double r24417 = a2;
        double r24418 = b2;
        double r24419 = r24417 / r24418;
        double r24420 = r24416 * r24419;
        double r24421 = -inf.0;
        bool r24422 = r24420 <= r24421;
        double r24423 = r24415 / r24417;
        double r24424 = r24414 / r24423;
        double r24425 = r24424 / r24418;
        double r24426 = -6.1104555447917e-315;
        bool r24427 = r24420 <= r24426;
        double r24428 = 1.2496181882040065e-305;
        bool r24429 = r24420 <= r24428;
        double r24430 = 1;
        double r24431 = r24415 * r24418;
        double r24432 = r24414 * r24417;
        double r24433 = r24431 / r24432;
        double r24434 = r24430 / r24433;
        double r24435 = 5.4538264489329024e+302;
        bool r24436 = r24420 <= r24435;
        double r24437 = r24436 ? r24420 : r24434;
        double r24438 = r24429 ? r24434 : r24437;
        double r24439 = r24427 ? r24420 : r24438;
        double r24440 = r24422 ? r24425 : r24439;
        return r24440;
}

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 r24441, r24442, r24443, r24444, r24445, r24446, r24447;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r24441);
        mpfr_init(r24442);
        mpfr_init(r24443);
        mpfr_init(r24444);
        mpfr_init(r24445);
        mpfr_init(r24446);
        mpfr_init(r24447);
}

double f_im(double a1, double a2, double b1, double b2) {
        mpfr_set_d(r24441, a1, MPFR_RNDN);
        mpfr_set_d(r24442, a2, MPFR_RNDN);
        mpfr_mul(r24443, r24441, r24442, MPFR_RNDN);
        mpfr_set_d(r24444, b1, MPFR_RNDN);
        mpfr_set_d(r24445, b2, MPFR_RNDN);
        mpfr_mul(r24446, r24444, r24445, MPFR_RNDN);
        mpfr_div(r24447, r24443, r24446, MPFR_RNDN);
        return mpfr_get_d(r24447, MPFR_RNDN);
}

static mpfr_t r24448, r24449, r24450, r24451, r24452, r24453, r24454, r24455, r24456, r24457, r24458, r24459, r24460, r24461, r24462, r24463, r24464, r24465, r24466, r24467, r24468, r24469, r24470, r24471, r24472, r24473, r24474;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r24448);
        mpfr_init(r24449);
        mpfr_init(r24450);
        mpfr_init(r24451);
        mpfr_init(r24452);
        mpfr_init(r24453);
        mpfr_init(r24454);
        mpfr_init_set_str(r24455, "-inf.0", 10, MPFR_RNDN);
        mpfr_init(r24456);
        mpfr_init(r24457);
        mpfr_init(r24458);
        mpfr_init(r24459);
        mpfr_init_set_str(r24460, "-6.1104555447917e-315", 10, MPFR_RNDN);
        mpfr_init(r24461);
        mpfr_init_set_str(r24462, "1.2496181882040065e-305", 10, MPFR_RNDN);
        mpfr_init(r24463);
        mpfr_init_set_str(r24464, "1", 10, MPFR_RNDN);
        mpfr_init(r24465);
        mpfr_init(r24466);
        mpfr_init(r24467);
        mpfr_init(r24468);
        mpfr_init_set_str(r24469, "5.4538264489329024e+302", 10, MPFR_RNDN);
        mpfr_init(r24470);
        mpfr_init(r24471);
        mpfr_init(r24472);
        mpfr_init(r24473);
        mpfr_init(r24474);
}

double f_fm(double a1, double a2, double b1, double b2) {
        mpfr_set_d(r24448, a1, MPFR_RNDN);
        mpfr_set_d(r24449, b1, MPFR_RNDN);
        mpfr_div(r24450, r24448, r24449, MPFR_RNDN);
        mpfr_set_d(r24451, a2, MPFR_RNDN);
        mpfr_set_d(r24452, b2, MPFR_RNDN);
        mpfr_div(r24453, r24451, r24452, MPFR_RNDN);
        mpfr_mul(r24454, r24450, r24453, MPFR_RNDN);
        ;
        mpfr_set_si(r24456, mpfr_cmp(r24454, r24455) <= 0, MPFR_RNDN);
        mpfr_div(r24457, r24449, r24451, MPFR_RNDN);
        mpfr_div(r24458, r24448, r24457, MPFR_RNDN);
        mpfr_div(r24459, r24458, r24452, MPFR_RNDN);
        ;
        mpfr_set_si(r24461, mpfr_cmp(r24454, r24460) <= 0, MPFR_RNDN);
        ;
        mpfr_set_si(r24463, mpfr_cmp(r24454, r24462) <= 0, MPFR_RNDN);
        ;
        mpfr_mul(r24465, r24449, r24452, MPFR_RNDN);
        mpfr_mul(r24466, r24448, r24451, MPFR_RNDN);
        mpfr_div(r24467, r24465, r24466, MPFR_RNDN);
        mpfr_div(r24468, r24464, r24467, MPFR_RNDN);
        ;
        mpfr_set_si(r24470, mpfr_cmp(r24454, r24469) <= 0, MPFR_RNDN);
        if (mpfr_get_si(r24470, MPFR_RNDN)) { mpfr_set(r24471, r24454, MPFR_RNDN); } else { mpfr_set(r24471, r24468, MPFR_RNDN); };
        if (mpfr_get_si(r24463, MPFR_RNDN)) { mpfr_set(r24472, r24468, MPFR_RNDN); } else { mpfr_set(r24472, r24471, MPFR_RNDN); };
        if (mpfr_get_si(r24461, MPFR_RNDN)) { mpfr_set(r24473, r24454, MPFR_RNDN); } else { mpfr_set(r24473, r24472, MPFR_RNDN); };
        if (mpfr_get_si(r24456, MPFR_RNDN)) { mpfr_set(r24474, r24459, MPFR_RNDN); } else { mpfr_set(r24474, r24473, MPFR_RNDN); };
        return mpfr_get_d(r24474, MPFR_RNDN);
}

static mpfr_t r24475, r24476, r24477, r24478, r24479, r24480, r24481, r24482, r24483, r24484, r24485, r24486, r24487, r24488, r24489, r24490, r24491, r24492, r24493, r24494, r24495, r24496, r24497, r24498, r24499, r24500, r24501;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r24475);
        mpfr_init(r24476);
        mpfr_init(r24477);
        mpfr_init(r24478);
        mpfr_init(r24479);
        mpfr_init(r24480);
        mpfr_init(r24481);
        mpfr_init_set_str(r24482, "-inf.0", 10, MPFR_RNDN);
        mpfr_init(r24483);
        mpfr_init(r24484);
        mpfr_init(r24485);
        mpfr_init(r24486);
        mpfr_init_set_str(r24487, "-6.1104555447917e-315", 10, MPFR_RNDN);
        mpfr_init(r24488);
        mpfr_init_set_str(r24489, "1.2496181882040065e-305", 10, MPFR_RNDN);
        mpfr_init(r24490);
        mpfr_init_set_str(r24491, "1", 10, MPFR_RNDN);
        mpfr_init(r24492);
        mpfr_init(r24493);
        mpfr_init(r24494);
        mpfr_init(r24495);
        mpfr_init_set_str(r24496, "5.4538264489329024e+302", 10, MPFR_RNDN);
        mpfr_init(r24497);
        mpfr_init(r24498);
        mpfr_init(r24499);
        mpfr_init(r24500);
        mpfr_init(r24501);
}

double f_dm(double a1, double a2, double b1, double b2) {
        mpfr_set_d(r24475, a1, MPFR_RNDN);
        mpfr_set_d(r24476, b1, MPFR_RNDN);
        mpfr_div(r24477, r24475, r24476, MPFR_RNDN);
        mpfr_set_d(r24478, a2, MPFR_RNDN);
        mpfr_set_d(r24479, b2, MPFR_RNDN);
        mpfr_div(r24480, r24478, r24479, MPFR_RNDN);
        mpfr_mul(r24481, r24477, r24480, MPFR_RNDN);
        ;
        mpfr_set_si(r24483, mpfr_cmp(r24481, r24482) <= 0, MPFR_RNDN);
        mpfr_div(r24484, r24476, r24478, MPFR_RNDN);
        mpfr_div(r24485, r24475, r24484, MPFR_RNDN);
        mpfr_div(r24486, r24485, r24479, MPFR_RNDN);
        ;
        mpfr_set_si(r24488, mpfr_cmp(r24481, r24487) <= 0, MPFR_RNDN);
        ;
        mpfr_set_si(r24490, mpfr_cmp(r24481, r24489) <= 0, MPFR_RNDN);
        ;
        mpfr_mul(r24492, r24476, r24479, MPFR_RNDN);
        mpfr_mul(r24493, r24475, r24478, MPFR_RNDN);
        mpfr_div(r24494, r24492, r24493, MPFR_RNDN);
        mpfr_div(r24495, r24491, r24494, MPFR_RNDN);
        ;
        mpfr_set_si(r24497, mpfr_cmp(r24481, r24496) <= 0, MPFR_RNDN);
        if (mpfr_get_si(r24497, MPFR_RNDN)) { mpfr_set(r24498, r24481, MPFR_RNDN); } else { mpfr_set(r24498, r24495, MPFR_RNDN); };
        if (mpfr_get_si(r24490, MPFR_RNDN)) { mpfr_set(r24499, r24495, MPFR_RNDN); } else { mpfr_set(r24499, r24498, MPFR_RNDN); };
        if (mpfr_get_si(r24488, MPFR_RNDN)) { mpfr_set(r24500, r24481, MPFR_RNDN); } else { mpfr_set(r24500, r24499, MPFR_RNDN); };
        if (mpfr_get_si(r24483, MPFR_RNDN)) { mpfr_set(r24501, r24486, MPFR_RNDN); } else { mpfr_set(r24501, r24500, MPFR_RNDN); };
        return mpfr_get_d(r24501, MPFR_RNDN);
}