#include #include #include #include #include char *name = "Bearing on a great circle"; double f_if(float lambda1, float lambda2, float phi1, float phi2) { float r24084 = lambda1; float r24085 = lambda2; float r24086 = r24084 - r24085; float r24087 = sin(r24086); float r24088 = phi2; float r24089 = cos(r24088); float r24090 = r24087 * r24089; float r24091 = phi1; float r24092 = cos(r24091); float r24093 = sin(r24088); float r24094 = r24092 * r24093; float r24095 = sin(r24091); float r24096 = r24095 * r24089; float r24097 = cos(r24086); float r24098 = r24096 * r24097; float r24099 = r24094 - r24098; float r24100 = atan2(r24090, r24099); return r24100; } double f_id(double lambda1, double lambda2, double phi1, double phi2) { double r24101 = lambda1; double r24102 = lambda2; double r24103 = r24101 - r24102; double r24104 = sin(r24103); double r24105 = phi2; double r24106 = cos(r24105); double r24107 = r24104 * r24106; double r24108 = phi1; double r24109 = cos(r24108); double r24110 = sin(r24105); double r24111 = r24109 * r24110; double r24112 = sin(r24108); double r24113 = r24112 * r24106; double r24114 = cos(r24103); double r24115 = r24113 * r24114; double r24116 = r24111 - r24115; double r24117 = atan2(r24107, r24116); return r24117; } double f_of(float lambda1, float lambda2, float phi1, float phi2) { float r24118 = lambda1; float r24119 = sin(r24118); float r24120 = lambda2; float r24121 = cos(r24120); float r24122 = r24119 * r24121; float r24123 = cos(r24118); float r24124 = cbrt(r24123); float r24125 = r24124 * r24124; float r24126 = -r24120; float r24127 = sin(r24126); float r24128 = r24124 * r24127; float r24129 = r24125 * r24128; float r24130 = r24122 + r24129; float r24131 = phi2; float r24132 = cos(r24131); float r24133 = r24130 * r24132; float r24134 = phi1; float r24135 = cos(r24134); float r24136 = sin(r24131); float r24137 = r24135 * r24136; float r24138 = sin(r24134); float r24139 = r24138 * r24132; float r24140 = r24123 * r24121; float r24141 = r24119 * r24127; float r24142 = r24140 - r24141; float r24143 = r24139 * r24142; float r24144 = r24137 - r24143; float r24145 = atan2(r24133, r24144); return r24145; } double f_od(double lambda1, double lambda2, double phi1, double phi2) { double r24146 = lambda1; double r24147 = sin(r24146); double r24148 = lambda2; double r24149 = cos(r24148); double r24150 = r24147 * r24149; double r24151 = cos(r24146); double r24152 = cbrt(r24151); double r24153 = r24152 * r24152; double r24154 = -r24148; double r24155 = sin(r24154); double r24156 = r24152 * r24155; double r24157 = r24153 * r24156; double r24158 = r24150 + r24157; double r24159 = phi2; double r24160 = cos(r24159); double r24161 = r24158 * r24160; double r24162 = phi1; double r24163 = cos(r24162); double r24164 = sin(r24159); double r24165 = r24163 * r24164; double r24166 = sin(r24162); double r24167 = r24166 * r24160; double r24168 = r24151 * r24149; double r24169 = r24147 * r24155; double r24170 = r24168 - r24169; double r24171 = r24167 * r24170; double r24172 = r24165 - r24171; double r24173 = atan2(r24161, r24172); return r24173; } 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 r24174, r24175, r24176, r24177, r24178, r24179, r24180, r24181, r24182, r24183, r24184, r24185, r24186, r24187, r24188, r24189, r24190; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r24174); mpfr_init(r24175); mpfr_init(r24176); mpfr_init(r24177); mpfr_init(r24178); mpfr_init(r24179); mpfr_init(r24180); mpfr_init(r24181); mpfr_init(r24182); mpfr_init(r24183); mpfr_init(r24184); mpfr_init(r24185); mpfr_init(r24186); mpfr_init(r24187); mpfr_init(r24188); mpfr_init(r24189); mpfr_init(r24190); } double f_im(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r24174, lambda1, MPFR_RNDN); mpfr_set_d(r24175, lambda2, MPFR_RNDN); mpfr_sub(r24176, r24174, r24175, MPFR_RNDN); mpfr_sin(r24177, r24176, MPFR_RNDN); mpfr_set_d(r24178, phi2, MPFR_RNDN); mpfr_cos(r24179, r24178, MPFR_RNDN); mpfr_mul(r24180, r24177, r24179, MPFR_RNDN); mpfr_set_d(r24181, phi1, MPFR_RNDN); mpfr_cos(r24182, r24181, MPFR_RNDN); mpfr_sin(r24183, r24178, MPFR_RNDN); mpfr_mul(r24184, r24182, r24183, MPFR_RNDN); mpfr_sin(r24185, r24181, MPFR_RNDN); mpfr_mul(r24186, r24185, r24179, MPFR_RNDN); mpfr_cos(r24187, r24176, MPFR_RNDN); mpfr_mul(r24188, r24186, r24187, MPFR_RNDN); mpfr_sub(r24189, r24184, r24188, MPFR_RNDN); mpfr_atan2(r24190, r24180, r24189, MPFR_RNDN); return mpfr_get_d(r24190, MPFR_RNDN); } static mpfr_t r24191, r24192, r24193, r24194, r24195, r24196, r24197, r24198, r24199, r24200, r24201, r24202, r24203, r24204, r24205, r24206, r24207, r24208, r24209, r24210, r24211, r24212, r24213, r24214, r24215, r24216, r24217, r24218; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r24191); mpfr_init(r24192); mpfr_init(r24193); mpfr_init(r24194); mpfr_init(r24195); mpfr_init(r24196); mpfr_init(r24197); mpfr_init(r24198); mpfr_init(r24199); mpfr_init(r24200); mpfr_init(r24201); mpfr_init(r24202); mpfr_init(r24203); mpfr_init(r24204); mpfr_init(r24205); mpfr_init(r24206); mpfr_init(r24207); mpfr_init(r24208); mpfr_init(r24209); mpfr_init(r24210); mpfr_init(r24211); mpfr_init(r24212); mpfr_init(r24213); mpfr_init(r24214); mpfr_init(r24215); mpfr_init(r24216); mpfr_init(r24217); mpfr_init(r24218); } double f_fm(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r24191, lambda1, MPFR_RNDN); mpfr_sin(r24192, r24191, MPFR_RNDN); mpfr_set_d(r24193, lambda2, MPFR_RNDN); mpfr_cos(r24194, r24193, MPFR_RNDN); mpfr_mul(r24195, r24192, r24194, MPFR_RNDN); mpfr_cos(r24196, r24191, MPFR_RNDN); mpfr_cbrt(r24197, r24196, MPFR_RNDN); mpfr_mul(r24198, r24197, r24197, MPFR_RNDN); mpfr_neg(r24199, r24193, MPFR_RNDN); mpfr_sin(r24200, r24199, MPFR_RNDN); mpfr_mul(r24201, r24197, r24200, MPFR_RNDN); mpfr_mul(r24202, r24198, r24201, MPFR_RNDN); mpfr_add(r24203, r24195, r24202, MPFR_RNDN); mpfr_set_d(r24204, phi2, MPFR_RNDN); mpfr_cos(r24205, r24204, MPFR_RNDN); mpfr_mul(r24206, r24203, r24205, MPFR_RNDN); mpfr_set_d(r24207, phi1, MPFR_RNDN); mpfr_cos(r24208, r24207, MPFR_RNDN); mpfr_sin(r24209, r24204, MPFR_RNDN); mpfr_mul(r24210, r24208, r24209, MPFR_RNDN); mpfr_sin(r24211, r24207, MPFR_RNDN); mpfr_mul(r24212, r24211, r24205, MPFR_RNDN); mpfr_mul(r24213, r24196, r24194, MPFR_RNDN); mpfr_mul(r24214, r24192, r24200, MPFR_RNDN); mpfr_sub(r24215, r24213, r24214, MPFR_RNDN); mpfr_mul(r24216, r24212, r24215, MPFR_RNDN); mpfr_sub(r24217, r24210, r24216, MPFR_RNDN); mpfr_atan2(r24218, r24206, r24217, MPFR_RNDN); return mpfr_get_d(r24218, MPFR_RNDN); } static mpfr_t r24219, r24220, r24221, r24222, r24223, r24224, r24225, r24226, r24227, r24228, r24229, r24230, r24231, r24232, r24233, r24234, r24235, r24236, r24237, r24238, r24239, r24240, r24241, r24242, r24243, r24244, r24245, r24246; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r24219); mpfr_init(r24220); mpfr_init(r24221); mpfr_init(r24222); mpfr_init(r24223); mpfr_init(r24224); mpfr_init(r24225); mpfr_init(r24226); mpfr_init(r24227); mpfr_init(r24228); mpfr_init(r24229); mpfr_init(r24230); mpfr_init(r24231); mpfr_init(r24232); mpfr_init(r24233); mpfr_init(r24234); mpfr_init(r24235); mpfr_init(r24236); mpfr_init(r24237); mpfr_init(r24238); mpfr_init(r24239); mpfr_init(r24240); mpfr_init(r24241); mpfr_init(r24242); mpfr_init(r24243); mpfr_init(r24244); mpfr_init(r24245); mpfr_init(r24246); } double f_dm(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r24219, lambda1, MPFR_RNDN); mpfr_sin(r24220, r24219, MPFR_RNDN); mpfr_set_d(r24221, lambda2, MPFR_RNDN); mpfr_cos(r24222, r24221, MPFR_RNDN); mpfr_mul(r24223, r24220, r24222, MPFR_RNDN); mpfr_cos(r24224, r24219, MPFR_RNDN); mpfr_cbrt(r24225, r24224, MPFR_RNDN); mpfr_mul(r24226, r24225, r24225, MPFR_RNDN); mpfr_neg(r24227, r24221, MPFR_RNDN); mpfr_sin(r24228, r24227, MPFR_RNDN); mpfr_mul(r24229, r24225, r24228, MPFR_RNDN); mpfr_mul(r24230, r24226, r24229, MPFR_RNDN); mpfr_add(r24231, r24223, r24230, MPFR_RNDN); mpfr_set_d(r24232, phi2, MPFR_RNDN); mpfr_cos(r24233, r24232, MPFR_RNDN); mpfr_mul(r24234, r24231, r24233, MPFR_RNDN); mpfr_set_d(r24235, phi1, MPFR_RNDN); mpfr_cos(r24236, r24235, MPFR_RNDN); mpfr_sin(r24237, r24232, MPFR_RNDN); mpfr_mul(r24238, r24236, r24237, MPFR_RNDN); mpfr_sin(r24239, r24235, MPFR_RNDN); mpfr_mul(r24240, r24239, r24233, MPFR_RNDN); mpfr_mul(r24241, r24224, r24222, MPFR_RNDN); mpfr_mul(r24242, r24220, r24228, MPFR_RNDN); mpfr_sub(r24243, r24241, r24242, MPFR_RNDN); mpfr_mul(r24244, r24240, r24243, MPFR_RNDN); mpfr_sub(r24245, r24238, r24244, MPFR_RNDN); mpfr_atan2(r24246, r24234, r24245, MPFR_RNDN); return mpfr_get_d(r24246, MPFR_RNDN); }