source: libcfa/src/iostream.cfa @ 0e398ad

ADTast-experimental
Last change on this file since 0e398ad was fb907d3, checked in by Peter A. Buhr <pabuhr@…>, 21 months ago

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