#include #include #include #include #include char *name = "Henrywood and Agarwal, Equation (12)"; double f_if(float d, float h, float l, float M, float D) { float r30464051 = d; float r30464052 = h; float r30464053 = r30464051 / r30464052; float r30464054 = 1; float r30464055 = 2; float r30464056 = r30464054 / r30464055; float r30464057 = pow(r30464053, r30464056); float r30464058 = l; float r30464059 = r30464051 / r30464058; float r30464060 = pow(r30464059, r30464056); float r30464061 = r30464057 * r30464060; float r30464062 = M; float r30464063 = D; float r30464064 = r30464062 * r30464063; float r30464065 = r30464055 * r30464051; float r30464066 = r30464064 / r30464065; float r30464067 = pow(r30464066, r30464055); float r30464068 = r30464056 * r30464067; float r30464069 = r30464052 / r30464058; float r30464070 = r30464068 * r30464069; float r30464071 = r30464054 - r30464070; float r30464072 = r30464061 * r30464071; return r30464072; } double f_id(double d, double h, double l, double M, double D) { double r30464073 = d; double r30464074 = h; double r30464075 = r30464073 / r30464074; double r30464076 = 1; double r30464077 = 2; double r30464078 = r30464076 / r30464077; double r30464079 = pow(r30464075, r30464078); double r30464080 = l; double r30464081 = r30464073 / r30464080; double r30464082 = pow(r30464081, r30464078); double r30464083 = r30464079 * r30464082; double r30464084 = M; double r30464085 = D; double r30464086 = r30464084 * r30464085; double r30464087 = r30464077 * r30464073; double r30464088 = r30464086 / r30464087; double r30464089 = pow(r30464088, r30464077); double r30464090 = r30464078 * r30464089; double r30464091 = r30464074 / r30464080; double r30464092 = r30464090 * r30464091; double r30464093 = r30464076 - r30464092; double r30464094 = r30464083 * r30464093; return r30464094; } double f_of(float d, float h, float l, float M, float D) { float r30464095 = l; float r30464096 = -4.5653656343947e-310; bool r30464097 = r30464095 <= r30464096; float r30464098 = d; float r30464099 = h; float r30464100 = r30464098 / r30464099; float r30464101 = 1; float r30464102 = 2; float r30464103 = r30464101 / r30464102; float r30464104 = pow(r30464100, r30464103); float r30464105 = r30464098 / r30464095; float r30464106 = pow(r30464105, r30464103); float r30464107 = r30464104 * r30464106; float r30464108 = M; float r30464109 = D; float r30464110 = r30464108 * r30464109; float r30464111 = r30464102 * r30464098; float r30464112 = r30464110 / r30464111; float r30464113 = pow(r30464112, r30464102); float r30464114 = r30464103 * r30464113; float r30464115 = r30464114 * r30464099; float r30464116 = r30464101 / r30464095; float r30464117 = r30464115 * r30464116; float r30464118 = r30464101 - r30464117; float r30464119 = r30464107 * r30464118; float r30464120 = 1.3970315963788558e-77; bool r30464121 = r30464095 <= r30464120; float r30464122 = pow(r30464098, r30464103); float r30464123 = pow(r30464116, r30464103); float r30464124 = r30464122 * r30464123; float r30464125 = r30464104 * r30464124; float r30464126 = r30464101 * r30464113; float r30464127 = r30464126 * r30464099; float r30464128 = r30464102 * r30464095; float r30464129 = r30464127 / r30464128; float r30464130 = r30464101 - r30464129; float r30464131 = r30464125 * r30464130; float r30464132 = 5.781959260454722e+222; bool r30464133 = r30464095 <= r30464132; float r30464134 = r30464101 / r30464099; float r30464135 = pow(r30464134, r30464103); float r30464136 = r30464122 * r30464135; float r30464137 = r30464136 * r30464106; float r30464138 = r30464137 * r30464130; float r30464139 = cbrt(r30464100); float r30464140 = r30464139 * r30464139; float r30464141 = pow(r30464140, r30464103); float r30464142 = pow(r30464139, r30464103); float r30464143 = r30464141 * r30464142; float r30464144 = r30464143 * r30464106; float r30464145 = r30464099 / r30464095; float r30464146 = r30464114 * r30464145; float r30464147 = r30464101 - r30464146; float r30464148 = r30464144 * r30464147; float r30464149 = r30464133 ? r30464138 : r30464148; float r30464150 = r30464121 ? r30464131 : r30464149; float r30464151 = r30464097 ? r30464119 : r30464150; return r30464151; } double f_od(double d, double h, double l, double M, double D) { double r30464152 = l; double r30464153 = -4.5653656343947e-310; bool r30464154 = r30464152 <= r30464153; double r30464155 = d; double r30464156 = h; double r30464157 = r30464155 / r30464156; double r30464158 = 1; double r30464159 = 2; double r30464160 = r30464158 / r30464159; double r30464161 = pow(r30464157, r30464160); double r30464162 = r30464155 / r30464152; double r30464163 = pow(r30464162, r30464160); double r30464164 = r30464161 * r30464163; double r30464165 = M; double r30464166 = D; double r30464167 = r30464165 * r30464166; double r30464168 = r30464159 * r30464155; double r30464169 = r30464167 / r30464168; double r30464170 = pow(r30464169, r30464159); double r30464171 = r30464160 * r30464170; double r30464172 = r30464171 * r30464156; double r30464173 = r30464158 / r30464152; double r30464174 = r30464172 * r30464173; double r30464175 = r30464158 - r30464174; double r30464176 = r30464164 * r30464175; double r30464177 = 1.3970315963788558e-77; bool r30464178 = r30464152 <= r30464177; double r30464179 = pow(r30464155, r30464160); double r30464180 = pow(r30464173, r30464160); double r30464181 = r30464179 * r30464180; double r30464182 = r30464161 * r30464181; double r30464183 = r30464158 * r30464170; double r30464184 = r30464183 * r30464156; double r30464185 = r30464159 * r30464152; double r30464186 = r30464184 / r30464185; double r30464187 = r30464158 - r30464186; double r30464188 = r30464182 * r30464187; double r30464189 = 5.781959260454722e+222; bool r30464190 = r30464152 <= r30464189; double r30464191 = r30464158 / r30464156; double r30464192 = pow(r30464191, r30464160); double r30464193 = r30464179 * r30464192; double r30464194 = r30464193 * r30464163; double r30464195 = r30464194 * r30464187; double r30464196 = cbrt(r30464157); double r30464197 = r30464196 * r30464196; double r30464198 = pow(r30464197, r30464160); double r30464199 = pow(r30464196, r30464160); double r30464200 = r30464198 * r30464199; double r30464201 = r30464200 * r30464163; double r30464202 = r30464156 / r30464152; double r30464203 = r30464171 * r30464202; double r30464204 = r30464158 - r30464203; double r30464205 = r30464201 * r30464204; double r30464206 = r30464190 ? r30464195 : r30464205; double r30464207 = r30464178 ? r30464188 : r30464206; double r30464208 = r30464154 ? r30464176 : r30464207; return r30464208; } 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 r30464209, r30464210, r30464211, r30464212, r30464213, r30464214, r30464215, r30464216, r30464217, r30464218, r30464219, r30464220, r30464221, r30464222, r30464223, r30464224, r30464225, r30464226, r30464227, r30464228, r30464229, r30464230; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r30464209); mpfr_init(r30464210); mpfr_init(r30464211); mpfr_init_set_str(r30464212, "1", 10, MPFR_RNDN); mpfr_init_set_str(r30464213, "2", 10, MPFR_RNDN); mpfr_init(r30464214); mpfr_init(r30464215); mpfr_init(r30464216); mpfr_init(r30464217); mpfr_init(r30464218); mpfr_init(r30464219); mpfr_init(r30464220); mpfr_init(r30464221); mpfr_init(r30464222); mpfr_init(r30464223); mpfr_init(r30464224); mpfr_init(r30464225); mpfr_init(r30464226); mpfr_init(r30464227); mpfr_init(r30464228); mpfr_init(r30464229); mpfr_init(r30464230); } double f_im(double d, double h, double l, double M, double D) { mpfr_set_d(r30464209, d, MPFR_RNDN); mpfr_set_d(r30464210, h, MPFR_RNDN); mpfr_div(r30464211, r30464209, r30464210, MPFR_RNDN); ; ; mpfr_div(r30464214, r30464212, r30464213, MPFR_RNDN); mpfr_pow(r30464215, r30464211, r30464214, MPFR_RNDN); mpfr_set_d(r30464216, l, MPFR_RNDN); mpfr_div(r30464217, r30464209, r30464216, MPFR_RNDN); mpfr_pow(r30464218, r30464217, r30464214, MPFR_RNDN); mpfr_mul(r30464219, r30464215, r30464218, MPFR_RNDN); mpfr_set_d(r30464220, M, MPFR_RNDN); mpfr_set_d(r30464221, D, MPFR_RNDN); mpfr_mul(r30464222, r30464220, r30464221, MPFR_RNDN); mpfr_mul(r30464223, r30464213, r30464209, MPFR_RNDN); mpfr_div(r30464224, r30464222, r30464223, MPFR_RNDN); mpfr_pow(r30464225, r30464224, r30464213, MPFR_RNDN); mpfr_mul(r30464226, r30464214, r30464225, MPFR_RNDN); mpfr_div(r30464227, r30464210, r30464216, MPFR_RNDN); mpfr_mul(r30464228, r30464226, r30464227, MPFR_RNDN); mpfr_sub(r30464229, r30464212, r30464228, MPFR_RNDN); mpfr_mul(r30464230, r30464219, r30464229, MPFR_RNDN); return mpfr_get_d(r30464230, MPFR_RNDN); } static mpfr_t r30464231, r30464232, r30464233, r30464234, r30464235, r30464236, r30464237, r30464238, r30464239, r30464240, r30464241, r30464242, r30464243, r30464244, r30464245, r30464246, r30464247, r30464248, r30464249, r30464250, r30464251, r30464252, r30464253, r30464254, r30464255, r30464256, r30464257, r30464258, r30464259, r30464260, r30464261, r30464262, r30464263, r30464264, r30464265, r30464266, r30464267, r30464268, r30464269, r30464270, r30464271, r30464272, r30464273, r30464274, r30464275, r30464276, r30464277, r30464278, r30464279, r30464280, r30464281, r30464282, r30464283, r30464284, r30464285, r30464286, r30464287; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r30464231); mpfr_init_set_str(r30464232, "-4.5653656343947e-310", 10, MPFR_RNDN); mpfr_init(r30464233); mpfr_init(r30464234); mpfr_init(r30464235); mpfr_init(r30464236); mpfr_init_set_str(r30464237, "1", 10, MPFR_RNDN); mpfr_init_set_str(r30464238, "2", 10, MPFR_RNDN); mpfr_init(r30464239); mpfr_init(r30464240); mpfr_init(r30464241); mpfr_init(r30464242); mpfr_init(r30464243); mpfr_init(r30464244); mpfr_init(r30464245); mpfr_init(r30464246); mpfr_init(r30464247); mpfr_init(r30464248); mpfr_init(r30464249); mpfr_init(r30464250); mpfr_init(r30464251); mpfr_init(r30464252); mpfr_init(r30464253); mpfr_init(r30464254); mpfr_init(r30464255); mpfr_init_set_str(r30464256, "1.3970315963788558e-77", 10, MPFR_RNDN); mpfr_init(r30464257); mpfr_init(r30464258); mpfr_init(r30464259); mpfr_init(r30464260); mpfr_init(r30464261); mpfr_init(r30464262); mpfr_init(r30464263); mpfr_init(r30464264); mpfr_init(r30464265); mpfr_init(r30464266); mpfr_init(r30464267); mpfr_init_set_str(r30464268, "5.781959260454722e+222", 10, MPFR_RNDN); mpfr_init(r30464269); mpfr_init(r30464270); mpfr_init(r30464271); mpfr_init(r30464272); mpfr_init(r30464273); mpfr_init(r30464274); mpfr_init(r30464275); mpfr_init(r30464276); mpfr_init(r30464277); mpfr_init(r30464278); mpfr_init(r30464279); mpfr_init(r30464280); mpfr_init(r30464281); mpfr_init(r30464282); mpfr_init(r30464283); mpfr_init(r30464284); mpfr_init(r30464285); mpfr_init(r30464286); mpfr_init(r30464287); } double f_fm(double d, double h, double l, double M, double D) { mpfr_set_d(r30464231, l, MPFR_RNDN); ; mpfr_set_si(r30464233, mpfr_cmp(r30464231, r30464232) <= 0, MPFR_RNDN); mpfr_set_d(r30464234, d, MPFR_RNDN); mpfr_set_d(r30464235, h, MPFR_RNDN); mpfr_div(r30464236, r30464234, r30464235, MPFR_RNDN); ; ; mpfr_div(r30464239, r30464237, r30464238, MPFR_RNDN); mpfr_pow(r30464240, r30464236, r30464239, MPFR_RNDN); mpfr_div(r30464241, r30464234, r30464231, MPFR_RNDN); mpfr_pow(r30464242, r30464241, r30464239, MPFR_RNDN); mpfr_mul(r30464243, r30464240, r30464242, MPFR_RNDN); mpfr_set_d(r30464244, M, MPFR_RNDN); mpfr_set_d(r30464245, D, MPFR_RNDN); mpfr_mul(r30464246, r30464244, r30464245, MPFR_RNDN); mpfr_mul(r30464247, r30464238, r30464234, MPFR_RNDN); mpfr_div(r30464248, r30464246, r30464247, MPFR_RNDN); mpfr_pow(r30464249, r30464248, r30464238, MPFR_RNDN); mpfr_mul(r30464250, r30464239, r30464249, MPFR_RNDN); mpfr_mul(r30464251, r30464250, r30464235, MPFR_RNDN); mpfr_div(r30464252, r30464237, r30464231, MPFR_RNDN); mpfr_mul(r30464253, r30464251, r30464252, MPFR_RNDN); mpfr_sub(r30464254, r30464237, r30464253, MPFR_RNDN); mpfr_mul(r30464255, r30464243, r30464254, MPFR_RNDN); ; mpfr_set_si(r30464257, mpfr_cmp(r30464231, r30464256) <= 0, MPFR_RNDN); mpfr_pow(r30464258, r30464234, r30464239, MPFR_RNDN); mpfr_pow(r30464259, r30464252, r30464239, MPFR_RNDN); mpfr_mul(r30464260, r30464258, r30464259, MPFR_RNDN); mpfr_mul(r30464261, r30464240, r30464260, MPFR_RNDN); mpfr_mul(r30464262, r30464237, r30464249, MPFR_RNDN); mpfr_mul(r30464263, r30464262, r30464235, MPFR_RNDN); mpfr_mul(r30464264, r30464238, r30464231, MPFR_RNDN); mpfr_div(r30464265, r30464263, r30464264, MPFR_RNDN); mpfr_sub(r30464266, r30464237, r30464265, MPFR_RNDN); mpfr_mul(r30464267, r30464261, r30464266, MPFR_RNDN); ; mpfr_set_si(r30464269, mpfr_cmp(r30464231, r30464268) <= 0, MPFR_RNDN); mpfr_div(r30464270, r30464237, r30464235, MPFR_RNDN); mpfr_pow(r30464271, r30464270, r30464239, MPFR_RNDN); mpfr_mul(r30464272, r30464258, r30464271, MPFR_RNDN); mpfr_mul(r30464273, r30464272, r30464242, MPFR_RNDN); mpfr_mul(r30464274, r30464273, r30464266, MPFR_RNDN); mpfr_cbrt(r30464275, r30464236, MPFR_RNDN); mpfr_mul(r30464276, r30464275, r30464275, MPFR_RNDN); mpfr_pow(r30464277, r30464276, r30464239, MPFR_RNDN); mpfr_pow(r30464278, r30464275, r30464239, MPFR_RNDN); mpfr_mul(r30464279, r30464277, r30464278, MPFR_RNDN); mpfr_mul(r30464280, r30464279, r30464242, MPFR_RNDN); mpfr_div(r30464281, r30464235, r30464231, MPFR_RNDN); mpfr_mul(r30464282, r30464250, r30464281, MPFR_RNDN); mpfr_sub(r30464283, r30464237, r30464282, MPFR_RNDN); mpfr_mul(r30464284, r30464280, r30464283, MPFR_RNDN); if (mpfr_get_si(r30464269, MPFR_RNDN)) { mpfr_set(r30464285, r30464274, MPFR_RNDN); } else { mpfr_set(r30464285, r30464284, MPFR_RNDN); }; if (mpfr_get_si(r30464257, MPFR_RNDN)) { mpfr_set(r30464286, r30464267, MPFR_RNDN); } else { mpfr_set(r30464286, r30464285, MPFR_RNDN); }; if (mpfr_get_si(r30464233, MPFR_RNDN)) { mpfr_set(r30464287, r30464255, MPFR_RNDN); } else { mpfr_set(r30464287, r30464286, MPFR_RNDN); }; return mpfr_get_d(r30464287, MPFR_RNDN); } static mpfr_t r30464288, r30464289, r30464290, r30464291, r30464292, r30464293, r30464294, r30464295, r30464296, r30464297, r30464298, r30464299, r30464300, r30464301, r30464302, r30464303, r30464304, r30464305, r30464306, r30464307, r30464308, r30464309, r30464310, r30464311, r30464312, r30464313, r30464314, r30464315, r30464316, r30464317, r30464318, r30464319, r30464320, r30464321, r30464322, r30464323, r30464324, r30464325, r30464326, r30464327, r30464328, r30464329, r30464330, r30464331, r30464332, r30464333, r30464334, r30464335, r30464336, r30464337, r30464338, r30464339, r30464340, r30464341, r30464342, r30464343, r30464344; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r30464288); mpfr_init_set_str(r30464289, "-4.5653656343947e-310", 10, MPFR_RNDN); mpfr_init(r30464290); mpfr_init(r30464291); mpfr_init(r30464292); mpfr_init(r30464293); mpfr_init_set_str(r30464294, "1", 10, MPFR_RNDN); mpfr_init_set_str(r30464295, "2", 10, MPFR_RNDN); mpfr_init(r30464296); mpfr_init(r30464297); mpfr_init(r30464298); mpfr_init(r30464299); mpfr_init(r30464300); mpfr_init(r30464301); mpfr_init(r30464302); mpfr_init(r30464303); mpfr_init(r30464304); mpfr_init(r30464305); mpfr_init(r30464306); mpfr_init(r30464307); mpfr_init(r30464308); mpfr_init(r30464309); mpfr_init(r30464310); mpfr_init(r30464311); mpfr_init(r30464312); mpfr_init_set_str(r30464313, "1.3970315963788558e-77", 10, MPFR_RNDN); mpfr_init(r30464314); mpfr_init(r30464315); mpfr_init(r30464316); mpfr_init(r30464317); mpfr_init(r30464318); mpfr_init(r30464319); mpfr_init(r30464320); mpfr_init(r30464321); mpfr_init(r30464322); mpfr_init(r30464323); mpfr_init(r30464324); mpfr_init_set_str(r30464325, "5.781959260454722e+222", 10, MPFR_RNDN); mpfr_init(r30464326); mpfr_init(r30464327); mpfr_init(r30464328); mpfr_init(r30464329); mpfr_init(r30464330); mpfr_init(r30464331); mpfr_init(r30464332); mpfr_init(r30464333); mpfr_init(r30464334); mpfr_init(r30464335); mpfr_init(r30464336); mpfr_init(r30464337); mpfr_init(r30464338); mpfr_init(r30464339); mpfr_init(r30464340); mpfr_init(r30464341); mpfr_init(r30464342); mpfr_init(r30464343); mpfr_init(r30464344); } double f_dm(double d, double h, double l, double M, double D) { mpfr_set_d(r30464288, l, MPFR_RNDN); ; mpfr_set_si(r30464290, mpfr_cmp(r30464288, r30464289) <= 0, MPFR_RNDN); mpfr_set_d(r30464291, d, MPFR_RNDN); mpfr_set_d(r30464292, h, MPFR_RNDN); mpfr_div(r30464293, r30464291, r30464292, MPFR_RNDN); ; ; mpfr_div(r30464296, r30464294, r30464295, MPFR_RNDN); mpfr_pow(r30464297, r30464293, r30464296, MPFR_RNDN); mpfr_div(r30464298, r30464291, r30464288, MPFR_RNDN); mpfr_pow(r30464299, r30464298, r30464296, MPFR_RNDN); mpfr_mul(r30464300, r30464297, r30464299, MPFR_RNDN); mpfr_set_d(r30464301, M, MPFR_RNDN); mpfr_set_d(r30464302, D, MPFR_RNDN); mpfr_mul(r30464303, r30464301, r30464302, MPFR_RNDN); mpfr_mul(r30464304, r30464295, r30464291, MPFR_RNDN); mpfr_div(r30464305, r30464303, r30464304, MPFR_RNDN); mpfr_pow(r30464306, r30464305, r30464295, MPFR_RNDN); mpfr_mul(r30464307, r30464296, r30464306, MPFR_RNDN); mpfr_mul(r30464308, r30464307, r30464292, MPFR_RNDN); mpfr_div(r30464309, r30464294, r30464288, MPFR_RNDN); mpfr_mul(r30464310, r30464308, r30464309, MPFR_RNDN); mpfr_sub(r30464311, r30464294, r30464310, MPFR_RNDN); mpfr_mul(r30464312, r30464300, r30464311, MPFR_RNDN); ; mpfr_set_si(r30464314, mpfr_cmp(r30464288, r30464313) <= 0, MPFR_RNDN); mpfr_pow(r30464315, r30464291, r30464296, MPFR_RNDN); mpfr_pow(r30464316, r30464309, r30464296, MPFR_RNDN); mpfr_mul(r30464317, r30464315, r30464316, MPFR_RNDN); mpfr_mul(r30464318, r30464297, r30464317, MPFR_RNDN); mpfr_mul(r30464319, r30464294, r30464306, MPFR_RNDN); mpfr_mul(r30464320, r30464319, r30464292, MPFR_RNDN); mpfr_mul(r30464321, r30464295, r30464288, MPFR_RNDN); mpfr_div(r30464322, r30464320, r30464321, MPFR_RNDN); mpfr_sub(r30464323, r30464294, r30464322, MPFR_RNDN); mpfr_mul(r30464324, r30464318, r30464323, MPFR_RNDN); ; mpfr_set_si(r30464326, mpfr_cmp(r30464288, r30464325) <= 0, MPFR_RNDN); mpfr_div(r30464327, r30464294, r30464292, MPFR_RNDN); mpfr_pow(r30464328, r30464327, r30464296, MPFR_RNDN); mpfr_mul(r30464329, r30464315, r30464328, MPFR_RNDN); mpfr_mul(r30464330, r30464329, r30464299, MPFR_RNDN); mpfr_mul(r30464331, r30464330, r30464323, MPFR_RNDN); mpfr_cbrt(r30464332, r30464293, MPFR_RNDN); mpfr_mul(r30464333, r30464332, r30464332, MPFR_RNDN); mpfr_pow(r30464334, r30464333, r30464296, MPFR_RNDN); mpfr_pow(r30464335, r30464332, r30464296, MPFR_RNDN); mpfr_mul(r30464336, r30464334, r30464335, MPFR_RNDN); mpfr_mul(r30464337, r30464336, r30464299, MPFR_RNDN); mpfr_div(r30464338, r30464292, r30464288, MPFR_RNDN); mpfr_mul(r30464339, r30464307, r30464338, MPFR_RNDN); mpfr_sub(r30464340, r30464294, r30464339, MPFR_RNDN); mpfr_mul(r30464341, r30464337, r30464340, MPFR_RNDN); if (mpfr_get_si(r30464326, MPFR_RNDN)) { mpfr_set(r30464342, r30464331, MPFR_RNDN); } else { mpfr_set(r30464342, r30464341, MPFR_RNDN); }; if (mpfr_get_si(r30464314, MPFR_RNDN)) { mpfr_set(r30464343, r30464324, MPFR_RNDN); } else { mpfr_set(r30464343, r30464342, MPFR_RNDN); }; if (mpfr_get_si(r30464290, MPFR_RNDN)) { mpfr_set(r30464344, r30464312, MPFR_RNDN); } else { mpfr_set(r30464344, r30464343, MPFR_RNDN); }; return mpfr_get_d(r30464344, MPFR_RNDN); }