Context Navigation

iostream.cfa @ 86530e7

ADTast-experimentalenumforall-pointer-decaypthread-emulationqualifiedEnum

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

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

Download in other formats:

Original Format