Context Navigation

iostream.cfa@ bd47f35

Visit:

ADT arm-eh ast-experimental enum forall-pointer-decay jacob/cs343-translation new-ast-unique-expr pthread-emulation qualifiedEnum

Last change on this file since bd47f35 was 7cfef0d, checked in by Peter A. Buhr <pabuhr@…>, 5 years ago
move routines floor, ceiling, ceiling_div from bitmanip.hfa to math.hfa
Property mode set to `100644`
File size: 40.2 KB

Line
1	//
2	// Cforall Version 1.0.0 Copyright (C) 2015 University of Waterloo
3	//
4	// The contents of this file are covered under the licence agreement in the
5	// file "LICENCE" distributed with Cforall.
6	//
7	// iostream.cfa --
8	//
9	// Author : Peter A. Buhr
10	// Created On : Wed May 27 17:56:53 2015
11	// Last Modified By : Peter A. Buhr
12	// Last Modified On : Mon Aug 24 08:31:35 2020
13	// Update Count : 1130
14	//
15
16	#include "iostream.hfa"
17
18	#include <stdio.h>
19	#include <stdbool.h> // true/false
20	#include <stdint.h> // UINT64_MAX
21	#include <float.h> // DBL_DIG, LDBL_DIG
22	#include <complex.h> // creal, cimag
23	//#include <string.h> // strlen, strcmp, memcpy
24	extern "C" {
25	extern size_t strlen (const char *__s) __attribute__ ((__nothrow__ , __leaf__)) __attribute__ ((__pure__)) __attribute__ ((__nonnull__ (1)));
26	extern int strcmp (const char __s1, const char __s2) __attribute__ ((__nothrow__ , __leaf__)) __attribute__ ((__pure__)) __attribute__ ((__nonnull__ (1, 2)));
27	extern char strcpy (char __restrict __dest, const char *__restrict __src) __attribute__ ((__nothrow__ , __leaf__)) __attribute__ ((__nonnull__ (1, 2)));
28	extern void memcpy (void __restrict __dest, const void *__restrict __src, size_t __n) __attribute__ ((__nothrow__ , __leaf__)) __attribute__ ((__nonnull__ (1, 2)));
29	} // extern "C"
30
31	#include "math.hfa" // isfinite, floor, ceiling_div
32	#include "bitmanip.hfa" // high1
33
34
35	// ********************************* ostream *********************************
36
37
38	forall( dtype ostype \| ostream( ostype ) ) {
39	ostype & ?\|?( ostype & os, bool b ) {
40	if ( $sepPrt( os ) ) fmt( os, "%s", $sepGetCur( os ) );
41	fmt( os, "%s", b ? "true" : "false" );
42	return os;
43	} // ?\|?
44	void ?\|?( ostype & os, bool b ) {
45	(ostype &)(os \| b); ends( os );
46	} // ?\|?
47
48	ostype & ?\|?( ostype & os, char c ) {
49	fmt( os, "%c", c );
50	if ( c == '\n' ) $setNL( os, true );
51	return sepOff( os );
52	} // ?\|?
53	void ?\|?( ostype & os, char c ) {
54	(ostype &)(os \| c); ends( os );
55	} // ?\|?
56
57	ostype & ?\|?( ostype & os, signed char sc ) {
58	if ( $sepPrt( os ) ) fmt( os, "%s", $sepGetCur( os ) );
59	fmt( os, "%hhd", sc );
60	return os;
61	} // ?\|?
62	void ?\|?( ostype & os, signed char sc ) {
63	(ostype &)(os \| sc); ends( os );
64	} // ?\|?
65
66	ostype & ?\|?( ostype & os, unsigned char usc ) {
67	if ( $sepPrt( os ) ) fmt( os, "%s", $sepGetCur( os ) );
68	fmt( os, "%hhu", usc );
69	return os;
70	} // ?\|?
71	void ?\|?( ostype & os, unsigned char usc ) {
72	(ostype &)(os \| usc); ends( os );
73	} // ?\|?
74
75	ostype & ?\|?( ostype & os, short int si ) {
76	if ( $sepPrt( os ) ) fmt( os, "%s", $sepGetCur( os ) );
77	fmt( os, "%hd", si );
78	return os;
79	} // ?\|?
80	void & ?\|?( ostype & os, short int si ) {
81	(ostype &)(os \| si); ends( os );
82	} // ?\|?
83
84	ostype & ?\|?( ostype & os, unsigned short int usi ) {
85	if ( $sepPrt( os ) ) fmt( os, "%s", $sepGetCur( os ) );
86	fmt( os, "%hu", usi );
87	return os;
88	} // ?\|?
89	void & ?\|?( ostype & os, unsigned short int usi ) {
90	(ostype &)(os \| usi); ends( os );
91	} // ?\|?
92
93	ostype & ?\|?( ostype & os, int i ) {
94	if ( $sepPrt( os ) ) fmt( os, "%s", $sepGetCur( os ) );
95	fmt( os, "%d", i );
96	return os;
97	} // ?\|?
98	void & ?\|?( ostype & os, int i ) {
99	(ostype &)(os \| i); ends( os );
100	} // ?\|?
101
102	ostype & ?\|?( ostype & os, unsigned int ui ) {
103	if ( $sepPrt( os ) ) fmt( os, "%s", $sepGetCur( os ) );
104	fmt( os, "%u", ui );
105	return os;
106	} // ?\|?
107	void & ?\|?( ostype & os, unsigned int ui ) {
108	(ostype &)(os \| ui); ends( os );
109	} // ?\|?
110
111	ostype & ?\|?( ostype & os, long int li ) {
112	if ( $sepPrt( os ) ) fmt( os, "%s", $sepGetCur( os ) );
113	fmt( os, "%ld", li );
114	return os;
115	} // ?\|?
116	void & ?\|?( ostype & os, long int li ) {
117	(ostype &)(os \| li); ends( os );
118	} // ?\|?
119
120	ostype & ?\|?( ostype & os, unsigned long int uli ) {
121	if ( $sepPrt( os ) ) fmt( os, "%s", $sepGetCur( os ) );
122	fmt( os, "%lu", uli );
123	return os;
124	} // ?\|?
125	void & ?\|?( ostype & os, unsigned long int uli ) {
126	(ostype &)(os \| uli); ends( os );
127	} // ?\|?
128
129	ostype & ?\|?( ostype & os, long long int lli ) {
130	if ( $sepPrt( os ) ) fmt( os, "%s", $sepGetCur( os ) );
131	fmt( os, "%lld", lli );
132	return os;
133	} // ?\|?
134	void & ?\|?( ostype & os, long long int lli ) {
135	(ostype &)(os \| lli); ends( os );
136	} // ?\|?
137
138	ostype & ?\|?( ostype & os, unsigned long long int ulli ) {
139	if ( $sepPrt( os ) ) fmt( os, "%s", $sepGetCur( os ) );
140	fmt( os, "%llu", ulli );
141	return os;
142	} // ?\|?
143	void & ?\|?( ostype & os, unsigned long long int ulli ) {
144	(ostype &)(os \| ulli); ends( os );
145	} // ?\|?
146
147	#if defined( __SIZEOF_INT128__ )
148	// UINT64_MAX 18_446_744_073_709_551_615_ULL
149	#define P10_UINT64 10_000_000_000_000_000_000_ULL // 19 zeroes
150
151	static inline void base10_128( ostype & os, unsigned int128 val ) {
152	#if defined(__GNUC__) && __GNUC_PREREQ(7,0) // gcc version >= 7
153	if ( val > P10_UINT64 ) {
154	#else
155	if ( (uint64_t)(val >> 64) != 0 \|\| (uint64_t)val > P10_UINT64 ) { // patch gcc 5 & 6 -O3 bug
156	#endif // __GNUC_PREREQ(7,0)
157	base10_128( os, val / P10_UINT64 ); // recursive
158	fmt( os, "%.19lu", (uint64_t)(val % P10_UINT64) );
159	} else {
160	fmt( os, "%lu", (uint64_t)val );
161	} // if
162	} // base10_128
163
164	static inline void base10_128( ostype & os, int128 val ) {
165	if ( val < 0 ) {
166	fmt( os, "-" ); // leading negative sign
167	val = -val;
168	} // if
169	base10_128( os, (unsigned int128)val ); // print zero/positive value
170	} // base10_128
171
172	ostype & ?\|?( ostype & os, int128 llli ) {
173	if ( $sepPrt( os ) ) fmt( os, "%s", $sepGetCur( os ) );
174	base10_128( os, llli );
175	return os;
176	} // ?\|?
177	void & ?\|?( ostype & os, int128 llli ) {
178	(ostype &)(os \| llli); ends( os );
179	} // ?\|?
180
181	ostype & ?\|?( ostype & os, unsigned int128 ullli ) {
182	if ( $sepPrt( os ) ) fmt( os, "%s", $sepGetCur( os ) );
183	base10_128( os, ullli );
184	return os;
185	} // ?\|?
186	void & ?\|?( ostype & os, unsigned int128 ullli ) {
187	(ostype &)(os \| ullli); ends( os );
188	} // ?\|?
189	#endif // __SIZEOF_INT128__
190
191	#define PrintWithDP( os, format, val, ... ) \
192	{ \
193	enum { size = 48 }; \
194	char buf[size]; \
195	int len = snprintf( buf, size, format, ##__VA_ARGS__, val ); \
196	fmt( os, "%s", buf ); \
197	if ( isfinite( val ) ) { /* if number, always print decimal point */ \
198	for ( int i = 0;; i += 1 ) { \
199	if ( i == len ) { fmt( os, "." ); break; } \
200	if ( buf[i] == '.' ) break; \
201	} /* for */ \
202	} /* if */ \
203	}
204
205	ostype & ?\|?( ostype & os, float f ) {
206	if ( $sepPrt( os ) ) fmt( os, "%s", $sepGetCur( os ) );
207	PrintWithDP( os, "%g", f );
208	return os;
209	} // ?\|?
210	void & ?\|?( ostype & os, float f ) {
211	(ostype &)(os \| f); ends( os );
212	} // ?\|?
213
214	ostype & ?\|?( ostype & os, double d ) {
215	if ( $sepPrt( os ) ) fmt( os, "%s", $sepGetCur( os ) );
216	PrintWithDP( os, "%.*lg", d, DBL_DIG );
217	return os;
218	} // ?\|?
219	void & ?\|?( ostype & os, double d ) {
220	(ostype &)(os \| d); ends( os );
221	} // ?\|?
222
223	ostype & ?\|?( ostype & os, long double ld ) {
224	if ( $sepPrt( os ) ) fmt( os, "%s", $sepGetCur( os ) );
225	PrintWithDP( os, "%.*Lg", ld, LDBL_DIG );
226	return os;
227	} // ?\|?
228	void & ?\|?( ostype & os, long double ld ) {
229	(ostype &)(os \| ld); ends( os );
230	} // ?\|?
231
232	ostype & ?\|?( ostype & os, float _Complex fc ) {
233	if ( $sepPrt( os ) ) fmt( os, "%s", $sepGetCur( os ) );
234	// os \| crealf( fc ) \| nonl;
235	PrintWithDP( os, "%g", crealf( fc ) );
236	PrintWithDP( os, "%+g", cimagf( fc ) );
237	fmt( os, "i" );
238	return os;
239	} // ?\|?
240	void & ?\|?( ostype & os, float _Complex fc ) {
241	(ostype &)(os \| fc); ends( os );
242	} // ?\|?
243
244	ostype & ?\|?( ostype & os, double _Complex dc ) {
245	if ( $sepPrt( os ) ) fmt( os, "%s", $sepGetCur( os ) );
246	// os \| creal( dc ) \| nonl;
247	PrintWithDP( os, "%.*lg", creal( dc ), DBL_DIG );
248	PrintWithDP( os, "%+.*lg", cimag( dc ), DBL_DIG );
249	fmt( os, "i" );
250	return os;
251	} // ?\|?
252	void & ?\|?( ostype & os, double _Complex dc ) {
253	(ostype &)(os \| dc); ends( os );
254	} // ?\|?
255
256	ostype & ?\|?( ostype & os, long double _Complex ldc ) {
257	if ( $sepPrt( os ) ) fmt( os, "%s", $sepGetCur( os ) );
258	// os \| creall( ldc ) \|\| nonl;
259	PrintWithDP( os, "%.*Lg", creall( ldc ), LDBL_DIG );
260	PrintWithDP( os, "%+.*Lg", cimagl( ldc ), LDBL_DIG );
261	fmt( os, "i" );
262	return os;
263	} // ?\|?
264	void & ?\|?( ostype & os, long double _Complex ldc ) {
265	(ostype &)(os \| ldc); ends( os );
266	} // ?\|?
267
268	ostype & ?\|?( ostype & os, const char str[] ) {
269	enum { Open = 1, Close, OpenClose };
270	static const unsigned char mask[256] @= {
271	// opening delimiters, no space after
272	['('] : Open, ['['] : Open, ['{'] : Open,
273	['='] : Open, ['$'] : Open, [(unsigned char)'£'] : Open, [(unsigned char)'¥'] : Open,
274	[(unsigned char)'¡'] : Open, [(unsigned char)'¿'] : Open, [(unsigned char)'«'] : Open,
275	// closing delimiters, no space before
276	[','] : Close, ['.'] : Close, [';'] : Close, ['!'] : Close, ['?'] : Close,
277	['%'] : Close, [(unsigned char)'¢'] : Close, [(unsigned char)'»'] : Close,
278	[')'] : Close, [']'] : Close, ['}'] : Close,
279	// opening-closing delimiters, no space before or after
280	['\''] : OpenClose, ['`'] : OpenClose, ['"'] : OpenClose, [':'] : OpenClose,
281	[' '] : OpenClose, ['\f'] : OpenClose, ['\n'] : OpenClose, ['\r'] : OpenClose, ['\t'] : OpenClose, ['\v'] : OpenClose, // isspace
282	}; // mask
283
284	if ( str[0] == '\0' ) { sepOff( os ); return os; } // null string => no separator
285
286	// first character IS NOT spacing or closing punctuation => add left separator
287	unsigned char ch = str[0]; // must make unsigned
288	if ( $sepPrt( os ) && mask[ ch ] != Close && mask[ ch ] != OpenClose ) {
289	fmt( os, "%s", $sepGetCur( os ) );
290	} // if
291
292	// if string starts line, must reset to determine open state because separator is off
293	$sepReset( os ); // reset separator
294
295	// last character IS spacing or opening punctuation => turn off separator for next item
296	size_t len = strlen( str );
297	ch = str[len - 1]; // must make unsigned
298	if ( $sepPrt( os ) && mask[ ch ] != Open && mask[ ch ] != OpenClose ) {
299	sepOn( os );
300	} else {
301	sepOff( os );
302	} // if
303	if ( ch == '\n' ) $setNL( os, true ); // check AFTER $sepPrt call above as it resets NL flag
304	return write( os, str, len );
305	} // ?\|?
306
307	void ?\|?( ostype & os, const char str[] ) {
308	(ostype &)(os \| str); ends( os );
309	} // ?\|?
310
311	// ostype & ?\|?( ostype & os, const char16_t * str ) {
312	// if ( $sepPrt( os ) ) fmt( os, "%s", $sepGetCur( os ) );
313	// fmt( os, "%ls", str );
314	// return os;
315	// } // ?\|?
316
317	// #if ! ( __ARM_ARCH_ISA_ARM == 1 && __ARM_32BIT_STATE == 1 ) // char32_t == wchar_t => ambiguous
318	// ostype & ?\|?( ostype & os, const char32_t * str ) {
319	// if ( $sepPrt( os ) ) fmt( os, "%s", $sepGetCur( os ) );
320	// fmt( os, "%ls", str );
321	// return os;
322	// } // ?\|?
323	// #endif // ! ( __ARM_ARCH_ISA_ARM == 1 && __ARM_32BIT_STATE == 1 )
324
325	// ostype & ?\|?( ostype & os, const wchar_t * str ) {
326	// if ( $sepPrt( os ) ) fmt( os, "%s", $sepGetCur( os ) );
327	// fmt( os, "%ls", str );
328	// return os;
329	// } // ?\|?
330
331	ostype & ?\|?( ostype & os, const void * p ) {
332	if ( $sepPrt( os ) ) fmt( os, "%s", $sepGetCur( os ) );
333	fmt( os, "%p", p );
334	return os;
335	} // ?\|?
336	void ?\|?( ostype & os, const void * p ) {
337	(ostype &)(os \| p); ends( os );
338	} // ?\|?
339
340	// manipulators
341	ostype & ?\|?( ostype & os, ostype & (* manip)( ostype & ) ) {
342	(ostype &)(manip( os ));
343	return os;
344	} // ?\|?
345	void ?\|?( ostype & os, ostype & (* manip)( ostype & ) ) {
346	(ostype &)(manip( os ));
347	if ( $getPrt( os ) ) ends( os ); // something printed ?
348	$setPrt( os, false ); // turn off
349	} // ?\|?
350
351	ostype & sep( ostype & os ) {
352	return (ostype &)(os \| sepGet( os ));
353	} // sep
354
355	ostype & sepTuple( ostype & os ) {
356	return os \| sepGetTuple( os );
357	} // sepTuple
358
359	ostype & nl( ostype & os ) {
360	(ostype &)(os \| '\n');
361	$setPrt( os, false ); // turn off
362	$setNL( os, true );
363	flush( os );
364	return sepOff( os ); // prepare for next line
365	} // nl
366
367	ostype & nonl( ostype & os ) {
368	$setPrt( os, false ); // turn off
369	return os;
370	} // nonl
371
372	ostype & sepOn( ostype & os ) {
373	sepOn( os ); // call void returning
374	return os;
375	} // sepOn
376
377	ostype & sepOff( ostype & os ) {
378	sepOff( os ); // call void returning
379	return os;
380	} // sepOff
381
382	ostype & sepEnable( ostype & os ) {
383	sepEnable( os ); // call void returning
384	return os;
385	} // sepEnable
386
387	ostype & sepDisable( ostype & os ) {
388	sepDisable( os ); // call void returning
389	return os;
390	} // sepDisable
391
392	ostype & nlOn( ostype & os ) {
393	nlOn( os ); // call void returning
394	return os;
395	} // nlOn
396
397	ostype & nlOff( ostype & os ) {
398	nlOff( os ); // call void returning
399	return os;
400	} // nlOff
401	} // distribution
402
403	// tuples
404	forall( dtype ostype, otype T, ttype Params \| writeable( T, ostype ) \| { ostype & ?\|?( ostype &, Params ); } ) {
405	ostype & ?\|?( ostype & os, T arg, Params rest ) {
406	(ostype &)(os \| arg); // print first argument
407	$sepSetCur( os, sepGetTuple( os ) ); // switch to tuple separator
408	(ostype &)(os \| rest); // print remaining arguments
409	$sepSetCur( os, sepGet( os ) ); // switch to regular separator
410	return os;
411	} // ?\|?
412	void ?\|?( ostype & os, T arg, Params rest ) {
413	// (ostype &)(?\|?( os, arg, rest )); ends( os );
414	(ostype &)(os \| arg); // print first argument
415	$sepSetCur( os, sepGetTuple( os ) ); // switch to tuple separator
416	(ostype &)(os \| rest); // print remaining arguments
417	$sepSetCur( os, sepGet( os ) ); // switch to regular separator
418	ends( os );
419	} // ?\|?
420	} // distribution
421
422	// writes the range [begin, end) to the given stream
423	forall( dtype ostype, otype elt_type \| writeable( elt_type, ostype ), otype iterator_type \| iterator( iterator_type, elt_type ) ) {
424	void write( iterator_type begin, iterator_type end, ostype & os ) {
425	void print( elt_type i ) { os \| i; }
426	for_each( begin, end, print );
427	} // ?\|?
428
429	void write_reverse( iterator_type begin, iterator_type end, ostype & os ) {
430	void print( elt_type i ) { os \| i; }
431	for_each_reverse( begin, end, print );
432	} // ?\|?
433	} // distribution
434
435	// ********************************* manipulators *********************************
436
437	// ********************************* integral *********************************
438
439	static const char * shortbin[] = { "0", "1", "10", "11", "100", "101", "110", "111", "1000", "1001", "1010", "1011", "1100", "1101", "1110", "1111" };
440	static const char * longbin[] = { "0000", "0001", "0010", "0011", "0100", "0101", "0110", "0111", "1000", "1001", "1010", "1011", "1100", "1101", "1110", "1111" };
441
442	// Default prefix for non-decimal prints is 0b, 0, 0x.
443	#define IntegralFMTImpl( T, IFMTNP, IFMTP ) \
444	forall( dtype ostype \| ostream( ostype ) ) { \
445	ostype & ?\|?( ostype & os, _Ostream_Manip(T) f ) { \
446	if ( $sepPrt( os ) ) fmt( os, "%s", $sepGetCur( os ) ); \
447	\
448	if ( f.base == 'b' \|\| f.base == 'B' ) { /* bespoke binary format */ \
449	int bits = high1( f.val ); /* position of most significant bit */ \
450	if ( bits == 0 ) bits = 1; /* 0 value => force one bit to print */ \
451	int spaces; \
452	if ( ! f.flags.left ) { /* right justified ? */ \
453	/* Note, base prefix then zero padding or spacing then prefix. */ \
454	if ( f.flags.pc ) { \
455	spaces = f.wd - f.pc; \
456	if ( ! f.flags.nobsdp ) { spaces -= 2; } /* base prefix takes space */ \
457	if ( spaces > 0 ) fmt( os, "%s", spaces, " " ); / space pad */ \
458	if ( ! f.flags.nobsdp ) { fmt( os, "0%c", f.base ); } \
459	spaces = f.pc - bits; \
460	if ( spaces > 0 ) fmt( os, "%0d", spaces, 0 ); / zero pad */ \
461	} else { \
462	spaces = f.wd - bits; \
463	if ( ! f.flags.nobsdp ) { spaces -= 2; } /* base prefix takes space */ \
464	if ( f.flags.pad0 ) { \
465	if ( ! f.flags.nobsdp ) { fmt( os, "0%c", f.base ); } \
466	if ( spaces > 0 ) fmt( os, "%0d", spaces, 0 ); / zero pad */ \
467	} else { \
468	if ( spaces > 0 ) fmt( os, "%s", spaces, " " ); / space pad */ \
469	if ( ! f.flags.nobsdp ) { fmt( os, "0%c", f.base ); } \
470	} /* if */ \
471	} /* if */ \
472	} else { \
473	if ( ! f.flags.nobsdp ) fmt( os, "0%c", f.base ); \
474	if ( f.flags.pc ) { \
475	spaces = f.pc - bits; \
476	if ( spaces > 0 ) fmt( os, "%0d", spaces, 0 ); / zero pad */ \
477	spaces = f.wd - f.pc; \
478	} else { /* pad0 flag ignored with left flag */ \
479	spaces = f.wd - bits; \
480	} /* if */ \
481	if ( ! f.flags.nobsdp ) { spaces -= 2; } /* base prefix takes space */ \
482	} /* if */ \
483	int shift = floor( bits - 1, 4 ); \
484	typeof( f.val ) temp = f.val; \
485	fmt( os, "%s", shortbin[(temp >> shift) & 0xf] ); \
486	for () { \
487	shift -= 4; \
488	if ( shift < 0 ) break; \
489	temp = f.val; \
490	fmt( os, "%s", longbin[(temp >> shift) & 0xf] ); \
491	} /* for */ \
492	if ( f.flags.left && spaces > 0 ) fmt( os, "%*s", spaces, " " ); \
493	return os; \
494	} /* if */ \
495	\
496	char fmtstr[sizeof(IFMTP)]; /* sizeof includes '\0' */ \
497	if ( ! f.flags.pc ) memcpy( &fmtstr, IFMTNP, sizeof(IFMTNP) ); \
498	else memcpy( &fmtstr, IFMTP, sizeof(IFMTP) ); \
499	int star = 4; /* position before first '' / \
500	\
501	/* Insert flags into spaces before '', from right to left. / \
502	if ( ! f.flags.nobsdp ) { fmtstr[star] = '#'; star -= 1; } \
503	if ( f.flags.left ) { fmtstr[star] = '-'; star -= 1; } \
504	if ( f.flags.sign ) { fmtstr[star] = '+'; star -= 1; } \
505	if ( f.flags.pad0 && ! f.flags.pc ) { fmtstr[star] = '0'; star -= 1; } \
506	fmtstr[star] = '%'; \
507	\
508	if ( ! f.flags.pc ) { /* no precision */ \
509	fmtstr[sizeof(IFMTNP)-2] = f.base; /* sizeof includes '\0' */ \
510	/* printf( "%s %c\n", &fmtstr[star], f.base ); */ \
511	fmt( os, &fmtstr[star], f.wd, f.val ); \
512	} else { /* precision */ \
513	fmtstr[sizeof(IFMTP)-2] = f.base; /* sizeof includes '\0' */ \
514	/* printf( "%s %c\n", &fmtstr[star], f.base ); */ \
515	fmt( os, &fmtstr[star], f.wd, f.pc, f.val ); \
516	} /* if */ \
517	return os; \
518	} /* ?\|? */ \
519	void ?\|?( ostype & os, _Ostream_Manip(T) f ) { (ostype &)(os \| f); ends( os ); } \
520	} // distribution
521
522	IntegralFMTImpl( signed char, "% hh ", "% .*hh " )
523	IntegralFMTImpl( unsigned char, "% hh ", "% .*hh " )
524	IntegralFMTImpl( signed short int, "% h ", "% .*h " )
525	IntegralFMTImpl( unsigned short int, "% h ", "% .*h " )
526	IntegralFMTImpl( signed int, "% * ", "% . " )
527	IntegralFMTImpl( unsigned int, "% * ", "% . " )
528	IntegralFMTImpl( signed long int, "% l ", "% .*l " )
529	IntegralFMTImpl( unsigned long int, "% l ", "% .*l " )
530	IntegralFMTImpl( signed long long int, "% ll ", "% .*ll " )
531	IntegralFMTImpl( unsigned long long int, "% ll ", "% .*ll " )
532
533	#if 0
534	#if defined( __SIZEOF_INT128__ )
535	// Default prefix for non-decimal prints is 0b, 0, 0x.
536	#define IntegralFMTImpl128( T, SIGNED, CODE, IFMTNP, IFMTP ) \
537	forall( dtype ostype \| ostream( ostype ) ) \
538	static void base10_128( ostype & os, _Ostream_Manip(T) f ) { \
539	if ( f.val > UINT64_MAX ) { \
540	unsigned long long int lsig = f.val % P10_UINT64; \
541	f.val /= P10_UINT64; /* msig */ \
542	base10_128( os, f ); /* recursion */ \
543	_Ostream_Manip(unsigned long long int) fmt @= { lsig, 0, 19, 'u', { .all : 0 } }; \
544	fmt.flags.nobsdp = true; \
545	/* printf( "fmt1 %c %lld %d\n", fmt.base, fmt.val, fmt.all ); */ \
546	sepOff( os ); \
547	(ostype &)(os \| fmt); \
548	} else { \
549	/* printf( "fmt2 %c %lld %d\n", f.base, (unsigned long long int)f.val, f.all ); */ \
550	_Ostream_Manip(SIGNED long long int) fmt @= { (SIGNED long long int)f.val, f.wd, f.pc, f.base, { .all : f.all } }; \
551	(ostype &)(os \| fmt); \
552	} /* if */ \
553	} /* base10_128 */ \
554	forall( dtype ostype \| ostream( ostype ) ) { \
555	ostype & ?\|?( ostype & os, _Ostream_Manip(T) f ) { \
556	if ( $sepPrt( os ) ) fmt( os, "%s", $sepGetCur( os ) ); \
557	\
558	if ( f.base == 'b' \| f.base == 'B' \| f.base == 'o' \| f.base == 'x' \| f.base == 'X' ) { \
559	unsigned long long int msig = (unsigned long long int)(f.val >> 64); \
560	unsigned long long int lsig = (unsigned long long int)(f.val); \
561	_Ostream_Manip(SIGNED long long int) fmt @= { msig, f.wd, f.pc, f.base, { .all : f.all } }; \
562	_Ostream_Manip(unsigned long long int) fmt2 @= { lsig, 0, 0, f.base, { .all : 0 } }; \
563	if ( msig == 0 ) { \
564	fmt.val = lsig; \
565	(ostype &)(os \| fmt); \
566	} else { \
567	fmt2.flags.pad0 = fmt2.flags.nobsdp = true; \
568	if ( f.base == 'b' \| f.base == 'B' ) { \
569	if ( fmt.flags.pc && fmt.pc > 64 ) fmt.pc -= 64; else { fmt.flags.pc = false; fmt.pc = 0; } \
570	if ( fmt.flags.left ) { \
571	fmt.flags.left = false; \
572	fmt.wd = 0; \
573	/* printf( "L %llo %llo %llo %d %d '%c' %x\n", msig, lsig, fmt.val, fmt.wd, fmt.pc, fmt.base, fmt.all ); */ \
574	fmt2.flags.left = true; \
575	int msigd = high1( msig ); \
576	fmt2.wd = f.wd - (fmt.pc > msigd ? fmt.pc : msigd); \
577	if ( ! fmt.flags.nobsdp ) fmt2.wd -= 2; /* compensate for 0b base specifier */ \
578	if ( (int)fmt2.wd < 64 ) fmt2.wd = 64; /* cast deals with negative value */ \
579	fmt2.flags.pc = true; fmt2.pc = 64; \
580	} else { \
581	if ( fmt.wd > 64 ) fmt.wd -= 64; \
582	else fmt.wd = 1; \
583	/* printf( "R %llo %llo %llo %d %d '%c' %x\n", msig, lsig, fmt.val, fmt.wd, fmt.pc, fmt.base, fmt.all ); */ \
584	fmt2.wd = 64; \
585	} /* if */ \
586	/* printf( "C %llo %d %d '%c' %x\n", fmt2.val, fmt2.wd, fmt2.pc, fmt2.base, fmt2.all ); */ \
587	(ostype &)(os \| fmt \| "" \| fmt2); \
588	} else if ( f.base == 'o' ) { \
589	if ( fmt.flags.pc && fmt.pc > 22 ) fmt.pc -= 22; else { fmt.flags.pc = false; fmt.pc = 0; } \
590	fmt.val = (unsigned long long int)fmt.val >> 2; \
591	fmt2.val = ((msig & 0x3) << 1) + ((lsig & 0x8000000000000000U) != 0); \
592	if ( fmt.flags.left ) { \
593	fmt.flags.left = false; \
594	fmt.wd = 0; \
595	/* printf( "L %llo %llo %llo %d %d '%c' %x %llo %d %d '%c' %x\n", msig, lsig, fmt.val, fmt.wd, fmt.pc, fmt.base, fmt.all, fmt2.val, fmt2.wd, fmt2.pc, fmt2.base, fmt2.all ); */ \
596	(ostype &)(os \| fmt \| "" \| fmt2); \
597	sepOff( os ); \
598	fmt2.flags.left = true; \
599	int msigd = ceiling_div( high1( fmt.val ), 3 ); \
600	fmt2.wd = f.wd - (fmt.pc > msigd ? fmt.pc : msigd); \
601	if ( ! fmt.flags.nobsdp ) fmt2.wd -= 1; /* compensate for 0 base specifier */ \
602	if ( (int)fmt2.wd < 21 ) fmt2.wd = 21; /* cast deals with negative value */ \
603	fmt2.flags.pc = true; fmt2.pc = 21; \
604	} else { \
605	if ( fmt.wd > 22 ) fmt.wd -= 22; \
606	else fmt.wd = 1; \
607	/* printf( "R %llo %llo %llo %d %d '%c' %x %llo %d %d '%c' %x\n", msig, lsig, fmt.val, fmt.wd, fmt.pc, fmt.base, fmt.all, fmt2.val, fmt2.wd, fmt2.pc, fmt2.base, fmt2.all ); */ \
608	(ostype &)(os \| fmt \| "" \| fmt2); \
609	sepOff( os ); \
610	fmt2.wd = 21; \
611	} /* if */ \
612	fmt2.val = lsig & 0x7fffffffffffffffU; \
613	/* printf( "\nC %llo %d %d '%c' %x\n", fmt2.val, fmt2.wd, fmt2.pc, fmt2.base, fmt2.all ); */ \
614	(ostype &)(os \| fmt2); \
615	} else { /* f.base == 'x' \| f.base == 'X' */ \
616	if ( fmt.flags.pc && fmt.pc > 16 ) fmt.pc -= 16; else { fmt.flags.pc = false; fmt.pc = 0; } \
617	if ( fmt.flags.left ) { \
618	fmt.flags.left = false; \
619	fmt.wd = 0; \
620	/* printf( "L %llo %llo %llo %d %d '%c' %x\n", msig, lsig, fmt.val, fmt.wd, fmt.pc, fmt.base, fmt.all ); */ \
621	fmt2.flags.left = true; \
622	int msigd = high1( msig ); \
623	fmt2.wd = f.wd - (fmt.pc > msigd ? fmt.pc : msigd); \
624	if ( ! fmt.flags.nobsdp ) fmt2.wd -= 2; /* compensate for 0x base specifier */ \
625	if ( (int)fmt2.wd < 16 ) fmt2.wd = 16; /* cast deals with negative value */ \
626	fmt2.flags.pc = true; fmt2.pc = 16; \
627	} else { \
628	if ( fmt.wd > 16 ) fmt.wd -= 16; \
629	else fmt.wd = 1; \
630	/* printf( "R %llo %llo %llo %d %d '%c' %x\n", msig, lsig, fmt.val, fmt.wd, fmt.pc, fmt.base, fmt.all ); */ \
631	fmt2.wd = 16; \
632	} /* if */ \
633	/* printf( "C %llo %d %d '%c' %x\n", fmt2.val, fmt2.wd, fmt2.pc, fmt2.base, fmt2.all ); */ \
634	(ostype &)(os \| fmt \| "" \| fmt2); \
635	} /* if */ \
636	} /* if */ \
637	} else { \
638	if ( CODE == 'd' ) { \
639	if ( f.val < 0 ) { fmt( os, "-" ); sepOff( os ); f.val = -f.val; f.flags.sign = false; } \
640	} /* if */ \
641	base10_128( os, f ); \
642	} /* if */ \
643	return os; \
644	} /* ?\|? */ \
645	void ?\|?( ostype & os, _Ostream_Manip(T) f ) { (ostype &)(os \| f); ends( os ); } \
646	} // distribution
647
648	IntegralFMTImpl128( int128, signed, 'd', "% ll ", "% .*ll " )
649	IntegralFMTImpl128( unsigned int128, unsigned, 'u', "% ll ", "% .*ll " )
650	#endif // __SIZEOF_INT128__
651	#endif // 0
652
653	#if 1
654	#if defined( __SIZEOF_INT128__ )
655	// Default prefix for non-decimal prints is 0b, 0, 0x.
656	forall( dtype ostype \| ostream( ostype ) )
657	static inline void base_128( ostype & os, unsigned int128 val, unsigned int128 power, _Ostream_Manip(uint64_t) & f, unsigned int maxdig, unsigned int bits, unsigned int cnt = 0 ) {
658	int wd = 1; // f.wd is never 0 because 0 implies left-pad
659	if ( val > power ) { // subdivide value into printable 64-bit values
660	base_128( os, val / power, power, f, maxdig, bits, cnt + 1 ); // recursive
661	f.val = val % power;
662	if ( cnt == 1 && f.flags.left ) { wd = f.wd; f.wd = maxdig; } // copy f.wd and reset for printing middle chunk
663	// printf( "R val:%#lx(%lu) wd:%u pc:%u base:%c neg:%d pc:%d left:%d nobsdp:%d sign:%d pad0:%d\n",
664	// f.val, f.val, f.wd, f.pc, f.base, f.flags.neg, f.flags.pc, f.flags.left, f.flags.nobsdp, f.flags.sign, f.flags.pad0 );
665	(ostype &)(os \| f);
666	if ( cnt == 1 ) {
667	if ( f.flags.left ) { wd -= maxdig; f.wd = wd < 0 ? 1 : wd; } // update and restore f.wd for printing end chunk
668	sepOff( os ); // no seperator between chunks
669	} // if
670	} else { // print start chunk
671	f.val = val;
672	// f.pc is unsigned => use wd
673	if ( f.flags.pc && f.pc > maxdig * cnt ) { wd = f.pc - maxdig * cnt; f.pc = wd < 0 ? 0 : wd; }
674	else { f.flags.pc = false; f.pc = 0; }
675
676	if ( ! f.flags.left ) { // right justify
677	wd = f.wd - maxdig * cnt;
678	f.wd = wd < 0 ? 1 : wd;
679	wd = maxdig;
680	} else { // left justify
681	if ( cnt != 0 ) { // value >= 2^64 ?
682	unsigned int dig, bs = 0;
683	// compute size of prefix digits and base
684	if ( f.base == 'd' \|\| f.base == 'u' ) { // no base prefix
685	dig = ceil( log10( f.val ) ); // use floating-point
686	if ( f.base == 'd' && (f.flags.neg \|\| f.flags.sign) ) bs = 1; // sign ?
687	} else {
688	dig = ceiling_div( high1( f.val ), bits );
689	if ( ! f.flags.nobsdp ) { // base prefix ?
690	if ( f.base == 'o' ) {
691	// 0 prefix for octal is not added for precision with leading zero
692	if ( f.pc <= dig ) bs = 1; // 1 character prefix
693	} else bs = 2; // 2 character prefix
694	} // if
695	} // if
696	wd = f.wd - (f.pc > dig ? f.pc : dig) - bs; // precision > leading digits ?
697	if ( wd < 0 ) wd = 1;
698	f.wd = 1;
699	} // if
700	// all manipulators handled implicitly for value < 2^64
701	} // if
702	// prior checks ensure wd not negative
703
704	if ( f.flags.neg ) f.val = -f.val;
705	// printf( "L val:%#lx(%lu) wd:%u pc:%u base:%c neg:%d pc:%d left:%d nobsdp:%d sign:%d pad0:%d\n",
706	// f.val, f.val, f.wd, f.pc, f.base, f.flags.neg, f.flags.pc, f.flags.left, f.flags.nobsdp, f.flags.sign, f.flags.pad0 );
707	(ostype &)(os \| f);
708
709	// remaining middle and end chunks are padded with 0s on the left
710	if ( ! f.flags.left ) { f.flags.pad0 = true; f.flags.pc = false; } // left pad with 0s
711	else { f.pc = maxdig; f.flags.pc = true; } // left pad with precision
712
713	if ( cnt != 0 ) sepOff( os ); // no seperator between chunks
714	f.wd = wd; // reset f.wd for next chunk
715	f.flags.sign = false; // no leading +/- sign
716	f.flags.nobsdp = true; // no leading base prefix
717	} // if
718	} // base_128
719
720	#define IntegralFMTImpl128( T ) \
721	forall( dtype ostype \| ostream( ostype ) ) { \
722	ostype & ?\|?( ostype & os, _Ostream_Manip(T) f ) { \
723	_Ostream_Manip(uint64_t) fmt; \
724	fmt.[wd, pc, base, all] = f.[wd, pc, base, all]; \
725	if ( f.base == 'b' \| f.base == 'B' ) { \
726	base_128( os, f.val, (unsigned int128)1 << 64, fmt, 64, 1 ); \
727	} else if ( f.base == 'o' ) { \
728	base_128( os, f.val, (unsigned int128)1 << 63, fmt, 21, 3 ); \
729	} else if ( f.base == 'd' \|\| f.base == 'u' ) { \
730	if ( f.base == 'd' && f.val < 0 ) { f.val = -f.val; fmt.flags.neg = true; } \
731	base_128( os, f.val, (unsigned int128)10_000_000_000_000_000_000UL, fmt, 19, 0 ); \
732	} else { \
733	base_128( os, f.val, (unsigned int128)1 << 64, fmt, 16, 4 ); \
734	} /* if */ \
735	return os; \
736	} /* ?\|? */ \
737	void ?\|?( ostype & os, _Ostream_Manip(T) f ) { (ostype &)(os \| f); ends( os ); } \
738	} // distribution
739
740	IntegralFMTImpl128( int128 )
741	IntegralFMTImpl128( unsigned int128 )
742	#endif // __SIZEOF_INT128__
743	#endif // 0
744
745	// ********************************* floating point *********************************
746
747	#define PrintWithDP2( os, format, val, ... ) \
748	{ \
749	enum { size = 48 }; \
750	char buf[size]; \
751	int bufbeg = 0, i, len = snprintf( buf, size, format, ##__VA_ARGS__, val ); \
752	if ( isfinite( val ) && (f.base != 'g' \|\| f.pc != 0) ) { /* if number, print decimal point */ \
753	for ( i = 0; i < len && buf[i] != '.' && buf[i] != 'e' && buf[i] != 'E'; i += 1 ); /* decimal point or scientific ? */ \
754	if ( i == len && ! f.flags.nobsdp ) { \
755	if ( ! f.flags.left ) { \
756	buf[i] = '.'; buf[i + 1] = '\0'; \
757	if ( buf[0] == ' ' ) bufbeg = 1; /* decimal point within width */ \
758	} else { \
759	for ( i = 0; i < len && buf[i] != ' '; i += 1 ); /* trailing blank ? */ \
760	buf[i] = '.'; \
761	if ( i == len ) buf[i + 1] = '\0'; \
762	} /* if */ \
763	} /* if */ \
764	} /* if */ \
765	fmt( os, "%s", &buf[bufbeg] ); \
766	}
767
768	#define FloatingPointFMTImpl( T, DFMTNP, DFMTP ) \
769	forall( dtype ostype \| ostream( ostype ) ) { \
770	ostype & ?\|?( ostype & os, _Ostream_Manip(T) f ) { \
771	if ( $sepPrt( os ) ) fmt( os, "%s", $sepGetCur( os ) ); \
772	char fmtstr[sizeof(DFMTP)]; /* sizeof includes '\0' */ \
773	if ( ! f.flags.pc ) memcpy( &fmtstr, DFMTNP, sizeof(DFMTNP) ); \
774	else memcpy( &fmtstr, DFMTP, sizeof(DFMTP) ); \
775	int star = 4; /* position before first '' / \
776	\
777	/* Insert flags into spaces before '', from right to left. / \
778	if ( f.flags.left ) { fmtstr[star] = '-'; star -= 1; } \
779	if ( f.flags.sign ) { fmtstr[star] = '+'; star -= 1; } \
780	if ( f.flags.pad0 ) { fmtstr[star] = '0'; star -= 1; } \
781	fmtstr[star] = '%'; \
782	\
783	if ( ! f.flags.pc ) { /* no precision */ \
784	fmtstr[sizeof(DFMTNP)-2] = f.base; /* sizeof includes '\0' */ \
785	/* printf( "%g %d %s\n", f.val, f.wd, &fmtstr[star]); */ \
786	PrintWithDP2( os, &fmtstr[star], f.val, f.wd ) \
787	} else { /* precision */ \
788	fmtstr[sizeof(DFMTP)-2] = f.base; /* sizeof includes '\0' */ \
789	/* printf( "%g %d %d %s\n", f.val, f.wd, f.pc, &fmtstr[star] ); */ \
790	PrintWithDP2( os, &fmtstr[star], f.val, f.wd, f.pc ) \
791	} /* if */ \
792	return os; \
793	} /* ?\|? */ \
794	\
795	void ?\|?( ostype & os, _Ostream_Manip(T) f ) { (ostype &)(os \| f); ends( os ); } \
796	} // distribution
797
798	FloatingPointFMTImpl( double, "% * ", "% . " )
799	FloatingPointFMTImpl( long double, "% L ", "% .*L " )
800
801	// ********************************* character *********************************
802
803	forall( dtype ostype \| ostream( ostype ) ) {
804	ostype & ?\|?( ostype & os, _Ostream_Manip(char) f ) {
805	if ( f.base != 'c' ) { // bespoke binary/octal/hex format
806	_Ostream_Manip(unsigned char) fmtuc @= { f.val, f.wd, f.pc, f.base, {'\0'} };
807	fmtuc.flags.pc = f.flags.pc;
808	fmtuc.flags.nobsdp = f.flags.nobsdp;
809	// os \| fmtuc \| nonl;
810	(ostype &)(os \| fmtuc);
811	return os;
812	} // if
813
814	if ( $sepPrt( os ) ) fmt( os, "%s", $sepGetCur( os ) );
815
816	#define CFMTNP "% * "
817	char fmtstr[sizeof(CFMTNP)]; // sizeof includes '\0'
818	memcpy( &fmtstr, CFMTNP, sizeof(CFMTNP) );
819	int star = 1; // position before first '*'
820
821	// Insert flags into spaces before '*', from right to left.
822	if ( f.flags.left ) { fmtstr[star] = '-'; star -= 1; }
823	fmtstr[star] = '%';
824
825	fmtstr[sizeof(CFMTNP)-2] = f.base; // sizeof includes '\0'
826	// printf( "%d %s\n", f.wd, &fmtstr[star] );
827	fmt( os, &fmtstr[star], f.wd, f.val );
828	return os;
829	} // ?\|?
830
831	void ?\|?( ostype & os, _Ostream_Manip(char) f ) { (ostype &)(os \| f); ends( os ); }
832	} // distribution
833
834	// ********************************* C string *********************************
835
836	forall( dtype ostype \| ostream( ostype ) ) {
837	ostype & ?\|?( ostype & os, _Ostream_Manip(const char *) f ) {
838	if ( ! f.val ) return os; // null pointer ?
839
840	if ( f.base != 's' ) { // bespoke binary/octal/hex format
841	_Ostream_Manip(unsigned char) fmtuc @= { 0, f.wd, f.pc, f.base, {'\0'} };
842	fmtuc.flags.pc = f.flags.pc;
843	fmtuc.flags.nobsdp = f.flags.nobsdp;
844	for ( unsigned int i = 0; f.val[i] != '\0'; i += 1 ) {
845	fmtuc.val = f.val[i];
846	// os \| fmtuc \| nonl;
847	(ostype &)(os \| fmtuc);
848	} // for
849	return os;
850	} // if
851
852	if ( $sepPrt( os ) ) fmt( os, "%s", $sepGetCur( os ) );
853
854	#define SFMTNP "% * "
855	#define SFMTP "% . "
856	char fmtstr[sizeof(SFMTP)]; // sizeof includes '\0'
857	if ( ! f.flags.pc ) memcpy( &fmtstr, SFMTNP, sizeof(SFMTNP) );
858	else memcpy( &fmtstr, SFMTP, sizeof(SFMTP) );
859	int star = 1; // position before first '*'
860
861	// Insert flags into spaces before '*', from right to left.
862	if ( f.flags.left ) { fmtstr[star] = '-'; star -= 1; }
863	fmtstr[star] = '%';
864
865	if ( ! f.flags.pc ) { // no precision
866	// printf( "%d %s\n", f.wd, &fmtstr[star] );
867	fmtstr[sizeof(SFMTNP)-2] = f.base; // sizeof includes '\0'
868	fmt( os, &fmtstr[star], f.wd, f.val );
869	} else { // precision
870	fmtstr[sizeof(SFMTP)-2] = f.base; // sizeof includes '\0'
871	// printf( "%d %d %s\n", f.wd, f.pc, &fmtstr[star] );
872	fmt( os, &fmtstr[star], f.wd, f.pc, f.val );
873	} // if
874	return os;
875	} // ?\|?
876
877	void ?\|?( ostype & os, _Ostream_Manip(const char *) f ) { (ostype &)(os \| f); ends( os ); }
878	} // distribution
879
880
881	// ********************************* istream *********************************
882
883
884	forall( dtype istype \| istream( istype ) ) {
885	istype & ?\|?( istype & is, bool & b ) {
886	char val[6];
887	fmt( is, "%5s", val );
888	if ( strcmp( val, "true" ) == 0 ) b = true;
889	else if ( strcmp( val, "false" ) == 0 ) b = false;
890	else {
891	fprintf( stderr, "invalid Boolean constant\n" );
892	abort(); // cannot use abort stream
893	} // if
894	return is;
895	} // ?\|?
896
897	istype & ?\|?( istype & is, char & c ) {
898	char temp;
899	for () {
900	fmt( is, "%c", &temp ); // must pass pointer through varg to fmt
901	// do not overwrite parameter with newline unless appropriate
902	if ( temp != '\n' \|\| getANL( is ) ) { c = temp; break; }
903	if ( eof( is ) ) break;
904	} // for
905	return is;
906	} // ?\|?
907
908	istype & ?\|?( istype & is, signed char & sc ) {
909	fmt( is, "%hhi", &sc );
910	return is;
911	} // ?\|?
912
913	istype & ?\|?( istype & is, unsigned char & usc ) {
914	fmt( is, "%hhi", &usc );
915	return is;
916	} // ?\|?
917
918	istype & ?\|?( istype & is, short int & si ) {
919	fmt( is, "%hi", &si );
920	return is;
921	} // ?\|?
922
923	istype & ?\|?( istype & is, unsigned short int & usi ) {
924	fmt( is, "%hi", &usi );
925	return is;
926	} // ?\|?
927
928	istype & ?\|?( istype & is, int & i ) {
929	fmt( is, "%i", &i );
930	return is;
931	} // ?\|?
932
933	istype & ?\|?( istype & is, unsigned int & ui ) {
934	fmt( is, "%i", &ui );
935	return is;
936	} // ?\|?
937
938	istype & ?\|?( istype & is, long int & li ) {
939	fmt( is, "%li", &li );
940	return is;
941	} // ?\|?
942
943	istype & ?\|?( istype & is, unsigned long int & ulli ) {
944	fmt( is, "%li", &ulli );
945	return is;
946	} // ?\|?
947
948	istype & ?\|?( istype & is, long long int & lli ) {
949	fmt( is, "%lli", &lli );
950	return is;
951	} // ?\|?
952
953	istype & ?\|?( istype & is, unsigned long long int & ulli ) {
954	fmt( is, "%lli", &ulli );
955	return is;
956	} // ?\|?
957
958	#if defined( __SIZEOF_INT128__ )
959	istype & ?\|?( istype & is, int128 & i128 ) {
960	return (istype &)(is \| (unsigned int128 &)i128);
961	} // ?\|?
962
963	istype & ?\|?( istype & is, unsigned int128 & ui128 ) {
964	char s[40];
965	bool sign = false;
966
967	if ( fmt( is, " %[-]", s ) == 1 ) sign = true; // skip whitespace, negative sign ?
968	// If the input is too large, the value returned is undefined. If there is no input, no value is returned
969	if ( fmt( is, "%39[0-9]%*[0-9]", s ) == 1 ) { // take first 39 characters, ignore remaining
970	ui128 = 0;
971	for ( unsigned int i = 0; s[i] != '\0'; i += 1 ) {
972	ui128 = ui128 * 10 + s[i] - '0';
973	} // for
974	if ( sign ) ui128 = -ui128;
975	} else if ( sign ) ungetc( is, '-' ); // return minus when no digits
976	return is;
977	} // ?\|?
978	#endif // __SIZEOF_INT128__
979
980	istype & ?\|?( istype & is, float & f ) {
981	fmt( is, "%f", &f );
982	return is;
983	} // ?\|?
984
985	istype & ?\|?( istype & is, double & d ) {
986	fmt( is, "%lf", &d );
987	return is;
988	} // ?\|?
989
990	istype & ?\|?( istype & is, long double & ld ) {
991	fmt( is, "%Lf", &ld );
992	return is;
993	} // ?\|?
994
995
996	istype & ?\|?( istype & is, float _Complex & fc ) {
997	float re, im;
998	fmt( is, "%f%fi", &re, &im );
999	fc = re + im * _Complex_I;
1000	return is;
1001	} // ?\|?
1002
1003	istype & ?\|?( istype & is, double _Complex & dc ) {
1004	double re, im;
1005	fmt( is, "%lf%lfi", &re, &im );
1006	dc = re + im * _Complex_I;
1007	return is;
1008	} // ?\|?
1009
1010	istype & ?\|?( istype & is, long double _Complex & ldc ) {
1011	long double re, im;
1012	fmt( is, "%Lf%Lfi", &re, &im );
1013	ldc = re + im * _Complex_I;
1014	return is;
1015	} // ?\|?
1016
1017	// istype & ?\|?( istype & is, const char fmt[] ) {
1018	// fmt( is, fmt, "" );
1019	// return is;
1020	// } // ?\|?
1021
1022	istype & ?\|?( istype & is, char * s ) {
1023	fmt( is, "%s", s );
1024	return is;
1025	} // ?\|?
1026
1027	// manipulators
1028	istype & ?\|?( istype & is, istype & (* manip)( istype & ) ) {
1029	return manip( is );
1030	} // ?\|?
1031
1032	istype & nl( istype & is ) {
1033	fmt( is, "%*[^\n]" ); // ignore characters to newline
1034	return is;
1035	} // nl
1036
1037	istype & nlOn( istype & is ) {
1038	nlOn( is ); // call void returning
1039	return is;
1040	} // nlOn
1041
1042	istype & nlOff( istype & is ) {
1043	nlOff( is ); // call void returning
1044	return is;
1045	} // nlOff
1046	} // distribution
1047
1048	// ********************************* manipulators *********************************
1049
1050	forall( dtype istype \| istream( istype ) )
1051	istype & ?\|?( istype & is, _Istream_Cstr f ) {
1052	// skip xxx
1053	if ( ! f.s ) {
1054	// printf( "skip %s %d\n", f.scanset, f.wd );
1055	if ( f.wd == -1 ) fmt( is, f.scanset, "" ); // no input arguments
1056	else for ( f.wd ) fmt( is, "%*c" );
1057	return is;
1058	} // if
1059	size_t len = 0;
1060	if ( f.scanset ) len = strlen( f.scanset );
1061	char fmtstr[len + 16];
1062	int start = 1;
1063	fmtstr[0] = '%';
1064	if ( f.flags.ignore ) { fmtstr[1] = '*'; start += 1; }
1065	if ( f.wd != -1 ) { start += sprintf( &fmtstr[start], "%d", f.wd ); }
1066	// cstr %s, %s, %ws, %ws
1067	if ( ! f.scanset ) {
1068	fmtstr[start] = 's'; fmtstr[start + 1] = '\0';
1069	// printf( "cstr %s\n", fmtstr );
1070	fmt( is, fmtstr, f.s );
1071	return is;
1072	} // if
1073	// incl %[xxx], %[xxx], %w[xxx], %w[xxx]
1074	// excl %[^xxx], %[^xxx], %w[^xxx], %w[^xxx]
1075	fmtstr[start] = '['; start += 1;
1076	if ( f.flags.inex ) { fmtstr[start] = '^'; start += 1; }
1077	strcpy( &fmtstr[start], f.scanset ); // copy includes '\0'
1078	len += start;
1079	fmtstr[len] = ']'; fmtstr[len + 1] = '\0';
1080	// printf( "incl/excl %s\n", fmtstr );
1081	fmt( is, fmtstr, f.s );
1082	return is;
1083	} // ?\|?
1084
1085	forall( dtype istype \| istream( istype ) )
1086	istype & ?\|?( istype & is, _Istream_Char f ) {
1087	fmt( is, "%*c" ); // argument variable unused
1088	return is;
1089	} // ?\|?
1090
1091	#define InputFMTImpl( T, CODE ) \
1092	forall( dtype istype \| istream( istype ) ) \
1093	istype & ?\|?( istype & is, _Istream_Manip(T) f ) { \
1094	enum { size = 16 }; \
1095	char fmtstr[size]; \
1096	if ( f.wd == -1 ) { \
1097	snprintf( fmtstr, size, "%%%s%s", f.ignore ? "*" : "", CODE ); \
1098	} else { \
1099	snprintf( fmtstr, size, "%%%s%d%s", f.ignore ? "*" : "", f.wd, CODE ); \
1100	} /* if */ \
1101	/* printf( "%d %s %p\n", f.wd, fmtstr, &f.val ); */ \
1102	fmt( is, fmtstr, &f.val ); \
1103	return is; \
1104	} // ?\|?
1105
1106	InputFMTImpl( signed char, "hhi" )
1107	InputFMTImpl( unsigned char, "hhi" )
1108	InputFMTImpl( signed short int, "hi" )
1109	InputFMTImpl( unsigned short int, "hi" )
1110	InputFMTImpl( signed int, "i" )
1111	InputFMTImpl( unsigned int, "i" )
1112	InputFMTImpl( signed long int, "li" )
1113	InputFMTImpl( unsigned long int, "li" )
1114	InputFMTImpl( signed long long int, "lli" )
1115	InputFMTImpl( unsigned long long int, "lli" )
1116
1117	InputFMTImpl( float, "f" )
1118	InputFMTImpl( double, "lf" )
1119	InputFMTImpl( long double, "Lf" )
1120
1121	forall( dtype istype \| istream( istype ) )
1122	istype & ?\|?( istype & is, _Istream_Manip(float _Complex) fc ) {
1123	float re, im;
1124	_Istream_Manip(float) fmtuc @= { re, fc.wd, fc.ignore };
1125	is \| fmtuc;
1126	&fmtuc.val = &im;
1127	is \| fmtuc;
1128	if ( ! fc.ignore ) fc.val = re + im * _Complex_I; // re/im are uninitialized for ignore
1129	return is;
1130	} // ?\|?
1131
1132	forall( dtype istype \| istream( istype ) )
1133	istype & ?\|?( istype & is, _Istream_Manip(double _Complex) dc ) {
1134	double re, im;
1135	_Istream_Manip(double) fmtuc @= { re, dc.wd, dc.ignore };
1136	is \| fmtuc;
1137	&fmtuc.val = &im;
1138	is \| fmtuc;
1139	if ( ! dc.ignore ) dc.val = re + im * _Complex_I; // re/im are uninitialized for ignore
1140	return is;
1141	} // ?\|?
1142
1143	forall( dtype istype \| istream( istype ) )
1144	istype & ?\|?( istype & is, _Istream_Manip(long double _Complex) ldc ) {
1145	long double re, im;
1146	_Istream_Manip(long double) fmtuc @= { re, ldc.wd, ldc.ignore };
1147	is \| fmtuc;
1148	&fmtuc.val = &im;
1149	is \| fmtuc;
1150	if ( ! ldc.ignore ) ldc.val = re + im * _Complex_I; // re/im are uninitialized for ignore
1151	return is;
1152	} // ?\|?
1153
1154	// Local Variables: //
1155	// tab-width: 4 //
1156	// compile-command: "cfa iostream.cfa" //
1157	// End: //

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

Original Format