Exp_msk1
word0(&d2) -= 31*Exp_msk1; /* adjust exponent */
word0(&eps) -= (P-1)*Exp_msk1;
&& word0(&d) & (Exp_mask & ~Exp_msk1)
word0(&eps) -= (P-1)*Exp_msk1;
#define Big0 (Frac_mask1 | Exp_msk1*(DBL_MAX_EXP+Bias-1))
z |= Exp_msk1;
word0(&db) += (k >> 2)*Exp_msk1;
word0(&da) += k*Exp_msk1;
word0(&db) += k*Exp_msk1;
word0(&da) += (k >> 2)*Exp_msk1;
word0(&rv) -= P*Exp_msk1;
if (y <= (P-1)*Exp_msk1 && aadj > 1.) {
word0(&rv0) = Exp_1 - 2*P*Exp_msk1;
word0(&rv) -= P*Exp_msk1;
> Exp_msk1*(DBL_MAX_EXP+Bias-1-P))
word0(&rv) += P*Exp_msk1;
word0(&rv) -= P*Exp_msk1;
> Exp_msk1*(DBL_MAX_EXP+Bias-P))
if (z > Exp_msk1*(DBL_MAX_EXP+Bias-1-P)) {
word0(&rv) += P*Exp_msk1;
word0(&rv) = (P+2)*Exp_msk1;
if (!scale || y > 2*P*Exp_msk1)
<= 2*P*Exp_msk1)
word0(&adj) += (2*P+1)*Exp_msk1 - y;
P*Exp_msk1) {
word0(&rv) += P*Exp_msk1;
word0(&rv) -= P*Exp_msk1;
if (scale && (y = word0(&rv) & Exp_mask) <= 2*P*Exp_msk1)
word0(&adj) += (2*P+1)*Exp_msk1 - y;
if ((word0(&rv) & Exp_mask) <= P*Exp_msk1) {
word0(&rv) += P*Exp_msk1;
word0(&rv) -= P*Exp_msk1;
|| (word0(&rv) & Exp_mask) <= (2*P+1)*Exp_msk1
|| (word0(&rv) & Exp_mask) <= Exp_msk1
(scale && (y = word0(&rv) & Exp_mask) <= 2*P*Exp_msk1)
+ Exp_msk1
| Exp_msk1 >> 4
if (L < Exp_msk1)
if (L <= (scale ? (2*P+1)*Exp_msk1 : Exp_msk1))
if (L <= Exp_msk1)
L -= Exp_msk1;
if (L <= (2*P+1)*Exp_msk1) {
if (L > (P+2)*Exp_msk1)
L = (word0(&rv) & Exp_mask) - Exp_msk1;
if (y == Exp_msk1*(DBL_MAX_EXP+Bias-1)) {
word0(&rv) -= P*Exp_msk1;
Exp_msk1*(DBL_MAX_EXP+Bias-P)) {
word0(&rv) += P*Exp_msk1;
if (scale && y <= 2*P*Exp_msk1) {
word0(&aadj1) += (2*P+1)*Exp_msk1 - y;
if ((word0(&rv) & Exp_mask) <= P*Exp_msk1) {
word0(&rv) += P*Exp_msk1;
if ((word0(&rv) & Exp_mask) < P*Exp_msk1)
if ((word0(&rv) & Exp_mask) <= P*Exp_msk1)
L = word0(d) | Exp_msk1;
word0(&adj) -= P*Exp_msk1;
> Exp_msk1*(DBL_MAX_EXP+Bias-1-P))
word0(&adj) += P*Exp_msk1;
L = (word0(x) & Exp_mask) - (P-1)*Exp_msk1;
L |= Exp_msk1 >> 4;