1 #include "stdio_impl.h"
12 /* Some useful macros */
14 #define MAX(a,b) ((a)>(b) ? (a) : (b))
15 #define MIN(a,b) ((a)<(b) ? (a) : (b))
17 /* Convenient bit representation for modifier flags, which all fall
18 * within 31 codepoints of the space character. */
20 #define ALT_FORM (1U<<'#'-' ')
21 #define ZERO_PAD (1U<<'0'-' ')
22 #define LEFT_ADJ (1U<<'-'-' ')
23 #define PAD_POS (1U<<' '-' ')
24 #define MARK_POS (1U<<'+'-' ')
25 #define GROUPED (1U<<'\''-' ')
27 #define FLAGMASK (ALT_FORM|ZERO_PAD|LEFT_ADJ|PAD_POS|MARK_POS|GROUPED)
29 #if UINT_MAX == ULONG_MAX
33 #if SIZE_MAX != ULONG_MAX || UINTMAX_MAX != ULLONG_MAX
37 /* State machine to accept length modifiers + conversion specifiers.
38 * Result is 0 on failure, or an argument type to pop on success. */
41 BARE, LPRE, LLPRE, HPRE, HHPRE, BIGLPRE,
44 PTR, INT, UINT, ULLONG,
51 SHORT, USHORT, CHAR, UCHAR,
53 LLONG, SIZET, IMAX, UMAX, PDIFF, UIPTR,
67 #define S(x) [(x)-'A']
69 static const unsigned char states[]['z'-'A'+1] = {
71 S('d') = INT, S('i') = INT,
72 S('o') = UINT, S('u') = UINT, S('x') = UINT, S('X') = UINT,
73 S('e') = DBL, S('f') = DBL, S('g') = DBL, S('a') = DBL,
74 S('E') = DBL, S('F') = DBL, S('G') = DBL, S('A') = DBL,
75 S('c') = CHAR, S('C') = INT,
76 S('s') = PTR, S('S') = PTR, S('p') = UIPTR, S('n') = PTR,
78 S('l') = LPRE, S('h') = HPRE, S('L') = BIGLPRE,
79 S('z') = ZTPRE, S('j') = JPRE, S('t') = ZTPRE,
80 }, { /* 1: l-prefixed */
81 S('d') = LONG, S('i') = LONG,
82 S('o') = ULONG, S('u') = ULONG, S('x') = ULONG, S('X') = ULONG,
83 S('e') = DBL, S('f') = DBL, S('g') = DBL, S('a') = DBL,
84 S('E') = DBL, S('F') = DBL, S('G') = DBL, S('A') = DBL,
85 S('c') = INT, S('s') = PTR, S('n') = PTR,
87 }, { /* 2: ll-prefixed */
88 S('d') = LLONG, S('i') = LLONG,
89 S('o') = ULLONG, S('u') = ULLONG,
90 S('x') = ULLONG, S('X') = ULLONG,
92 }, { /* 3: h-prefixed */
93 S('d') = SHORT, S('i') = SHORT,
94 S('o') = USHORT, S('u') = USHORT,
95 S('x') = USHORT, S('X') = USHORT,
98 }, { /* 4: hh-prefixed */
99 S('d') = CHAR, S('i') = CHAR,
100 S('o') = UCHAR, S('u') = UCHAR,
101 S('x') = UCHAR, S('X') = UCHAR,
103 }, { /* 5: L-prefixed */
104 S('e') = LDBL, S('f') = LDBL, S('g') = LDBL, S('a') = LDBL,
105 S('E') = LDBL, S('F') = LDBL, S('G') = LDBL, S('A') = LDBL,
107 }, { /* 6: z- or t-prefixed (assumed to be same size) */
108 S('d') = PDIFF, S('i') = PDIFF,
109 S('o') = SIZET, S('u') = SIZET,
110 S('x') = SIZET, S('X') = SIZET,
112 }, { /* 7: j-prefixed */
113 S('d') = IMAX, S('i') = IMAX,
114 S('o') = UMAX, S('u') = UMAX,
115 S('x') = UMAX, S('X') = UMAX,
120 #define OOB(x) ((unsigned)(x)-'A' > 'z'-'A')
129 static void pop_arg(union arg *arg, int type, va_list *ap)
131 /* Give the compiler a hint for optimizing the switch. */
132 if ((unsigned)type > MAXSTATE) return;
134 case PTR: arg->p = va_arg(*ap, void *);
135 break; case INT: arg->i = va_arg(*ap, int);
136 break; case UINT: arg->i = va_arg(*ap, unsigned int);
138 break; case LONG: arg->i = va_arg(*ap, long);
139 break; case ULONG: arg->i = va_arg(*ap, unsigned long);
141 break; case ULLONG: arg->i = va_arg(*ap, unsigned long long);
142 break; case SHORT: arg->i = (short)va_arg(*ap, int);
143 break; case USHORT: arg->i = (unsigned short)va_arg(*ap, int);
144 break; case CHAR: arg->i = (signed char)va_arg(*ap, int);
145 break; case UCHAR: arg->i = (unsigned char)va_arg(*ap, int);
147 break; case LLONG: arg->i = va_arg(*ap, long long);
148 break; case SIZET: arg->i = va_arg(*ap, size_t);
149 break; case IMAX: arg->i = va_arg(*ap, intmax_t);
150 break; case UMAX: arg->i = va_arg(*ap, uintmax_t);
151 break; case PDIFF: arg->i = va_arg(*ap, ptrdiff_t);
152 break; case UIPTR: arg->i = (uintptr_t)va_arg(*ap, void *);
154 break; case DBL: arg->f = va_arg(*ap, double);
155 break; case LDBL: arg->f = va_arg(*ap, long double);
159 static void out(FILE *f, const char *s, size_t l)
161 if (!(f->flags & F_ERR)) __fwritex((void *)s, l, f);
164 static void pad(FILE *f, char c, int w, int l, int fl)
167 if (fl & (LEFT_ADJ | ZERO_PAD) || l >= w) return;
169 memset(pad, c, l>sizeof pad ? sizeof pad : l);
170 for (; l >= sizeof pad; l -= sizeof pad)
171 out(f, pad, sizeof pad);
175 static const char xdigits[16] = {
179 static char *fmt_x(uintmax_t x, char *s, int lower)
181 for (; x; x>>=4) *--s = xdigits[(x&15)]|lower;
185 static char *fmt_o(uintmax_t x, char *s)
187 for (; x; x>>=3) *--s = '0' + (x&7);
191 static char *fmt_u(uintmax_t x, char *s)
194 for ( ; x>ULONG_MAX; x/=10) *--s = '0' + x%10;
195 for (y=x; y; y/=10) *--s = '0' + y%10;
199 /* Do not override this check. The floating point printing code below
200 * depends on the float.h constants being right. If they are wrong, it
201 * may overflow the stack. */
202 #if LDBL_MANT_DIG == 53
203 typedef char compiler_defines_long_double_incorrectly[9-(int)sizeof(long double)];
206 static int fmt_fp(FILE *f, long double y, int w, int p, int fl, int t)
208 uint32_t big[(LDBL_MANT_DIG+28)/29 + 1 // mantissa expansion
209 + (LDBL_MAX_EXP+LDBL_MANT_DIG+28+8)/9]; // exponent expansion
210 uint32_t *a, *d, *r, *z;
211 int e2=0, e, i, j, l;
212 char buf[9+LDBL_MANT_DIG/4], *s;
213 const char *prefix="-0X+0X 0X-0x+0x 0x";
215 char ebuf0[3*sizeof(int)], *ebuf=&ebuf0[3*sizeof(int)], *estr;
220 } else if (fl & MARK_POS) {
222 } else if (fl & PAD_POS) {
224 } else prefix++, pl=0;
227 char *s = (t&32)?"inf":"INF";
228 if (y!=y) s=(t&32)?"nan":"NAN";
229 pad(f, ' ', w, 3+pl, fl&~ZERO_PAD);
232 pad(f, ' ', w, 3+pl, fl^LEFT_ADJ);
236 y = frexpl(y, &e2) * 2;
240 long double round = 8.0;
243 if (t&32) prefix += 9;
246 if (p<0 || p>=LDBL_MANT_DIG/4-1) re=0;
247 else re=LDBL_MANT_DIG/4-1-p;
250 while (re--) round*=16;
262 estr=fmt_u(e2<0 ? -e2 : e2, ebuf);
263 if (estr==ebuf) *--estr='0';
264 *--estr = (e2<0 ? '-' : '+');
265 *--estr = t+('p'-'a');
270 *s++=xdigits[x]|(t&32);
272 if (s-buf==1 && (y||p>0||(fl&ALT_FORM))) *s++='.';
275 if (p > INT_MAX-2-(ebuf-estr)-pl)
277 if (p && s-buf-2 < p)
278 l = (p+2) + (ebuf-estr);
280 l = (s-buf) + (ebuf-estr);
282 pad(f, ' ', w, pl+l, fl);
284 pad(f, '0', w, pl+l, fl^ZERO_PAD);
286 pad(f, '0', l-(ebuf-estr)-(s-buf), 0, 0);
287 out(f, estr, ebuf-estr);
288 pad(f, ' ', w, pl+l, fl^LEFT_ADJ);
293 if (y) y *= 0x1p28, e2-=28;
296 else a=r=z=big+sizeof(big)/sizeof(*big) - LDBL_MANT_DIG - 1;
300 y = 1000000000*(y-*z++);
306 for (d=z-1; d>=a; d--) {
307 uint64_t x = ((uint64_t)*d<<sh)+carry;
309 carry = x / 1000000000;
311 if (carry) *--a = carry;
312 while (z>a && !z[-1]) z--;
316 uint32_t carry=0, *b;
317 int sh=MIN(9,-e2), need=1+(p+LDBL_MANT_DIG/3U+8)/9;
318 for (d=a; d<z; d++) {
319 uint32_t rm = *d & (1<<sh)-1;
320 *d = (*d>>sh) + carry;
321 carry = (1000000000>>sh) * rm;
324 if (carry) *z++ = carry;
325 /* Avoid (slow!) computation past requested precision */
326 b = (t|32)=='f' ? r : a;
327 if (z-b > need) z = b+need;
331 if (a<z) for (i=10, e=9*(r-a); *a>=i; i*=10, e++);
334 /* Perform rounding: j is precision after the radix (possibly neg) */
335 j = p - ((t|32)!='f')*e - ((t|32)=='g' && p);
338 /* We avoid C's broken division of negative numbers */
339 d = r + 1 + ((j+9*LDBL_MAX_EXP)/9 - LDBL_MAX_EXP);
342 for (i=10, j++; j<9; i*=10, j++);
344 /* Are there any significant digits past j? */
346 long double round = 2/LDBL_EPSILON;
348 if ((*d/i & 1) || (i==1000000000 && d>a && (d[-1]&1)))
350 if (x<i/2) small=0x0.8p0;
351 else if (x==i/2 && d+1==z) small=0x1.0p0;
353 if (pl && *prefix=='-') round*=-1, small*=-1;
355 /* Decide whether to round by probing round+small */
356 if (round+small != round) {
358 while (*d > 999999999) {
363 for (i=10, e=9*(r-a); *a>=i; i*=10, e++);
368 for (; z>a && !z[-1]; z--);
379 if (!(fl&ALT_FORM)) {
380 /* Count trailing zeros in last place */
381 if (z>a && z[-1]) for (i=10, j=0; z[-1]%i==0; i*=10, j++);
384 p = MIN(p,MAX(0,9*(z-r-1)-j));
386 p = MIN(p,MAX(0,9*(z-r-1)+e-j));
389 if (p > INT_MAX-1-(p || (fl&ALT_FORM)))
391 l = 1 + p + (p || (fl&ALT_FORM));
393 if (e > INT_MAX-l) return -1;
396 estr=fmt_u(e<0 ? -e : e, ebuf);
397 while(ebuf-estr<2) *--estr='0';
398 *--estr = (e<0 ? '-' : '+');
400 if (ebuf-estr > INT_MAX-l) return -1;
404 if (l > INT_MAX-pl) return -1;
405 pad(f, ' ', w, pl+l, fl);
407 pad(f, '0', w, pl+l, fl^ZERO_PAD);
411 for (d=a; d<=r; d++) {
412 char *s = fmt_u(*d, buf+9);
413 if (d!=a) while (s>buf) *--s='0';
414 else if (s==buf+9) *--s='0';
417 if (p || (fl&ALT_FORM)) out(f, ".", 1);
418 for (; d<z && p>0; d++, p-=9) {
419 char *s = fmt_u(*d, buf+9);
420 while (s>buf) *--s='0';
423 pad(f, '0', p+9, 9, 0);
426 for (d=a; d<z && p>=0; d++) {
427 char *s = fmt_u(*d, buf+9);
428 if (s==buf+9) *--s='0';
429 if (d!=a) while (s>buf) *--s='0';
432 if (p>0||(fl&ALT_FORM)) out(f, ".", 1);
434 out(f, s, MIN(buf+9-s, p));
437 pad(f, '0', p+18, 18, 0);
438 out(f, estr, ebuf-estr);
441 pad(f, ' ', w, pl+l, fl^LEFT_ADJ);
446 static int getint(char **s) {
448 for (i=0; isdigit(**s); (*s)++) {
449 if (i > INT_MAX/10U || **s-'0' > INT_MAX-10*i) i = -1;
450 else i = 10*i + (**s-'0');
455 static int printf_core(FILE *f, const char *fmt, va_list *ap, union arg *nl_arg, int *nl_type)
457 char *a, *z, *s=(char *)fmt;
465 char buf[sizeof(uintmax_t)*3+3+LDBL_MANT_DIG/4];
472 /* This error is only specified for snprintf, but since it's
473 * unspecified for other forms, do the same. Stop immediately
474 * on overflow; otherwise %n could produce wrong results. */
475 if (l > INT_MAX - cnt) goto overflow;
477 /* Update output count, end loop when fmt is exhausted */
481 /* Handle literal text and %% format specifiers */
482 for (a=s; *s && *s!='%'; s++);
483 for (z=s; s[0]=='%' && s[1]=='%'; z++, s+=2);
484 if (z-a > INT_MAX-cnt) goto overflow;
489 if (isdigit(s[1]) && s[2]=='$') {
498 /* Read modifier flags */
499 for (fl=0; (unsigned)*s-' '<32 && (FLAGMASK&(1U<<*s-' ')); s++)
502 /* Read field width */
504 if (isdigit(s[1]) && s[2]=='$') {
506 nl_type[s[1]-'0'] = INT;
507 w = nl_arg[s[1]-'0'].i;
510 w = f ? va_arg(*ap, int) : 0;
513 if (w<0) fl|=LEFT_ADJ, w=-w;
514 } else if ((w=getint(&s))<0) goto overflow;
517 if (*s=='.' && s[1]=='*') {
518 if (isdigit(s[2]) && s[3]=='$') {
519 nl_type[s[2]-'0'] = INT;
520 p = nl_arg[s[2]-'0'].i;
523 p = f ? va_arg(*ap, int) : 0;
527 } else if (*s=='.') {
536 /* Format specifier state machine */
539 if (OOB(*s)) goto inval;
541 st=states[st]S(*s++);
545 /* Check validity of argument type (nl/normal) */
547 if (argpos>=0) goto inval;
549 if (argpos>=0) nl_type[argpos]=st, arg=nl_arg[argpos];
550 else if (f) pop_arg(&arg, st, ap);
556 z = buf + sizeof(buf);
561 /* Transform ls,lc -> S,C */
562 if (ps && (t&15)==3) t&=~32;
564 /* - and 0 flags are mutually exclusive */
565 if (fl & LEFT_ADJ) fl &= ~ZERO_PAD;
570 case BARE: *(int *)arg.p = cnt; break;
571 case LPRE: *(long *)arg.p = cnt; break;
572 case LLPRE: *(long long *)arg.p = cnt; break;
573 case HPRE: *(unsigned short *)arg.p = cnt; break;
574 case HHPRE: *(unsigned char *)arg.p = cnt; break;
575 case ZTPRE: *(size_t *)arg.p = cnt; break;
576 case JPRE: *(uintmax_t *)arg.p = cnt; break;
580 p = MAX(p, 2*sizeof(void*));
584 a = fmt_x(arg.i, z, t&32);
585 if (arg.i && (fl & ALT_FORM)) prefix+=(t>>4), pl=2;
589 if ((fl&ALT_FORM) && p<z-a+1) prefix+=5, pl=1;
593 if (arg.i>INTMAX_MAX) {
595 } else if (fl & MARK_POS) {
597 } else if (fl & PAD_POS) {
603 if (xp && p<0) goto overflow;
604 if (p>=0) fl &= ~ZERO_PAD;
609 p = MAX(p, z-a + !arg.i);
616 if (1) a = strerror(errno); else
618 a = arg.p ? arg.p : "(null)";
619 z = a + strnlen(a, p<0 ? INT_MAX : p);
620 if (p<0 && *z) goto overflow;
631 for (i=l=0; i<p && *ws && (l=wctomb(mb, *ws++))>=0 && l<=p-i; i+=l);
633 if (i > INT_MAX) goto overflow;
635 pad(f, ' ', w, p, fl);
637 for (i=0; i<0U+p && *ws && i+(l=wctomb(mb, *ws++))<=p; i+=l)
639 pad(f, ' ', w, p, fl^LEFT_ADJ);
642 case 'e': case 'f': case 'g': case 'a':
643 case 'E': case 'F': case 'G': case 'A':
644 if (xp && p<0) goto overflow;
645 l = fmt_fp(f, arg.f, w, p, fl, t);
646 if (l<0) goto overflow;
650 if (p < z-a) p = z-a;
651 if (p > INT_MAX-pl) goto overflow;
652 if (w < pl+p) w = pl+p;
653 if (w > INT_MAX-cnt) goto overflow;
655 pad(f, ' ', w, pl+p, fl);
657 pad(f, '0', w, pl+p, fl^ZERO_PAD);
658 pad(f, '0', p, z-a, 0);
660 pad(f, ' ', w, pl+p, fl^LEFT_ADJ);
668 for (i=1; i<=NL_ARGMAX && nl_type[i]; i++)
669 pop_arg(nl_arg+i, nl_type[i], ap);
670 for (; i<=NL_ARGMAX && !nl_type[i]; i++);
671 if (i<=NL_ARGMAX) goto inval;
682 int vfprintf(FILE *restrict f, const char *restrict fmt, va_list ap)
685 int nl_type[NL_ARGMAX+1] = {0};
686 union arg nl_arg[NL_ARGMAX+1];
687 unsigned char internal_buf[80], *saved_buf = 0;
691 /* the copy allows passing va_list* even if va_list is an array */
693 if (printf_core(0, fmt, &ap2, nl_arg, nl_type) < 0) {
699 olderr = f->flags & F_ERR;
700 if (f->mode < 1) f->flags &= ~F_ERR;
703 f->wpos = f->wbase = f->buf = internal_buf;
704 f->buf_size = sizeof internal_buf;
705 f->wend = internal_buf + sizeof internal_buf;
707 ret = printf_core(f, fmt, &ap2, nl_arg, nl_type);
710 if (!f->wpos) ret = -1;
713 f->wpos = f->wbase = f->wend = 0;
715 if (f->flags & F_ERR) ret = -1;