#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 r22816 = R;
        float r22817 = lambda1;
        float r22818 = lambda2;
        float r22819 = r22817 - r22818;
        float r22820 = phi1;
        float r22821 = phi2;
        float r22822 = r22820 + r22821;
        float r22823 = 2;
        float r22824 = r22822 / r22823;
        float r22825 = cos(r22824);
        float r22826 = r22819 * r22825;
        float r22827 = r22826 * r22826;
        float r22828 = r22820 - r22821;
        float r22829 = r22828 * r22828;
        float r22830 = r22827 + r22829;
        float r22831 = sqrt(r22830);
        float r22832 = r22816 * r22831;
        return r22832;
}

double f_id(double R, double lambda1, double lambda2, double phi1, double phi2) {
        double r22833 = R;
        double r22834 = lambda1;
        double r22835 = lambda2;
        double r22836 = r22834 - r22835;
        double r22837 = phi1;
        double r22838 = phi2;
        double r22839 = r22837 + r22838;
        double r22840 = 2;
        double r22841 = r22839 / r22840;
        double r22842 = cos(r22841);
        double r22843 = r22836 * r22842;
        double r22844 = r22843 * r22843;
        double r22845 = r22837 - r22838;
        double r22846 = r22845 * r22845;
        double r22847 = r22844 + r22846;
        double r22848 = sqrt(r22847);
        double r22849 = r22833 * r22848;
        return r22849;
}


double f_of(float R, float lambda1, float lambda2, float phi1, float phi2) {
        float r22850 = phi1;
        float r22851 = phi2;
        float r22852 = r22850 + r22851;
        float r22853 = 2;
        float r22854 = r22852 / r22853;
        float r22855 = cos(r22854);
        float r22856 = log1p(r22855);
        float r22857 = expm1(r22856);
        float r22858 = lambda1;
        float r22859 = lambda2;
        float r22860 = r22858 - r22859;
        float r22861 = r22857 * r22860;
        float r22862 = r22850 - r22851;
        float r22863 = hypot(r22861, r22862);
        float r22864 = R;
        float r22865 = r22863 * r22864;
        return r22865;
}

double f_od(double R, double lambda1, double lambda2, double phi1, double phi2) {
        double r22866 = phi1;
        double r22867 = phi2;
        double r22868 = r22866 + r22867;
        double r22869 = 2;
        double r22870 = r22868 / r22869;
        double r22871 = cos(r22870);
        double r22872 = log1p(r22871);
        double r22873 = expm1(r22872);
        double r22874 = lambda1;
        double r22875 = lambda2;
        double r22876 = r22874 - r22875;
        double r22877 = r22873 * r22876;
        double r22878 = r22866 - r22867;
        double r22879 = hypot(r22877, r22878);
        double r22880 = R;
        double r22881 = r22879 * r22880;
        return r22881;
}

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 r22882, r22883, r22884, r22885, r22886, r22887, r22888, r22889, r22890, r22891, r22892, r22893, r22894, r22895, r22896, r22897, r22898;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init(r22882);
        mpfr_init(r22883);
        mpfr_init(r22884);
        mpfr_init(r22885);
        mpfr_init(r22886);
        mpfr_init(r22887);
        mpfr_init(r22888);
        mpfr_init_set_str(r22889, "2", 10, MPFR_RNDN);
        mpfr_init(r22890);
        mpfr_init(r22891);
        mpfr_init(r22892);
        mpfr_init(r22893);
        mpfr_init(r22894);
        mpfr_init(r22895);
        mpfr_init(r22896);
        mpfr_init(r22897);
        mpfr_init(r22898);
}

double f_im(double R, double lambda1, double lambda2, double phi1, double phi2) {
        mpfr_set_d(r22882, R, MPFR_RNDN);
        mpfr_set_d(r22883, lambda1, MPFR_RNDN);
        mpfr_set_d(r22884, lambda2, MPFR_RNDN);
        mpfr_sub(r22885, r22883, r22884, MPFR_RNDN);
        mpfr_set_d(r22886, phi1, MPFR_RNDN);
        mpfr_set_d(r22887, phi2, MPFR_RNDN);
        mpfr_add(r22888, r22886, r22887, MPFR_RNDN);
        ;
        mpfr_div(r22890, r22888, r22889, MPFR_RNDN);
        mpfr_cos(r22891, r22890, MPFR_RNDN);
        mpfr_mul(r22892, r22885, r22891, MPFR_RNDN);
        mpfr_mul(r22893, r22892, r22892, MPFR_RNDN);
        mpfr_sub(r22894, r22886, r22887, MPFR_RNDN);
        mpfr_mul(r22895, r22894, r22894, MPFR_RNDN);
        mpfr_add(r22896, r22893, r22895, MPFR_RNDN);
        mpfr_sqrt(r22897, r22896, MPFR_RNDN);
        mpfr_mul(r22898, r22882, r22897, MPFR_RNDN);
        return mpfr_get_d(r22898, MPFR_RNDN);
}

