#include #include #include #include #include char *name = "Spherical law of cosines"; double f_if(float R, float lambda1, float lambda2, float phi1, float phi2) { float r24213 = phi1; float r24214 = sin(r24213); float r24215 = phi2; float r24216 = sin(r24215); float r24217 = r24214 * r24216; float r24218 = cos(r24213); float r24219 = cos(r24215); float r24220 = r24218 * r24219; float r24221 = lambda1; float r24222 = lambda2; float r24223 = r24221 - r24222; float r24224 = cos(r24223); float r24225 = r24220 * r24224; float r24226 = r24217 + r24225; float r24227 = acos(r24226); float r24228 = R; float r24229 = r24227 * r24228; return r24229; } double f_id(double R, double lambda1, double lambda2, double phi1, double phi2) { double r24230 = phi1; double r24231 = sin(r24230); double r24232 = phi2; double r24233 = sin(r24232); double r24234 = r24231 * r24233; double r24235 = cos(r24230); double r24236 = cos(r24232); double r24237 = r24235 * r24236; double r24238 = lambda1; double r24239 = lambda2; double r24240 = r24238 - r24239; double r24241 = cos(r24240); double r24242 = r24237 * r24241; double r24243 = r24234 + r24242; double r24244 = acos(r24243); double r24245 = R; double r24246 = r24244 * r24245; return r24246; } double f_of(float R, float lambda1, float lambda2, float phi1, float phi2) { float r24247 = phi2; float r24248 = cos(r24247); float r24249 = phi1; float r24250 = cos(r24249); float r24251 = r24248 * r24250; float r24252 = lambda1; float r24253 = cos(r24252); float r24254 = lambda2; float r24255 = cos(r24254); float r24256 = r24253 * r24255; float r24257 = sin(r24252); float r24258 = -r24254; float r24259 = sin(r24258); float r24260 = r24257 * r24259; float r24261 = log1p(r24260); float r24262 = expm1(r24261); float r24263 = r24256 - r24262; float r24264 = sin(r24247); float r24265 = sin(r24249); float r24266 = r24264 * r24265; float r24267 = fma(r24251, r24263, r24266); float r24268 = acos(r24267); float r24269 = R; float r24270 = r24268 * r24269; return r24270; } double f_od(double R, double lambda1, double lambda2, double phi1, double phi2) { double r24271 = phi2; double r24272 = cos(r24271); double r24273 = phi1; double r24274 = cos(r24273); double r24275 = r24272 * r24274; double r24276 = lambda1; double r24277 = cos(r24276); double r24278 = lambda2; double r24279 = cos(r24278); double r24280 = r24277 * r24279; double r24281 = sin(r24276); double r24282 = -r24278; double r24283 = sin(r24282); double r24284 = r24281 * r24283; double r24285 = log1p(r24284); double r24286 = expm1(r24285); double r24287 = r24280 - r24286; double r24288 = sin(r24271); double r24289 = sin(r24273); double r24290 = r24288 * r24289; double r24291 = fma(r24275, r24287, r24290); double r24292 = acos(r24291); double r24293 = R; double r24294 = r24292 * r24293; return r24294; } 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 r24295, r24296, r24297, r24298, r24299, r24300, r24301, r24302, r24303, r24304, r24305, r24306, r24307, r24308, r24309, r24310, r24311; void setup_mpfr_f_im() { mpfr_set_default_prec(2128); mpfr_init(r24295); mpfr_init(r24296); mpfr_init(r24297); mpfr_init(r24298); mpfr_init(r24299); mpfr_init(r24300); mpfr_init(r24301); mpfr_init(r24302); mpfr_init(r24303); mpfr_init(r24304); mpfr_init(r24305); mpfr_init(r24306); mpfr_init(r24307); mpfr_init(r24308); mpfr_init(r24309); mpfr_init(r24310); mpfr_init(r24311); } double f_im(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r24295, phi1, MPFR_RNDN); mpfr_sin(r24296, r24295, MPFR_RNDN); mpfr_set_d(r24297, phi2, MPFR_RNDN); mpfr_sin(r24298, r24297, MPFR_RNDN); mpfr_mul(r24299, r24296, r24298, MPFR_RNDN); mpfr_cos(r24300, r24295, MPFR_RNDN); mpfr_cos(r24301, r24297, MPFR_RNDN); mpfr_mul(r24302, r24300, r24301, MPFR_RNDN); mpfr_set_d(r24303, lambda1, MPFR_RNDN); mpfr_set_d(r24304, lambda2, MPFR_RNDN); mpfr_sub(r24305, r24303, r24304, MPFR_RNDN); mpfr_cos(r24306, r24305, MPFR_RNDN); mpfr_mul(r24307, r24302, r24306, MPFR_RNDN); mpfr_add(r24308, r24299, r24307, MPFR_RNDN); mpfr_acos(r24309, r24308, MPFR_RNDN); mpfr_set_d(r24310, R, MPFR_RNDN); mpfr_mul(r24311, r24309, r24310, MPFR_RNDN); return mpfr_get_d(r24311, MPFR_RNDN); } static mpfr_t r24312, r24313, r24314, r24315, r24316, r24317, r24318, r24319, r24320, r24321, r24322, r24323, r24324, r24325, r24326, r24327, r24328, r24329, r24330, r24331, r24332, r24333, r24334, r24335; void setup_mpfr_f_fm() { mpfr_set_default_prec(2128); mpfr_init(r24312); mpfr_init(r24313); mpfr_init(r24314); mpfr_init(r24315); mpfr_init(r24316); mpfr_init(r24317); mpfr_init(r24318); mpfr_init(r24319); mpfr_init(r24320); mpfr_init(r24321); mpfr_init(r24322); mpfr_init(r24323); mpfr_init(r24324); mpfr_init(r24325); mpfr_init(r24326); mpfr_init(r24327); mpfr_init(r24328); mpfr_init(r24329); mpfr_init(r24330); mpfr_init(r24331); mpfr_init(r24332); mpfr_init(r24333); mpfr_init(r24334); mpfr_init(r24335); } double f_fm(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r24312, phi2, MPFR_RNDN); mpfr_cos(r24313, r24312, MPFR_RNDN); mpfr_set_d(r24314, phi1, MPFR_RNDN); mpfr_cos(r24315, r24314, MPFR_RNDN); mpfr_mul(r24316, r24313, r24315, MPFR_RNDN); mpfr_set_d(r24317, lambda1, MPFR_RNDN); mpfr_cos(r24318, r24317, MPFR_RNDN); mpfr_set_d(r24319, lambda2, MPFR_RNDN); mpfr_cos(r24320, r24319, MPFR_RNDN); mpfr_mul(r24321, r24318, r24320, MPFR_RNDN); mpfr_sin(r24322, r24317, MPFR_RNDN); mpfr_neg(r24323, r24319, MPFR_RNDN); mpfr_sin(r24324, r24323, MPFR_RNDN); mpfr_mul(r24325, r24322, r24324, MPFR_RNDN); mpfr_log1p(r24326, r24325, MPFR_RNDN); mpfr_expm1(r24327, r24326, MPFR_RNDN); mpfr_sub(r24328, r24321, r24327, MPFR_RNDN); mpfr_sin(r24329, r24312, MPFR_RNDN); mpfr_sin(r24330, r24314, MPFR_RNDN); mpfr_mul(r24331, r24329, r24330, MPFR_RNDN); mpfr_fma(r24332, r24316, r24328, r24331, MPFR_RNDN); mpfr_acos(r24333, r24332, MPFR_RNDN); mpfr_set_d(r24334, R, MPFR_RNDN); mpfr_mul(r24335, r24333, r24334, MPFR_RNDN); return mpfr_get_d(r24335, MPFR_RNDN); } static mpfr_t r24336, r24337, r24338, r24339, r24340, r24341, r24342, r24343, r24344, r24345, r24346, r24347, r24348, r24349, r24350, r24351, r24352, r24353, r24354, r24355, r24356, r24357, r24358, r24359; void setup_mpfr_f_dm() { mpfr_set_default_prec(2128); mpfr_init(r24336); mpfr_init(r24337); mpfr_init(r24338); mpfr_init(r24339); mpfr_init(r24340); mpfr_init(r24341); mpfr_init(r24342); mpfr_init(r24343); mpfr_init(r24344); mpfr_init(r24345); mpfr_init(r24346); mpfr_init(r24347); mpfr_init(r24348); mpfr_init(r24349); mpfr_init(r24350); mpfr_init(r24351); mpfr_init(r24352); mpfr_init(r24353); mpfr_init(r24354); mpfr_init(r24355); mpfr_init(r24356); mpfr_init(r24357); mpfr_init(r24358); mpfr_init(r24359); } double f_dm(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r24336, phi2, MPFR_RNDN); mpfr_cos(r24337, r24336, MPFR_RNDN); mpfr_set_d(r24338, phi1, MPFR_RNDN); mpfr_cos(r24339, r24338, MPFR_RNDN); mpfr_mul(r24340, r24337, r24339, MPFR_RNDN); mpfr_set_d(r24341, lambda1, MPFR_RNDN); mpfr_cos(r24342, r24341, MPFR_RNDN); mpfr_set_d(r24343, lambda2, MPFR_RNDN); mpfr_cos(r24344, r24343, MPFR_RNDN); mpfr_mul(r24345, r24342, r24344, MPFR_RNDN); mpfr_sin(r24346, r24341, MPFR_RNDN); mpfr_neg(r24347, r24343, MPFR_RNDN); mpfr_sin(r24348, r24347, MPFR_RNDN); mpfr_mul(r24349, r24346, r24348, MPFR_RNDN); mpfr_log1p(r24350, r24349, MPFR_RNDN); mpfr_expm1(r24351, r24350, MPFR_RNDN); mpfr_sub(r24352, r24345, r24351, MPFR_RNDN); mpfr_sin(r24353, r24336, MPFR_RNDN); mpfr_sin(r24354, r24338, MPFR_RNDN); mpfr_mul(r24355, r24353, r24354, MPFR_RNDN); mpfr_fma(r24356, r24340, r24352, r24355, MPFR_RNDN); mpfr_acos(r24357, r24356, MPFR_RNDN); mpfr_set_d(r24358, R, MPFR_RNDN); mpfr_mul(r24359, r24357, r24358, MPFR_RNDN); return mpfr_get_d(r24359, MPFR_RNDN); }