#include #include #include #include #include char *name = "From Warwick Tucker's Validated Numerics"; double f_if() { float r16280058 = 333.75; float r16280059 = 33096; float r16280060 = 6; float r16280061 = pow(r16280059, r16280060); float r16280062 = r16280058 * r16280061; float r16280063 = 77617; float r16280064 = r16280063 * r16280063; float r16280065 = 11; float r16280066 = r16280065 * r16280064; float r16280067 = r16280059 * r16280059; float r16280068 = r16280066 * r16280067; float r16280069 = -r16280061; float r16280070 = r16280068 + r16280069; float r16280071 = -121; float r16280072 = 4; float r16280073 = pow(r16280059, r16280072); float r16280074 = r16280071 * r16280073; float r16280075 = r16280070 + r16280074; float r16280076 = -2; float r16280077 = r16280075 + r16280076; float r16280078 = r16280064 * r16280077; float r16280079 = r16280062 + r16280078; float r16280080 = 5.5; float r16280081 = 8; float r16280082 = pow(r16280059, r16280081); float r16280083 = r16280080 * r16280082; float r16280084 = r16280079 + r16280083; float r16280085 = 2; float r16280086 = r16280085 * r16280059; float r16280087 = r16280063 / r16280086; float r16280088 = r16280084 + r16280087; return r16280088; } double f_id() { double r16280089 = 333.75; double r16280090 = 33096; double r16280091 = 6; double r16280092 = pow(r16280090, r16280091); double r16280093 = r16280089 * r16280092; double r16280094 = 77617; double r16280095 = r16280094 * r16280094; double r16280096 = 11; double r16280097 = r16280096 * r16280095; double r16280098 = r16280090 * r16280090; double r16280099 = r16280097 * r16280098; double r16280100 = -r16280092; double r16280101 = r16280099 + r16280100; double r16280102 = -121; double r16280103 = 4; double r16280104 = pow(r16280090, r16280103); double r16280105 = r16280102 * r16280104; double r16280106 = r16280101 + r16280105; double r16280107 = -2; double r16280108 = r16280106 + r16280107; double r16280109 = r16280095 * r16280108; double r16280110 = r16280093 + r16280109; double r16280111 = 5.5; double r16280112 = 8; double r16280113 = pow(r16280090, r16280112); double r16280114 = r16280111 * r16280113; double r16280115 = r16280110 + r16280114; double r16280116 = 2; double r16280117 = r16280116 * r16280090; double r16280118 = r16280094 / r16280117; double r16280119 = r16280115 + r16280118; return r16280119; } double f_of() { float r16280120 = 333.75; float r16280121 = 33096; float r16280122 = 6; float r16280123 = pow(r16280121, r16280122); float r16280124 = r16280120 * r16280123; float r16280125 = 77617; float r16280126 = r16280125 * r16280125; float r16280127 = 11; float r16280128 = r16280127 * r16280126; float r16280129 = r16280121 * r16280121; float r16280130 = r16280128 * r16280129; float r16280131 = -r16280123; float r16280132 = r16280130 + r16280131; float r16280133 = -121; float r16280134 = 4; float r16280135 = pow(r16280121, r16280134); float r16280136 = r16280133 * r16280135; float r16280137 = r16280132 + r16280136; float r16280138 = -2; float r16280139 = r16280137 + r16280138; float r16280140 = r16280126 * r16280139; float r16280141 = r16280124 + r16280140; float r16280142 = r16280141 * r16280141; float r16280143 = 5.5; float r16280144 = 8; float r16280145 = pow(r16280121, r16280144); float r16280146 = r16280143 * r16280145; float r16280147 = r16280141 - r16280146; float r16280148 = r16280142 / r16280147; float r16280149 = r16280148 * r16280148; float r16280150 = r16280146 * r16280146; float r16280151 = r16280150 / r16280147; float r16280152 = 2; float r16280153 = r16280152 * r16280121; float r16280154 = r16280125 / r16280153; float r16280155 = r16280151 - r16280154; float r16280156 = r16280155 * r16280155; float r16280157 = r16280149 - r16280156; float r16280158 = r16280142 - r16280150; float r16280159 = r16280150 / r16280158; float r16280160 = r16280141 + r16280146; float r16280161 = r16280159 * r16280160; float r16280162 = r16280161 - r16280154; float r16280163 = r16280148 + r16280162; float r16280164 = r16280157 / r16280163; return r16280164; } double f_od() { double r16280165 = 333.75; double r16280166 = 33096; double r16280167 = 6; double r16280168 = pow(r16280166, r16280167); double r16280169 = r16280165 * r16280168; double r16280170 = 77617; double r16280171 = r16280170 * r16280170; double r16280172 = 11; double r16280173 = r16280172 * r16280171; double r16280174 = r16280166 * r16280166; double r16280175 = r16280173 * r16280174; double r16280176 = -r16280168; double r16280177 = r16280175 + r16280176; double r16280178 = -121; double r16280179 = 4; double r16280180 = pow(r16280166, r16280179); double r16280181 = r16280178 * r16280180; double r16280182 = r16280177 + r16280181; double r16280183 = -2; double r16280184 = r16280182 + r16280183; double r16280185 = r16280171 * r16280184; double r16280186 = r16280169 + r16280185; double r16280187 = r16280186 * r16280186; double r16280188 = 5.5; double r16280189 = 8; double r16280190 = pow(r16280166, r16280189); double r16280191 = r16280188 * r16280190; double r16280192 = r16280186 - r16280191; double r16280193 = r16280187 / r16280192; double r16280194 = r16280193 * r16280193; double r16280195 = r16280191 * r16280191; double r16280196 = r16280195 / r16280192; double r16280197 = 2; double r16280198 = r16280197 * r16280166; double r16280199 = r16280170 / r16280198; double r16280200 = r16280196 - r16280199; double r16280201 = r16280200 * r16280200; double r16280202 = r16280194 - r16280201; double r16280203 = r16280187 - r16280195; double r16280204 = r16280195 / r16280203; double r16280205 = r16280186 + r16280191; double r16280206 = r16280204 * r16280205; double r16280207 = r16280206 - r16280199; double r16280208 = r16280193 + r16280207; double r16280209 = r16280202 / r16280208; return r16280209; } 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 r16280210, r16280211, r16280212, r16280213, r16280214, r16280215, r16280216, r16280217, r16280218, r16280219, r16280220, r16280221, r16280222, r16280223, r16280224, r16280225, r16280226, r16280227, r16280228, r16280229, r16280230, r16280231, r16280232, r16280233, r16280234, r16280235, r16280236, r16280237, r16280238, r16280239, r16280240; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r16280210, "333.75", 10, MPFR_RNDN); mpfr_init_set_str(r16280211, "33096", 10, MPFR_RNDN); mpfr_init_set_str(r16280212, "6", 10, MPFR_RNDN); mpfr_init(r16280213); mpfr_init(r16280214); mpfr_init_set_str(r16280215, "77617", 10, MPFR_RNDN); mpfr_init(r16280216); mpfr_init_set_str(r16280217, "11", 10, MPFR_RNDN); mpfr_init(r16280218); mpfr_init(r16280219); mpfr_init(r16280220); mpfr_init(r16280221); mpfr_init(r16280222); mpfr_init_set_str(r16280223, "-121", 10, MPFR_RNDN); mpfr_init_set_str(r16280224, "4", 10, MPFR_RNDN); mpfr_init(r16280225); mpfr_init(r16280226); mpfr_init(r16280227); mpfr_init_set_str(r16280228, "-2", 10, MPFR_RNDN); mpfr_init(r16280229); mpfr_init(r16280230); mpfr_init(r16280231); mpfr_init_set_str(r16280232, "5.5", 10, MPFR_RNDN); mpfr_init_set_str(r16280233, "8", 10, MPFR_RNDN); mpfr_init(r16280234); mpfr_init(r16280235); mpfr_init(r16280236); mpfr_init_set_str(r16280237, "2", 10, MPFR_RNDN); mpfr_init(r16280238); mpfr_init(r16280239); mpfr_init(r16280240); } double f_im() { ; ; ; mpfr_pow(r16280213, r16280211, r16280212, MPFR_RNDN); mpfr_mul(r16280214, r16280210, r16280213, MPFR_RNDN); ; mpfr_mul(r16280216, r16280215, r16280215, MPFR_RNDN); ; mpfr_mul(r16280218, r16280217, r16280216, MPFR_RNDN); mpfr_mul(r16280219, r16280211, r16280211, MPFR_RNDN); mpfr_mul(r16280220, r16280218, r16280219, MPFR_RNDN); mpfr_neg(r16280221, r16280213, MPFR_RNDN); mpfr_add(r16280222, r16280220, r16280221, MPFR_RNDN); ; ; mpfr_pow(r16280225, r16280211, r16280224, MPFR_RNDN); mpfr_mul(r16280226, r16280223, r16280225, MPFR_RNDN); mpfr_add(r16280227, r16280222, r16280226, MPFR_RNDN); ; mpfr_add(r16280229, r16280227, r16280228, MPFR_RNDN); mpfr_mul(r16280230, r16280216, r16280229, MPFR_RNDN); mpfr_add(r16280231, r16280214, r16280230, MPFR_RNDN); ; ; mpfr_pow(r16280234, r16280211, r16280233, MPFR_RNDN); mpfr_mul(r16280235, r16280232, r16280234, MPFR_RNDN); mpfr_add(r16280236, r16280231, r16280235, MPFR_RNDN); ; mpfr_mul(r16280238, r16280237, r16280211, MPFR_RNDN); mpfr_div(r16280239, r16280215, r16280238, MPFR_RNDN); mpfr_add(r16280240, r16280236, r16280239, MPFR_RNDN); return mpfr_get_d(r16280240, MPFR_RNDN); } static mpfr_t r16280241, r16280242, r16280243, r16280244, r16280245, r16280246, r16280247, r16280248, r16280249, r16280250, r16280251, r16280252, r16280253, r16280254, r16280255, r16280256, r16280257, r16280258, r16280259, r16280260, r16280261, r16280262, r16280263, r16280264, r16280265, r16280266, r16280267, r16280268, r16280269, r16280270, r16280271, r16280272, r16280273, r16280274, r16280275, r16280276, r16280277, r16280278, r16280279, r16280280, r16280281, r16280282, r16280283, r16280284, r16280285; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r16280241, "333.75", 10, MPFR_RNDN); mpfr_init_set_str(r16280242, "33096", 10, MPFR_RNDN); mpfr_init_set_str(r16280243, "6", 10, MPFR_RNDN); mpfr_init(r16280244); mpfr_init(r16280245); mpfr_init_set_str(r16280246, "77617", 10, MPFR_RNDN); mpfr_init(r16280247); mpfr_init_set_str(r16280248, "11", 10, MPFR_RNDN); mpfr_init(r16280249); mpfr_init(r16280250); mpfr_init(r16280251); mpfr_init(r16280252); mpfr_init(r16280253); mpfr_init_set_str(r16280254, "-121", 10, MPFR_RNDN); mpfr_init_set_str(r16280255, "4", 10, MPFR_RNDN); mpfr_init(r16280256); mpfr_init(r16280257); mpfr_init(r16280258); mpfr_init_set_str(r16280259, "-2", 10, MPFR_RNDN); mpfr_init(r16280260); mpfr_init(r16280261); mpfr_init(r16280262); mpfr_init(r16280263); mpfr_init_set_str(r16280264, "5.5", 10, MPFR_RNDN); mpfr_init_set_str(r16280265, "8", 10, MPFR_RNDN); mpfr_init(r16280266); mpfr_init(r16280267); mpfr_init(r16280268); mpfr_init(r16280269); mpfr_init(r16280270); mpfr_init(r16280271); mpfr_init(r16280272); mpfr_init_set_str(r16280273, "2", 10, MPFR_RNDN); mpfr_init(r16280274); mpfr_init(r16280275); mpfr_init(r16280276); mpfr_init(r16280277); mpfr_init(r16280278); mpfr_init(r16280279); mpfr_init(r16280280); mpfr_init(r16280281); mpfr_init(r16280282); mpfr_init(r16280283); mpfr_init(r16280284); mpfr_init(r16280285); } double f_fm() { ; ; ; mpfr_pow(r16280244, r16280242, r16280243, MPFR_RNDN); mpfr_mul(r16280245, r16280241, r16280244, MPFR_RNDN); ; mpfr_mul(r16280247, r16280246, r16280246, MPFR_RNDN); ; mpfr_mul(r16280249, r16280248, r16280247, MPFR_RNDN); mpfr_mul(r16280250, r16280242, r16280242, MPFR_RNDN); mpfr_mul(r16280251, r16280249, r16280250, MPFR_RNDN); mpfr_neg(r16280252, r16280244, MPFR_RNDN); mpfr_add(r16280253, r16280251, r16280252, MPFR_RNDN); ; ; mpfr_pow(r16280256, r16280242, r16280255, MPFR_RNDN); mpfr_mul(r16280257, r16280254, r16280256, MPFR_RNDN); mpfr_add(r16280258, r16280253, r16280257, MPFR_RNDN); ; mpfr_add(r16280260, r16280258, r16280259, MPFR_RNDN); mpfr_mul(r16280261, r16280247, r16280260, MPFR_RNDN); mpfr_add(r16280262, r16280245, r16280261, MPFR_RNDN); mpfr_mul(r16280263, r16280262, r16280262, MPFR_RNDN); ; ; mpfr_pow(r16280266, r16280242, r16280265, MPFR_RNDN); mpfr_mul(r16280267, r16280264, r16280266, MPFR_RNDN); mpfr_sub(r16280268, r16280262, r16280267, MPFR_RNDN); mpfr_div(r16280269, r16280263, r16280268, MPFR_RNDN); mpfr_mul(r16280270, r16280269, r16280269, MPFR_RNDN); mpfr_mul(r16280271, r16280267, r16280267, MPFR_RNDN); mpfr_div(r16280272, r16280271, r16280268, MPFR_RNDN); ; mpfr_mul(r16280274, r16280273, r16280242, MPFR_RNDN); mpfr_div(r16280275, r16280246, r16280274, MPFR_RNDN); mpfr_sub(r16280276, r16280272, r16280275, MPFR_RNDN); mpfr_mul(r16280277, r16280276, r16280276, MPFR_RNDN); mpfr_sub(r16280278, r16280270, r16280277, MPFR_RNDN); mpfr_sub(r16280279, r16280263, r16280271, MPFR_RNDN); mpfr_div(r16280280, r16280271, r16280279, MPFR_RNDN); mpfr_add(r16280281, r16280262, r16280267, MPFR_RNDN); mpfr_mul(r16280282, r16280280, r16280281, MPFR_RNDN); mpfr_sub(r16280283, r16280282, r16280275, MPFR_RNDN); mpfr_add(r16280284, r16280269, r16280283, MPFR_RNDN); mpfr_div(r16280285, r16280278, r16280284, MPFR_RNDN); return mpfr_get_d(r16280285, MPFR_RNDN); } static mpfr_t r16280286, r16280287, r16280288, r16280289, r16280290, r16280291, r16280292, r16280293, r16280294, r16280295, r16280296, r16280297, r16280298, r16280299, r16280300, r16280301, r16280302, r16280303, r16280304, r16280305, r16280306, r16280307, r16280308, r16280309, r16280310, r16280311, r16280312, r16280313, r16280314, r16280315, r16280316, r16280317, r16280318, r16280319, r16280320, r16280321, r16280322, r16280323, r16280324, r16280325, r16280326, r16280327, r16280328, r16280329, r16280330; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r16280286, "333.75", 10, MPFR_RNDN); mpfr_init_set_str(r16280287, "33096", 10, MPFR_RNDN); mpfr_init_set_str(r16280288, "6", 10, MPFR_RNDN); mpfr_init(r16280289); mpfr_init(r16280290); mpfr_init_set_str(r16280291, "77617", 10, MPFR_RNDN); mpfr_init(r16280292); mpfr_init_set_str(r16280293, "11", 10, MPFR_RNDN); mpfr_init(r16280294); mpfr_init(r16280295); mpfr_init(r16280296); mpfr_init(r16280297); mpfr_init(r16280298); mpfr_init_set_str(r16280299, "-121", 10, MPFR_RNDN); mpfr_init_set_str(r16280300, "4", 10, MPFR_RNDN); mpfr_init(r16280301); mpfr_init(r16280302); mpfr_init(r16280303); mpfr_init_set_str(r16280304, "-2", 10, MPFR_RNDN); mpfr_init(r16280305); mpfr_init(r16280306); mpfr_init(r16280307); mpfr_init(r16280308); mpfr_init_set_str(r16280309, "5.5", 10, MPFR_RNDN); mpfr_init_set_str(r16280310, "8", 10, MPFR_RNDN); mpfr_init(r16280311); mpfr_init(r16280312); mpfr_init(r16280313); mpfr_init(r16280314); mpfr_init(r16280315); mpfr_init(r16280316); mpfr_init(r16280317); mpfr_init_set_str(r16280318, "2", 10, MPFR_RNDN); mpfr_init(r16280319); mpfr_init(r16280320); mpfr_init(r16280321); mpfr_init(r16280322); mpfr_init(r16280323); mpfr_init(r16280324); mpfr_init(r16280325); mpfr_init(r16280326); mpfr_init(r16280327); mpfr_init(r16280328); mpfr_init(r16280329); mpfr_init(r16280330); } double f_dm() { ; ; ; mpfr_pow(r16280289, r16280287, r16280288, MPFR_RNDN); mpfr_mul(r16280290, r16280286, r16280289, MPFR_RNDN); ; mpfr_mul(r16280292, r16280291, r16280291, MPFR_RNDN); ; mpfr_mul(r16280294, r16280293, r16280292, MPFR_RNDN); mpfr_mul(r16280295, r16280287, r16280287, MPFR_RNDN); mpfr_mul(r16280296, r16280294, r16280295, MPFR_RNDN); mpfr_neg(r16280297, r16280289, MPFR_RNDN); mpfr_add(r16280298, r16280296, r16280297, MPFR_RNDN); ; ; mpfr_pow(r16280301, r16280287, r16280300, MPFR_RNDN); mpfr_mul(r16280302, r16280299, r16280301, MPFR_RNDN); mpfr_add(r16280303, r16280298, r16280302, MPFR_RNDN); ; mpfr_add(r16280305, r16280303, r16280304, MPFR_RNDN); mpfr_mul(r16280306, r16280292, r16280305, MPFR_RNDN); mpfr_add(r16280307, r16280290, r16280306, MPFR_RNDN); mpfr_mul(r16280308, r16280307, r16280307, MPFR_RNDN); ; ; mpfr_pow(r16280311, r16280287, r16280310, MPFR_RNDN); mpfr_mul(r16280312, r16280309, r16280311, MPFR_RNDN); mpfr_sub(r16280313, r16280307, r16280312, MPFR_RNDN); mpfr_div(r16280314, r16280308, r16280313, MPFR_RNDN); mpfr_mul(r16280315, r16280314, r16280314, MPFR_RNDN); mpfr_mul(r16280316, r16280312, r16280312, MPFR_RNDN); mpfr_div(r16280317, r16280316, r16280313, MPFR_RNDN); ; mpfr_mul(r16280319, r16280318, r16280287, MPFR_RNDN); mpfr_div(r16280320, r16280291, r16280319, MPFR_RNDN); mpfr_sub(r16280321, r16280317, r16280320, MPFR_RNDN); mpfr_mul(r16280322, r16280321, r16280321, MPFR_RNDN); mpfr_sub(r16280323, r16280315, r16280322, MPFR_RNDN); mpfr_sub(r16280324, r16280308, r16280316, MPFR_RNDN); mpfr_div(r16280325, r16280316, r16280324, MPFR_RNDN); mpfr_add(r16280326, r16280307, r16280312, MPFR_RNDN); mpfr_mul(r16280327, r16280325, r16280326, MPFR_RNDN); mpfr_sub(r16280328, r16280327, r16280320, MPFR_RNDN); mpfr_add(r16280329, r16280314, r16280328, MPFR_RNDN); mpfr_div(r16280330, r16280323, r16280329, MPFR_RNDN); return mpfr_get_d(r16280330, MPFR_RNDN); }