#include #include #include #include #include char *name = "Bearing on a great circle"; double f_if(float lambda1, float lambda2, float phi1, float phi2) { float r26268 = lambda1; float r26269 = lambda2; float r26270 = r26268 - r26269; float r26271 = sin(r26270); float r26272 = phi2; float r26273 = cos(r26272); float r26274 = r26271 * r26273; float r26275 = phi1; float r26276 = cos(r26275); float r26277 = sin(r26272); float r26278 = r26276 * r26277; float r26279 = sin(r26275); float r26280 = r26279 * r26273; float r26281 = cos(r26270); float r26282 = r26280 * r26281; float r26283 = r26278 - r26282; float r26284 = atan2(r26274, r26283); return r26284; } double f_id(double lambda1, double lambda2, double phi1, double phi2) { double r26285 = lambda1; double r26286 = lambda2; double r26287 = r26285 - r26286; double r26288 = sin(r26287); double r26289 = phi2; double r26290 = cos(r26289); double r26291 = r26288 * r26290; double r26292 = phi1; double r26293 = cos(r26292); double r26294 = sin(r26289); double r26295 = r26293 * r26294; double r26296 = sin(r26292); double r26297 = r26296 * r26290; double r26298 = cos(r26287); double r26299 = r26297 * r26298; double r26300 = r26295 - r26299; double r26301 = atan2(r26291, r26300); return r26301; } double f_of(float lambda1, float lambda2, float phi1, float phi2) { float r26302 = lambda1; float r26303 = sin(r26302); float r26304 = lambda2; float r26305 = cos(r26304); float r26306 = r26303 * r26305; float r26307 = cos(r26302); float r26308 = -r26304; float r26309 = sin(r26308); float r26310 = r26307 * r26309; float r26311 = r26306 + r26310; float r26312 = phi2; float r26313 = cos(r26312); float r26314 = r26311 * r26313; float r26315 = phi1; float r26316 = cos(r26315); float r26317 = sin(r26312); float r26318 = r26316 * r26317; float r26319 = sin(r26315); float r26320 = r26319 * r26313; float r26321 = r26307 * r26305; float r26322 = r26320 * r26321; float r26323 = r26319 * r26303; float r26324 = sin(r26304); float r26325 = r26313 * r26324; float r26326 = r26323 * r26325; float r26327 = r26322 + r26326; float r26328 = r26318 - r26327; float r26329 = atan2(r26314, r26328); return r26329; } double f_od(double lambda1, double lambda2, double phi1, double phi2) { double r26330 = lambda1; double r26331 = sin(r26330); double r26332 = lambda2; double r26333 = cos(r26332); double r26334 = r26331 * r26333; double r26335 = cos(r26330); double r26336 = -r26332; double r26337 = sin(r26336); double r26338 = r26335 * r26337; double r26339 = r26334 + r26338; double r26340 = phi2; double r26341 = cos(r26340); double r26342 = r26339 * r26341; double r26343 = phi1; double r26344 = cos(r26343); double r26345 = sin(r26340); double r26346 = r26344 * r26345; double r26347 = sin(r26343); double r26348 = r26347 * r26341; double r26349 = r26335 * r26333; double r26350 = r26348 * r26349; double r26351 = r26347 * r26331; double r26352 = sin(r26332); double r26353 = r26341 * r26352; double r26354 = r26351 * r26353; double r26355 = r26350 + r26354; double r26356 = r26346 - r26355; double r26357 = atan2(r26342, r26356); return r26357; } 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 r26358, r26359, r26360, r26361, r26362, r26363, r26364, r26365, r26366, r26367, r26368, r26369, r26370, r26371, r26372, r26373, r26374; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r26358); mpfr_init(r26359); mpfr_init(r26360); mpfr_init(r26361); mpfr_init(r26362); mpfr_init(r26363); mpfr_init(r26364); mpfr_init(r26365); mpfr_init(r26366); mpfr_init(r26367); mpfr_init(r26368); mpfr_init(r26369); mpfr_init(r26370); mpfr_init(r26371); mpfr_init(r26372); mpfr_init(r26373); mpfr_init(r26374); } double f_im(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r26358, lambda1, MPFR_RNDN); mpfr_set_d(r26359, lambda2, MPFR_RNDN); mpfr_sub(r26360, r26358, r26359, MPFR_RNDN); mpfr_sin(r26361, r26360, MPFR_RNDN); mpfr_set_d(r26362, phi2, MPFR_RNDN); mpfr_cos(r26363, r26362, MPFR_RNDN); mpfr_mul(r26364, r26361, r26363, MPFR_RNDN); mpfr_set_d(r26365, phi1, MPFR_RNDN); mpfr_cos(r26366, r26365, MPFR_RNDN); mpfr_sin(r26367, r26362, MPFR_RNDN); mpfr_mul(r26368, r26366, r26367, MPFR_RNDN); mpfr_sin(r26369, r26365, MPFR_RNDN); mpfr_mul(r26370, r26369, r26363, MPFR_RNDN); mpfr_cos(r26371, r26360, MPFR_RNDN); mpfr_mul(r26372, r26370, r26371, MPFR_RNDN); mpfr_sub(r26373, r26368, r26372, MPFR_RNDN); mpfr_atan2(r26374, r26364, r26373, MPFR_RNDN); return mpfr_get_d(r26374, MPFR_RNDN); } static mpfr_t r26375, r26376, r26377, r26378, r26379, r26380, r26381, r26382, r26383, r26384, r26385, r26386, r26387, r26388, r26389, r26390, r26391, r26392, r26393, r26394, r26395, r26396, r26397, r26398, r26399, r26400, r26401, r26402; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r26375); mpfr_init(r26376); mpfr_init(r26377); mpfr_init(r26378); mpfr_init(r26379); mpfr_init(r26380); mpfr_init(r26381); mpfr_init(r26382); mpfr_init(r26383); mpfr_init(r26384); mpfr_init(r26385); mpfr_init(r26386); mpfr_init(r26387); mpfr_init(r26388); mpfr_init(r26389); mpfr_init(r26390); mpfr_init(r26391); mpfr_init(r26392); mpfr_init(r26393); mpfr_init(r26394); mpfr_init(r26395); mpfr_init(r26396); mpfr_init(r26397); mpfr_init(r26398); mpfr_init(r26399); mpfr_init(r26400); mpfr_init(r26401); mpfr_init(r26402); } double f_fm(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r26375, lambda1, MPFR_RNDN); mpfr_sin(r26376, r26375, MPFR_RNDN); mpfr_set_d(r26377, lambda2, MPFR_RNDN); mpfr_cos(r26378, r26377, MPFR_RNDN); mpfr_mul(r26379, r26376, r26378, MPFR_RNDN); mpfr_cos(r26380, r26375, MPFR_RNDN); mpfr_neg(r26381, r26377, MPFR_RNDN); mpfr_sin(r26382, r26381, MPFR_RNDN); mpfr_mul(r26383, r26380, r26382, MPFR_RNDN); mpfr_add(r26384, r26379, r26383, MPFR_RNDN); mpfr_set_d(r26385, phi2, MPFR_RNDN); mpfr_cos(r26386, r26385, MPFR_RNDN); mpfr_mul(r26387, r26384, r26386, MPFR_RNDN); mpfr_set_d(r26388, phi1, MPFR_RNDN); mpfr_cos(r26389, r26388, MPFR_RNDN); mpfr_sin(r26390, r26385, MPFR_RNDN); mpfr_mul(r26391, r26389, r26390, MPFR_RNDN); mpfr_sin(r26392, r26388, MPFR_RNDN); mpfr_mul(r26393, r26392, r26386, MPFR_RNDN); mpfr_mul(r26394, r26380, r26378, MPFR_RNDN); mpfr_mul(r26395, r26393, r26394, MPFR_RNDN); mpfr_mul(r26396, r26392, r26376, MPFR_RNDN); mpfr_sin(r26397, r26377, MPFR_RNDN); mpfr_mul(r26398, r26386, r26397, MPFR_RNDN); mpfr_mul(r26399, r26396, r26398, MPFR_RNDN); mpfr_add(r26400, r26395, r26399, MPFR_RNDN); mpfr_sub(r26401, r26391, r26400, MPFR_RNDN); mpfr_atan2(r26402, r26387, r26401, MPFR_RNDN); return mpfr_get_d(r26402, MPFR_RNDN); } static mpfr_t r26403, r26404, r26405, r26406, r26407, r26408, r26409, r26410, r26411, r26412, r26413, r26414, r26415, r26416, r26417, r26418, r26419, r26420, r26421, r26422, r26423, r26424, r26425, r26426, r26427, r26428, r26429, r26430; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r26403); mpfr_init(r26404); mpfr_init(r26405); mpfr_init(r26406); mpfr_init(r26407); mpfr_init(r26408); mpfr_init(r26409); mpfr_init(r26410); mpfr_init(r26411); mpfr_init(r26412); mpfr_init(r26413); mpfr_init(r26414); mpfr_init(r26415); mpfr_init(r26416); mpfr_init(r26417); mpfr_init(r26418); mpfr_init(r26419); mpfr_init(r26420); mpfr_init(r26421); mpfr_init(r26422); mpfr_init(r26423); mpfr_init(r26424); mpfr_init(r26425); mpfr_init(r26426); mpfr_init(r26427); mpfr_init(r26428); mpfr_init(r26429); mpfr_init(r26430); } double f_dm(double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r26403, lambda1, MPFR_RNDN); mpfr_sin(r26404, r26403, MPFR_RNDN); mpfr_set_d(r26405, lambda2, MPFR_RNDN); mpfr_cos(r26406, r26405, MPFR_RNDN); mpfr_mul(r26407, r26404, r26406, MPFR_RNDN); mpfr_cos(r26408, r26403, MPFR_RNDN); mpfr_neg(r26409, r26405, MPFR_RNDN); mpfr_sin(r26410, r26409, MPFR_RNDN); mpfr_mul(r26411, r26408, r26410, MPFR_RNDN); mpfr_add(r26412, r26407, r26411, MPFR_RNDN); mpfr_set_d(r26413, phi2, MPFR_RNDN); mpfr_cos(r26414, r26413, MPFR_RNDN); mpfr_mul(r26415, r26412, r26414, MPFR_RNDN); mpfr_set_d(r26416, phi1, MPFR_RNDN); mpfr_cos(r26417, r26416, MPFR_RNDN); mpfr_sin(r26418, r26413, MPFR_RNDN); mpfr_mul(r26419, r26417, r26418, MPFR_RNDN); mpfr_sin(r26420, r26416, MPFR_RNDN); mpfr_mul(r26421, r26420, r26414, MPFR_RNDN); mpfr_mul(r26422, r26408, r26406, MPFR_RNDN); mpfr_mul(r26423, r26421, r26422, MPFR_RNDN); mpfr_mul(r26424, r26420, r26404, MPFR_RNDN); mpfr_sin(r26425, r26405, MPFR_RNDN); mpfr_mul(r26426, r26414, r26425, MPFR_RNDN); mpfr_mul(r26427, r26424, r26426, MPFR_RNDN); mpfr_add(r26428, r26423, r26427, MPFR_RNDN); mpfr_sub(r26429, r26419, r26428, MPFR_RNDN); mpfr_atan2(r26430, r26415, r26429, MPFR_RNDN); return mpfr_get_d(r26430, MPFR_RNDN); }