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

char *name = "triangle3";

double f_if(float a, float b, float c) {
        float r19847 = a;
        float r19848 = b;
        float r19849 = r19847 + r19848;
        float r19850 = c;
        float r19851 = r19849 + r19850;
        float r19852 = 2.0f;
        float r19853 = r19851 / r19852;
        float r19854 = r19853 - r19847;
        float r19855 = r19853 * r19854;
        float r19856 = r19853 - r19848;
        float r19857 = r19855 * r19856;
        float r19858 = r19853 - r19850;
        float r19859 = r19857 * r19858;
        float r19860 = sqrt(r19859);
        return r19860;
}

double f_id(double a, double b, double c) {
        double r19861 = a;
        double r19862 = b;
        double r19863 = r19861 + r19862;
        double r19864 = c;
        double r19865 = r19863 + r19864;
        double r19866 = 2.0;
        double r19867 = r19865 / r19866;
        double r19868 = r19867 - r19861;
        double r19869 = r19867 * r19868;
        double r19870 = r19867 - r19862;
        double r19871 = r19869 * r19870;
        double r19872 = r19867 - r19864;
        double r19873 = r19871 * r19872;
        double r19874 = sqrt(r19873);
        return r19874;
}


double f_of(float a, float b, float c) {
        float r19875 = 0.5f;
        float r19876 = r19875 * r19875;
        float r19877 = c;
        float r19878 = a;
        float r19879 = r19877 - r19878;
        float r19880 = b;
        float r19881 = r19879 + r19880;
        float r19882 = r19876 * r19881;
        float r19883 = r19878 - r19877;
        float r19884 = r19880 + r19883;
        float r19885 = r19878 - r19880;
        float r19886 = r19885 + r19877;
        float r19887 = r19884 * r19886;
        float r19888 = r19878 + r19877;
        float r19889 = r19880 + r19888;
        float r19890 = 2.0f;
        float r19891 = r19890 / r19875;
        float r19892 = r19889 / r19891;
        float r19893 = r19887 * r19892;
        float r19894 = r19882 * r19893;
        float r19895 = sqrt(r19894);
        return r19895;
}

double f_od(double a, double b, double c) {
        double r19896 = 0.5;
        double r19897 = r19896 * r19896;
        double r19898 = c;
        double r19899 = a;
        double r19900 = r19898 - r19899;
        double r19901 = b;
        double r19902 = r19900 + r19901;
        double r19903 = r19897 * r19902;
        double r19904 = r19899 - r19898;
        double r19905 = r19901 + r19904;
        double r19906 = r19899 - r19901;
        double r19907 = r19906 + r19898;
        double r19908 = r19905 * r19907;
        double r19909 = r19899 + r19898;
        double r19910 = r19901 + r19909;
        double r19911 = 2.0;
        double r19912 = r19911 / r19896;
        double r19913 = r19910 / r19912;
        double r19914 = r19908 * r19913;
        double r19915 = r19903 * r19914;
        double r19916 = sqrt(r19915);
        return r19916;
}

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 r19917, r19918, r19919, r19920, r19921, r19922, r19923, r19924, r19925, r19926, r19927, r19928, r19929, r19930;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r19917);
        mpfr_init(r19918);
        mpfr_init(r19919);
        mpfr_init(r19920);
        mpfr_init(r19921);
        mpfr_init_set_str(r19922, "2", 10, MPFR_RNDN);
        mpfr_init(r19923);
        mpfr_init(r19924);
        mpfr_init(r19925);
        mpfr_init(r19926);
        mpfr_init(r19927);
        mpfr_init(r19928);
        mpfr_init(r19929);
        mpfr_init(r19930);
}

double f_im(double a, double b, double c) {
        mpfr_set_d(r19917, a, MPFR_RNDN);
        mpfr_set_d(r19918, b, MPFR_RNDN);
        mpfr_add(r19919, r19917, r19918, MPFR_RNDN);
        mpfr_set_d(r19920, c, MPFR_RNDN);
        mpfr_add(r19921, r19919, r19920, MPFR_RNDN);
        ;
        mpfr_div(r19923, r19921, r19922, MPFR_RNDN);
        mpfr_sub(r19924, r19923, r19917, MPFR_RNDN);
        mpfr_mul(r19925, r19923, r19924, MPFR_RNDN);
        mpfr_sub(r19926, r19923, r19918, MPFR_RNDN);
        mpfr_mul(r19927, r19925, r19926, MPFR_RNDN);
        mpfr_sub(r19928, r19923, r19920, MPFR_RNDN);
        mpfr_mul(r19929, r19927, r19928, MPFR_RNDN);
        mpfr_sqrt(r19930, r19929, MPFR_RNDN);
        return mpfr_get_d(r19930, MPFR_RNDN);
}

