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

char *name = "Linear.V4:$cdot from linear-1.19.1.3";

double f_if(float x, float y, float z, float t, float a, float b, float c, float i) {
        float r26538 = x;
        float r26539 = y;
        float r26540 = r26538 * r26539;
        float r26541 = z;
        float r26542 = t;
        float r26543 = r26541 * r26542;
        float r26544 = r26540 + r26543;
        float r26545 = a;
        float r26546 = b;
        float r26547 = r26545 * r26546;
        float r26548 = r26544 + r26547;
        float r26549 = c;
        float r26550 = i;
        float r26551 = r26549 * r26550;
        float r26552 = r26548 + r26551;
        return r26552;
}

double f_id(double x, double y, double z, double t, double a, double b, double c, double i) {
        double r26553 = x;
        double r26554 = y;
        double r26555 = r26553 * r26554;
        double r26556 = z;
        double r26557 = t;
        double r26558 = r26556 * r26557;
        double r26559 = r26555 + r26558;
        double r26560 = a;
        double r26561 = b;
        double r26562 = r26560 * r26561;
        double r26563 = r26559 + r26562;
        double r26564 = c;
        double r26565 = i;
        double r26566 = r26564 * r26565;
        double r26567 = r26563 + r26566;
        return r26567;
}


double f_of(float x, float y, float z, float t, float a, float b, float c, float i) {
        float r26568 = x;
        float r26569 = y;
        float r26570 = r26568 * r26569;
        float r26571 = z;
        float r26572 = t;
        float r26573 = r26571 * r26572;
        float r26574 = r26570 + r26573;
        float r26575 = a;
        float r26576 = b;
        float r26577 = r26575 * r26576;
        float r26578 = r26574 + r26577;
        float r26579 = c;
        float r26580 = i;
        float r26581 = r26579 * r26580;
        float r26582 = r26578 + r26581;
        return r26582;
}

double f_od(double x, double y, double z, double t, double a, double b, double c, double i) {
        double r26583 = x;
        double r26584 = y;
        double r26585 = r26583 * r26584;
        double r26586 = z;
        double r26587 = t;
        double r26588 = r26586 * r26587;
        double r26589 = r26585 + r26588;
        double r26590 = a;
        double r26591 = b;
        double r26592 = r26590 * r26591;
        double r26593 = r26589 + r26592;
        double r26594 = c;
        double r26595 = i;
        double r26596 = r26594 * r26595;
        double r26597 = r26593 + r26596;
        return r26597;
}

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 r26598, r26599, r26600, r26601, r26602, r26603, r26604, r26605, r26606, r26607, r26608, r26609, r26610, r26611, r26612;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r26598);
        mpfr_init(r26599);
        mpfr_init(r26600);
        mpfr_init(r26601);
        mpfr_init(r26602);
        mpfr_init(r26603);
        mpfr_init(r26604);
        mpfr_init(r26605);
        mpfr_init(r26606);
        mpfr_init(r26607);
        mpfr_init(r26608);
        mpfr_init(r26609);
        mpfr_init(r26610);
        mpfr_init(r26611);
        mpfr_init(r26612);
}

double f_im(double x, double y, double z, double t, double a, double b, double c, double i) {
        mpfr_set_d(r26598, x, MPFR_RNDN);
        mpfr_set_d(r26599, y, MPFR_RNDN);
        mpfr_mul(r26600, r26598, r26599, MPFR_RNDN);
        mpfr_set_d(r26601, z, MPFR_RNDN);
        mpfr_set_d(r26602, t, MPFR_RNDN);
        mpfr_mul(r26603, r26601, r26602, MPFR_RNDN);
        mpfr_add(r26604, r26600, r26603, MPFR_RNDN);
        mpfr_set_d(r26605, a, MPFR_RNDN);
        mpfr_set_d(r26606, b, MPFR_RNDN);
        mpfr_mul(r26607, r26605, r26606, MPFR_RNDN);
        mpfr_add(r26608, r26604, r26607, MPFR_RNDN);
        mpfr_set_d(r26609, c, MPFR_RNDN);
        mpfr_set_d(r26610, i, MPFR_RNDN);
        mpfr_mul(r26611, r26609, r26610, MPFR_RNDN);
        mpfr_add(r26612, r26608, r26611, MPFR_RNDN);
        return mpfr_get_d(r26612, MPFR_RNDN);
}

