12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061626364656667686970717273747576777879808182838485868788899091929394959697989910010110210310410510610710810911011111211311411511611711811912012112212312412512612712812913013113213313413513613713813914014114214314414514614714814915015115215315415515615715815916016116216316416516616716816917017117217317417517617717817918018118218318418518618718818919019119219319419519619719819920020120220320420520620720820921021121221321421521621721821922022122222322422522622722822923023123223323423523623723823924024124224324424524624724824925025125225325425525625725825926026126226326426526626726826927027127227327427527627727827928028128228328428528628728828929029129229329429529629729829930030130230330430530630730830931031131231331431531631731831932032132232332432532632732832933033133233333433533633733833934034134234334434534634734834935035135235335435535635735835936036136236336436536636736836937037137237337437537637737837938038138238338438538638738838939039139239339439539639739839940040140240340440540640740840941041141241341441541641741841942042142242342442542642742842943043143243343443543643743843944044144244344444544644744844945045145245345445545645745845946046146246346446546646746846947047147247347447547647747847948048148248348448548648748848949049149249349449549649749849950050150250350450550650750850951051151251351451551651751851952052152252352452552652752852953053153253353453553653753853954054154254354454554654754854955055155255355455555655755855956056156256356456556656756856957057157257357457557657757857958058158258358458558658758858959059159259359459559659759859960060160260360460560660760860961061161261361461561661761861962062162262362462562662762862963063163263363463563663763863964064164264364464564664764864965065165265365465565665765865966066166266366466566666766866967067167267367467567667767867968068168268368468568668768868969069169269369469569669769869970070170270370470570670770870971071171271371471571671771871972072172272372472572672772872973073173273373473573673773873974074174274374474574674774874975075175275375475575675775875976076176276376476576676776876977077177277377477577677777877978078178278378478578678778878979079179279379479579679779879980080180280380480580680780880981081181281381481581681781881982082182282382482582682782882983083183283383483583683783883984084184284384484584684784884985085185285385485585685785885986086186286386486586686786886987087187287387487587687787887988088188288388488588688788888989089189289389489589689789889990090190290390490590690790890991091191291391491591691791891992092192292392492592692792892993093193293393493593693793893994094194294394494594694794894995095195295395495595695795895996096196296396496596696796896997097197297397497597697797897998098198298398498598698798898999099199299399499599699799899910001001100210031004100510061007100810091010101110121013101410151016101710181019102010211022102310241025102610271028102910301031103210331034103510361037103810391040104110421043104410451046104710481049105010511052105310541055105610571058105910601061106210631064106510661067106810691070107110721073107410751076107710781079108010811082108310841085108610871088108910901091109210931094109510961097109810991100110111021103110411051106110711081109111011111112111311141115111611171118111911201121112211231124112511261127112811291130113111321133113411351136113711381139114011411142114311441145114611471148114911501151115211531154115511561157115811591160116111621163116411651166116711681169117011711172117311741175117611771178117911801181118211831184118511861187118811891190119111921193119411951196119711981199120012011202120312041205120612071208120912101211121212131214121512161217121812191220122112221223122412251226122712281229123012311232123312341235123612371238123912401241124212431244124512461247124812491250125112521253125412551256125712581259126012611262126312641265126612671268126912701271127212731274127512761277127812791280128112821283128412851286128712881289129012911292129312941295129612971298129913001301130213031304130513061307130813091310131113121313131413151316131713181319 |
- #include <config.h>
- #include <float.h>
- #include <limits.h>
- #include <math.h>
- #include <string.h>
- #include <unistd.h>
- #include <stdlib.h>
- #include <stdbool.h>
- #define NDEBUG
- #include <assert.h>
- #ifdef HAVE_ERRNO_H
- #include <errno.h>
- #endif
- #include <stdio.h>
- #include <stdarg.h>
- #ifdef HAVE_FENV_H
- #include "quadmath-rounding-mode.h"
- #endif
- #include "quadmath-printf.h"
- #include "fpioconst.h"
- #ifdef USE_I18N_NUMBER_H
- #include "_i18n_number.h"
- #endif
- #define outchar(ch) \
- do \
- { \
- register const int outc = (ch); \
- if (PUTC (outc, fp) == EOF) \
- { \
- if (buffer_malloced) \
- free (wbuffer); \
- return -1; \
- } \
- ++done; \
- } while (0)
- #define PRINT(ptr, wptr, len) \
- do \
- { \
- register size_t outlen = (len); \
- if (len > 20) \
- { \
- if (PUT (fp, wide ? (const char *) wptr : ptr, outlen) != outlen) \
- { \
- if (buffer_malloced) \
- free (wbuffer); \
- return -1; \
- } \
- ptr += outlen; \
- done += outlen; \
- } \
- else \
- { \
- if (wide) \
- while (outlen-- > 0) \
- outchar (*wptr++); \
- else \
- while (outlen-- > 0) \
- outchar (*ptr++); \
- } \
- } while (0)
- #define PADN(ch, len) \
- do \
- { \
- if (PAD (fp, ch, len) != len) \
- { \
- if (buffer_malloced) \
- free (wbuffer); \
- return -1; \
- } \
- done += len; \
- } \
- while (0)
- #define MPN_VAR(name) mp_limb_t *name; mp_size_t name##size
- #define MPN_ASSIGN(dst,src) \
- memcpy (dst, src, (dst##size = src##size) * sizeof (mp_limb_t))
- #define MPN_GE(u,v) \
- (u##size > v##size || (u##size == v##size && mpn_cmp (u, v, u##size) >= 0))
- extern mp_size_t mpn_extract_flt128 (mp_ptr res_ptr, mp_size_t size,
- int *expt, int *is_neg,
- __float128 value) attribute_hidden;
- static unsigned int guess_grouping (unsigned int intdig_max,
- const char *grouping);
- static wchar_t *group_number (wchar_t *buf, wchar_t *bufend,
- unsigned int intdig_no, const char *grouping,
- wchar_t thousands_sep, int ngroups);
- int
- __quadmath_printf_fp (struct __quadmath_printf_file *fp,
- const struct printf_info *info,
- const void *const *args)
- {
-
- __float128 fpnum;
-
- const char *decimal;
- wchar_t decimalwc;
-
- const char *thousands_sep = NULL;
- wchar_t thousands_sepwc = L_('\0');
- const char *grouping;
-
- const char *special = NULL;
- const wchar_t *wspecial = NULL;
-
- mp_limb_t fp_input[(FLT128_MANT_DIG + BITS_PER_MP_LIMB - 1) / BITS_PER_MP_LIMB];
-
- int to_shift = 0;
-
- MPN_VAR(frac);
-
- int exponent;
-
- int expsign = 0;
-
- int is_neg = 0;
-
- MPN_VAR(scale);
-
- MPN_VAR(tmp);
-
- wchar_t last_digit, next_digit;
- bool more_bits;
-
- int type;
-
- int done = 0;
-
- mp_limb_t cy;
-
- int wide = info->wide;
-
- wchar_t *wbuffer = NULL;
-
- int buffer_malloced = 0;
- auto wchar_t hack_digit (void);
- wchar_t hack_digit (void)
- {
- mp_limb_t hi;
- if (expsign != 0 && type == 'f' && exponent-- > 0)
- hi = 0;
- else if (scalesize == 0)
- {
- hi = frac[fracsize - 1];
- frac[fracsize - 1] = mpn_mul_1 (frac, frac, fracsize - 1, 10);
- }
- else
- {
- if (fracsize < scalesize)
- hi = 0;
- else
- {
- hi = mpn_divmod (tmp, frac, fracsize, scale, scalesize);
- tmp[fracsize - scalesize] = hi;
- hi = tmp[0];
- fracsize = scalesize;
- while (fracsize != 0 && frac[fracsize - 1] == 0)
- --fracsize;
- if (fracsize == 0)
- {
-
- fracsize = 1;
- return L_('0') + hi;
- }
- }
- mp_limb_t _cy = mpn_mul_1 (frac, frac, fracsize, 10);
- if (_cy != 0)
- frac[fracsize++] = _cy;
- }
- return L_('0') + hi;
- }
-
- #ifdef USE_NL_LANGINFO
- if (info->extra == 0)
- decimal = nl_langinfo (DECIMAL_POINT);
- else
- {
- decimal = nl_langinfo (MON_DECIMAL_POINT);
- if (*decimal == '\0')
- decimal = nl_langinfo (DECIMAL_POINT);
- }
-
- assert (*decimal != '\0');
- #elif defined USE_LOCALECONV
- const struct lconv *lc = localeconv ();
- if (info->extra == 0)
- decimal = lc->decimal_point;
- else
- {
- decimal = lc->mon_decimal_point;
- if (decimal == NULL || *decimal == '\0')
- decimal = lc->decimal_point;
- }
- if (decimal == NULL || *decimal == '\0')
- decimal = ".";
- #else
- decimal = ".";
- #endif
- #ifdef USE_NL_LANGINFO_WC
- if (info->extra == 0)
- decimalwc = nl_langinfo_wc (_NL_NUMERIC_DECIMAL_POINT_WC);
- else
- {
- decimalwc = nl_langinfo_wc (_NL_MONETARY_DECIMAL_POINT_WC);
- if (decimalwc == L_('\0'))
- decimalwc = nl_langinfo_wc (_NL_NUMERIC_DECIMAL_POINT_WC);
- }
-
- assert (decimalwc != L_('\0'));
- #else
- decimalwc = L_('.');
- #endif
- #if defined USE_NL_LANGINFO && defined USE_NL_LANGINFO_WC
- if (info->group)
- {
- if (info->extra == 0)
- grouping = nl_langinfo (GROUPING);
- else
- grouping = nl_langinfo (MON_GROUPING);
- if (*grouping <= 0 || *grouping == CHAR_MAX)
- grouping = NULL;
- else
- {
-
- if (wide)
- {
- if (info->extra == 0)
- thousands_sepwc = nl_langinfo_wc (_NL_NUMERIC_THOUSANDS_SEP_WC);
- else
- thousands_sepwc = nl_langinfo_wc (_NL_MONETARY_THOUSANDS_SEP_WC);
- if (thousands_sepwc == L_('\0'))
- grouping = NULL;
- }
- else
- {
- if (info->extra == 0)
- thousands_sep = nl_langinfo (THOUSANDS_SEP);
- else
- thousands_sep = nl_langinfo (MON_THOUSANDS_SEP);
- if (*thousands_sep == '\0')
- grouping = NULL;
- }
- }
- }
- else
- #elif defined USE_NL_LANGINFO
- if (info->group && !wide)
- {
- if (info->extra == 0)
- grouping = nl_langinfo (GROUPING);
- else
- grouping = nl_langinfo (MON_GROUPING);
- if (*grouping <= 0 || *grouping == CHAR_MAX)
- grouping = NULL;
- else
- {
-
- if (info->extra == 0)
- thousands_sep = nl_langinfo (THOUSANDS_SEP);
- else
- thousands_sep = nl_langinfo (MON_THOUSANDS_SEP);
- if (*thousands_sep == '\0')
- grouping = NULL;
- }
- }
- else
- #elif defined USE_LOCALECONV
- if (info->group && !wide)
- {
- if (info->extra == 0)
- grouping = lc->grouping;
- else
- grouping = lc->mon_grouping;
- if (grouping == NULL || *grouping <= 0 || *grouping == CHAR_MAX)
- grouping = NULL;
- else
- {
-
- if (info->extra == 0)
- thousands_sep = lc->thousands_sep;
- else
- thousands_sep = lc->mon_thousands_sep;
- if (thousands_sep == NULL || *thousands_sep == '\0')
- grouping = NULL;
- }
- }
- else
- #endif
- grouping = NULL;
- if (grouping != NULL && !wide)
-
- thousands_sepwc = (wchar_t) 0xfffffffe;
-
- {
- fpnum = **(const __float128 **) args[0];
-
- if (isnanq (fpnum))
- {
- ieee854_float128 u = { .value = fpnum };
- is_neg = u.ieee.negative != 0;
- if (isupper (info->spec))
- {
- special = "NAN";
- wspecial = L_("NAN");
- }
- else
- {
- special = "nan";
- wspecial = L_("nan");
- }
- }
- else if (isinfq (fpnum))
- {
- is_neg = fpnum < 0;
- if (isupper (info->spec))
- {
- special = "INF";
- wspecial = L_("INF");
- }
- else
- {
- special = "inf";
- wspecial = L_("inf");
- }
- }
- else
- {
- fracsize = mpn_extract_flt128 (fp_input,
- (sizeof (fp_input) /
- sizeof (fp_input[0])),
- &exponent, &is_neg, fpnum);
- to_shift = 1 + fracsize * BITS_PER_MP_LIMB - FLT128_MANT_DIG;
- }
- }
- if (special)
- {
- int width = info->width;
- if (is_neg || info->showsign || info->space)
- --width;
- width -= 3;
- if (!info->left && width > 0)
- PADN (' ', width);
- if (is_neg)
- outchar ('-');
- else if (info->showsign)
- outchar ('+');
- else if (info->space)
- outchar (' ');
- PRINT (special, wspecial, 3);
- if (info->left && width > 0)
- PADN (' ', width);
- return done;
- }
-
- {
- mp_size_t bignum_size = ((ABS (exponent) + BITS_PER_MP_LIMB - 1)
- / BITS_PER_MP_LIMB
- + (FLT128_MANT_DIG / BITS_PER_MP_LIMB > 2 ? 8 : 4))
- * sizeof (mp_limb_t);
- frac = (mp_limb_t *) alloca (bignum_size);
- tmp = (mp_limb_t *) alloca (bignum_size);
- scale = (mp_limb_t *) alloca (bignum_size);
- }
-
- scalesize = 0;
- if (exponent > 2)
- {
-
- int scaleexpo = 0;
- int explog = FLT128_MAX_10_EXP_LOG;
- int exp10 = 0;
- const struct mp_power *powers = &_fpioconst_pow10[explog + 1];
- int cnt_h, cnt_l, i;
- if ((exponent + to_shift) % BITS_PER_MP_LIMB == 0)
- {
- MPN_COPY_DECR (frac + (exponent + to_shift) / BITS_PER_MP_LIMB,
- fp_input, fracsize);
- fracsize += (exponent + to_shift) / BITS_PER_MP_LIMB;
- }
- else
- {
- cy = mpn_lshift (frac + (exponent + to_shift) / BITS_PER_MP_LIMB,
- fp_input, fracsize,
- (exponent + to_shift) % BITS_PER_MP_LIMB);
- fracsize += (exponent + to_shift) / BITS_PER_MP_LIMB;
- if (cy)
- frac[fracsize++] = cy;
- }
- MPN_ZERO (frac, (exponent + to_shift) / BITS_PER_MP_LIMB);
- assert (powers > &_fpioconst_pow10[0]);
- do
- {
- --powers;
-
- if (exponent >= scaleexpo + powers->p_expo - 1)
- {
- if (scalesize == 0)
- {
- if (FLT128_MANT_DIG > _FPIO_CONST_OFFSET * BITS_PER_MP_LIMB)
- {
- #define _FPIO_CONST_SHIFT \
- (((FLT128_MANT_DIG + BITS_PER_MP_LIMB - 1) / BITS_PER_MP_LIMB) \
- - _FPIO_CONST_OFFSET)
-
- tmpsize = powers->arraysize + _FPIO_CONST_SHIFT;
- memcpy (tmp + _FPIO_CONST_SHIFT,
- &__tens[powers->arrayoff],
- tmpsize * sizeof (mp_limb_t));
- MPN_ZERO (tmp, _FPIO_CONST_SHIFT);
-
- exponent += _FPIO_CONST_SHIFT * BITS_PER_MP_LIMB;
- }
- else
- {
- tmpsize = powers->arraysize;
- memcpy (tmp, &__tens[powers->arrayoff],
- tmpsize * sizeof (mp_limb_t));
- }
- }
- else
- {
- cy = mpn_mul (tmp, scale, scalesize,
- &__tens[powers->arrayoff
- + _FPIO_CONST_OFFSET],
- powers->arraysize - _FPIO_CONST_OFFSET);
- tmpsize = scalesize + powers->arraysize - _FPIO_CONST_OFFSET;
- if (cy == 0)
- --tmpsize;
- }
- if (MPN_GE (frac, tmp))
- {
- int cnt;
- MPN_ASSIGN (scale, tmp);
- count_leading_zeros (cnt, scale[scalesize - 1]);
- scaleexpo = (scalesize - 2) * BITS_PER_MP_LIMB - cnt - 1;
- exp10 |= 1 << explog;
- }
- }
- --explog;
- }
- while (powers > &_fpioconst_pow10[0]);
- exponent = exp10;
-
- if (scalesize > 0)
- {
-
- for (i = 0; scale[i] == 0 && frac[i] == 0; i++)
- ;
-
- count_leading_zeros (cnt_h, scale[scalesize - 1]);
- if (cnt_h == 0)
- {
-
- if (i > 0)
- {
- MPN_COPY_INCR (scale, scale + i, scalesize - i);
- scalesize -= i;
- MPN_COPY_INCR (frac, frac + i, fracsize - i);
- fracsize -= i;
- }
- }
- else
- {
- if (scale[i] != 0)
- {
- count_trailing_zeros (cnt_l, scale[i]);
- if (frac[i] != 0)
- {
- int cnt_l2;
- count_trailing_zeros (cnt_l2, frac[i]);
- if (cnt_l2 < cnt_l)
- cnt_l = cnt_l2;
- }
- }
- else
- count_trailing_zeros (cnt_l, frac[i]);
-
- if (i == 0 && BITS_PER_MP_LIMB - cnt_h > cnt_l)
- {
-
- (void) mpn_lshift (scale, scale, scalesize, cnt_h);
- cy = mpn_lshift (frac, frac, fracsize, cnt_h);
- if (cy != 0)
- frac[fracsize++] = cy;
- }
- else if (BITS_PER_MP_LIMB - cnt_h <= cnt_l)
- {
-
- (void) mpn_rshift (scale, scale + i, scalesize - i,
- BITS_PER_MP_LIMB - cnt_h);
- scalesize -= i + 1;
- (void) mpn_rshift (frac, frac + i, fracsize - i,
- BITS_PER_MP_LIMB - cnt_h);
- fracsize -= frac[fracsize - i - 1] == 0 ? i + 1 : i;
- }
- else
- {
-
- (void) mpn_rshift (scale, scale + (i - 1),
- scalesize - (i - 1),
- BITS_PER_MP_LIMB - cnt_h);
- scalesize -= i;
- (void) mpn_rshift (frac, frac + (i - 1),
- fracsize - (i - 1),
- BITS_PER_MP_LIMB - cnt_h);
- fracsize -= frac[fracsize - (i - 1) - 1] == 0 ? i : i - 1;
- }
- }
- }
- }
- else if (exponent < 0)
- {
-
- int exp10 = 0;
- int explog = FLT128_MAX_10_EXP_LOG;
- const struct mp_power *powers = &_fpioconst_pow10[explog + 1];
-
- cy = mpn_lshift (frac, fp_input, fracsize, to_shift);
- frac[fracsize++] = cy;
- assert (cy == 1 || (frac[fracsize - 2] == 0 && frac[0] == 0));
- expsign = 1;
- exponent = -exponent;
- assert (powers != &_fpioconst_pow10[0]);
- do
- {
- --powers;
- if (exponent >= powers->m_expo)
- {
- int i, incr, cnt_h, cnt_l;
- mp_limb_t topval[2];
-
- if (fracsize < powers->arraysize - _FPIO_CONST_OFFSET)
- cy = mpn_mul (tmp, &__tens[powers->arrayoff
- + _FPIO_CONST_OFFSET],
- powers->arraysize - _FPIO_CONST_OFFSET,
- frac, fracsize);
- else
- cy = mpn_mul (tmp, frac, fracsize,
- &__tens[powers->arrayoff + _FPIO_CONST_OFFSET],
- powers->arraysize - _FPIO_CONST_OFFSET);
- tmpsize = fracsize + powers->arraysize - _FPIO_CONST_OFFSET;
- if (cy == 0)
- --tmpsize;
- count_leading_zeros (cnt_h, tmp[tmpsize - 1]);
- incr = (tmpsize - fracsize) * BITS_PER_MP_LIMB
- + BITS_PER_MP_LIMB - 1 - cnt_h;
- assert (incr <= powers->p_expo);
-
- if (incr == exponent + 3)
- {
- if (cnt_h <= BITS_PER_MP_LIMB - 4)
- {
- topval[0] = 0;
- topval[1]
- = ((mp_limb_t) 10) << (BITS_PER_MP_LIMB - 4 - cnt_h);
- }
- else
- {
- topval[0] = ((mp_limb_t) 10) << (BITS_PER_MP_LIMB - 4);
- topval[1] = 0;
- (void) mpn_lshift (topval, topval, 2,
- BITS_PER_MP_LIMB - cnt_h);
- }
- }
-
- if (incr < exponent + 3
- || (incr == exponent + 3 &&
- (tmp[tmpsize - 1] < topval[1]
- || (tmp[tmpsize - 1] == topval[1]
- && tmp[tmpsize - 2] < topval[0]))))
- {
-
- exponent -= incr;
- exp10 |= 1 << explog;
-
- if (exponent < 0)
- cnt_h += -exponent;
-
- for (i = 0; tmp[i] == 0; ++i);
- if (cnt_h == BITS_PER_MP_LIMB - 1)
- {
- MPN_COPY (frac, tmp + i, tmpsize - i);
- fracsize = tmpsize - i;
- }
- else
- {
- count_trailing_zeros (cnt_l, tmp[i]);
-
- if (i == 0 && BITS_PER_MP_LIMB - 1 - cnt_h > cnt_l)
- {
-
- cy = mpn_lshift (frac, tmp, tmpsize, cnt_h + 1);
- fracsize = tmpsize + 1;
- frac[fracsize - 1] = cy;
- }
- else if (BITS_PER_MP_LIMB - 1 - cnt_h <= cnt_l)
- {
- (void) mpn_rshift (frac, tmp + i, tmpsize - i,
- BITS_PER_MP_LIMB - 1 - cnt_h);
- fracsize = tmpsize - i;
- }
- else
- {
-
- (void) mpn_rshift (frac, tmp + (i - 1),
- tmpsize - (i - 1),
- BITS_PER_MP_LIMB - 1 - cnt_h);
- fracsize = tmpsize - (i - 1);
- }
- }
- }
- }
- --explog;
- }
- while (powers != &_fpioconst_pow10[1] && exponent > 0);
-
- if (exponent > 0)
- {
- int cnt_l;
- cy = mpn_mul_1 (tmp, frac, fracsize, 10);
- tmpsize = fracsize;
- assert (cy == 0 || tmp[tmpsize - 1] < 20);
- count_trailing_zeros (cnt_l, tmp[0]);
- if (cnt_l < MIN (4, exponent))
- {
- cy = mpn_lshift (frac, tmp, tmpsize,
- BITS_PER_MP_LIMB - MIN (4, exponent));
- if (cy != 0)
- frac[tmpsize++] = cy;
- }
- else
- (void) mpn_rshift (frac, tmp, tmpsize, MIN (4, exponent));
- fracsize = tmpsize;
- exp10 |= 1;
- assert (frac[fracsize - 1] < 10);
- }
- exponent = exp10;
- }
- else
- {
-
- assert (0 <= exponent && exponent < 3 &&
- exponent + to_shift < BITS_PER_MP_LIMB);
-
- cy = mpn_lshift (frac, fp_input, fracsize, (exponent + to_shift));
- frac[fracsize++] = cy;
- exponent = 0;
- }
- {
- int width = info->width;
- wchar_t *wstartp, *wcp;
- size_t chars_needed;
- int expscale;
- int intdig_max, intdig_no = 0;
- int fracdig_min;
- int fracdig_max;
- int dig_max;
- int significant;
- int ngroups = 0;
- char spec = tolower (info->spec);
- if (spec == 'e')
- {
- type = info->spec;
- intdig_max = 1;
- fracdig_min = fracdig_max = info->prec < 0 ? 6 : info->prec;
- chars_needed = 1 + 1 + (size_t) fracdig_max + 1 + 1 + 4;
-
- dig_max = __INT_MAX__;
- significant = 1;
- }
- else if (spec == 'f')
- {
- type = 'f';
- fracdig_min = fracdig_max = info->prec < 0 ? 6 : info->prec;
- dig_max = __INT_MAX__;
- significant = 1;
- if (expsign == 0)
- {
- intdig_max = exponent + 1;
-
- chars_needed = (size_t) exponent + 1 + 1 + (size_t) fracdig_max;
- }
- else
- {
- intdig_max = 1;
- chars_needed = 1 + 1 + (size_t) fracdig_max;
- }
- }
- else
- {
- dig_max = info->prec < 0 ? 6 : (info->prec == 0 ? 1 : info->prec);
- if ((expsign == 0 && exponent >= dig_max)
- || (expsign != 0 && exponent > 4))
- {
- if ('g' - 'G' == 'e' - 'E')
- type = 'E' + (info->spec - 'G');
- else
- type = isupper (info->spec) ? 'E' : 'e';
- fracdig_max = dig_max - 1;
- intdig_max = 1;
- chars_needed = 1 + 1 + (size_t) fracdig_max + 1 + 1 + 4;
- }
- else
- {
- type = 'f';
- intdig_max = expsign == 0 ? exponent + 1 : 0;
- fracdig_max = dig_max - intdig_max;
-
- chars_needed = (size_t) dig_max + 1 + 4;
- }
- fracdig_min = info->alt ? fracdig_max : 0;
- significant = 0;
- }
- if (grouping)
- {
-
- ngroups = guess_grouping (intdig_max, grouping);
-
- chars_needed += ngroups + 1;
- }
-
- if (__builtin_expect (chars_needed >= (size_t) -1 / sizeof (wchar_t) - 2
- || chars_needed < fracdig_max, 0))
- {
-
- #if defined HAVE_ERRNO_H && defined ERANGE
- errno = ERANGE;
- #endif
- return -1;
- }
- size_t wbuffer_to_alloc = (2 + chars_needed) * sizeof (wchar_t);
- buffer_malloced = wbuffer_to_alloc >= 4096;
- if (__builtin_expect (buffer_malloced, 0))
- {
- wbuffer = (wchar_t *) malloc (wbuffer_to_alloc);
- if (wbuffer == NULL)
-
- return -1;
- }
- else
- wbuffer = (wchar_t *) alloca (wbuffer_to_alloc);
- wcp = wstartp = wbuffer + 2;
-
- if (expsign == 0 || type != 'f')
- {
- assert (expsign == 0 || intdig_max == 1);
- while (intdig_no < intdig_max)
- {
- ++intdig_no;
- *wcp++ = hack_digit ();
- }
- significant = 1;
- if (info->alt
- || fracdig_min > 0
- || (fracdig_max > 0 && (fracsize > 1 || frac[0] != 0)))
- *wcp++ = decimalwc;
- }
- else
- {
-
- *wcp++ = L_('0');
- --exponent;
- *wcp++ = decimalwc;
- }
-
- int fracdig_no = 0;
- int added_zeros = 0;
- while (fracdig_no < fracdig_min + added_zeros
- || (fracdig_no < fracdig_max && (fracsize > 1 || frac[0] != 0)))
- {
- ++fracdig_no;
- *wcp = hack_digit ();
- if (*wcp++ != L_('0'))
- significant = 1;
- else if (significant == 0)
- {
- ++fracdig_max;
- if (fracdig_min > 0)
- ++added_zeros;
- }
- }
-
- last_digit = wcp[-1] != decimalwc ? wcp[-1] : wcp[-2];
- next_digit =hack_digit ();
- if (next_digit != L_('0') && next_digit != L_('5'))
- more_bits = true;
- else if (fracsize == 1 && frac[0] == 0)
-
- more_bits = false;
- else if (scalesize == 0)
- {
-
- size_t lcnt = fracsize;
- while (lcnt >= 1 && frac[lcnt - 1] == 0)
- --lcnt;
- more_bits = lcnt > 0;
- }
- else
- more_bits = true;
- #ifdef HAVE_FENV_H
- int rounding_mode = get_rounding_mode ();
- if (round_away (is_neg, (last_digit - L_('0')) & 1, next_digit >= L_('5'),
- more_bits, rounding_mode))
- {
- wchar_t *wtp = wcp;
- if (fracdig_no > 0)
- {
-
- int removed = 0;
- while (*--wtp != decimalwc && *wtp == L_('9'))
- {
- *wtp = L_('0');
- ++removed;
- }
- if (removed == fracdig_min && added_zeros > 0)
- --added_zeros;
- if (*wtp != decimalwc)
-
- (*wtp)++;
- else if (__builtin_expect (spec == 'g' && type == 'f' && info->alt
- && wtp == wstartp + 1
- && wstartp[0] == L_('0'),
- 0))
-
- --added_zeros;
- }
- if (fracdig_no == 0 || *wtp == decimalwc)
- {
-
- if (*(wtp - 1) == decimalwc)
- --wtp;
- while (--wtp >= wstartp && *wtp == L_('9'))
- *wtp = L_('0');
- if (wtp >= wstartp)
-
- (*wtp)++;
- else
-
- {
- if (type != 'f')
- {
- *wstartp = '1';
- exponent += expsign == 0 ? 1 : -1;
-
- if (exponent == 0)
- expsign = 0;
- }
- else if (intdig_no == dig_max)
- {
-
- *--wstartp = decimalwc;
- *--wstartp = L_('1');
- if (info->alt || fracdig_no > 0)
- {
-
- wstartp[intdig_no + 2] = L_('0');
- ++fracdig_no;
- }
- fracdig_no += intdig_no;
- intdig_no = 1;
- fracdig_max = intdig_max - intdig_no;
- ++exponent;
-
- type = isupper (info->spec) ? 'E' : 'e';
- }
- else
- {
-
- *--wstartp = L_('1');
- ++intdig_no;
- }
-
- if (intdig_no + fracdig_no > dig_max)
- {
- wcp -= intdig_no + fracdig_no - dig_max;
- fracdig_no -= intdig_no + fracdig_no - dig_max;
- }
- }
- }
- }
- #endif
-
- while (fracdig_no > fracdig_min + added_zeros && *(wcp - 1) == L_('0'))
- {
- --wcp;
- --fracdig_no;
- }
-
- if (fracdig_no == 0 && !info->alt && *(wcp - 1) == decimalwc)
- --wcp;
- if (grouping)
- {
-
- if (intdig_no != intdig_max)
- ngroups = guess_grouping (intdig_no, grouping);
-
- wcp = group_number (wstartp, wcp, intdig_no, grouping, thousands_sepwc,
- ngroups);
- }
-
- if (type != 'f')
- {
- if (__builtin_expect (expsign != 0 && exponent == 4 && spec == 'g', 0))
- {
-
- assert (wcp >= wstartp + 1);
- assert (wstartp[0] == L_('1'));
- memcpy (wstartp, L_("0.0001"), 6 * sizeof (wchar_t));
- wstartp[1] = decimalwc;
- if (wcp >= wstartp + 2)
- {
- size_t cnt;
- for (cnt = 0; cnt < wcp - (wstartp + 2); cnt++)
- wstartp[6 + cnt] = L_('0');
- wcp += 4;
- }
- else
- wcp += 5;
- }
- else
- {
- *wcp++ = (wchar_t) type;
- *wcp++ = expsign ? L_('-') : L_('+');
-
- expscale = 10;
- while (expscale <= exponent)
- expscale *= 10;
- if (exponent < 10)
-
- *wcp++ = L_('0');
- else
- do
- {
- expscale /= 10;
- *wcp++ = L_('0') + (exponent / expscale);
- exponent %= expscale;
- }
- while (expscale > 10);
- *wcp++ = L_('0') + exponent;
- }
- }
-
- if (is_neg || info->showsign || info->space)
- --width;
- width -= wcp - wstartp;
- if (!info->left && info->pad != '0' && width > 0)
- PADN (info->pad, width);
- if (is_neg)
- outchar ('-');
- else if (info->showsign)
- outchar ('+');
- else if (info->space)
- outchar (' ');
- if (!info->left && info->pad == '0' && width > 0)
- PADN ('0', width);
- {
- char *buffer = NULL;
- char *buffer_end __attribute__((__unused__)) = NULL;
- char *cp = NULL;
- char *tmpptr;
- if (! wide)
- {
-
- size_t decimal_len;
- size_t thousands_sep_len;
- wchar_t *copywc;
- #ifdef USE_I18N_NUMBER_H
- size_t factor = (info->i18n
- ? nl_langinfo_wc (_NL_CTYPE_MB_CUR_MAX)
- : 1);
- #else
- size_t factor = 1;
- #endif
- decimal_len = strlen (decimal);
- if (thousands_sep == NULL)
- thousands_sep_len = 0;
- else
- thousands_sep_len = strlen (thousands_sep);
- size_t nbuffer = (2 + chars_needed * factor + decimal_len
- + ngroups * thousands_sep_len);
- if (__builtin_expect (buffer_malloced, 0))
- {
- buffer = (char *) malloc (nbuffer);
- if (buffer == NULL)
- {
-
- free (wbuffer);
- return -1;
- }
- }
- else
- buffer = (char *) alloca (nbuffer);
- buffer_end = buffer + nbuffer;
-
- for (cp = buffer, copywc = wstartp; copywc < wcp; ++copywc)
- if (*copywc == decimalwc)
- memcpy (cp, decimal, decimal_len), cp += decimal_len;
- else if (*copywc == thousands_sepwc)
- memcpy (cp, thousands_sep, thousands_sep_len), cp += thousands_sep_len;
- else
- *cp++ = (char) *copywc;
- }
- tmpptr = buffer;
- #if USE_I18N_NUMBER_H
- if (__builtin_expect (info->i18n, 0))
- {
- tmpptr = _i18n_number_rewrite (tmpptr, cp, buffer_end);
- cp = buffer_end;
- assert ((uintptr_t) buffer <= (uintptr_t) tmpptr);
- assert ((uintptr_t) tmpptr < (uintptr_t) buffer_end);
- }
- #endif
- PRINT (tmpptr, wstartp, wide ? wcp - wstartp : cp - tmpptr);
-
- if (__builtin_expect (buffer_malloced, 0))
- {
- free (buffer);
- free (wbuffer);
- }
- }
- if (info->left && width > 0)
- PADN (info->pad, width);
- }
- return done;
- }
- static unsigned int
- guess_grouping (unsigned int intdig_max, const char *grouping)
- {
- unsigned int groups;
-
- if (*grouping == CHAR_MAX || *grouping <= 0)
-
- return 0;
- groups = 0;
- while (intdig_max > (unsigned int) *grouping)
- {
- ++groups;
- intdig_max -= *grouping++;
- if (*grouping == CHAR_MAX
- #if CHAR_MIN < 0
- || *grouping < 0
- #endif
- )
-
- break;
- else if (*grouping == 0)
- {
-
- groups += (intdig_max - 1) / grouping[-1];
- break;
- }
- }
- return groups;
- }
- static wchar_t *
- group_number (wchar_t *buf, wchar_t *bufend, unsigned int intdig_no,
- const char *grouping, wchar_t thousands_sep, int ngroups)
- {
- wchar_t *p;
- if (ngroups == 0)
- return bufend;
-
- memmove (buf + intdig_no + ngroups, buf + intdig_no,
- (bufend - (buf + intdig_no)) * sizeof (wchar_t));
- p = buf + intdig_no + ngroups - 1;
- do
- {
- unsigned int len = *grouping++;
- do
- *p-- = buf[--intdig_no];
- while (--len > 0);
- *p-- = thousands_sep;
- if (*grouping == CHAR_MAX
- #if CHAR_MIN < 0
- || *grouping < 0
- #endif
- )
-
- break;
- else if (*grouping == 0)
-
- --grouping;
- } while (intdig_no > (unsigned int) *grouping);
-
- do
- *p-- = buf[--intdig_no];
- while (p > buf);
- return bufend + ngroups;
- }
|