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

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

double f_if(float x, float y, float z, float t) {
        float r29338 = x;
        float r29339 = y;
        float r29340 = r29338 * r29339;
        float r29341 = z;
        float r29342 = t;
        float r29343 = r29341 * r29342;
        float r29344 = r29340 + r29343;
        return r29344;
}

double f_id(double x, double y, double z, double t) {
        double r29345 = x;
        double r29346 = y;
        double r29347 = r29345 * r29346;
        double r29348 = z;
        double r29349 = t;
        double r29350 = r29348 * r29349;
        double r29351 = r29347 + r29350;
        return r29351;
}


double f_of(float x, float y, float z, float t) {
        float r29352 = x;
        float r29353 = y;
        float r29354 = r29352 * r29353;
        float r29355 = z;
        float r29356 = t;
        float r29357 = r29355 * r29356;
        float r29358 = r29354 + r29357;
        return r29358;
}

double f_od(double x, double y, double z, double t) {
        double r29359 = x;
        double r29360 = y;
        double r29361 = r29359 * r29360;
        double r29362 = z;
        double r29363 = t;
        double r29364 = r29362 * r29363;
        double r29365 = r29361 + r29364;
        return r29365;
}

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

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

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

static mpfr_t r29373, r29374, r29375, r29376, r29377, r29378, r29379;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r29373);
        mpfr_init(r29374);
        mpfr_init(r29375);
        mpfr_init(r29376);
        mpfr_init(r29377);
        mpfr_init(r29378);
        mpfr_init(r29379);
}

double f_fm(double x, double y, double z, double t) {
        mpfr_set_d(r29373, x, MPFR_RNDN);
        mpfr_set_d(r29374, y, MPFR_RNDN);
        mpfr_mul(r29375, r29373, r29374, MPFR_RNDN);
        mpfr_set_d(r29376, z, MPFR_RNDN);
        mpfr_set_d(r29377, t, MPFR_RNDN);
        mpfr_mul(r29378, r29376, r29377, MPFR_RNDN);
        mpfr_add(r29379, r29375, r29378, MPFR_RNDN);
        return mpfr_get_d(r29379, MPFR_RNDN);
}

static mpfr_t r29380, r29381, r29382, r29383, r29384, r29385, r29386;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r29380);
        mpfr_init(r29381);
        mpfr_init(r29382);
        mpfr_init(r29383);
        mpfr_init(r29384);
        mpfr_init(r29385);
        mpfr_init(r29386);
}

double f_dm(double x, double y, double z, double t) {
        mpfr_set_d(r29380, x, MPFR_RNDN);
        mpfr_set_d(r29381, y, MPFR_RNDN);
        mpfr_mul(r29382, r29380, r29381, MPFR_RNDN);
        mpfr_set_d(r29383, z, MPFR_RNDN);
        mpfr_set_d(r29384, t, MPFR_RNDN);
        mpfr_mul(r29385, r29383, r29384, MPFR_RNDN);
        mpfr_add(r29386, r29382, r29385, MPFR_RNDN);
        return mpfr_get_d(r29386, MPFR_RNDN);
}