Context Navigation

iostream.cfa@ feb999f

Visit:

Last change on this file since feb999f was 6b228cae, checked in by Peter A. Buhr <pabuhr@…>, 19 months ago
more testing for missing data when reading
Property mode set to `100644`
File size: 42.9 KB

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

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

Download in other formats:

Original Format