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