4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync/* $NetBSD: strtod.c,v 1.4.14.1 2008/04/08 21:10:55 jdc Exp $ */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync/****************************************************************
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsyncThe author of this software is David M. Gay.
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsyncCopyright (C) 1998-2001 by Lucent Technologies
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsyncAll Rights Reserved
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsyncPermission to use, copy, modify, and distribute this software and
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsyncits documentation for any purpose and without fee is hereby
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsyncgranted, provided that the above copyright notice appear in all
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsynccopies and that both that the copyright notice and this
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsyncpermission notice and warranty disclaimer appear in supporting
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsyncdocumentation, and that the name of Lucent or any of its entities
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsyncnot be used in advertising or publicity pertaining to
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsyncdistribution of the software without specific, written prior
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsyncpermission.
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsyncLUCENT DISCLAIMS ALL WARRANTIES WITH REGARD TO THIS SOFTWARE,
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsyncINCLUDING ALL IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS.
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsyncIN NO EVENT SHALL LUCENT OR ANY OF ITS ENTITIES BE LIABLE FOR ANY
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsyncSPECIAL, INDIRECT OR CONSEQUENTIAL DAMAGES OR ANY DAMAGES
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsyncWHATSOEVER RESULTING FROM LOSS OF USE, DATA OR PROFITS, WHETHER
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsyncIN AN ACTION OF CONTRACT, NEGLIGENCE OR OTHER TORTIOUS ACTION,
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsyncARISING OUT OF OR IN CONNECTION WITH THE USE OR PERFORMANCE OF
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsyncTHIS SOFTWARE.
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync****************************************************************/
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync/* Please send bug reports to David M. Gay (dmg at acm dot org,
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync * with " at " changed at "@" and " dot " changed to "."). */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#include <LibConfig.h>
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#include "gdtoaimp.h"
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifndef NO_FENV_H
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#include <fenv.h>
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef USE_LOCALE
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#include "locale.h"
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef IEEE_Arith
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifndef NO_IEEE_Scale
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#define Avoid_Underflow
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#undef tinytens
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync/* The factor of 2^53 in tinytens[4] helps us avoid setting the underflow */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync/* flag unnecessarily. It leads to a song and dance at the end of strtod. */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsyncstatic CONST double tinytens[] = { 1e-16, 1e-32, 1e-64, 1e-128,
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync 9007199254740992.e-256
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync };
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef Honor_FLT_ROUNDS
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#define Rounding rounding
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#undef Check_FLT_ROUNDS
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#define Check_FLT_ROUNDS
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#else
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#define Rounding Flt_Rounds
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync//#ifndef __HAVE_LONG_DOUBLE
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync//__strong_alias(_strtold, strtod)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync//__weak_alias(strtold, _strtold)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync//#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#if defined(_MSC_VER) /* Handle Microsoft VC++ compiler specifics. */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync// Disable: warning C4700: uninitialized local variable 'xx' used
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#pragma warning ( disable : 4700 )
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif /* defined(_MSC_VER) */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsyncdouble
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsyncstrtod(CONST char *s00, char **se)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync{
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef Avoid_Underflow
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync int scale;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync #endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync int bb2, bb5, bbe, bd2, bd5, bbbits, bs2, c, decpt, dsign,
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync e, e1, esign, i, j, k, nd, nd0, nf, nz, nz0, sign;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync CONST char *s, *s0, *s1;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync double aadj, aadj1, adj, rv, rv0;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync Long L;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync ULong y, z;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync Bigint *bb = NULL, *bb1, *bd0;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync Bigint *bd = NULL, *bs = NULL, *delta = NULL; /* pacify gcc */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef SET_INEXACT
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync int inexact, oldinexact;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef Honor_FLT_ROUNDS
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync int rounding;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync sign = nz0 = nz = decpt = 0;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync dval(rv) = 0.;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync for(s = s00;;s++) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync switch(*s) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync case '-':
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync sign = 1;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* FALLTHROUGH */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync case '+':
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (*++s)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto break2;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* FALLTHROUGH */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync case 0:
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto ret0;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync case '\t':
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync case '\n':
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync case '\v':
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync case '\f':
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync case '\r':
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync case ' ':
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync continue;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync default:
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto break2;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync break2:
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (*s == '0') {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifndef NO_HEX_FP
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync static FPI fpi = { 53, 1-1023-53+1, 2046-1023-53+1, 1, SI };
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync Long expt;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync ULong bits[2];
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync switch(s[1]) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync case 'x':
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync case 'X':
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#if defined(FE_DOWNWARD) && defined(FE_TONEAREST) && defined(FE_TOWARDZERO) && defined(FE_UPWARD)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync FPI fpi1 = fpi;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync switch(fegetround()) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync case FE_TOWARDZERO: fpi1.rounding = 0; break;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync case FE_UPWARD: fpi1.rounding = 2; break;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync case FE_DOWNWARD: fpi1.rounding = 3;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#else
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync switch((i = gethex(&s, &fpi, &expt, &bb, sign)) & STRTOG_Retmask) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync case STRTOG_NoNumber:
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync s = s00;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync sign = 0;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* FALLTHROUGH */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync case STRTOG_Zero:
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync break;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync default:
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (bb) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync copybits(bits, fpi.nbits, bb);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync Bfree(bb);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync ULtod((/* LINTED */(U*)&rv)->L, bits, expt, i);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }}
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto ret;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync nz0 = 1;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync while(*++s == '0') ;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (!*s)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto ret;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync s0 = s;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync y = z = 0;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync for(nd = nf = 0; (c = *s) >= '0' && c <= '9'; nd++, s++)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (nd < 9)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync y = 10*y + c - '0';
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync else if (nd < 16)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync z = 10*z + c - '0';
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync nd0 = nd;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef USE_LOCALE
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (c == *localeconv()->decimal_point)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#else
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (c == '.')
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync decpt = 1;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync c = *++s;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (!nd) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync for(; c == '0'; c = *++s)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync nz++;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (c > '0' && c <= '9') {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync s0 = s;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync nf += nz;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync nz = 0;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto have_dig;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto dig_done;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync for(; c >= '0' && c <= '9'; c = *++s) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync have_dig:
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync nz++;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (c -= '0') {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync nf += nz;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync for(i = 1; i < nz; i++)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (nd++ < 9)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync y *= 10;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync else if (nd <= DBL_DIG + 1)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync z *= 10;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (nd++ < 9)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync y = 10*y + c;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync else if (nd <= DBL_DIG + 1)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync z = 10*z + c;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync nz = 0;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync dig_done:
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync e = 0;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (c == 'e' || c == 'E') {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (!nd && !nz && !nz0) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto ret0;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync s00 = s;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync esign = 0;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync switch(c = *++s) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync case '-':
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync esign = 1;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* FALLTHROUGH */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync case '+':
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync c = *++s;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (c >= '0' && c <= '9') {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync while(c == '0')
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync c = *++s;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (c > '0' && c <= '9') {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync L = c - '0';
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync s1 = s;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync while((c = *++s) >= '0' && c <= '9')
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync L = 10*L + c - '0';
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (s - s1 > 8 || L > 19999)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* Avoid confusion from exponents
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync * so large that e might overflow.
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync e = 19999; /* safe for 16 bit ints */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync else
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync e = (int)L;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (esign)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync e = -e;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync else
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync e = 0;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync else
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync s = s00;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (!nd) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (!nz && !nz0) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef INFNAN_CHECK
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* Check for Nan and Infinity */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifndef No_Hex_NaN
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync ULong bits[2];
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync static FPI fpinan = /* only 52 explicit bits */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync { 52, 1-1023-53+1, 2046-1023-53+1, 1, SI };
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif // No_Hex_NaN
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (!decpt)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync switch(c) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync case 'i':
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync case 'I':
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (match(&s,"nf")) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync --s;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (!match(&s,"inity"))
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync ++s;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word0(rv) = 0x7ff00000;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word1(rv) = 0;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto ret;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync break;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync case 'n':
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync case 'N':
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (match(&s, "an")) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifndef No_Hex_NaN
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (*s == '(' /*)*/
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync && hexnan(&s, &fpinan, bits)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync == STRTOG_NaNbits) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word0(rv) = (UINT32)(0x7ff00000U | bits[1]);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word1(rv) = (UINT32)bits[0];
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync else {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word0(rv) = NAN_WORD0;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word1(rv) = NAN_WORD1;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifndef No_Hex_NaN
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto ret;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif /* INFNAN_CHECK */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync ret0:
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync s = s00;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync sign = 0;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto ret;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync e1 = e -= nf;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* Now we have nd0 digits, starting at s0, followed by a
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync * decimal point, followed by nd-nd0 digits. The number we're
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync * after is the integer represented by those digits times
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync * 10**e */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (!nd0)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync nd0 = nd;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync k = nd < DBL_DIG + 1 ? nd : DBL_DIG + 1;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync dval(rv) = (double)y;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (k > 9) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef SET_INEXACT
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (k > DBL_DIG)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync oldinexact = get_inexact();
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync dval(rv) = tens[k - 9] * dval(rv) + z;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync bd0 = 0;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (nd <= DBL_DIG
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifndef RND_PRODQUOT
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifndef Honor_FLT_ROUNDS
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync && Flt_Rounds == 1
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync ) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (!e)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto ret;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (e > 0) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (e <= Ten_pmax) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef VAX
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto vax_ovfl_check;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#else
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef Honor_FLT_ROUNDS
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* round correctly FLT_ROUNDS = 2 or 3 */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (sign) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync rv = -rv;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync sign = 0;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* rv = */ rounded_product(dval(rv), tens[e]);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto ret;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync i = DBL_DIG - nd;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (e <= Ten_pmax + i) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* A fancier test would sometimes let us do
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync * this for larger i values.
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef Honor_FLT_ROUNDS
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* round correctly FLT_ROUNDS = 2 or 3 */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (sign) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync rv = -rv;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync sign = 0;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync e -= i;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync dval(rv) *= tens[i];
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef VAX
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* VAX exponent range is so narrow we must
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync * worry about overflow here...
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync vax_ovfl_check:
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word0(rv) -= P*Exp_msk1;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* rv = */ rounded_product(dval(rv), tens[e]);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if ((word0(rv) & Exp_mask)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync > Exp_msk1*(DBL_MAX_EXP+Bias-1-P))
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto ovfl;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word0(rv) += P*Exp_msk1;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#else
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* rv = */ rounded_product(dval(rv), tens[e]);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto ret;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifndef Inaccurate_Divide
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync else if (e >= -Ten_pmax) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef Honor_FLT_ROUNDS
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* round correctly FLT_ROUNDS = 2 or 3 */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (sign) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync rv = -rv;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync sign = 0;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* rv = */ rounded_quotient(dval(rv), tens[-e]);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto ret;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync e1 += nd - k;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef IEEE_Arith
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef SET_INEXACT
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync inexact = 1;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (k <= DBL_DIG)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync oldinexact = get_inexact();
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef Avoid_Underflow
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync scale = 0;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef Honor_FLT_ROUNDS
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if ((rounding = Flt_Rounds) >= 2) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (sign)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync rounding = rounding == 2 ? 0 : 2;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync else
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (rounding != 2)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync rounding = 0;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif /*IEEE_Arith*/
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* Get starting approximation = rv * 10**e1 */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (e1 > 0) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if ( (i = e1 & 15) !=0)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync dval(rv) *= tens[i];
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (e1 &= ~15) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (e1 > DBL_MAX_10_EXP) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync ovfl:
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifndef NO_ERRNO
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync errno = ERANGE;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* Can't trust HUGE_VAL */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef IEEE_Arith
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef Honor_FLT_ROUNDS
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync switch(rounding) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync case 0: /* toward 0 */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync case 3: /* toward -infinity */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word0(rv) = Big0;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word1(rv) = Big1;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync break;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync default:
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word0(rv) = Exp_mask;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word1(rv) = 0;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#else /*Honor_FLT_ROUNDS*/
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word0(rv) = Exp_mask;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word1(rv) = 0;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif /*Honor_FLT_ROUNDS*/
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef SET_INEXACT
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* set overflow bit */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync dval(rv0) = 1e300;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync dval(rv0) *= dval(rv0);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#else /*IEEE_Arith*/
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word0(rv) = Big0;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word1(rv) = Big1;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif /*IEEE_Arith*/
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (bd0)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto retfree;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto ret;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync e1 = (unsigned int)e1 >> 4;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync for(j = 0; e1 > 1; j++, e1 = (unsigned int)e1 >> 1)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (e1 & 1)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync dval(rv) *= bigtens[j];
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* The last multiplication could overflow. */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word0(rv) -= P*Exp_msk1;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync dval(rv) *= bigtens[j];
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if ((z = word0(rv) & Exp_mask)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync > Exp_msk1*(DBL_MAX_EXP+Bias-P))
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto ovfl;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (z > Exp_msk1*(DBL_MAX_EXP+Bias-1-P)) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* set to largest number */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* (Can't trust DBL_MAX) */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word0(rv) = Big0;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word1(rv) = Big1;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync else
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word0(rv) += P*Exp_msk1;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync else if (e1 < 0) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync e1 = -e1;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if ( (i = e1 & 15) !=0)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync dval(rv) /= tens[i];
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (e1 >>= 4) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (e1 >= 1 << n_bigtens)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto undfl;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef Avoid_Underflow
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (e1 & Scale_Bit)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync scale = 2*P;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync for(j = 0; e1 > 0; j++, e1 = (unsigned int)e1 >> 1)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (e1 & 1)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync dval(rv) *= tinytens[j];
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (scale && (j = 2*P + 1 - ((word0(rv) & Exp_mask)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync >> Exp_shift)) > 0) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* scaled rv is denormal; zap j low bits */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (j >= 32) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word1(rv) = 0;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (j >= 53)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word0(rv) = (P+2)*Exp_msk1;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync else
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word0(rv) &= 0xffffffff << (j-32);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync else
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word1(rv) &= 0xffffffff << j;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#else
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync for(j = 0; e1 > 1; j++, e1 >>= 1)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (e1 & 1)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync dval(rv) *= tinytens[j];
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* The last multiplication could underflow. */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync dval(rv0) = dval(rv);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync dval(rv) *= tinytens[j];
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (!dval(rv)) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync dval(rv) = 2.*dval(rv0);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync dval(rv) *= tinytens[j];
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (!dval(rv)) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync undfl:
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync dval(rv) = 0.;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifndef NO_ERRNO
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync errno = ERANGE;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (bd0)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto retfree;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto ret;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifndef Avoid_Underflow
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word0(rv) = Tiny0;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word1(rv) = Tiny1;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* The refinement below will clean
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync * this approximation up.
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* Now the hard part -- adjusting rv to the correct value.*/
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* Put digits into bd: true value = bd * 10^e */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync bd0 = s2b(s0, nd0, nd, y);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (bd0 == NULL)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto ovfl;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync for(;;) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync bd = Balloc(bd0->k);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (bd == NULL)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto ovfl;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync Bcopy(bd, bd0);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync bb = d2b(dval(rv), &bbe, &bbbits); /* rv = bb * 2^bbe */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (bb == NULL)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto ovfl;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync bs = i2b(1);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (bs == NULL)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto ovfl;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (e >= 0) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync bb2 = bb5 = 0;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync bd2 = bd5 = e;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync else {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync bb2 = bb5 = -e;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync bd2 = bd5 = 0;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (bbe >= 0)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync bb2 += bbe;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync else
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync bd2 -= bbe;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync bs2 = bb2;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef Honor_FLT_ROUNDS
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (rounding != 1)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync bs2++;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef Avoid_Underflow
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync j = bbe - scale;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync i = j + bbbits - 1; /* logb(rv) */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (i < Emin) /* denormal */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync j += P - Emin;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync else
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync j = P + 1 - bbbits;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#else /*Avoid_Underflow*/
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef Sudden_Underflow
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef IBM
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync j = 1 + 4*P - 3 - bbbits + ((bbe + bbbits - 1) & 3);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#else
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync j = P + 1 - bbbits;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#else /*Sudden_Underflow*/
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync j = bbe;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync i = j + bbbits - 1; /* logb(rv) */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (i < Emin) /* denormal */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync j += P - Emin;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync else
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync j = P + 1 - bbbits;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif /*Sudden_Underflow*/
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif /*Avoid_Underflow*/
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync bb2 += j;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync bd2 += j;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef Avoid_Underflow
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync bd2 += scale;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync i = bb2 < bd2 ? bb2 : bd2;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (i > bs2)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync i = bs2;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (i > 0) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync bb2 -= i;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync bd2 -= i;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync bs2 -= i;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (bb5 > 0) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync bs = pow5mult(bs, bb5);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (bs == NULL)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto ovfl;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync bb1 = mult(bs, bb);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (bb1 == NULL)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto ovfl;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync Bfree(bb);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync bb = bb1;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (bb2 > 0) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync bb = lshift(bb, bb2);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (bb == NULL)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto ovfl;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (bd5 > 0) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync bd = pow5mult(bd, bd5);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (bd == NULL)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto ovfl;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (bd2 > 0) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync bd = lshift(bd, bd2);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (bd == NULL)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto ovfl;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (bs2 > 0) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync bs = lshift(bs, bs2);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (bs == NULL)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto ovfl;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync delta = diff(bb, bd);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (delta == NULL)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto ovfl;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync dsign = delta->sign;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync delta->sign = 0;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync i = cmp(delta, bs);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef Honor_FLT_ROUNDS
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (rounding != 1) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (i < 0) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* Error is less than an ulp */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (!delta->x[0] && delta->wds <= 1) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* exact */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef SET_INEXACT
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync inexact = 0;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync break;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (rounding) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (dsign) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync adj = 1.;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto apply_adj;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync else if (!dsign) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync adj = -1.;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (!word1(rv)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync && !(word0(rv) & Frac_mask)) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync y = word0(rv) & Exp_mask;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef Avoid_Underflow
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (!scale || y > 2*P*Exp_msk1)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#else
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (y)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync delta = lshift(delta,Log2P);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (cmp(delta, bs) <= 0)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync adj = -0.5;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync apply_adj:
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef Avoid_Underflow
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (scale && (y = word0(rv) & Exp_mask)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync <= 2*P*Exp_msk1)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word0(adj) += (2*P+1)*Exp_msk1 - y;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#else
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef Sudden_Underflow
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if ((word0(rv) & Exp_mask) <=
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync P*Exp_msk1) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word0(rv) += P*Exp_msk1;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync dval(rv) += adj*ulp(dval(rv));
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word0(rv) -= P*Exp_msk1;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync else
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif /*Sudden_Underflow*/
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif /*Avoid_Underflow*/
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync dval(rv) += adj*ulp(dval(rv));
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync break;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync adj = ratio(delta, bs);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (adj < 1.)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync adj = 1.;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (adj <= 0x7ffffffe) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* adj = rounding ? ceil(adj) : floor(adj); */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync y = adj;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (y != adj) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (!((rounding>>1) ^ dsign))
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync y++;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync adj = y;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef Avoid_Underflow
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (scale && (y = word0(rv) & Exp_mask) <= 2*P*Exp_msk1)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word0(adj) += (2*P+1)*Exp_msk1 - y;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#else
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef Sudden_Underflow
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if ((word0(rv) & Exp_mask) <= P*Exp_msk1) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word0(rv) += P*Exp_msk1;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync adj *= ulp(dval(rv));
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (dsign)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync dval(rv) += adj;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync else
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync dval(rv) -= adj;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word0(rv) -= P*Exp_msk1;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto cont;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif /*Sudden_Underflow*/
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif /*Avoid_Underflow*/
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync adj *= ulp(dval(rv));
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (dsign)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync dval(rv) += adj;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync else
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync dval(rv) -= adj;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto cont;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif /*Honor_FLT_ROUNDS*/
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (i < 0) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* Error is less than half an ulp -- check for
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync * special case of mantissa a power of two.
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (dsign || word1(rv) || word0(rv) & Bndry_mask
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef IEEE_Arith
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef Avoid_Underflow
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync || (word0(rv) & Exp_mask) <= (2*P+1)*Exp_msk1
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#else
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync || (word0(rv) & Exp_mask) <= Exp_msk1
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync ) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef SET_INEXACT
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (!delta->x[0] && delta->wds <= 1)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync inexact = 0;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync break;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (!delta->x[0] && delta->wds <= 1) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* exact result */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef SET_INEXACT
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync inexact = 0;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync break;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync delta = lshift(delta,Log2P);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (cmp(delta, bs) > 0)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto drop_down;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync break;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (i == 0) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* exactly half-way between */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (dsign) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if ((word0(rv) & Bndry_mask1) == Bndry_mask1
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync && word1(rv) == (
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef Avoid_Underflow
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync (scale && (y = word0(rv) & Exp_mask) <= 2*P*Exp_msk1)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync ? (0xffffffff & (0xffffffff << (2*P+1-(y>>Exp_shift)))) :
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync 0xffffffff)) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /*boundary case -- increment exponent*/
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word0(rv) = (word0(rv) & Exp_mask)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync + Exp_msk1
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef IBM
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync | Exp_msk1 >> 4
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync ;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word1(rv) = 0;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef Avoid_Underflow
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync dsign = 0;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync break;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync else if (!(word0(rv) & Bndry_mask) && !word1(rv)) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync drop_down:
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* boundary case -- decrement exponent */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef Sudden_Underflow /*{{*/
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync L = word0(rv) & Exp_mask;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef IBM
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (L < Exp_msk1)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#else
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef Avoid_Underflow
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (L <= (scale ? (2*P+1)*Exp_msk1 : Exp_msk1))
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#else
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (L <= Exp_msk1)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif /*Avoid_Underflow*/
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif /*IBM*/
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto undfl;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync L -= Exp_msk1;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#else /*Sudden_Underflow}{*/
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef Avoid_Underflow
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (scale) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync L = word0(rv) & Exp_mask;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (L <= (2*P+1)*Exp_msk1) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (L > (P+2)*Exp_msk1)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* round even ==> */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* accept rv */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync break;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* rv = smallest denormal */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto undfl;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif /*Avoid_Underflow*/
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync L = (word0(rv) & Exp_mask) - Exp_msk1;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif /*Sudden_Underflow}*/
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word0(rv) = (UINT32)(L | Bndry_mask1);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word1(rv) = 0xffffffffU;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef IBM
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto cont;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#else
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync break;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifndef ROUND_BIASED
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (!(word1(rv) & LSB))
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync break;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (dsign)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync dval(rv) += ulp(dval(rv));
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifndef ROUND_BIASED
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync else {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync dval(rv) -= ulp(dval(rv));
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifndef Sudden_Underflow
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (!dval(rv))
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto undfl;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef Avoid_Underflow
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync dsign = 1 - dsign;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync break;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if ((aadj = ratio(delta, bs)) <= 2.) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (dsign)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync aadj = aadj1 = 1.;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync else if (word1(rv) || word0(rv) & Bndry_mask) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifndef Sudden_Underflow
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (word1(rv) == Tiny1 && !word0(rv))
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto undfl;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync aadj = 1.;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync aadj1 = -1.;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync else {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* special case -- power of FLT_RADIX to be */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* rounded down... */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (aadj < 2./FLT_RADIX)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync aadj = 1./FLT_RADIX;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync else
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync aadj *= 0.5;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync aadj1 = -aadj;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync else {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync aadj *= 0.5;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync aadj1 = dsign ? aadj : -aadj;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef Check_FLT_ROUNDS
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync switch(Rounding) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync case 2: /* towards +infinity */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync aadj1 -= 0.5;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync break;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync case 0: /* towards 0 */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync case 3: /* towards -infinity */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync aadj1 += 0.5;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#else
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (Flt_Rounds == 0)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync aadj1 += 0.5;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif /*Check_FLT_ROUNDS*/
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync y = word0(rv) & Exp_mask;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* Check for overflow */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (y == Exp_msk1*(DBL_MAX_EXP+Bias-1)) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync dval(rv0) = dval(rv);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word0(rv) -= P*Exp_msk1;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync adj = aadj1 * ulp(dval(rv));
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync dval(rv) += adj;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if ((word0(rv) & Exp_mask) >=
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync Exp_msk1*(DBL_MAX_EXP+Bias-P)) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (word0(rv0) == Big0 && word1(rv0) == Big1)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto ovfl;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word0(rv) = Big0;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word1(rv) = Big1;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto cont;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync else
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word0(rv) += P*Exp_msk1;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync else {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef Avoid_Underflow
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (scale && y <= 2*P*Exp_msk1) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (aadj <= 0x7fffffff) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if ((z = (uint32_t)aadj) == 0)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync z = 1;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync aadj = (double)z;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync aadj1 = dsign ? aadj : -aadj;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word0(aadj1) += (UINT32)((2*P+1)*Exp_msk1 - y);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync adj = aadj1 * ulp(dval(rv));
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync dval(rv) += adj;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#else
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef Sudden_Underflow
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if ((word0(rv) & Exp_mask) <= P*Exp_msk1) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync dval(rv0) = dval(rv);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word0(rv) += P*Exp_msk1;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync adj = aadj1 * ulp(dval(rv));
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync dval(rv) += adj;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef IBM
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if ((word0(rv) & Exp_mask) < P*Exp_msk1)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#else
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if ((word0(rv) & Exp_mask) <= P*Exp_msk1)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (word0(rv0) == Tiny0
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync && word1(rv0) == Tiny1)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto undfl;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word0(rv) = Tiny0;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word1(rv) = Tiny1;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync goto cont;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync else
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word0(rv) -= P*Exp_msk1;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync else {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync adj = aadj1 * ulp(dval(rv));
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync dval(rv) += adj;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#else /*Sudden_Underflow*/
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* Compute adj so that the IEEE rounding rules will
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync * correctly round rv + adj in some half-way cases.
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync * If rv * ulp(rv) is denormalized (i.e.,
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync * y <= (P-1)*Exp_msk1), we must adjust aadj to avoid
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync * trouble from bits lost to denormalization;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync * example: 1.2e-307 .
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (y <= (P-1)*Exp_msk1 && aadj > 1.) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync aadj1 = (double)(int)(aadj + 0.5);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (!dsign)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync aadj1 = -aadj1;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync adj = aadj1 * ulp(dval(rv));
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync dval(rv) += adj;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif /*Sudden_Underflow*/
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif /*Avoid_Underflow*/
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync z = word0(rv) & Exp_mask;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifndef SET_INEXACT
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef Avoid_Underflow
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (!scale)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (y == z) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* Can we stop now? */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync L = (Long)aadj;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync aadj -= L;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* The tolerances below are conservative. */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (dsign || word1(rv) || word0(rv) & Bndry_mask) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (aadj < .4999999 || aadj > .5000001)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync break;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync else if (aadj < .4999999/FLT_RADIX)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync break;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync cont:
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync Bfree(bb);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync Bfree(bd);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync Bfree(bs);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync Bfree(delta);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef SET_INEXACT
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (inexact) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (!oldinexact) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word0(rv0) = Exp_1 + (70 << Exp_shift);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word1(rv0) = 0;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync dval(rv0) += 1.;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync else if (!oldinexact)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync clear_inexact();
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef Avoid_Underflow
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (scale) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word0(rv0) = Exp_1 - 2*P*Exp_msk1;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync word1(rv0) = 0;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync dval(rv) *= dval(rv0);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifndef NO_ERRNO
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* try to avoid the bug of testing an 8087 register value */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (word0(rv) == 0 && word1(rv) == 0)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync errno = ERANGE;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif /* Avoid_Underflow */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#ifdef SET_INEXACT
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (inexact && !(word0(rv) & Exp_mask)) {
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync /* set underflow bit */
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync dval(rv0) = 1e-300;
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync dval(rv0) *= dval(rv0);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync }
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync#endif
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync retfree:
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync Bfree(bb);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync Bfree(bd);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync Bfree(bs);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync Bfree(bd0);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync Bfree(delta);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync ret:
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync if (se)
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync *se = __UNCONST(s);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync return sign ? -dval(rv) : dval(rv);
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync}
4fd606d1f5abe38e1f42c38de1d2e895166bd0f4vboxsync