#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 r26069 = x;
        float r26070 = y;
        float r26071 = r26069 * r26070;
        float r26072 = z;
        float r26073 = t;
        float r26074 = r26072 * r26073;
        float r26075 = r26071 + r26074;
        float r26076 = a;
        float r26077 = b;
        float r26078 = r26076 * r26077;
        float r26079 = r26075 + r26078;
        float r26080 = c;
        float r26081 = i;
        float r26082 = r26080 * r26081;
        float r26083 = r26079 + r26082;
        return r26083;
}

double f_id(double x, double y, double z, double t, double a, double b, double c, double i) {
        double r26084 = x;
        double r26085 = y;
        double r26086 = r26084 * r26085;
        double r26087 = z;
        double r26088 = t;
        double r26089 = r26087 * r26088;
        double r26090 = r26086 + r26089;
        double r26091 = a;
        double r26092 = b;
        double r26093 = r26091 * r26092;
        double r26094 = r26090 + r26093;
        double r26095 = c;
        double r26096 = i;
        double r26097 = r26095 * r26096;
        double r26098 = r26094 + r26097;
        return r26098;
}


double f_of(float x, float y, float z, float t, float a, float b, float c, float i) {
        float r26099 = x;
        float r26100 = y;
        float r26101 = r26099 * r26100;
        float r26102 = z;
        float r26103 = t;
        float r26104 = r26102 * r26103;
        float r26105 = r26101 + r26104;
        float r26106 = a;
        float r26107 = b;
        float r26108 = r26106 * r26107;
        float r26109 = r26105 + r26108;
        float r26110 = c;
        float r26111 = i;
        float r26112 = r26110 * r26111;
        float r26113 = r26109 + r26112;
        return r26113;
}

double f_od(double x, double y, double z, double t, double a, double b, double c, double i) {
        double r26114 = x;
        double r26115 = y;
        double r26116 = r26114 * r26115;
        double r26117 = z;
        double r26118 = t;
        double r26119 = r26117 * r26118;
        double r26120 = r26116 + r26119;
        double r26121 = a;
        double r26122 = b;
        double r26123 = r26121 * r26122;
        double r26124 = r26120 + r26123;
        double r26125 = c;
        double r26126 = i;
        double r26127 = r26125 * r26126;
        double r26128 = r26124 + r26127;
        return r26128;
}

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 r26129, r26130, r26131, r26132, r26133, r26134, r26135, r26136, r26137, r26138, r26139, r26140, r26141, r26142, r26143;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r26129);
        mpfr_init(r26130);
        mpfr_init(r26131);
        mpfr_init(r26132);
        mpfr_init(r26133);
        mpfr_init(r26134);
        mpfr_init(r26135);
        mpfr_init(r26136);
        mpfr_init(r26137);
        mpfr_init(r26138);
        mpfr_init(r26139);
        mpfr_init(r26140);
        mpfr_init(r26141);
        mpfr_init(r26142);
        mpfr_init(r26143);
}

double f_im(double x, double y, double z, double t, double a, double b, double c, double i) {
        mpfr_set_d(r26129, x, MPFR_RNDN);
        mpfr_set_d(r26130, y, MPFR_RNDN);
        mpfr_mul(r26131, r26129, r26130, MPFR_RNDN);
        mpfr_set_d(r26132, z, MPFR_RNDN);
        mpfr_set_d(r26133, t, MPFR_RNDN);
        mpfr_mul(r26134, r26132, r26133, MPFR_RNDN);
        mpfr_add(r26135, r26131, r26134, MPFR_RNDN);
        mpfr_set_d(r26136, a, MPFR_RNDN);
        mpfr_set_d(r26137, b, MPFR_RNDN);
        mpfr_mul(r26138, r26136, r26137, MPFR_RNDN);
        mpfr_add(r26139, r26135, r26138, MPFR_RNDN);
        mpfr_set_d(r26140, c, MPFR_RNDN);
        mpfr_set_d(r26141, i, MPFR_RNDN);
        mpfr_mul(r26142, r26140, r26141, MPFR_RNDN);
        mpfr_add(r26143, r26139, r26142, MPFR_RNDN);
        return mpfr_get_d(r26143, MPFR_RNDN);
}

