dval 1593 util.c dval(x) = x_; dval 1623 util.c return dval(a); dval 1688 util.c return dval(d); dval 1704 util.c dval(d) = d_; dval 1834 util.c dval(da) = b2d(a, &ka); dval 1835 util.c dval(db) = b2d(b, &kb); dval 1845 util.c dval(da) *= 1 << k; dval 1851 util.c dval(db) *= 1 << k; dval 1861 util.c return dval(da) / dval(db); dval 2007 util.c dval(rv) = 0.; dval 2095 util.c dval(rv) = ldexp(adj, nd0); dval 2247 util.c dval(rv) = y; dval 2253 util.c dval(rv) = tens[k - 9] * dval(rv) + z; dval 2273 util.c dval(rv) = -dval(rv); dval 2277 util.c /* rv = */ rounded_product(dval(rv), tens[e]); dval 2289 util.c dval(rv) = -dval(rv); dval 2294 util.c dval(rv) *= tens[i]; dval 2301 util.c /* rv = */ rounded_product(dval(rv), tens[e]); dval 2307 util.c /* rv = */ rounded_product(dval(rv), tens[e]); dval 2317 util.c dval(rv) = -dval(rv); dval 2321 util.c /* rv = */ rounded_quotient(dval(rv), tens[-e]); dval 2352 util.c dval(rv) *= tens[i]; dval 2378 util.c dval(rv0) = 1e300; dval 2379 util.c dval(rv0) *= dval(rv0); dval 2392 util.c dval(rv) *= bigtens[j]; dval 2395 util.c dval(rv) *= bigtens[j]; dval 2412 util.c dval(rv) /= tens[i]; dval 2421 util.c dval(rv) *= tinytens[j]; dval 2438 util.c dval(rv) *= tinytens[j]; dval 2440 util.c dval(rv0) = dval(rv); dval 2441 util.c dval(rv) *= tinytens[j]; dval 2442 util.c if (!dval(rv)) { dval 2443 util.c dval(rv) = 2.*dval(rv0); dval 2444 util.c dval(rv) *= tinytens[j]; dval 2446 util.c if (!dval(rv)) { dval 2448 util.c dval(rv) = 0.; dval 2476 util.c bb = d2b(dval(rv), &bbe, &bbbits); /* rv = bb * 2^bbe */ dval 2593 util.c dval(rv) += adj*ulp(dval(rv)); dval 2599 util.c dval(rv) += adj*ulp(dval(rv)); dval 2622 util.c adj *= ulp(dval(rv)); dval 2624 util.c dval(rv) += adj; dval 2626 util.c dval(rv) -= adj; dval 2632 util.c adj *= ulp(dval(rv)); dval 2634 util.c dval(rv) += adj; dval 2636 util.c dval(rv) -= adj; dval 2741 util.c dval(rv) += ulp(dval(rv)); dval 2744 util.c dval(rv) -= ulp(dval(rv)); dval 2746 util.c if (!dval(rv)) dval 2758 util.c aadj = dval(aadj1) = 1.; dval 2765 util.c dval(aadj1) = -1.; dval 2775 util.c dval(aadj1) = -aadj; dval 2780 util.c dval(aadj1) = dsign ? aadj : -aadj; dval 2784 util.c dval(aadj1) -= 0.5; dval 2788 util.c dval(aadj1) += 0.5; dval 2792 util.c dval(aadj1) += 0.5; dval 2800 util.c dval(rv0) = dval(rv); dval 2802 util.c adj = dval(aadj1) * ulp(dval(rv)); dval 2803 util.c dval(rv) += adj; dval 2822 util.c dval(aadj1) = dsign ? aadj : -aadj; dval 2826 util.c adj = dval(aadj1) * ulp(dval(rv)); dval 2827 util.c dval(rv) += adj; dval 2831 util.c dval(rv0) = dval(rv); dval 2833 util.c adj = dval(aadj1) * ulp(dval(rv)); dval 2834 util.c dval(rv) += adj; dval 2851 util.c adj = dval(aadj1) * ulp(dval(rv)); dval 2852 util.c dval(rv) += adj; dval 2863 util.c dval(aadj1) = (double)(int)(aadj + 0.5); dval 2865 util.c dval(aadj1) = -dval(aadj1); dval 2867 util.c adj = dval(aadj1) * ulp(dval(rv)); dval 2868 util.c dval(rv) += adj; dval 2901 util.c dval(rv0) += 1.; dval 2911 util.c dval(rv) *= dval(rv0); dval 2922 util.c dval(rv0) = 1e-300; dval 2923 util.c dval(rv0) *= dval(rv0); dval 2935 util.c return sign ? -dval(rv) : dval(rv); dval 3185 util.c dval(d) = d_; dval 3219 util.c dval(d) += 0; /* normalize */ dval 3221 util.c if (!dval(d)) { dval 3240 util.c b = d2b(dval(d), &be, &bbits); dval 3246 util.c dval(d2) = dval(d); dval 3251 util.c dval(d2) /= 1 << j; dval 3290 util.c dval(d2) = x; dval 3296 util.c ds = (dval(d2)-1.5)*0.289529654602168 + 0.1760912590558 + i*0.301029995663981; dval 3302 util.c if (dval(d) < tens[k]) dval 3378 util.c dval(d2) = dval(d); dval 3388 util.c dval(d) /= bigtens[n_bigtens-1]; dval 3396 util.c dval(d) /= ds; dval 3399 util.c dval(d) *= tens[j1 & 0xf]; dval 3403 util.c dval(d) *= bigtens[i]; dval 3406 util.c if (k_check && dval(d) < 1. && ilim > 0) { dval 3411 util.c dval(d) *= 10.; dval 3414 util.c dval(eps) = ieps*dval(d) + 7.; dval 3418 util.c dval(d) -= 5.; dval 3419 util.c if (dval(d) > dval(eps)) dval 3421 util.c if (dval(d) < -dval(eps)) dval 3430 util.c dval(eps) = 0.5/tens[ilim-1] - dval(eps); dval 3432 util.c L = (int)dval(d); dval 3433 util.c dval(d) -= L; dval 3435 util.c if (dval(d) < dval(eps)) dval 3437 util.c if (1. - dval(d) < dval(eps)) dval 3441 util.c dval(eps) *= 10.; dval 3442 util.c dval(d) *= 10.; dval 3448 util.c dval(eps) *= tens[ilim-1]; dval 3449 util.c for (i = 1;; i++, dval(d) *= 10.) { dval 3450 util.c L = (Long)(dval(d)); dval 3451 util.c if (!(dval(d) -= L)) dval 3455 util.c if (dval(d) > 0.5 + dval(eps)) dval 3457 util.c else if (dval(d) < 0.5 - dval(eps)) { dval 3470 util.c dval(d) = dval(d2); dval 3482 util.c if (ilim < 0 || dval(d) <= 5*ds) dval 3486 util.c for (i = 1;; i++, dval(d) *= 10.) { dval 3487 util.c L = (Long)(dval(d) / ds); dval 3488 util.c dval(d) -= L*ds; dval 3491 util.c if (dval(d) < 0) { dval 3493 util.c dval(d) += ds; dval 3497 util.c if (!dval(d)) { dval 3511 util.c dval(d) += dval(d); dval 3512 util.c if (dval(d) > ds || (dval(d) == ds && (L & 1))) { dval 3788 util.c dval(d) += 1.;