Exp_msk1
word0(da) += (k >> 2)*Exp_msk1;
word0(db) += (k >> 2)*Exp_msk1;
word0(da) += k*Exp_msk1;
word0(db) += k*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;
+ Exp_msk1
| Exp_msk1 >> 4
if (L < Exp_msk1)
if (L <= Exp_msk1)
L -= 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 ((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)
word0(rv) -= P*Exp_msk1;
if (y <= (P-1)*Exp_msk1 && aadj >= 1.) {
word0(d2) -= 31*Exp_msk1; /* adjust exponent */
word0(eps) -= (P-1)*Exp_msk1;
#define Big0 (Frac_mask1 | Exp_msk1*(DBL_MAX_EXP+Bias-1))
L = (word0(x) & Exp_mask) - (P-1)*Exp_msk1;
L |= Exp_msk1 >> 4;
z |= Exp_msk1;