static mpfr_t r22899, r22900, r22901, r22902, r22903, r22904, r22905, r22906, r22907, r22908, r22909, r22910, r22911, r22912, r22913, r22914;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r22899);
        mpfr_init(r22900);
        mpfr_init(r22901);
        mpfr_init_set_str(r22902, "2", 10, MPFR_RNDN);
        mpfr_init(r22903);
        mpfr_init(r22904);
        mpfr_init(r22905);
        mpfr_init(r22906);
        mpfr_init(r22907);
        mpfr_init(r22908);
        mpfr_init(r22909);
        mpfr_init(r22910);
        mpfr_init(r22911);
        mpfr_init(r22912);
        mpfr_init(r22913);
        mpfr_init(r22914);
}

double f_fm(double R, double lambda1, double lambda2, double phi1, double phi2) {
        mpfr_set_d(r22899, phi1, MPFR_RNDN);
        mpfr_set_d(r22900, phi2, MPFR_RNDN);
        mpfr_add(r22901, r22899, r22900, MPFR_RNDN);
        ;
        mpfr_div(r22903, r22901, r22902, MPFR_RNDN);
        mpfr_cos(r22904, r22903, MPFR_RNDN);
        mpfr_log1p(r22905, r22904, MPFR_RNDN);
        mpfr_expm1(r22906, r22905, MPFR_RNDN);
        mpfr_set_d(r22907, lambda1, MPFR_RNDN);
        mpfr_set_d(r22908, lambda2, MPFR_RNDN);
        mpfr_sub(r22909, r22907, r22908, MPFR_RNDN);
        mpfr_mul(r22910, r22906, r22909, MPFR_RNDN);
        mpfr_sub(r22911, r22899, r22900, MPFR_RNDN);
        mpfr_hypot(r22912, r22910, r22911, MPFR_RNDN);
        mpfr_set_d(r22913, R, MPFR_RNDN);
        mpfr_mul(r22914, r22912, r22913, MPFR_RNDN);
        return mpfr_get_d(r22914, MPFR_RNDN);
}

static mpfr_t r22915, r22916, r22917, r22918, r22919, r22920, r22921, r22922, r22923, r22924, r22925, r22926, r22927, r22928, r22929, r22930;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r22915);
        mpfr_init(r22916);
        mpfr_init(r22917);
        mpfr_init_set_str(r22918, "2", 10, MPFR_RNDN);
        mpfr_init(r22919);
        mpfr_init(r22920);
        mpfr_init(r22921);
        mpfr_init(r22922);
        mpfr_init(r22923);
        mpfr_init(r22924);
        mpfr_init(r22925);
        mpfr_init(r22926);
        mpfr_init(r22927);
        mpfr_init(r22928);
        mpfr_init(r22929);
        mpfr_init(r22930);
}

double f_dm(double R, double lambda1, double lambda2, double phi1, double phi2) {
        mpfr_set_d(r22915, phi1, MPFR_RNDN);
        mpfr_set_d(r22916, phi2, MPFR_RNDN);
        mpfr_add(r22917, r22915, r22916, MPFR_RNDN);
        ;
        mpfr_div(r22919, r22917, r22918, MPFR_RNDN);
        mpfr_cos(r22920, r22919, MPFR_RNDN);
        mpfr_log1p(r22921, r22920, MPFR_RNDN);
        mpfr_expm1(r22922, r22921, MPFR_RNDN);
        mpfr_set_d(r22923, lambda1, MPFR_RNDN);
        mpfr_set_d(r22924, lambda2, MPFR_RNDN);
        mpfr_sub(r22925, r22923, r22924, MPFR_RNDN);
        mpfr_mul(r22926, r22922, r22925, MPFR_RNDN);
        mpfr_sub(r22927, r22915, r22916, MPFR_RNDN);
        mpfr_hypot(r22928, r22926, r22927, MPFR_RNDN);
        mpfr_set_d(r22929, R, MPFR_RNDN);
        mpfr_mul(r22930, r22928, r22929, MPFR_RNDN);
        return mpfr_get_d(r22930, MPFR_RNDN);
}