Context Navigation

iostream.cfa@ 53e4562

Visit:

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

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

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

Download in other formats:

Original Format