static mpfr_t r26144, r26145, r26146, r26147, r26148, r26149, r26150, r26151, r26152, r26153, r26154, r26155, r26156, r26157, r26158;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26144);
        mpfr_init(r26145);
        mpfr_init(r26146);
        mpfr_init(r26147);
        mpfr_init(r26148);
        mpfr_init(r26149);
        mpfr_init(r26150);
        mpfr_init(r26151);
        mpfr_init(r26152);
        mpfr_init(r26153);
        mpfr_init(r26154);
        mpfr_init(r26155);
        mpfr_init(r26156);
        mpfr_init(r26157);
        mpfr_init(r26158);
}

double f_fm(double x, double y, double z, double t, double a, double b, double c, double i) {
        mpfr_set_d(r26144, x, MPFR_RNDN);
        mpfr_set_d(r26145, y, MPFR_RNDN);
        mpfr_mul(r26146, r26144, r26145, MPFR_RNDN);
        mpfr_set_d(r26147, z, MPFR_RNDN);
        mpfr_set_d(r26148, t, MPFR_RNDN);
        mpfr_mul(r26149, r26147, r26148, MPFR_RNDN);
        mpfr_add(r26150, r26146, r26149, MPFR_RNDN);
        mpfr_set_d(r26151, a, MPFR_RNDN);
        mpfr_set_d(r26152, b, MPFR_RNDN);
        mpfr_mul(r26153, r26151, r26152, MPFR_RNDN);
        mpfr_add(r26154, r26150, r26153, MPFR_RNDN);
        mpfr_set_d(r26155, c, MPFR_RNDN);
        mpfr_set_d(r26156, i, MPFR_RNDN);
        mpfr_mul(r26157, r26155, r26156, MPFR_RNDN);
        mpfr_add(r26158, r26154, r26157, MPFR_RNDN);
        return mpfr_get_d(r26158, MPFR_RNDN);
}

static mpfr_t r26159, r26160, r26161, r26162, r26163, r26164, r26165, r26166, r26167, r26168, r26169, r26170, r26171, r26172, r26173;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26159);
        mpfr_init(r26160);
        mpfr_init(r26161);
        mpfr_init(r26162);
        mpfr_init(r26163);
        mpfr_init(r26164);
        mpfr_init(r26165);
        mpfr_init(r26166);
        mpfr_init(r26167);
        mpfr_init(r26168);
        mpfr_init(r26169);
        mpfr_init(r26170);
        mpfr_init(r26171);
        mpfr_init(r26172);
        mpfr_init(r26173);
}

double f_dm(double x, double y, double z, double t, double a, double b, double c, double i) {
        mpfr_set_d(r26159, x, MPFR_RNDN);
        mpfr_set_d(r26160, y, MPFR_RNDN);
        mpfr_mul(r26161, r26159, r26160, MPFR_RNDN);
        mpfr_set_d(r26162, z, MPFR_RNDN);
        mpfr_set_d(r26163, t, MPFR_RNDN);
        mpfr_mul(r26164, r26162, r26163, MPFR_RNDN);
        mpfr_add(r26165, r26161, r26164, MPFR_RNDN);
        mpfr_set_d(r26166, a, MPFR_RNDN);
        mpfr_set_d(r26167, b, MPFR_RNDN);
        mpfr_mul(r26168, r26166, r26167, MPFR_RNDN);
        mpfr_add(r26169, r26165, r26168, MPFR_RNDN);
        mpfr_set_d(r26170, c, MPFR_RNDN);
        mpfr_set_d(r26171, i, MPFR_RNDN);
        mpfr_mul(r26172, r26170, r26171, MPFR_RNDN);
        mpfr_add(r26173, r26169, r26172, MPFR_RNDN);
        return mpfr_get_d(r26173, MPFR_RNDN);
}