static mpfr_t r19931, r19932, r19933, r19934, r19935, r19936, r19937, r19938, r19939, r19940, r19941, r19942, r19943, r19944, r19945, r19946, r19947, r19948, r19949, r19950, r19951;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r19931, "1/2", 10, MPFR_RNDN);
        mpfr_init(r19932);
        mpfr_init(r19933);
        mpfr_init(r19934);
        mpfr_init(r19935);
        mpfr_init(r19936);
        mpfr_init(r19937);
        mpfr_init(r19938);
        mpfr_init(r19939);
        mpfr_init(r19940);
        mpfr_init(r19941);
        mpfr_init(r19942);
        mpfr_init(r19943);
        mpfr_init(r19944);
        mpfr_init(r19945);
        mpfr_init_set_str(r19946, "2", 10, MPFR_RNDN);
        mpfr_init(r19947);
        mpfr_init(r19948);
        mpfr_init(r19949);
        mpfr_init(r19950);
        mpfr_init(r19951);
}

double f_fm(double a, double b, double c) {
        ;
        mpfr_mul(r19932, r19931, r19931, MPFR_RNDN);
        mpfr_set_d(r19933, c, MPFR_RNDN);
        mpfr_set_d(r19934, a, MPFR_RNDN);
        mpfr_sub(r19935, r19933, r19934, MPFR_RNDN);
        mpfr_set_d(r19936, b, MPFR_RNDN);
        mpfr_add(r19937, r19935, r19936, MPFR_RNDN);
        mpfr_mul(r19938, r19932, r19937, MPFR_RNDN);
        mpfr_sub(r19939, r19934, r19933, MPFR_RNDN);
        mpfr_add(r19940, r19936, r19939, MPFR_RNDN);
        mpfr_sub(r19941, r19934, r19936, MPFR_RNDN);
        mpfr_add(r19942, r19941, r19933, MPFR_RNDN);
        mpfr_mul(r19943, r19940, r19942, MPFR_RNDN);
        mpfr_add(r19944, r19934, r19933, MPFR_RNDN);
        mpfr_add(r19945, r19936, r19944, MPFR_RNDN);
        ;
        mpfr_div(r19947, r19946, r19931, MPFR_RNDN);
        mpfr_div(r19948, r19945, r19947, MPFR_RNDN);
        mpfr_mul(r19949, r19943, r19948, MPFR_RNDN);
        mpfr_mul(r19950, r19938, r19949, MPFR_RNDN);
        mpfr_sqrt(r19951, r19950, MPFR_RNDN);
        return mpfr_get_d(r19951, MPFR_RNDN);
}

static mpfr_t r19952, r19953, r19954, r19955, r19956, r19957, r19958, r19959, r19960, r19961, r19962, r19963, r19964, r19965, r19966, r19967, r19968, r19969, r19970, r19971, r19972;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r19952, "1/2", 10, MPFR_RNDN);
        mpfr_init(r19953);
        mpfr_init(r19954);
        mpfr_init(r19955);
        mpfr_init(r19956);
        mpfr_init(r19957);
        mpfr_init(r19958);
        mpfr_init(r19959);
        mpfr_init(r19960);
        mpfr_init(r19961);
        mpfr_init(r19962);
        mpfr_init(r19963);
        mpfr_init(r19964);
        mpfr_init(r19965);
        mpfr_init(r19966);
        mpfr_init_set_str(r19967, "2", 10, MPFR_RNDN);
        mpfr_init(r19968);
        mpfr_init(r19969);
        mpfr_init(r19970);
        mpfr_init(r19971);
        mpfr_init(r19972);
}

double f_dm(double a, double b, double c) {
        ;
        mpfr_mul(r19953, r19952, r19952, MPFR_RNDN);
        mpfr_set_d(r19954, c, MPFR_RNDN);
        mpfr_set_d(r19955, a, MPFR_RNDN);
        mpfr_sub(r19956, r19954, r19955, MPFR_RNDN);
        mpfr_set_d(r19957, b, MPFR_RNDN);
        mpfr_add(r19958, r19956, r19957, MPFR_RNDN);
        mpfr_mul(r19959, r19953, r19958, MPFR_RNDN);
        mpfr_sub(r19960, r19955, r19954, MPFR_RNDN);
        mpfr_add(r19961, r19957, r19960, MPFR_RNDN);
        mpfr_sub(r19962, r19955, r19957, MPFR_RNDN);
        mpfr_add(r19963, r19962, r19954, MPFR_RNDN);
        mpfr_mul(r19964, r19961, r19963, MPFR_RNDN);
        mpfr_add(r19965, r19955, r19954, MPFR_RNDN);
        mpfr_add(r19966, r19957, r19965, MPFR_RNDN);
        ;
        mpfr_div(r19968, r19967, r19952, MPFR_RNDN);
        mpfr_div(r19969, r19966, r19968, MPFR_RNDN);
        mpfr_mul(r19970, r19964, r19969, MPFR_RNDN);
        mpfr_mul(r19971, r19959, r19970, MPFR_RNDN);
        mpfr_sqrt(r19972, r19971, MPFR_RNDN);
        return mpfr_get_d(r19972, MPFR_RNDN);
}