static mpfr_t r26613, r26614, r26615, r26616, r26617, r26618, r26619, r26620, r26621, r26622, r26623, r26624, r26625, r26626, r26627;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26613);
        mpfr_init(r26614);
        mpfr_init(r26615);
        mpfr_init(r26616);
        mpfr_init(r26617);
        mpfr_init(r26618);
        mpfr_init(r26619);
        mpfr_init(r26620);
        mpfr_init(r26621);
        mpfr_init(r26622);
        mpfr_init(r26623);
        mpfr_init(r26624);
        mpfr_init(r26625);
        mpfr_init(r26626);
        mpfr_init(r26627);
}

double f_fm(double x, double y, double z, double t, double a, double b, double c, double i) {
        mpfr_set_d(r26613, x, MPFR_RNDN);
        mpfr_set_d(r26614, y, MPFR_RNDN);
        mpfr_mul(r26615, r26613, r26614, MPFR_RNDN);
        mpfr_set_d(r26616, z, MPFR_RNDN);
        mpfr_set_d(r26617, t, MPFR_RNDN);
        mpfr_mul(r26618, r26616, r26617, MPFR_RNDN);
        mpfr_add(r26619, r26615, r26618, MPFR_RNDN);
        mpfr_set_d(r26620, a, MPFR_RNDN);
        mpfr_set_d(r26621, b, MPFR_RNDN);
        mpfr_mul(r26622, r26620, r26621, MPFR_RNDN);
        mpfr_add(r26623, r26619, r26622, MPFR_RNDN);
        mpfr_set_d(r26624, c, MPFR_RNDN);
        mpfr_set_d(r26625, i, MPFR_RNDN);
        mpfr_mul(r26626, r26624, r26625, MPFR_RNDN);
        mpfr_add(r26627, r26623, r26626, MPFR_RNDN);
        return mpfr_get_d(r26627, MPFR_RNDN);
}

static mpfr_t r26628, r26629, r26630, r26631, r26632, r26633, r26634, r26635, r26636, r26637, r26638, r26639, r26640, r26641, r26642;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26628);
        mpfr_init(r26629);
        mpfr_init(r26630);
        mpfr_init(r26631);
        mpfr_init(r26632);
        mpfr_init(r26633);
        mpfr_init(r26634);
        mpfr_init(r26635);
        mpfr_init(r26636);
        mpfr_init(r26637);
        mpfr_init(r26638);
        mpfr_init(r26639);
        mpfr_init(r26640);
        mpfr_init(r26641);
        mpfr_init(r26642);
}

double f_dm(double x, double y, double z, double t, double a, double b, double c, double i) {
        mpfr_set_d(r26628, x, MPFR_RNDN);
        mpfr_set_d(r26629, y, MPFR_RNDN);
        mpfr_mul(r26630, r26628, r26629, MPFR_RNDN);
        mpfr_set_d(r26631, z, MPFR_RNDN);
        mpfr_set_d(r26632, t, MPFR_RNDN);
        mpfr_mul(r26633, r26631, r26632, MPFR_RNDN);
        mpfr_add(r26634, r26630, r26633, MPFR_RNDN);
        mpfr_set_d(r26635, a, MPFR_RNDN);
        mpfr_set_d(r26636, b, MPFR_RNDN);
        mpfr_mul(r26637, r26635, r26636, MPFR_RNDN);
        mpfr_add(r26638, r26634, r26637, MPFR_RNDN);
        mpfr_set_d(r26639, c, MPFR_RNDN);
        mpfr_set_d(r26640, i, MPFR_RNDN);
        mpfr_mul(r26641, r26639, r26640, MPFR_RNDN);
        mpfr_add(r26642, r26638, r26641, MPFR_RNDN);
        return mpfr_get_d(r26642, MPFR_RNDN);
}

