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

char *name = "Optimisation.CirclePacking:place from circle-packing-0.1.0.4, B";

double f_if(float x, float y, float z, float t) {
        float r26242 = x;
        float r26243 = y;
        float r26244 = r26242 + r26243;
        float r26245 = z;
        float r26246 = r26244 - r26245;
        float r26247 = t;
        float r26248 = 2.0;
        float r26249 = r26247 * r26248;
        float r26250 = r26246 / r26249;
        return r26250;
}

double f_id(double x, double y, double z, double t) {
        double r26251 = x;
        double r26252 = y;
        double r26253 = r26251 + r26252;
        double r26254 = z;
        double r26255 = r26253 - r26254;
        double r26256 = t;
        double r26257 = 2.0;
        double r26258 = r26256 * r26257;
        double r26259 = r26255 / r26258;
        return r26259;
}


double f_of(float x, float y, float z, float t) {
        float r26260 = x;
        float r26261 = y;
        float r26262 = r26260 + r26261;
        float r26263 = z;
        float r26264 = r26262 - r26263;
        float r26265 = t;
        float r26266 = 2.0;
        float r26267 = r26265 * r26266;
        float r26268 = r26264 / r26267;
        return r26268;
}

double f_od(double x, double y, double z, double t) {
        double r26269 = x;
        double r26270 = y;
        double r26271 = r26269 + r26270;
        double r26272 = z;
        double r26273 = r26271 - r26272;
        double r26274 = t;
        double r26275 = 2.0;
        double r26276 = r26274 * r26275;
        double r26277 = r26273 / r26276;
        return r26277;
}

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 r26278, r26279, r26280, r26281, r26282, r26283, r26284, r26285, r26286;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r26278);
        mpfr_init(r26279);
        mpfr_init(r26280);
        mpfr_init(r26281);
        mpfr_init(r26282);
        mpfr_init(r26283);
        mpfr_init_set_str(r26284, "2.0", 10, MPFR_RNDN);
        mpfr_init(r26285);
        mpfr_init(r26286);
}

double f_im(double x, double y, double z, double t) {
        mpfr_set_d(r26278, x, MPFR_RNDN);
        mpfr_set_d(r26279, y, MPFR_RNDN);
        mpfr_add(r26280, r26278, r26279, MPFR_RNDN);
        mpfr_set_d(r26281, z, MPFR_RNDN);
        mpfr_sub(r26282, r26280, r26281, MPFR_RNDN);
        mpfr_set_d(r26283, t, MPFR_RNDN);
        ;
        mpfr_mul(r26285, r26283, r26284, MPFR_RNDN);
        mpfr_div(r26286, r26282, r26285, MPFR_RNDN);
        return mpfr_get_d(r26286, MPFR_RNDN);
}

static mpfr_t r26287, r26288, r26289, r26290, r26291, r26292, r26293, r26294, r26295;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r26287);
        mpfr_init(r26288);
        mpfr_init(r26289);
        mpfr_init(r26290);
        mpfr_init(r26291);
        mpfr_init(r26292);
        mpfr_init_set_str(r26293, "2.0", 10, MPFR_RNDN);
        mpfr_init(r26294);
        mpfr_init(r26295);
}

double f_fm(double x, double y, double z, double t) {
        mpfr_set_d(r26287, x, MPFR_RNDN);
        mpfr_set_d(r26288, y, MPFR_RNDN);
        mpfr_add(r26289, r26287, r26288, MPFR_RNDN);
        mpfr_set_d(r26290, z, MPFR_RNDN);
        mpfr_sub(r26291, r26289, r26290, MPFR_RNDN);
        mpfr_set_d(r26292, t, MPFR_RNDN);
        ;
        mpfr_mul(r26294, r26292, r26293, MPFR_RNDN);
        mpfr_div(r26295, r26291, r26294, MPFR_RNDN);
        return mpfr_get_d(r26295, MPFR_RNDN);
}

static mpfr_t r26296, r26297, r26298, r26299, r26300, r26301, r26302, r26303, r26304;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r26296);
        mpfr_init(r26297);
        mpfr_init(r26298);
        mpfr_init(r26299);
        mpfr_init(r26300);
        mpfr_init(r26301);
        mpfr_init_set_str(r26302, "2.0", 10, MPFR_RNDN);
        mpfr_init(r26303);
        mpfr_init(r26304);
}

double f_dm(double x, double y, double z, double t) {
        mpfr_set_d(r26296, x, MPFR_RNDN);
        mpfr_set_d(r26297, y, MPFR_RNDN);
        mpfr_add(r26298, r26296, r26297, MPFR_RNDN);
        mpfr_set_d(r26299, z, MPFR_RNDN);
        mpfr_sub(r26300, r26298, r26299, MPFR_RNDN);
        mpfr_set_d(r26301, t, MPFR_RNDN);
        ;
        mpfr_mul(r26303, r26301, r26302, MPFR_RNDN);
        mpfr_div(r26304, r26300, r26303, MPFR_RNDN);
        return mpfr_get_d(r26304, MPFR_RNDN);
}