Lines Matching refs:rv
187 double rv; in sulp() local
190 rv = ulp(dval(x)); in sulp()
192 return rv; /* Is there an example where i <= 0 ? */ in sulp()
197 return rv * u.d; in sulp()
252 U aadj1, rv, rv0; local
277 dval(rv) = 0.;
334 ULtod(rv.i, bits, exp, i);
337 if ((dword0(rv)&Exp_mask) == 0)
445 dword0(rv) = 0x7ff00000;
447 dword1(rv) = 0;
459 dword0(rv) = 0x7ff00000 | bits[1];
461 dword1(rv) = bits[0];
466 dval(rv) = nan ("");
490 dval(rv) = y;
496 dval(rv) = tens[k - 9] * dval(rv) + z;
516 dval(rv) = -dval(rv);
520 /* rv = */ rounded_product(dval(rv), tens[e]);
532 dval(rv) = -dval(rv);
537 dval(rv) *= tens[i];
543 dword0(rv) -= P*Exp_msk1;
544 /* rv = */ rounded_product(dval(rv), tens[e]);
545 if ((dword0(rv) & Exp_mask)
548 dword0(rv) += P*Exp_msk1;
550 /* rv = */ rounded_product(dval(rv), tens[e]);
560 dval(rv) = -dval(rv);
564 /* rv = */ rounded_quotient(dval(rv), tens[-e]);
595 dval(rv) *= tens[i];
608 dword0(rv) = Big0;
610 dword1(rv) = Big1;
614 dword0(rv) = Exp_mask;
616 dword1(rv) = 0;
620 dword0(rv) = Exp_mask;
622 dword1(rv) = 0;
631 dword0(rv) = Big0;
633 dword1(rv) = Big1;
643 dval(rv) *= bigtens[j];
645 dword0(rv) -= P*Exp_msk1;
646 dval(rv) *= bigtens[j];
647 if ((z = dword0(rv) & Exp_mask)
653 dword0(rv) = Big0;
655 dword1(rv) = Big1;
659 dword0(rv) += P*Exp_msk1;
665 dval(rv) /= tens[i];
674 dval(rv) *= tinytens[j];
675 if (scale && (j = 2*P + 1 - ((dword0(rv) & Exp_mask)
680 dword1(rv) = 0;
683 dword0(rv) = (P+2)*Exp_msk1;
685 dword0(rv) &= 0xffffffff << (j-32);
689 dword1(rv) &= 0xffffffff << j;
695 dval(rv) *= tinytens[j];
697 dval(rv0) = dval(rv);
698 dval(rv) *= tinytens[j];
699 if (!dval(rv)) {
700 dval(rv) = 2.*dval(rv0);
701 dval(rv) *= tinytens[j];
703 if (!dval(rv)) {
705 dval(rv) = 0.;
715 dword0(rv) = Tiny0;
716 dword1(rv) = Tiny1;
718 dword0(rv) = Tiny1;
741 bb = d2b(dval(rv), &bbe, &bbbits); /* rv = bb * 2^bbe */
863 if (!dword1(rv)
864 && !(dword0(rv) & Frac_mask)) {
865 y = dword0(rv) & Exp_mask;
879 if (scale && (y = dword0(rv) & Exp_mask)
884 if ((dword0(rv) & Exp_mask) <=
886 dword0(rv) += P*Exp_msk1;
887 dval(rv) += adj*ulp(dval(rv));
888 dword0(rv) -= P*Exp_msk1;
893 dval(rv) += adj*ulp(dval(rv));
910 if (scale && (y = dword0(rv) & Exp_mask) <= 2*P*Exp_msk1)
914 if ((dword0(rv) & Exp_mask) <= P*Exp_msk1) {
915 dword0(rv) += P*Exp_msk1;
916 adj *= ulp(dval(rv));
918 dval(rv) += adj;
920 dval(rv) -= adj;
921 dword0(rv) -= P*Exp_msk1;
926 adj *= ulp(dval(rv));
928 if (dword0(rv) == Big0 && dword1(rv) == Big1)
930 dval(rv) += adj;
932 dval(rv) -= adj;
941 if (dsign || dword1(rv) || dword0(rv) & Bndry_mask
944 || (dword0(rv) & Exp_mask) <= (2*P+1)*Exp_msk1
946 || (dword0(rv) & Exp_mask) <= Exp_msk1
971 if ((dword0(rv) & Bndry_mask1) == Bndry_mask1
972 && dword1(rv) == (
974 (scale && (y = dword0(rv) & Exp_mask) <= 2*P*Exp_msk1)
979 if (dword0(rv) == Big0 && dword1(rv) == Big1)
981 dword0(rv) = (dword0(rv) & Exp_mask)
988 dword1(rv) = 0;
996 else if (!(dword0(rv) & Bndry_mask) && !dword1(rv)) {
1000 L = dword0(rv) & Exp_mask;
1015 L = dword0(rv) & Exp_mask;
1026 L = (dword0(rv) & Exp_mask) - Exp_msk1;
1028 dword0(rv) = L | Bndry_mask1;
1030 dword1(rv) = 0xffffffff;
1041 if (!(dword0(rv) & Lsb1))
1044 else if (!(dword1(rv) & Lsb))
1047 if (!(dword1(rv) & LSB))
1053 dval(rv) += sulp(rv, scale);
1055 dval(rv) += ulp(dval(rv));
1060 dval(rv) -= sulp(rv, scale);
1062 dval(rv) -= ulp(dval(rv));
1065 if (!dval(rv))
1078 else if (dword1(rv) || dword0(rv) & Bndry_mask) {
1080 if (dword1(rv) == Tiny1 && !dword0(rv))
1114 y = dword0(rv) & Exp_mask;
1119 dval(rv0) = dval(rv);
1120 dword0(rv) -= P*Exp_msk1;
1121 adj = dval(aadj1) * ulp(dval(rv));
1122 dval(rv) += adj;
1123 if ((dword0(rv) & Exp_mask) >=
1127 dword0(rv) = Big0;
1129 dword1(rv) = Big1;
1134 dword0(rv) += P*Exp_msk1;
1147 adj = dval(aadj1) * ulp(dval(rv));
1148 dval(rv) += adj;
1151 if ((dword0(rv) & Exp_mask) <= P*Exp_msk1) {
1152 dval(rv0) = dval(rv);
1153 dword0(rv) += P*Exp_msk1;
1154 adj = dval(aadj1) * ulp(dval(rv));
1155 dval(rv) += adj;
1157 if ((dword0(rv) & Exp_mask) < P*Exp_msk1)
1159 if ((dword0(rv) & Exp_mask) <= P*Exp_msk1)
1166 dword0(rv) = Tiny0;
1167 dword1(rv) = Tiny1;
1169 dword0(rv) = Tiny1;
1174 dword0(rv) -= P*Exp_msk1;
1177 adj = dval(aadj1) * ulp(dval(rv));
1178 dval(rv) += adj;
1193 adj = dval(aadj1) * ulp(dval(rv));
1194 dval(rv) += adj;
1198 z = dword0(rv) & Exp_mask;
1217 if (dsign || dword1(rv) || dword0(rv) & Bndry_mask) {
1250 dval(rv) *= dval(rv0);
1253 if ((dword0(rv) & Exp_mask) == 0)
1259 if (inexact && !(dword0(rv) & Exp_mask)) {
1274 return sign ? -dval(rv) : dval(rv);