#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 r8022 = phi1; float r8023 = sin(r8022); float r8024 = phi2; float r8025 = sin(r8024); float r8026 = r8023 * r8025; float r8027 = cos(r8022); float r8028 = cos(r8024); float r8029 = r8027 * r8028; float r8030 = lambda1; float r8031 = lambda2; float r8032 = r8030 - r8031; float r8033 = cos(r8032); float r8034 = r8029 * r8033; float r8035 = r8026 + r8034; float r8036 = acos(r8035); float r8037 = R; float r8038 = r8036 * r8037; return r8038; } double f_id(double R, double lambda1, double lambda2, double phi1, double phi2) { double r8039 = phi1; double r8040 = sin(r8039); double r8041 = phi2; double r8042 = sin(r8041); double r8043 = r8040 * r8042; double r8044 = cos(r8039); double r8045 = cos(r8041); double r8046 = r8044 * r8045; double r8047 = lambda1; double r8048 = lambda2; double r8049 = r8047 - r8048; double r8050 = cos(r8049); double r8051 = r8046 * r8050; double r8052 = r8043 + r8051; double r8053 = acos(r8052); double r8054 = R; double r8055 = r8053 * r8054; return r8055; } double f_of(float R, float lambda1, float lambda2, float phi1, float phi2) { float r8056 = lambda2; float r8057 = sin(r8056); float r8058 = lambda1; float r8059 = sin(r8058); float r8060 = cos(r8056); float r8061 = cos(r8058); float r8062 = r8060 * r8061; float r8063 = fma(r8057, r8059, r8062); float r8064 = phi1; float r8065 = cos(r8064); float r8066 = phi2; float r8067 = cos(r8066); float r8068 = r8065 * r8067; float r8069 = sin(r8066); float r8070 = sin(r8064); float r8071 = r8069 * r8070; float r8072 = fma(r8063, r8068, r8071); float r8073 = expm1(r8072); float r8074 = log1p(r8073); float r8075 = acos(r8074); float r8076 = exp(r8075); float r8077 = log(r8076); float r8078 = R; float r8079 = r8077 * r8078; return r8079; } double f_od(double R, double lambda1, double lambda2, double phi1, double phi2) { double r8080 = lambda2; double r8081 = sin(r8080); double r8082 = lambda1; double r8083 = sin(r8082); double r8084 = cos(r8080); double r8085 = cos(r8082); double r8086 = r8084 * r8085; double r8087 = fma(r8081, r8083, r8086); double r8088 = phi1; double r8089 = cos(r8088); double r8090 = phi2; double r8091 = cos(r8090); double r8092 = r8089 * r8091; double r8093 = sin(r8090); double r8094 = sin(r8088); double r8095 = r8093 * r8094; double r8096 = fma(r8087, r8092, r8095); double r8097 = expm1(r8096); double r8098 = log1p(r8097); double r8099 = acos(r8098); double r8100 = exp(r8099); double r8101 = log(r8100); double r8102 = R; double r8103 = r8101 * r8102; return r8103; } 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 r8104, r8105, r8106, r8107, r8108, r8109, r8110, r8111, r8112, r8113, r8114, r8115, r8116, r8117, r8118, r8119, r8120; void setup_mpfr_f_im() { mpfr_set_default_prec(2128); mpfr_init(r8104); mpfr_init(r8105); mpfr_init(r8106); mpfr_init(r8107); mpfr_init(r8108); mpfr_init(r8109); mpfr_init(r8110); mpfr_init(r8111); mpfr_init(r8112); mpfr_init(r8113); mpfr_init(r8114); mpfr_init(r8115); mpfr_init(r8116); mpfr_init(r8117); mpfr_init(r8118); mpfr_init(r8119); mpfr_init(r8120); } double f_im(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r8104, phi1, MPFR_RNDN); mpfr_sin(r8105, r8104, MPFR_RNDN); mpfr_set_d(r8106, phi2, MPFR_RNDN); mpfr_sin(r8107, r8106, MPFR_RNDN); mpfr_mul(r8108, r8105, r8107, MPFR_RNDN); mpfr_cos(r8109, r8104, MPFR_RNDN); mpfr_cos(r8110, r8106, MPFR_RNDN); mpfr_mul(r8111, r8109, r8110, MPFR_RNDN); mpfr_set_d(r8112, lambda1, MPFR_RNDN); mpfr_set_d(r8113, lambda2, MPFR_RNDN); mpfr_sub(r8114, r8112, r8113, MPFR_RNDN); mpfr_cos(r8115, r8114, MPFR_RNDN); mpfr_mul(r8116, r8111, r8115, MPFR_RNDN); mpfr_add(r8117, r8108, r8116, MPFR_RNDN); mpfr_acos(r8118, r8117, MPFR_RNDN); mpfr_set_d(r8119, R, MPFR_RNDN); mpfr_mul(r8120, r8118, r8119, MPFR_RNDN); return mpfr_get_d(r8120, MPFR_RNDN); } static mpfr_t r8121, r8122, r8123, r8124, r8125, r8126, r8127, r8128, r8129, r8130, r8131, r8132, r8133, r8134, r8135, r8136, r8137, r8138, r8139, r8140, r8141, r8142, r8143, r8144; void setup_mpfr_f_fm() { mpfr_set_default_prec(2128); mpfr_init(r8121); mpfr_init(r8122); mpfr_init(r8123); mpfr_init(r8124); mpfr_init(r8125); mpfr_init(r8126); mpfr_init(r8127); mpfr_init(r8128); mpfr_init(r8129); mpfr_init(r8130); mpfr_init(r8131); mpfr_init(r8132); mpfr_init(r8133); mpfr_init(r8134); mpfr_init(r8135); mpfr_init(r8136); mpfr_init(r8137); mpfr_init(r8138); mpfr_init(r8139); mpfr_init(r8140); mpfr_init(r8141); mpfr_init(r8142); mpfr_init(r8143); mpfr_init(r8144); } double f_fm(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r8121, lambda2, MPFR_RNDN); mpfr_sin(r8122, r8121, MPFR_RNDN); mpfr_set_d(r8123, lambda1, MPFR_RNDN); mpfr_sin(r8124, r8123, MPFR_RNDN); mpfr_cos(r8125, r8121, MPFR_RNDN); mpfr_cos(r8126, r8123, MPFR_RNDN); mpfr_mul(r8127, r8125, r8126, MPFR_RNDN); mpfr_fma(r8128, r8122, r8124, r8127, MPFR_RNDN); mpfr_set_d(r8129, phi1, MPFR_RNDN); mpfr_cos(r8130, r8129, MPFR_RNDN); mpfr_set_d(r8131, phi2, MPFR_RNDN); mpfr_cos(r8132, r8131, MPFR_RNDN); mpfr_mul(r8133, r8130, r8132, MPFR_RNDN); mpfr_sin(r8134, r8131, MPFR_RNDN); mpfr_sin(r8135, r8129, MPFR_RNDN); mpfr_mul(r8136, r8134, r8135, MPFR_RNDN); mpfr_fma(r8137, r8128, r8133, r8136, MPFR_RNDN); mpfr_expm1(r8138, r8137, MPFR_RNDN); mpfr_log1p(r8139, r8138, MPFR_RNDN); mpfr_acos(r8140, r8139, MPFR_RNDN); mpfr_exp(r8141, r8140, MPFR_RNDN); mpfr_log(r8142, r8141, MPFR_RNDN); mpfr_set_d(r8143, R, MPFR_RNDN); mpfr_mul(r8144, r8142, r8143, MPFR_RNDN); return mpfr_get_d(r8144, MPFR_RNDN); } static mpfr_t r8145, r8146, r8147, r8148, r8149, r8150, r8151, r8152, r8153, r8154, r8155, r8156, r8157, r8158, r8159, r8160, r8161, r8162, r8163, r8164, r8165, r8166, r8167, r8168; void setup_mpfr_f_dm() { mpfr_set_default_prec(2128); mpfr_init(r8145); mpfr_init(r8146); mpfr_init(r8147); mpfr_init(r8148); mpfr_init(r8149); mpfr_init(r8150); mpfr_init(r8151); mpfr_init(r8152); mpfr_init(r8153); mpfr_init(r8154); mpfr_init(r8155); mpfr_init(r8156); mpfr_init(r8157); mpfr_init(r8158); mpfr_init(r8159); mpfr_init(r8160); mpfr_init(r8161); mpfr_init(r8162); mpfr_init(r8163); mpfr_init(r8164); mpfr_init(r8165); mpfr_init(r8166); mpfr_init(r8167); mpfr_init(r8168); } double f_dm(double R, double lambda1, double lambda2, double phi1, double phi2) { mpfr_set_d(r8145, lambda2, MPFR_RNDN); mpfr_sin(r8146, r8145, MPFR_RNDN); mpfr_set_d(r8147, lambda1, MPFR_RNDN); mpfr_sin(r8148, r8147, MPFR_RNDN); mpfr_cos(r8149, r8145, MPFR_RNDN); mpfr_cos(r8150, r8147, MPFR_RNDN); mpfr_mul(r8151, r8149, r8150, MPFR_RNDN); mpfr_fma(r8152, r8146, r8148, r8151, MPFR_RNDN); mpfr_set_d(r8153, phi1, MPFR_RNDN); mpfr_cos(r8154, r8153, MPFR_RNDN); mpfr_set_d(r8155, phi2, MPFR_RNDN); mpfr_cos(r8156, r8155, MPFR_RNDN); mpfr_mul(r8157, r8154, r8156, MPFR_RNDN); mpfr_sin(r8158, r8155, MPFR_RNDN); mpfr_sin(r8159, r8153, MPFR_RNDN); mpfr_mul(r8160, r8158, r8159, MPFR_RNDN); mpfr_fma(r8161, r8152, r8157, r8160, MPFR_RNDN); mpfr_expm1(r8162, r8161, MPFR_RNDN); mpfr_log1p(r8163, r8162, MPFR_RNDN); mpfr_acos(r8164, r8163, MPFR_RNDN); mpfr_exp(r8165, r8164, MPFR_RNDN); mpfr_log(r8166, r8165, MPFR_RNDN); mpfr_set_d(r8167, R, MPFR_RNDN); mpfr_mul(r8168, r8166, r8167, MPFR_RNDN); return mpfr_get_d(r8168, MPFR_RNDN); }