#include #include #include #include #include char *name = "Distance on a great circle"; double f_if(float R, float lambda1, float lambda2, float phi1, float phi2) { float r23212 = R; float r23213 = 2; float r23214 = phi1; float r23215 = phi2; float r23216 = r23214 - r23215; float r23217 = r23216 / r23213; float r23218 = sin(r23217); float r23219 = pow(r23218, r23213); float r23220 = cos(r23214); float r23221 = cos(r23215); float r23222 = r23220 * r23221; float r23223 = lambda1; float r23224 = lambda2; float r23225 = r23223 - r23224; float r23226 = r23225 / r23213; float r23227 = sin(r23226); float r23228 = r23222 * r23227; float r23229 = r23228 * r23227; float r23230 = r23219 + r23229; float r23231 = sqrt(r23230); float r23232 = 1; float r23233 = r23232 - r23230; float r23234 = sqrt(r23233); float r23235 = atan2(r23231, r23234); float r23236 = r23213 * r23235; float r23237 = r23212 * r23236; return r23237; } double f_id(double R, double lambda1, double lambda2, double phi1, double phi2) { double r23238 = R; double r23239 = 2; double r23240 = phi1; double r23241 = phi2; double r23242 = r23240 - r23241; double r23243 = r23242 / r23239; double r23244 = sin(r23243); double r23245 = pow(r23244, r23239); double r23246 = cos(r23240); double r23247 = cos(r23241); double r23248 = r23246 * r23247; double r23249 = lambda1; double r23250 = lambda2; double r23251 = r23249 - r23250; double r23252 = r23251 / r23239; double r23253 = sin(r23252); double r23254 = r23248 * r23253; double r23255 = r23254 * r23253; double r23256 = r23245 + r23255; double r23257 = sqrt(r23256); double r23258 = 1; double r23259 = r23258 - r23256; double r23260 = sqrt(r23259); double r23261 = atan2(r23257, r23260); double r23262 = r23239 * r23261; double r23263 = r23238 * r23262; return r23263; } double f_of(float R, float lambda1, float lambda2, float phi1, float phi2) { float r23264 = R; float r23265 = 2; float r23266 = phi1; float r23267 = phi2; float r23268 = r23266 - r23267; float r23269 = r23268 / r23265; float r23270 = sin(r23269); float r23271 = pow(r23270, r23265); float r23272 = cos(r23266); float r23273 = cos(r23267); float r23274 = r23272 * r23273; float r23275 = lambda1; float r23276 = lambda2; float r23277 = r23275 - r23276; float r23278 = r23277 / r23265; float r23279 = sin(r23278); float r23280 = expm1(r23279); float r23281 = log1p(r23280); float r23282 = r23274 * r23281; float r23283 = r23282 * r23279; float r23284 = r23271 + r23283; float r23285 = sqrt(r23284); float r23286 = 1; float r23287 = r23286 - r23284; float r23288 = sqrt(r23287); float r23289 = atan2(r23285, r23288); float r23290 = r23265 * r23289; float r23291 = r23264 * r23290; return r23291; } double f_od(double R, double lambda1, double lambda2, double phi1, double phi2) { double r23292 = R; double r23293 = 2; double r23294 = phi1; double r23295 = phi2; double r23296 = r23294 - r23295; double r23297 = r23296 / r23293; double r23298 = sin(r23297); double r23299 = pow(r23298, r23293); double r23300 = cos(r23294); double r23301 = cos(r23295); double r23302 = r23300 * r23301; double r23303 = lambda1; double r23304 = lambda2; double r23305 = r23303 - r23304; double r23306 = r23305 / r23293; double r23307 = sin(r23306); double r23308 = expm1(r23307); double r23309 = log1p(r23308); double r23310 = r23302 * r23309; double r23311 = r23310 * r23307; double r23312 = r23299 + r23311; double r23313 = sqrt(r23312); double r23314 = 1; double r23315 = r23314 - r23312; double r23316 = sqrt(r23315); double r23317 = atan2(r23313, r23316); double r23318 = r23293 * r23317; double r23319 = r23292 * r23318; return r23319; } 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 r23320, r23321, r23322, r23323, r23324, r23325, r23326, r23327, r23328, r23329, r23330, r23331, r23332, r23333, r23334, r23335, r23336, r23337, r23338, r23339, r23340, r23341, r23342, r23343, r23344, r23345; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r23320); mpfr_init_set_str(r23321, "2", 10, MPFR_RNDN); mpfr_init(r23322); mpfr_init(r23323); mpfr_init(r23324); mpfr_init(r23325); mpfr_init(r23326); mpfr_init(r23327); mpfr_init(r23328); mpfr_init(r23329); mpfr_init(r23330); mpfr_init(r23331); mpfr_init(r23332); mpfr_init(r23333); mpfr_init(r23334); mpfr_init(r23335); mpfr_init(r23336); mpfr_init(r23337); mpfr_init(r23338); mpfr_init(r23339); mpfr_init_set_str(r23340, "1", 10, MPFR_RNDN); mpfr_init(r23341); mpfr_init(r23342); mpfr_init(r23343); mpfr_init(r23344); mpfr_init(r23345); } double f_im(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r23320, R, MPFR_RNDN); ; mpfr_set_d(r23322, phi1, MPFR_RNDN); mpfr_set_d(r23323, phi2, MPFR_RNDN); mpfr_sub(r23324, r23322, r23323, MPFR_RNDN); mpfr_div(r23325, r23324, r23321, MPFR_RNDN); mpfr_sin(r23326, r23325, MPFR_RNDN); mpfr_pow(r23327, r23326, r23321, MPFR_RNDN); mpfr_cos(r23328, r23322, MPFR_RNDN); mpfr_cos(r23329, r23323, MPFR_RNDN); mpfr_mul(r23330, r23328, r23329, MPFR_RNDN); mpfr_set_d(r23331, lambda1, MPFR_RNDN); mpfr_set_d(r23332, lambda2, MPFR_RNDN); mpfr_sub(r23333, r23331, r23332, MPFR_RNDN); mpfr_div(r23334, r23333, r23321, MPFR_RNDN); mpfr_sin(r23335, r23334, MPFR_RNDN); mpfr_mul(r23336, r23330, r23335, MPFR_RNDN); mpfr_mul(r23337, r23336, r23335, MPFR_RNDN); mpfr_add(r23338, r23327, r23337, MPFR_RNDN); mpfr_sqrt(r23339, r23338, MPFR_RNDN); ; mpfr_sub(r23341, r23340, r23338, MPFR_RNDN); mpfr_sqrt(r23342, r23341, MPFR_RNDN); mpfr_atan2(r23343, r23339, r23342, MPFR_RNDN); mpfr_mul(r23344, r23321, r23343, MPFR_RNDN); mpfr_mul(r23345, r23320, r23344, MPFR_RNDN); return mpfr_get_d(r23345, MPFR_RNDN); } static mpfr_t r23346, r23347, r23348, r23349, r23350, r23351, r23352, r23353, r23354, r23355, r23356, r23357, r23358, r23359, r23360, r23361, r23362, r23363, r23364, r23365, r23366, r23367, r23368, r23369, r23370, r23371, r23372, r23373; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r23346); mpfr_init_set_str(r23347, "2", 10, MPFR_RNDN); mpfr_init(r23348); mpfr_init(r23349); mpfr_init(r23350); mpfr_init(r23351); mpfr_init(r23352); mpfr_init(r23353); mpfr_init(r23354); mpfr_init(r23355); mpfr_init(r23356); mpfr_init(r23357); mpfr_init(r23358); mpfr_init(r23359); mpfr_init(r23360); mpfr_init(r23361); mpfr_init(r23362); mpfr_init(r23363); mpfr_init(r23364); mpfr_init(r23365); mpfr_init(r23366); mpfr_init(r23367); mpfr_init_set_str(r23368, "1", 10, MPFR_RNDN); mpfr_init(r23369); mpfr_init(r23370); mpfr_init(r23371); mpfr_init(r23372); mpfr_init(r23373); } double f_fm(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r23346, R, MPFR_RNDN); ; mpfr_set_d(r23348, phi1, MPFR_RNDN); mpfr_set_d(r23349, phi2, MPFR_RNDN); mpfr_sub(r23350, r23348, r23349, MPFR_RNDN); mpfr_div(r23351, r23350, r23347, MPFR_RNDN); mpfr_sin(r23352, r23351, MPFR_RNDN); mpfr_pow(r23353, r23352, r23347, MPFR_RNDN); mpfr_cos(r23354, r23348, MPFR_RNDN); mpfr_cos(r23355, r23349, MPFR_RNDN); mpfr_mul(r23356, r23354, r23355, MPFR_RNDN); mpfr_set_d(r23357, lambda1, MPFR_RNDN); mpfr_set_d(r23358, lambda2, MPFR_RNDN); mpfr_sub(r23359, r23357, r23358, MPFR_RNDN); mpfr_div(r23360, r23359, r23347, MPFR_RNDN); mpfr_sin(r23361, r23360, MPFR_RNDN); mpfr_expm1(r23362, r23361, MPFR_RNDN); mpfr_log1p(r23363, r23362, MPFR_RNDN); mpfr_mul(r23364, r23356, r23363, MPFR_RNDN); mpfr_mul(r23365, r23364, r23361, MPFR_RNDN); mpfr_add(r23366, r23353, r23365, MPFR_RNDN); mpfr_sqrt(r23367, r23366, MPFR_RNDN); ; mpfr_sub(r23369, r23368, r23366, MPFR_RNDN); mpfr_sqrt(r23370, r23369, MPFR_RNDN); mpfr_atan2(r23371, r23367, r23370, MPFR_RNDN); mpfr_mul(r23372, r23347, r23371, MPFR_RNDN); mpfr_mul(r23373, r23346, r23372, MPFR_RNDN); return mpfr_get_d(r23373, MPFR_RNDN); } static mpfr_t r23374, r23375, r23376, r23377, r23378, r23379, r23380, r23381, r23382, r23383, r23384, r23385, r23386, r23387, r23388, r23389, r23390, r23391, r23392, r23393, r23394, r23395, r23396, r23397, r23398, r23399, r23400, r23401; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r23374); mpfr_init_set_str(r23375, "2", 10, MPFR_RNDN); mpfr_init(r23376); mpfr_init(r23377); mpfr_init(r23378); mpfr_init(r23379); mpfr_init(r23380); mpfr_init(r23381); mpfr_init(r23382); mpfr_init(r23383); mpfr_init(r23384); mpfr_init(r23385); mpfr_init(r23386); mpfr_init(r23387); mpfr_init(r23388); mpfr_init(r23389); mpfr_init(r23390); mpfr_init(r23391); mpfr_init(r23392); mpfr_init(r23393); mpfr_init(r23394); mpfr_init(r23395); mpfr_init_set_str(r23396, "1", 10, MPFR_RNDN); mpfr_init(r23397); mpfr_init(r23398); mpfr_init(r23399); mpfr_init(r23400); mpfr_init(r23401); } double f_dm(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r23374, R, MPFR_RNDN); ; mpfr_set_d(r23376, phi1, MPFR_RNDN); mpfr_set_d(r23377, phi2, MPFR_RNDN); mpfr_sub(r23378, r23376, r23377, MPFR_RNDN); mpfr_div(r23379, r23378, r23375, MPFR_RNDN); mpfr_sin(r23380, r23379, MPFR_RNDN); mpfr_pow(r23381, r23380, r23375, MPFR_RNDN); mpfr_cos(r23382, r23376, MPFR_RNDN); mpfr_cos(r23383, r23377, MPFR_RNDN); mpfr_mul(r23384, r23382, r23383, MPFR_RNDN); mpfr_set_d(r23385, lambda1, MPFR_RNDN); mpfr_set_d(r23386, lambda2, MPFR_RNDN); mpfr_sub(r23387, r23385, r23386, MPFR_RNDN); mpfr_div(r23388, r23387, r23375, MPFR_RNDN); mpfr_sin(r23389, r23388, MPFR_RNDN); mpfr_expm1(r23390, r23389, MPFR_RNDN); mpfr_log1p(r23391, r23390, MPFR_RNDN); mpfr_mul(r23392, r23384, r23391, MPFR_RNDN); mpfr_mul(r23393, r23392, r23389, MPFR_RNDN); mpfr_add(r23394, r23381, r23393, MPFR_RNDN); mpfr_sqrt(r23395, r23394, MPFR_RNDN); ; mpfr_sub(r23397, r23396, r23394, MPFR_RNDN); mpfr_sqrt(r23398, r23397, MPFR_RNDN); mpfr_atan2(r23399, r23395, r23398, MPFR_RNDN); mpfr_mul(r23400, r23375, r23399, MPFR_RNDN); mpfr_mul(r23401, r23374, r23400, MPFR_RNDN); return mpfr_get_d(r23401, MPFR_RNDN); }