#include #include #include #include #include char *name = "Toniolo and Linder, Equation (7)"; double f_if(float x, float l, float t) { float r8158 = 2; float r8159 = sqrt(r8158); float r8160 = t; float r8161 = r8159 * r8160; float r8162 = x; float r8163 = 1; float r8164 = r8162 + r8163; float r8165 = r8162 - r8163; float r8166 = r8164 / r8165; float r8167 = l; float r8168 = r8167 * r8167; float r8169 = r8160 * r8160; float r8170 = r8158 * r8169; float r8171 = r8168 + r8170; float r8172 = r8166 * r8171; float r8173 = r8172 - r8168; float r8174 = sqrt(r8173); float r8175 = r8161 / r8174; return r8175; } double f_id(double x, double l, double t) { double r8176 = 2; double r8177 = sqrt(r8176); double r8178 = t; double r8179 = r8177 * r8178; double r8180 = x; double r8181 = 1; double r8182 = r8180 + r8181; double r8183 = r8180 - r8181; double r8184 = r8182 / r8183; double r8185 = l; double r8186 = r8185 * r8185; double r8187 = r8178 * r8178; double r8188 = r8176 * r8187; double r8189 = r8186 + r8188; double r8190 = r8184 * r8189; double r8191 = r8190 - r8186; double r8192 = sqrt(r8191); double r8193 = r8179 / r8192; return r8193; } double f_of(float x, float l, float t) { float r8194 = t; float r8195 = -3.135428033554305e+38; bool r8196 = r8194 <= r8195; float r8197 = 2; float r8198 = sqrt(r8197); float r8199 = r8194 * r8198; float r8200 = x; float r8201 = r8200 * r8200; float r8202 = r8198 * r8201; float r8203 = r8194 / r8202; float r8204 = 1; float r8205 = r8204 - r8197; float r8206 = r8203 * r8205; float r8207 = r8197 / r8200; float r8208 = r8207 / r8198; float r8209 = r8208 + r8198; float r8210 = r8194 * r8209; float r8211 = r8206 - r8210; float r8212 = r8199 / r8211; float r8213 = -7.96461310098498e-207; bool r8214 = r8194 <= r8213; float r8215 = sqrt(r8198); float r8216 = r8194 * r8215; float r8217 = r8216 * r8215; float r8218 = l; float r8219 = r8197 * r8218; float r8220 = r8218 / r8200; float r8221 = r8219 * r8220; float r8222 = 4; float r8223 = r8222 / r8200; float r8224 = r8197 + r8223; float r8225 = r8194 * r8194; float r8226 = r8224 * r8225; float r8227 = r8221 + r8226; float r8228 = sqrt(r8227); float r8229 = r8217 / r8228; float r8230 = 2.986334750317526e-159; bool r8231 = r8194 <= r8230; float r8232 = r8197 - r8223; float r8233 = r8194 * r8200; float r8234 = r8233 * r8194; float r8235 = r8224 * r8234; float r8236 = r8218 * r8218; float r8237 = r8197 * r8236; float r8238 = r8235 + r8237; float r8239 = r8232 * r8238; float r8240 = sqrt(r8239); float r8241 = r8200 * r8197; float r8242 = r8241 - r8222; float r8243 = sqrt(r8242); float r8244 = r8240 / r8243; float r8245 = r8199 / r8244; float r8246 = 6.8042531379524426e+140; bool r8247 = r8194 <= r8246; float r8248 = r8194 / r8198; float r8249 = r8248 / r8201; float r8250 = r8197 - r8204; float r8251 = r8249 * r8250; float r8252 = r8251 + r8210; float r8253 = r8199 / r8252; float r8254 = r8247 ? r8229 : r8253; float r8255 = r8231 ? r8245 : r8254; float r8256 = r8214 ? r8229 : r8255; float r8257 = r8196 ? r8212 : r8256; return r8257; } double f_od(double x, double l, double t) { double r8258 = t; double r8259 = -3.135428033554305e+38; bool r8260 = r8258 <= r8259; double r8261 = 2; double r8262 = sqrt(r8261); double r8263 = r8258 * r8262; double r8264 = x; double r8265 = r8264 * r8264; double r8266 = r8262 * r8265; double r8267 = r8258 / r8266; double r8268 = 1; double r8269 = r8268 - r8261; double r8270 = r8267 * r8269; double r8271 = r8261 / r8264; double r8272 = r8271 / r8262; double r8273 = r8272 + r8262; double r8274 = r8258 * r8273; double r8275 = r8270 - r8274; double r8276 = r8263 / r8275; double r8277 = -7.96461310098498e-207; bool r8278 = r8258 <= r8277; double r8279 = sqrt(r8262); double r8280 = r8258 * r8279; double r8281 = r8280 * r8279; double r8282 = l; double r8283 = r8261 * r8282; double r8284 = r8282 / r8264; double r8285 = r8283 * r8284; double r8286 = 4; double r8287 = r8286 / r8264; double r8288 = r8261 + r8287; double r8289 = r8258 * r8258; double r8290 = r8288 * r8289; double r8291 = r8285 + r8290; double r8292 = sqrt(r8291); double r8293 = r8281 / r8292; double r8294 = 2.986334750317526e-159; bool r8295 = r8258 <= r8294; double r8296 = r8261 - r8287; double r8297 = r8258 * r8264; double r8298 = r8297 * r8258; double r8299 = r8288 * r8298; double r8300 = r8282 * r8282; double r8301 = r8261 * r8300; double r8302 = r8299 + r8301; double r8303 = r8296 * r8302; double r8304 = sqrt(r8303); double r8305 = r8264 * r8261; double r8306 = r8305 - r8286; double r8307 = sqrt(r8306); double r8308 = r8304 / r8307; double r8309 = r8263 / r8308; double r8310 = 6.8042531379524426e+140; bool r8311 = r8258 <= r8310; double r8312 = r8258 / r8262; double r8313 = r8312 / r8265; double r8314 = r8261 - r8268; double r8315 = r8313 * r8314; double r8316 = r8315 + r8274; double r8317 = r8263 / r8316; double r8318 = r8311 ? r8293 : r8317; double r8319 = r8295 ? r8309 : r8318; double r8320 = r8278 ? r8293 : r8319; double r8321 = r8260 ? r8276 : r8320; return r8321; } 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 r8322, r8323, r8324, r8325, r8326, r8327, r8328, r8329, r8330, r8331, r8332, r8333, r8334, r8335, r8336, r8337, r8338, r8339; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init_set_str(r8322, "2", 10, MPFR_RNDN); mpfr_init(r8323); mpfr_init(r8324); mpfr_init(r8325); mpfr_init(r8326); mpfr_init_set_str(r8327, "1", 10, MPFR_RNDN); mpfr_init(r8328); mpfr_init(r8329); mpfr_init(r8330); mpfr_init(r8331); mpfr_init(r8332); mpfr_init(r8333); mpfr_init(r8334); mpfr_init(r8335); mpfr_init(r8336); mpfr_init(r8337); mpfr_init(r8338); mpfr_init(r8339); } double f_im(double x, double l, double t) { ; mpfr_sqrt(r8323, r8322, MPFR_RNDN); mpfr_set_d(r8324, t, MPFR_RNDN); mpfr_mul(r8325, r8323, r8324, MPFR_RNDN); mpfr_set_d(r8326, x, MPFR_RNDN); ; mpfr_add(r8328, r8326, r8327, MPFR_RNDN); mpfr_sub(r8329, r8326, r8327, MPFR_RNDN); mpfr_div(r8330, r8328, r8329, MPFR_RNDN); mpfr_set_d(r8331, l, MPFR_RNDN); mpfr_mul(r8332, r8331, r8331, MPFR_RNDN); mpfr_mul(r8333, r8324, r8324, MPFR_RNDN); mpfr_mul(r8334, r8322, r8333, MPFR_RNDN); mpfr_add(r8335, r8332, r8334, MPFR_RNDN); mpfr_mul(r8336, r8330, r8335, MPFR_RNDN); mpfr_sub(r8337, r8336, r8332, MPFR_RNDN); mpfr_sqrt(r8338, r8337, MPFR_RNDN); mpfr_div(r8339, r8325, r8338, MPFR_RNDN); return mpfr_get_d(r8339, MPFR_RNDN); } static mpfr_t r8340, r8341, r8342, r8343, r8344, r8345, r8346, r8347, r8348, r8349, r8350, r8351, r8352, r8353, r8354, r8355, r8356, r8357, r8358, r8359, r8360, r8361, r8362, r8363, r8364, r8365, r8366, r8367, r8368, r8369, r8370, r8371, r8372, r8373, r8374, r8375, r8376, r8377, r8378, r8379, r8380, r8381, r8382, r8383, r8384, r8385, r8386, r8387, r8388, r8389, r8390, r8391, r8392, r8393, r8394, r8395, r8396, r8397, r8398, r8399, r8400, r8401, r8402, r8403; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r8340); mpfr_init_set_str(r8341, "-3.135428033554305e+38", 10, MPFR_RNDN); mpfr_init(r8342); mpfr_init_set_str(r8343, "2", 10, MPFR_RNDN); mpfr_init(r8344); mpfr_init(r8345); mpfr_init(r8346); mpfr_init(r8347); mpfr_init(r8348); mpfr_init(r8349); mpfr_init_set_str(r8350, "1", 10, MPFR_RNDN); mpfr_init(r8351); mpfr_init(r8352); mpfr_init(r8353); mpfr_init(r8354); mpfr_init(r8355); mpfr_init(r8356); mpfr_init(r8357); mpfr_init(r8358); mpfr_init_set_str(r8359, "-7.96461310098498e-207", 10, MPFR_RNDN); mpfr_init(r8360); mpfr_init(r8361); mpfr_init(r8362); mpfr_init(r8363); mpfr_init(r8364); mpfr_init(r8365); mpfr_init(r8366); mpfr_init(r8367); mpfr_init_set_str(r8368, "4", 10, MPFR_RNDN); mpfr_init(r8369); mpfr_init(r8370); mpfr_init(r8371); mpfr_init(r8372); mpfr_init(r8373); mpfr_init(r8374); mpfr_init(r8375); mpfr_init_set_str(r8376, "2.986334750317526e-159", 10, MPFR_RNDN); mpfr_init(r8377); mpfr_init(r8378); mpfr_init(r8379); mpfr_init(r8380); mpfr_init(r8381); mpfr_init(r8382); mpfr_init(r8383); mpfr_init(r8384); mpfr_init(r8385); mpfr_init(r8386); mpfr_init(r8387); mpfr_init(r8388); mpfr_init(r8389); mpfr_init(r8390); mpfr_init(r8391); mpfr_init_set_str(r8392, "6.8042531379524426e+140", 10, MPFR_RNDN); mpfr_init(r8393); mpfr_init(r8394); mpfr_init(r8395); mpfr_init(r8396); mpfr_init(r8397); mpfr_init(r8398); mpfr_init(r8399); mpfr_init(r8400); mpfr_init(r8401); mpfr_init(r8402); mpfr_init(r8403); } double f_fm(double x, double l, double t) { mpfr_set_d(r8340, t, MPFR_RNDN); ; mpfr_set_si(r8342, mpfr_cmp(r8340, r8341) <= 0, MPFR_RNDN); ; mpfr_sqrt(r8344, r8343, MPFR_RNDN); mpfr_mul(r8345, r8340, r8344, MPFR_RNDN); mpfr_set_d(r8346, x, MPFR_RNDN); mpfr_mul(r8347, r8346, r8346, MPFR_RNDN); mpfr_mul(r8348, r8344, r8347, MPFR_RNDN); mpfr_div(r8349, r8340, r8348, MPFR_RNDN); ; mpfr_sub(r8351, r8350, r8343, MPFR_RNDN); mpfr_mul(r8352, r8349, r8351, MPFR_RNDN); mpfr_div(r8353, r8343, r8346, MPFR_RNDN); mpfr_div(r8354, r8353, r8344, MPFR_RNDN); mpfr_add(r8355, r8354, r8344, MPFR_RNDN); mpfr_mul(r8356, r8340, r8355, MPFR_RNDN); mpfr_sub(r8357, r8352, r8356, MPFR_RNDN); mpfr_div(r8358, r8345, r8357, MPFR_RNDN); ; mpfr_set_si(r8360, mpfr_cmp(r8340, r8359) <= 0, MPFR_RNDN); mpfr_sqrt(r8361, r8344, MPFR_RNDN); mpfr_mul(r8362, r8340, r8361, MPFR_RNDN); mpfr_mul(r8363, r8362, r8361, MPFR_RNDN); mpfr_set_d(r8364, l, MPFR_RNDN); mpfr_mul(r8365, r8343, r8364, MPFR_RNDN); mpfr_div(r8366, r8364, r8346, MPFR_RNDN); mpfr_mul(r8367, r8365, r8366, MPFR_RNDN); ; mpfr_div(r8369, r8368, r8346, MPFR_RNDN); mpfr_add(r8370, r8343, r8369, MPFR_RNDN); mpfr_mul(r8371, r8340, r8340, MPFR_RNDN); mpfr_mul(r8372, r8370, r8371, MPFR_RNDN); mpfr_add(r8373, r8367, r8372, MPFR_RNDN); mpfr_sqrt(r8374, r8373, MPFR_RNDN); mpfr_div(r8375, r8363, r8374, MPFR_RNDN); ; mpfr_set_si(r8377, mpfr_cmp(r8340, r8376) <= 0, MPFR_RNDN); mpfr_sub(r8378, r8343, r8369, MPFR_RNDN); mpfr_mul(r8379, r8340, r8346, MPFR_RNDN); mpfr_mul(r8380, r8379, r8340, MPFR_RNDN); mpfr_mul(r8381, r8370, r8380, MPFR_RNDN); mpfr_mul(r8382, r8364, r8364, MPFR_RNDN); mpfr_mul(r8383, r8343, r8382, MPFR_RNDN); mpfr_add(r8384, r8381, r8383, MPFR_RNDN); mpfr_mul(r8385, r8378, r8384, MPFR_RNDN); mpfr_sqrt(r8386, r8385, MPFR_RNDN); mpfr_mul(r8387, r8346, r8343, MPFR_RNDN); mpfr_sub(r8388, r8387, r8368, MPFR_RNDN); mpfr_sqrt(r8389, r8388, MPFR_RNDN); mpfr_div(r8390, r8386, r8389, MPFR_RNDN); mpfr_div(r8391, r8345, r8390, MPFR_RNDN); ; mpfr_set_si(r8393, mpfr_cmp(r8340, r8392) <= 0, MPFR_RNDN); mpfr_div(r8394, r8340, r8344, MPFR_RNDN); mpfr_div(r8395, r8394, r8347, MPFR_RNDN); mpfr_sub(r8396, r8343, r8350, MPFR_RNDN); mpfr_mul(r8397, r8395, r8396, MPFR_RNDN); mpfr_add(r8398, r8397, r8356, MPFR_RNDN); mpfr_div(r8399, r8345, r8398, MPFR_RNDN); if (mpfr_get_si(r8393, MPFR_RNDN)) { mpfr_set(r8400, r8375, MPFR_RNDN); } else { mpfr_set(r8400, r8399, MPFR_RNDN); }; if (mpfr_get_si(r8377, MPFR_RNDN)) { mpfr_set(r8401, r8391, MPFR_RNDN); } else { mpfr_set(r8401, r8400, MPFR_RNDN); }; if (mpfr_get_si(r8360, MPFR_RNDN)) { mpfr_set(r8402, r8375, MPFR_RNDN); } else { mpfr_set(r8402, r8401, MPFR_RNDN); }; if (mpfr_get_si(r8342, MPFR_RNDN)) { mpfr_set(r8403, r8358, MPFR_RNDN); } else { mpfr_set(r8403, r8402, MPFR_RNDN); }; return mpfr_get_d(r8403, MPFR_RNDN); } static mpfr_t r8404, r8405, r8406, r8407, r8408, r8409, r8410, r8411, r8412, r8413, r8414, r8415, r8416, r8417, r8418, r8419, r8420, r8421, r8422, r8423, r8424, r8425, r8426, r8427, r8428, r8429, r8430, r8431, r8432, r8433, r8434, r8435, r8436, r8437, r8438, r8439, r8440, r8441, r8442, r8443, r8444, r8445, r8446, r8447, r8448, r8449, r8450, r8451, r8452, r8453, r8454, r8455, r8456, r8457, r8458, r8459, r8460, r8461, r8462, r8463, r8464, r8465, r8466, r8467; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r8404); mpfr_init_set_str(r8405, "-3.135428033554305e+38", 10, MPFR_RNDN); mpfr_init(r8406); mpfr_init_set_str(r8407, "2", 10, MPFR_RNDN); mpfr_init(r8408); mpfr_init(r8409); mpfr_init(r8410); mpfr_init(r8411); mpfr_init(r8412); mpfr_init(r8413); mpfr_init_set_str(r8414, "1", 10, MPFR_RNDN); mpfr_init(r8415); mpfr_init(r8416); mpfr_init(r8417); mpfr_init(r8418); mpfr_init(r8419); mpfr_init(r8420); mpfr_init(r8421); mpfr_init(r8422); mpfr_init_set_str(r8423, "-7.96461310098498e-207", 10, MPFR_RNDN); mpfr_init(r8424); mpfr_init(r8425); mpfr_init(r8426); mpfr_init(r8427); mpfr_init(r8428); mpfr_init(r8429); mpfr_init(r8430); mpfr_init(r8431); mpfr_init_set_str(r8432, "4", 10, MPFR_RNDN); mpfr_init(r8433); mpfr_init(r8434); mpfr_init(r8435); mpfr_init(r8436); mpfr_init(r8437); mpfr_init(r8438); mpfr_init(r8439); mpfr_init_set_str(r8440, "2.986334750317526e-159", 10, MPFR_RNDN); mpfr_init(r8441); mpfr_init(r8442); mpfr_init(r8443); mpfr_init(r8444); mpfr_init(r8445); mpfr_init(r8446); mpfr_init(r8447); mpfr_init(r8448); mpfr_init(r8449); mpfr_init(r8450); mpfr_init(r8451); mpfr_init(r8452); mpfr_init(r8453); mpfr_init(r8454); mpfr_init(r8455); mpfr_init_set_str(r8456, "6.8042531379524426e+140", 10, MPFR_RNDN); mpfr_init(r8457); mpfr_init(r8458); mpfr_init(r8459); mpfr_init(r8460); mpfr_init(r8461); mpfr_init(r8462); mpfr_init(r8463); mpfr_init(r8464); mpfr_init(r8465); mpfr_init(r8466); mpfr_init(r8467); } double f_dm(double x, double l, double t) { mpfr_set_d(r8404, t, MPFR_RNDN); ; mpfr_set_si(r8406, mpfr_cmp(r8404, r8405) <= 0, MPFR_RNDN); ; mpfr_sqrt(r8408, r8407, MPFR_RNDN); mpfr_mul(r8409, r8404, r8408, MPFR_RNDN); mpfr_set_d(r8410, x, MPFR_RNDN); mpfr_mul(r8411, r8410, r8410, MPFR_RNDN); mpfr_mul(r8412, r8408, r8411, MPFR_RNDN); mpfr_div(r8413, r8404, r8412, MPFR_RNDN); ; mpfr_sub(r8415, r8414, r8407, MPFR_RNDN); mpfr_mul(r8416, r8413, r8415, MPFR_RNDN); mpfr_div(r8417, r8407, r8410, MPFR_RNDN); mpfr_div(r8418, r8417, r8408, MPFR_RNDN); mpfr_add(r8419, r8418, r8408, MPFR_RNDN); mpfr_mul(r8420, r8404, r8419, MPFR_RNDN); mpfr_sub(r8421, r8416, r8420, MPFR_RNDN); mpfr_div(r8422, r8409, r8421, MPFR_RNDN); ; mpfr_set_si(r8424, mpfr_cmp(r8404, r8423) <= 0, MPFR_RNDN); mpfr_sqrt(r8425, r8408, MPFR_RNDN); mpfr_mul(r8426, r8404, r8425, MPFR_RNDN); mpfr_mul(r8427, r8426, r8425, MPFR_RNDN); mpfr_set_d(r8428, l, MPFR_RNDN); mpfr_mul(r8429, r8407, r8428, MPFR_RNDN); mpfr_div(r8430, r8428, r8410, MPFR_RNDN); mpfr_mul(r8431, r8429, r8430, MPFR_RNDN); ; mpfr_div(r8433, r8432, r8410, MPFR_RNDN); mpfr_add(r8434, r8407, r8433, MPFR_RNDN); mpfr_mul(r8435, r8404, r8404, MPFR_RNDN); mpfr_mul(r8436, r8434, r8435, MPFR_RNDN); mpfr_add(r8437, r8431, r8436, MPFR_RNDN); mpfr_sqrt(r8438, r8437, MPFR_RNDN); mpfr_div(r8439, r8427, r8438, MPFR_RNDN); ; mpfr_set_si(r8441, mpfr_cmp(r8404, r8440) <= 0, MPFR_RNDN); mpfr_sub(r8442, r8407, r8433, MPFR_RNDN); mpfr_mul(r8443, r8404, r8410, MPFR_RNDN); mpfr_mul(r8444, r8443, r8404, MPFR_RNDN); mpfr_mul(r8445, r8434, r8444, MPFR_RNDN); mpfr_mul(r8446, r8428, r8428, MPFR_RNDN); mpfr_mul(r8447, r8407, r8446, MPFR_RNDN); mpfr_add(r8448, r8445, r8447, MPFR_RNDN); mpfr_mul(r8449, r8442, r8448, MPFR_RNDN); mpfr_sqrt(r8450, r8449, MPFR_RNDN); mpfr_mul(r8451, r8410, r8407, MPFR_RNDN); mpfr_sub(r8452, r8451, r8432, MPFR_RNDN); mpfr_sqrt(r8453, r8452, MPFR_RNDN); mpfr_div(r8454, r8450, r8453, MPFR_RNDN); mpfr_div(r8455, r8409, r8454, MPFR_RNDN); ; mpfr_set_si(r8457, mpfr_cmp(r8404, r8456) <= 0, MPFR_RNDN); mpfr_div(r8458, r8404, r8408, MPFR_RNDN); mpfr_div(r8459, r8458, r8411, MPFR_RNDN); mpfr_sub(r8460, r8407, r8414, MPFR_RNDN); mpfr_mul(r8461, r8459, r8460, MPFR_RNDN); mpfr_add(r8462, r8461, r8420, MPFR_RNDN); mpfr_div(r8463, r8409, r8462, MPFR_RNDN); if (mpfr_get_si(r8457, MPFR_RNDN)) { mpfr_set(r8464, r8439, MPFR_RNDN); } else { mpfr_set(r8464, r8463, MPFR_RNDN); }; if (mpfr_get_si(r8441, MPFR_RNDN)) { mpfr_set(r8465, r8455, MPFR_RNDN); } else { mpfr_set(r8465, r8464, MPFR_RNDN); }; if (mpfr_get_si(r8424, MPFR_RNDN)) { mpfr_set(r8466, r8439, MPFR_RNDN); } else { mpfr_set(r8466, r8465, MPFR_RNDN); }; if (mpfr_get_si(r8406, MPFR_RNDN)) { mpfr_set(r8467, r8422, MPFR_RNDN); } else { mpfr_set(r8467, r8466, MPFR_RNDN); }; return mpfr_get_d(r8467, MPFR_RNDN); }