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

char *name = "Linear.V3:cross from linear-1.19.1.3";

double f_if(float x, float y, float z, float t) {
        float r43736 = x;
        float r43737 = y;
        float r43738 = r43736 * r43737;
        float r43739 = z;
        float r43740 = t;
        float r43741 = r43739 * r43740;
        float r43742 = r43738 - r43741;
        return r43742;
}

double f_id(double x, double y, double z, double t) {
        double r43743 = x;
        double r43744 = y;
        double r43745 = r43743 * r43744;
        double r43746 = z;
        double r43747 = t;
        double r43748 = r43746 * r43747;
        double r43749 = r43745 - r43748;
        return r43749;
}


double f_of(float x, float y, float z, float t) {
        float r43750 = x;
        float r43751 = y;
        float r43752 = r43750 * r43751;
        float r43753 = z;
        float r43754 = t;
        float r43755 = r43753 * r43754;
        float r43756 = r43752 - r43755;
        return r43756;
}

double f_od(double x, double y, double z, double t) {
        double r43757 = x;
        double r43758 = y;
        double r43759 = r43757 * r43758;
        double r43760 = z;
        double r43761 = t;
        double r43762 = r43760 * r43761;
        double r43763 = r43759 - r43762;
        return r43763;
}

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 r43764, r43765, r43766, r43767, r43768, r43769, r43770;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r43764);
        mpfr_init(r43765);
        mpfr_init(r43766);
        mpfr_init(r43767);
        mpfr_init(r43768);
        mpfr_init(r43769);
        mpfr_init(r43770);
}

double f_im(double x, double y, double z, double t) {
        mpfr_set_d(r43764, x, MPFR_RNDN);
        mpfr_set_d(r43765, y, MPFR_RNDN);
        mpfr_mul(r43766, r43764, r43765, MPFR_RNDN);
        mpfr_set_d(r43767, z, MPFR_RNDN);
        mpfr_set_d(r43768, t, MPFR_RNDN);
        mpfr_mul(r43769, r43767, r43768, MPFR_RNDN);
        mpfr_sub(r43770, r43766, r43769, MPFR_RNDN);
        return mpfr_get_d(r43770, MPFR_RNDN);
}

static mpfr_t r43771, r43772, r43773, r43774, r43775, r43776, r43777;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r43771);
        mpfr_init(r43772);
        mpfr_init(r43773);
        mpfr_init(r43774);
        mpfr_init(r43775);
        mpfr_init(r43776);
        mpfr_init(r43777);
}

double f_fm(double x, double y, double z, double t) {
        mpfr_set_d(r43771, x, MPFR_RNDN);
        mpfr_set_d(r43772, y, MPFR_RNDN);
        mpfr_mul(r43773, r43771, r43772, MPFR_RNDN);
        mpfr_set_d(r43774, z, MPFR_RNDN);
        mpfr_set_d(r43775, t, MPFR_RNDN);
        mpfr_mul(r43776, r43774, r43775, MPFR_RNDN);
        mpfr_sub(r43777, r43773, r43776, MPFR_RNDN);
        return mpfr_get_d(r43777, MPFR_RNDN);
}

static mpfr_t r43778, r43779, r43780, r43781, r43782, r43783, r43784;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r43778);
        mpfr_init(r43779);
        mpfr_init(r43780);
        mpfr_init(r43781);
        mpfr_init(r43782);
        mpfr_init(r43783);
        mpfr_init(r43784);
}

double f_dm(double x, double y, double z, double t) {
        mpfr_set_d(r43778, x, MPFR_RNDN);
        mpfr_set_d(r43779, y, MPFR_RNDN);
        mpfr_mul(r43780, r43778, r43779, MPFR_RNDN);
        mpfr_set_d(r43781, z, MPFR_RNDN);
        mpfr_set_d(r43782, t, MPFR_RNDN);
        mpfr_mul(r43783, r43781, r43782, MPFR_RNDN);
        mpfr_sub(r43784, r43780, r43783, MPFR_RNDN);
        return mpfr_get_d(r43784, MPFR_RNDN);
}