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

char *name = "Equirectangular approximation to distance on a great circle";

double f_if(float R, float lambda1, float lambda2, float phi1, float phi2) {
        float r22861 = R;
        float r22862 = lambda1;
        float r22863 = lambda2;
        float r22864 = r22862 - r22863;
        float r22865 = phi1;
        float r22866 = phi2;
        float r22867 = r22865 + r22866;
        float r22868 = 2;
        float r22869 = r22867 / r22868;
        float r22870 = cos(r22869);
        float r22871 = r22864 * r22870;
        float r22872 = r22871 * r22871;
        float r22873 = r22865 - r22866;
        float r22874 = r22873 * r22873;
        float r22875 = r22872 + r22874;
        float r22876 = sqrt(r22875);
        float r22877 = r22861 * r22876;
        return r22877;
}

double f_id(double R, double lambda1, double lambda2, double phi1, double phi2) {
        double r22878 = R;
        double r22879 = lambda1;
        double r22880 = lambda2;
        double r22881 = r22879 - r22880;
        double r22882 = phi1;
        double r22883 = phi2;
        double r22884 = r22882 + r22883;
        double r22885 = 2;
        double r22886 = r22884 / r22885;
        double r22887 = cos(r22886);
        double r22888 = r22881 * r22887;
        double r22889 = r22888 * r22888;
        double r22890 = r22882 - r22883;
        double r22891 = r22890 * r22890;
        double r22892 = r22889 + r22891;
        double r22893 = sqrt(r22892);
        double r22894 = r22878 * r22893;
        return r22894;
}


double f_of(float R, float lambda1, float lambda2, float phi1, float phi2) {
        float r22895 = phi1;
        float r22896 = phi2;
        float r22897 = r22895 + r22896;
        float r22898 = 2;
        float r22899 = r22897 / r22898;
        float r22900 = cos(r22899);
        float r22901 = log1p(r22900);
        float r22902 = expm1(r22901);
        float r22903 = lambda1;
        float r22904 = lambda2;
        float r22905 = r22903 - r22904;
        float r22906 = r22902 * r22905;
        float r22907 = r22895 - r22896;
        float r22908 = hypot(r22906, r22907);
        float r22909 = R;
        float r22910 = r22908 * r22909;
        return r22910;
}

double f_od(double R, double lambda1, double lambda2, double phi1, double phi2) {
        double r22911 = phi1;
        double r22912 = phi2;
        double r22913 = r22911 + r22912;
        double r22914 = 2;
        double r22915 = r22913 / r22914;
        double r22916 = cos(r22915);
        double r22917 = log1p(r22916);
        double r22918 = expm1(r22917);
        double r22919 = lambda1;
        double r22920 = lambda2;
        double r22921 = r22919 - r22920;
        double r22922 = r22918 * r22921;
        double r22923 = r22911 - r22912;
        double r22924 = hypot(r22922, r22923);
        double r22925 = R;
        double r22926 = r22924 * r22925;
        return r22926;
}

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 r22927, r22928, r22929, r22930, r22931, r22932, r22933, r22934, r22935, r22936, r22937, r22938, r22939, r22940, r22941, r22942, r22943;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init(r22927);
        mpfr_init(r22928);
        mpfr_init(r22929);
        mpfr_init(r22930);
        mpfr_init(r22931);
        mpfr_init(r22932);
        mpfr_init(r22933);
        mpfr_init_set_str(r22934, "2", 10, MPFR_RNDN);
        mpfr_init(r22935);
        mpfr_init(r22936);
        mpfr_init(r22937);
        mpfr_init(r22938);
        mpfr_init(r22939);
        mpfr_init(r22940);
        mpfr_init(r22941);
        mpfr_init(r22942);
        mpfr_init(r22943);
}

double f_im(double R, double lambda1, double lambda2, double phi1, double phi2) {
        mpfr_set_d(r22927, R, MPFR_RNDN);
        mpfr_set_d(r22928, lambda1, MPFR_RNDN);
        mpfr_set_d(r22929, lambda2, MPFR_RNDN);
        mpfr_sub(r22930, r22928, r22929, MPFR_RNDN);
        mpfr_set_d(r22931, phi1, MPFR_RNDN);
        mpfr_set_d(r22932, phi2, MPFR_RNDN);
        mpfr_add(r22933, r22931, r22932, MPFR_RNDN);
        ;
        mpfr_div(r22935, r22933, r22934, MPFR_RNDN);
        mpfr_cos(r22936, r22935, MPFR_RNDN);
        mpfr_mul(r22937, r22930, r22936, MPFR_RNDN);
        mpfr_mul(r22938, r22937, r22937, MPFR_RNDN);
        mpfr_sub(r22939, r22931, r22932, MPFR_RNDN);
        mpfr_mul(r22940, r22939, r22939, MPFR_RNDN);
        mpfr_add(r22941, r22938, r22940, MPFR_RNDN);
        mpfr_sqrt(r22942, r22941, MPFR_RNDN);
        mpfr_mul(r22943, r22927, r22942, MPFR_RNDN);
        return mpfr_get_d(r22943, MPFR_RNDN);
}

