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

char *name = "Linear.V3:$cdot from linear-1.19.1.3, B";

double f_if(float x, float y, float z, float t, float a, float b) {
        float r29321 = x;
        float r29322 = y;
        float r29323 = r29321 * r29322;
        float r29324 = z;
        float r29325 = t;
        float r29326 = r29324 * r29325;
        float r29327 = r29323 + r29326;
        float r29328 = a;
        float r29329 = b;
        float r29330 = r29328 * r29329;
        float r29331 = r29327 + r29330;
        return r29331;
}

double f_id(double x, double y, double z, double t, double a, double b) {
        double r29332 = x;
        double r29333 = y;
        double r29334 = r29332 * r29333;
        double r29335 = z;
        double r29336 = t;
        double r29337 = r29335 * r29336;
        double r29338 = r29334 + r29337;
        double r29339 = a;
        double r29340 = b;
        double r29341 = r29339 * r29340;
        double r29342 = r29338 + r29341;
        return r29342;
}


double f_of(float x, float y, float z, float t, float a, float b) {
        float r29343 = x;
        float r29344 = y;
        float r29345 = r29343 * r29344;
        float r29346 = z;
        float r29347 = t;
        float r29348 = r29346 * r29347;
        float r29349 = r29345 + r29348;
        float r29350 = a;
        float r29351 = b;
        float r29352 = r29350 * r29351;
        float r29353 = r29349 + r29352;
        return r29353;
}

double f_od(double x, double y, double z, double t, double a, double b) {
        double r29354 = x;
        double r29355 = y;
        double r29356 = r29354 * r29355;
        double r29357 = z;
        double r29358 = t;
        double r29359 = r29357 * r29358;
        double r29360 = r29356 + r29359;
        double r29361 = a;
        double r29362 = b;
        double r29363 = r29361 * r29362;
        double r29364 = r29360 + r29363;
        return r29364;
}

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 r29365, r29366, r29367, r29368, r29369, r29370, r29371, r29372, r29373, r29374, r29375;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r29365);
        mpfr_init(r29366);
        mpfr_init(r29367);
        mpfr_init(r29368);
        mpfr_init(r29369);
        mpfr_init(r29370);
        mpfr_init(r29371);
        mpfr_init(r29372);
        mpfr_init(r29373);
        mpfr_init(r29374);
        mpfr_init(r29375);
}

double f_im(double x, double y, double z, double t, double a, double b) {
        mpfr_set_d(r29365, x, MPFR_RNDN);
        mpfr_set_d(r29366, y, MPFR_RNDN);
        mpfr_mul(r29367, r29365, r29366, MPFR_RNDN);
        mpfr_set_d(r29368, z, MPFR_RNDN);
        mpfr_set_d(r29369, t, MPFR_RNDN);
        mpfr_mul(r29370, r29368, r29369, MPFR_RNDN);
        mpfr_add(r29371, r29367, r29370, MPFR_RNDN);
        mpfr_set_d(r29372, a, MPFR_RNDN);
        mpfr_set_d(r29373, b, MPFR_RNDN);
        mpfr_mul(r29374, r29372, r29373, MPFR_RNDN);
        mpfr_add(r29375, r29371, r29374, MPFR_RNDN);
        return mpfr_get_d(r29375, MPFR_RNDN);
}

static mpfr_t r29376, r29377, r29378, r29379, r29380, r29381, r29382, r29383, r29384, r29385, r29386;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r29376);
        mpfr_init(r29377);
        mpfr_init(r29378);
        mpfr_init(r29379);
        mpfr_init(r29380);
        mpfr_init(r29381);
        mpfr_init(r29382);
        mpfr_init(r29383);
        mpfr_init(r29384);
        mpfr_init(r29385);
        mpfr_init(r29386);
}

double f_fm(double x, double y, double z, double t, double a, double b) {
        mpfr_set_d(r29376, x, MPFR_RNDN);
        mpfr_set_d(r29377, y, MPFR_RNDN);
        mpfr_mul(r29378, r29376, r29377, MPFR_RNDN);
        mpfr_set_d(r29379, z, MPFR_RNDN);
        mpfr_set_d(r29380, t, MPFR_RNDN);
        mpfr_mul(r29381, r29379, r29380, MPFR_RNDN);
        mpfr_add(r29382, r29378, r29381, MPFR_RNDN);
        mpfr_set_d(r29383, a, MPFR_RNDN);
        mpfr_set_d(r29384, b, MPFR_RNDN);
        mpfr_mul(r29385, r29383, r29384, MPFR_RNDN);
        mpfr_add(r29386, r29382, r29385, MPFR_RNDN);
        return mpfr_get_d(r29386, MPFR_RNDN);
}

static mpfr_t r29387, r29388, r29389, r29390, r29391, r29392, r29393, r29394, r29395, r29396, r29397;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r29387);
        mpfr_init(r29388);
        mpfr_init(r29389);
        mpfr_init(r29390);
        mpfr_init(r29391);
        mpfr_init(r29392);
        mpfr_init(r29393);
        mpfr_init(r29394);
        mpfr_init(r29395);
        mpfr_init(r29396);
        mpfr_init(r29397);
}

double f_dm(double x, double y, double z, double t, double a, double b) {
        mpfr_set_d(r29387, x, MPFR_RNDN);
        mpfr_set_d(r29388, y, MPFR_RNDN);
        mpfr_mul(r29389, r29387, r29388, MPFR_RNDN);
        mpfr_set_d(r29390, z, MPFR_RNDN);
        mpfr_set_d(r29391, t, MPFR_RNDN);
        mpfr_mul(r29392, r29390, r29391, MPFR_RNDN);
        mpfr_add(r29393, r29389, r29392, MPFR_RNDN);
        mpfr_set_d(r29394, a, MPFR_RNDN);
        mpfr_set_d(r29395, b, MPFR_RNDN);
        mpfr_mul(r29396, r29394, r29395, MPFR_RNDN);
        mpfr_add(r29397, r29393, r29396, MPFR_RNDN);
        return mpfr_get_d(r29397, MPFR_RNDN);
}