static mpfr_t r22944, r22945, r22946, r22947, r22948, r22949, r22950, r22951, r22952, r22953, r22954, r22955, r22956, r22957, r22958, r22959;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r22944);
        mpfr_init(r22945);
        mpfr_init(r22946);
        mpfr_init_set_str(r22947, "2", 10, MPFR_RNDN);
        mpfr_init(r22948);
        mpfr_init(r22949);
        mpfr_init(r22950);
        mpfr_init(r22951);
        mpfr_init(r22952);
        mpfr_init(r22953);
        mpfr_init(r22954);
        mpfr_init(r22955);
        mpfr_init(r22956);
        mpfr_init(r22957);
        mpfr_init(r22958);
        mpfr_init(r22959);
}

double f_fm(double R, double lambda1, double lambda2, double phi1, double phi2) {
        mpfr_set_d(r22944, phi1, MPFR_RNDN);
        mpfr_set_d(r22945, phi2, MPFR_RNDN);
        mpfr_add(r22946, r22944, r22945, MPFR_RNDN);
        ;
        mpfr_div(r22948, r22946, r22947, MPFR_RNDN);
        mpfr_cos(r22949, r22948, MPFR_RNDN);
        mpfr_log1p(r22950, r22949, MPFR_RNDN);
        mpfr_expm1(r22951, r22950, MPFR_RNDN);
        mpfr_set_d(r22952, lambda1, MPFR_RNDN);
        mpfr_set_d(r22953, lambda2, MPFR_RNDN);
        mpfr_sub(r22954, r22952, r22953, MPFR_RNDN);
        mpfr_mul(r22955, r22951, r22954, MPFR_RNDN);
        mpfr_sub(r22956, r22944, r22945, MPFR_RNDN);
        mpfr_hypot(r22957, r22955, r22956, MPFR_RNDN);
        mpfr_set_d(r22958, R, MPFR_RNDN);
        mpfr_mul(r22959, r22957, r22958, MPFR_RNDN);
        return mpfr_get_d(r22959, MPFR_RNDN);
}

static mpfr_t r22960, r22961, r22962, r22963, r22964, r22965, r22966, r22967, r22968, r22969, r22970, r22971, r22972, r22973, r22974, r22975;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r22960);
        mpfr_init(r22961);
        mpfr_init(r22962);
        mpfr_init_set_str(r22963, "2", 10, MPFR_RNDN);
        mpfr_init(r22964);
        mpfr_init(r22965);
        mpfr_init(r22966);
        mpfr_init(r22967);
        mpfr_init(r22968);
        mpfr_init(r22969);
        mpfr_init(r22970);
        mpfr_init(r22971);
        mpfr_init(r22972);
        mpfr_init(r22973);
        mpfr_init(r22974);
        mpfr_init(r22975);
}

double f_dm(double R, double lambda1, double lambda2, double phi1, double phi2) {
        mpfr_set_d(r22960, phi1, MPFR_RNDN);
        mpfr_set_d(r22961, phi2, MPFR_RNDN);
        mpfr_add(r22962, r22960, r22961, MPFR_RNDN);
        ;
        mpfr_div(r22964, r22962, r22963, MPFR_RNDN);
        mpfr_cos(r22965, r22964, MPFR_RNDN);
        mpfr_log1p(r22966, r22965, MPFR_RNDN);
        mpfr_expm1(r22967, r22966, MPFR_RNDN);
        mpfr_set_d(r22968, lambda1, MPFR_RNDN);
        mpfr_set_d(r22969, lambda2, MPFR_RNDN);
        mpfr_sub(r22970, r22968, r22969, MPFR_RNDN);
        mpfr_mul(r22971, r22967, r22970, MPFR_RNDN);
        mpfr_sub(r22972, r22960, r22961, MPFR_RNDN);
        mpfr_hypot(r22973, r22971, r22972, MPFR_RNDN);
        mpfr_set_d(r22974, R, MPFR_RNDN);
        mpfr_mul(r22975, r22973, r22974, MPFR_RNDN);
        return mpfr_get_d(r22975, MPFR_RNDN);
}