Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed
  • libcfa/src/iostream.cfa

    r7cfef0d re59e663  
    1010// Created On       : Wed May 27 17:56:53 2015
    1111// Last Modified By : Peter A. Buhr
    12 // Last Modified On : Mon Aug 24 08:31:35 2020
    13 // Update Count     : 1130
     12// Last Modified On : Thu Feb 20 15:53:23 2020
     13// Update Count     : 829
    1414//
    1515
    1616#include "iostream.hfa"
    1717
     18extern "C" {
    1819#include <stdio.h>
    1920#include <stdbool.h>                                                                    // true/false
    2021#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
    24 extern "C" {
     22//#include <string.h>                                                                   // strlen, strcmp
    2523extern size_t strlen (const char *__s) __attribute__ ((__nothrow__ , __leaf__)) __attribute__ ((__pure__)) __attribute__ ((__nonnull__ (1)));
    2624extern int strcmp (const char *__s1, const char *__s2) __attribute__ ((__nothrow__ , __leaf__)) __attribute__ ((__pure__)) __attribute__ ((__nonnull__ (1, 2)));
    2725extern char *strcpy (char *__restrict __dest, const char *__restrict __src) __attribute__ ((__nothrow__ , __leaf__)) __attribute__ ((__nonnull__ (1, 2)));
    2826extern void *memcpy (void *__restrict __dest, const void *__restrict __src, size_t __n) __attribute__ ((__nothrow__ , __leaf__)) __attribute__ ((__nonnull__ (1, 2)));
     27#include <float.h>                                                                              // DBL_DIG, LDBL_DIG
     28#include <math.h>                                                                               // isfinite
     29#include <complex.h>                                                                    // creal, cimag
    2930} // extern "C"
    3031
    31 #include "math.hfa"                                                                             // isfinite, floor, ceiling_div
    32 #include "bitmanip.hfa"                                                                 // high1
    33 
    34 
    35 // *********************************** ostream ***********************************
     32
     33//*********************************** ostream ***********************************
    3634
    3735
    3836forall( dtype ostype | ostream( ostype ) ) {
     37        ostype & ?|?( ostype & os, zero_t ) {
     38                if ( $sepPrt( os ) ) fmt( os, "%s", $sepGetCur( os ) );
     39                fmt( os, "%d", 0n );
     40                return os;
     41        } // ?|?
     42        void ?|?( ostype & os, zero_t z ) {
     43                (ostype &)(os | z); ends( os );
     44        } // ?|?
     45
     46        ostype & ?|?( ostype & os, one_t ) {
     47                if ( $sepPrt( os ) ) fmt( os, "%s", $sepGetCur( os ) );
     48                fmt( os, "%d", 1n );
     49                return os;
     50        } // ?|?
     51        void ?|?( ostype & os, one_t o ) {
     52                (ostype &)(os | o); ends( os );
     53        } // ?|?
     54
    3955        ostype & ?|?( ostype & os, bool b ) {
    4056                if ( $sepPrt( os ) ) fmt( os, "%s", $sepGetCur( os ) );
     
    149165        #define P10_UINT64 10_000_000_000_000_000_000_ULL       // 19 zeroes
    150166
    151         static inline void base10_128( ostype & os, unsigned int128 val ) {
    152 #if defined(__GNUC__) && __GNUC_PREREQ(7,0)                             // gcc version >= 7
    153                 if ( val > P10_UINT64 ) {
    154 #else
    155                 if ( (uint64_t)(val >> 64) != 0 || (uint64_t)val > P10_UINT64 ) { // patch gcc 5 & 6 -O3 bug
    156 #endif // __GNUC_PREREQ(7,0)
     167        static void base10_128( ostype & os, unsigned int128 val ) {
     168                if ( val > UINT64_MAX ) {
    157169                        base10_128( os, val / P10_UINT64 );                     // recursive
    158170                        fmt( os, "%.19lu", (uint64_t)(val % P10_UINT64) );
     
    162174        } // base10_128
    163175
    164         static inline void base10_128( ostype & os, int128 val ) {
     176        static void base10_128( ostype & os, int128 val ) {
    165177                if ( val < 0 ) {
    166178                        fmt( os, "-" );                                                         // leading negative sign
     
    433445} // distribution
    434446
    435 // *********************************** manipulators ***********************************
    436 
    437 // *********************************** integral ***********************************
     447//*********************************** manipulators ***********************************
     448
     449//*********************************** integral ***********************************
    438450
    439451static const char * shortbin[] = { "0", "1", "10", "11", "100", "101", "110", "111", "1000", "1001", "1010", "1011", "1100", "1101", "1110", "1111" };
     
    441453
    442454// Default prefix for non-decimal prints is 0b, 0, 0x.
    443 #define IntegralFMTImpl( T, IFMTNP, IFMTP ) \
     455#define IntegralFMTImpl( T, CODE, IFMTNP, IFMTP ) \
    444456forall( dtype ostype | ostream( ostype ) ) { \
    445457        ostype & ?|?( ostype & os, _Ostream_Manip(T) f ) { \
     
    447459\
    448460                if ( f.base == 'b' || f.base == 'B' ) {                 /* bespoke binary format */ \
    449                         int bits = high1( f.val );                                      /* position of most significant bit */ \
    450                         if ( bits == 0 ) bits = 1;                                      /* 0 value => force one bit to print */ \
    451                         int spaces; \
     461                        int bits;                                                                                                       \
     462                        if ( f.val == (T){0} ) bits = 1;                        /* force at least one bit to print */ \
     463                        else bits = sizeof(long long int) * 8 - __builtin_clzll( f.val ); /* position of most significant bit */ \
     464                        bits = bits > sizeof(f.val) * 8 ? sizeof(f.val) * 8 : bits; \
     465                        int spaces = f.wd - bits;                                       /* can be negative */ \
     466                        if ( ! f.flags.nobsdp ) { spaces -= 2; }        /* base prefix takes space */ \
     467                        /* printf( "%d %d\n", bits, spaces ); */ \
    452468                        if ( ! f.flags.left ) {                                         /* right justified ? */ \
    453469                                /* Note, base prefix then zero padding or spacing then prefix. */ \
    454                                 if ( f.flags.pc ) { \
    455                                         spaces = f.wd - f.pc; \
    456                                         if ( ! f.flags.nobsdp ) { spaces -= 2; } /* base prefix takes space */ \
     470                                if ( f.flags.pad0 || f.flags.pc ) { \
     471                                        if ( ! f.flags.nobsdp ) { fmt( os, "0%c", f.base ); } \
     472                                        if ( f.flags.pc ) spaces = f.pc - bits; \
     473                                        if ( spaces > 0 ) fmt( os, "%0*d", spaces, 0 ); /* zero pad */ \
     474                                } else { \
    457475                                        if ( spaces > 0 ) fmt( os, "%*s", spaces, " " ); /* space pad */ \
    458476                                        if ( ! f.flags.nobsdp ) { fmt( os, "0%c", f.base ); } \
    459                                         spaces = f.pc - bits; \
    460                                         if ( spaces > 0 ) fmt( os, "%0*d", spaces, 0 ); /* zero pad */ \
    461                                 } else { \
    462                                         spaces = f.wd - bits; \
    463                                         if ( ! f.flags.nobsdp ) { spaces -= 2; } /* base prefix takes space */ \
    464                                         if ( f.flags.pad0 ) { \
    465                                                 if ( ! f.flags.nobsdp ) { fmt( os, "0%c", f.base ); } \
    466                                                 if ( spaces > 0 ) fmt( os, "%0*d", spaces, 0 ); /* zero pad */ \
    467                                         } else { \
    468                                                 if ( spaces > 0 ) fmt( os, "%*s", spaces, " " ); /* space pad */ \
    469                                                 if ( ! f.flags.nobsdp ) { fmt( os, "0%c", f.base ); } \
    470                                         } /* if */ \
    471477                                } /* if */ \
    472                         } else { \
    473                                 if ( ! f.flags.nobsdp ) fmt( os, "0%c", f.base ); \
    474                                 if ( f.flags.pc ) { \
    475                                         spaces = f.pc - bits; \
    476                                         if ( spaces > 0 ) fmt( os, "%0*d", spaces, 0 ); /* zero pad */ \
    477                                         spaces = f.wd - f.pc; \
    478                                 } else { /* pad0 flag ignored with left flag */ \
    479                                         spaces = f.wd - bits; \
    480                                 } /* if */ \
    481                                 if ( ! f.flags.nobsdp ) { spaces -= 2; } /* base prefix takes space */ \
     478                        } else if ( ! f.flags.nobsdp ) { \
     479                                fmt( os, "0%c", f.base ); \
    482480                        } /* if */ \
    483                         int shift = floor( bits - 1, 4 ); \
     481                        int shift = (bits - 1) / 4 * 4; /* floor( bits - 1, 4 ) */ \
    484482                        typeof( f.val ) temp = f.val; \
    485483                        fmt( os, "%s", shortbin[(temp >> shift) & 0xf] ); \
     
    492490                        if ( f.flags.left && spaces > 0 ) fmt( os, "%*s", spaces, " " ); \
    493491                        return os; \
    494                 } /* if */ \
     492                } /* if  */ \
    495493\
    496494                char fmtstr[sizeof(IFMTP)];                                             /* sizeof includes '\0' */ \
     
    502500                if ( ! f.flags.nobsdp ) { fmtstr[star] = '#'; star -= 1; } \
    503501                if ( f.flags.left ) { fmtstr[star] = '-'; star -= 1; } \
    504                 if ( f.flags.sign ) { fmtstr[star] = '+'; star -= 1; } \
     502                if ( f.flags.sign && f.base == CODE ) { fmtstr[star] = '+'; star -= 1; } \
    505503                if ( f.flags.pad0 && ! f.flags.pc ) { fmtstr[star] = '0'; star -= 1; } \
    506504                fmtstr[star] = '%'; \
     
    508506                if ( ! f.flags.pc ) {                                                   /* no precision */ \
    509507                        fmtstr[sizeof(IFMTNP)-2] = f.base;                      /* sizeof includes '\0' */ \
    510                         /* printf( "%s %c\n", &fmtstr[star], f.base ); */ \
     508                        /* printf( "%s %c %c\n", &fmtstr[star], f.base, CODE ); */ \
    511509                        fmt( os, &fmtstr[star], f.wd, f.val ); \
    512510                } else {                                                                                /* precision */ \
    513511                        fmtstr[sizeof(IFMTP)-2] = f.base;                       /* sizeof includes '\0' */ \
    514                         /* printf( "%s %c\n", &fmtstr[star], f.base ); */ \
     512                        /* printf( "%s %c %c\n", &fmtstr[star], f.base, CODE ); */ \
    515513                        fmt( os, &fmtstr[star], f.wd, f.pc, f.val ); \
    516514                } /* if */ \
     
    520518} // distribution
    521519
    522 IntegralFMTImpl( signed char, "%    *hh ", "%    *.*hh " )
    523 IntegralFMTImpl( unsigned char, "%    *hh ", "%    *.*hh " )
    524 IntegralFMTImpl( signed short int, "%    *h ", "%    *.*h " )
    525 IntegralFMTImpl( unsigned short int, "%    *h ", "%    *.*h " )
    526 IntegralFMTImpl( signed int, "%    * ", "%    *.* " )
    527 IntegralFMTImpl( unsigned int, "%    * ", "%    *.* " )
    528 IntegralFMTImpl( signed long int, "%    *l ", "%    *.*l " )
    529 IntegralFMTImpl( unsigned long int, "%    *l ", "%    *.*l " )
    530 IntegralFMTImpl( signed long long int, "%    *ll ", "%    *.*ll " )
    531 IntegralFMTImpl( unsigned long long int, "%    *ll ", "%    *.*ll " )
    532 
    533 #if 0
     520IntegralFMTImpl( signed char, 'd', "%    *hh ", "%    *.*hh " )
     521IntegralFMTImpl( unsigned char, 'u', "%    *hh ", "%    *.*hh " )
     522IntegralFMTImpl( signed short int, 'd', "%    *h ", "%    *.*h " )
     523IntegralFMTImpl( unsigned short int, 'u', "%    *h ", "%    *.*h " )
     524IntegralFMTImpl( signed int, 'd', "%    * ", "%    *.* " )
     525IntegralFMTImpl( unsigned int, 'u', "%    * ", "%    *.* " )
     526IntegralFMTImpl( signed long int, 'd', "%    *l ", "%    *.*l " )
     527IntegralFMTImpl( unsigned long int, 'u', "%    *l ", "%    *.*l " )
     528IntegralFMTImpl( signed long long int, 'd', "%    *ll ", "%    *.*ll " )
     529IntegralFMTImpl( unsigned long long int, 'u', "%    *ll ", "%    *.*ll " )
     530
     531
    534532#if defined( __SIZEOF_INT128__ )
    535533// Default prefix for non-decimal prints is 0b, 0, 0x.
    536534#define IntegralFMTImpl128( T, SIGNED, CODE, IFMTNP, IFMTP ) \
    537535forall( dtype ostype | ostream( ostype ) ) \
    538 static void base10_128( ostype & os, _Ostream_Manip(T) f ) { \
    539         if ( f.val > UINT64_MAX ) { \
    540                 unsigned long long int lsig = f.val % P10_UINT64; \
    541                 f.val /= P10_UINT64; /* msig */ \
    542                 base10_128( os, f ); /* recursion */ \
    543                 _Ostream_Manip(unsigned long long int) fmt @= { lsig, 0, 19, 'u', { .all : 0 } }; \
    544                 fmt.flags.nobsdp = true; \
    545                 /* printf( "fmt1 %c %lld %d\n", fmt.base, fmt.val, fmt.all ); */ \
     536static void base10_128( ostype & os, _Ostream_Manip(T) fmt ) { \
     537        if ( fmt.val > UINT64_MAX ) { \
     538                fmt.val /= P10_UINT64; \
     539                base10_128( os, fmt ); /* recursive */ \
     540                _Ostream_Manip(unsigned long long int) fmt2 @= { (uint64_t)(fmt.val % P10_UINT64), 0, 19, 'u', { .all : 0 } }; \
     541                fmt2.flags.nobsdp = true; \
     542                printf( "fmt2 %c %lld %d\n", fmt2.base, fmt2.val, fmt2.all );   \
    546543                sepOff( os ); \
    547                 (ostype &)(os | fmt); \
     544                (ostype &)(os | fmt2); \
    548545        } else { \
    549                 /* printf( "fmt2 %c %lld %d\n", f.base, (unsigned long long int)f.val, f.all ); */ \
    550                 _Ostream_Manip(SIGNED long long int) fmt @= { (SIGNED long long int)f.val, f.wd, f.pc, f.base, { .all : f.all } }; \
     546                printf( "fmt %c %lld %d\n", fmt.base, fmt.val, fmt.all ); \
    551547                (ostype &)(os | fmt); \
    552548        } /* if */ \
    553 } /* base10_128 */ \
     549} /* base10_128 */                                                \
    554550forall( dtype ostype | ostream( ostype ) ) { \
    555551        ostype & ?|?( ostype & os, _Ostream_Manip(T) f ) { \
    556552                if ( $sepPrt( os ) ) fmt( os, "%s", $sepGetCur( os ) ); \
    557553\
    558                 if ( f.base == 'b' | f.base == 'B' | f.base == 'o' | f.base == 'x' | f.base == 'X' ) { \
     554                if ( f.base == 'b' | f.base == 'o' | f.base == 'x' | f.base == 'X' ) { \
    559555                        unsigned long long int msig = (unsigned long long int)(f.val >> 64); \
    560556                        unsigned long long int lsig = (unsigned long long int)(f.val); \
     
    566562                        } else { \
    567563                                fmt2.flags.pad0 = fmt2.flags.nobsdp = true;     \
    568                                 if ( f.base == 'b' | f.base == 'B' ) { \
    569                                         if ( fmt.flags.pc && fmt.pc > 64 ) fmt.pc -= 64; else { fmt.flags.pc = false; fmt.pc = 0; } \
    570                                         if ( fmt.flags.left ) { \
    571                                                 fmt.flags.left = false; \
    572                                                 fmt.wd = 0; \
    573                                                 /* printf( "L %llo %llo %llo %d %d '%c' %x\n", msig, lsig, fmt.val, fmt.wd, fmt.pc, fmt.base, fmt.all ); */ \
    574                                                 fmt2.flags.left = true; \
    575                                                 int msigd = high1( msig ); \
    576                                                 fmt2.wd = f.wd - (fmt.pc > msigd ? fmt.pc : msigd); \
    577                                                 if ( ! fmt.flags.nobsdp ) fmt2.wd -= 2; /* compensate for 0b base specifier */ \
    578                                                 if ( (int)fmt2.wd < 64 ) fmt2.wd = 64; /* cast deals with negative value */ \
    579                                                 fmt2.flags.pc = true; fmt2.pc = 64; \
    580                                         } else { \
    581                                                 if ( fmt.wd > 64 ) fmt.wd -= 64; \
    582                                                 else fmt.wd = 1; \
    583                                                 /* printf( "R %llo %llo %llo %d %d '%c' %x\n", msig, lsig, fmt.val, fmt.wd, fmt.pc, fmt.base, fmt.all ); */ \
    584                                                 fmt2.wd = 64; \
    585                                         } /* if */ \
    586                                         /* printf( "C %llo %d %d '%c' %x\n", fmt2.val, fmt2.wd, fmt2.pc, fmt2.base, fmt2.all ); */ \
     564                                if ( f.base == 'b' ) { \
     565                                        if ( f.wd > 64 ) fmt.wd = f.wd - 64; \
     566                                        fmt2.wd = 64; \
    587567                                        (ostype &)(os | fmt | "" | fmt2); \
    588568                                } else if ( f.base == 'o' ) { \
    589                                         if ( fmt.flags.pc && fmt.pc > 22 ) fmt.pc -= 22; else { fmt.flags.pc = false; fmt.pc = 0; } \
    590569                                        fmt.val = (unsigned long long int)fmt.val >> 2; \
    591                                         fmt2.val = ((msig & 0x3) << 1) + ((lsig & 0x8000000000000000U) != 0); \
    592                                         if ( fmt.flags.left ) { \
    593                                                 fmt.flags.left = false; \
    594                                                 fmt.wd = 0; \
    595                                                 /* 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 ); */ \
    596                                                 (ostype &)(os | fmt | "" | fmt2); \
    597                                                 sepOff( os ); \
    598                                                 fmt2.flags.left = true; \
    599                                                 int msigd = ceiling_div( high1( fmt.val ), 3 ); \
    600                                                 fmt2.wd = f.wd - (fmt.pc > msigd ? fmt.pc : msigd); \
    601                                                 if ( ! fmt.flags.nobsdp ) fmt2.wd -= 1; /* compensate for 0 base specifier */ \
    602                                                 if ( (int)fmt2.wd < 21 ) fmt2.wd = 21; /* cast deals with negative value */ \
    603                                                 fmt2.flags.pc = true; fmt2.pc = 21; \
     570                                        if ( f.wd > 21 ) fmt.wd = f.wd - 21; \
     571                                        fmt2.wd = 1; \
     572                                        fmt2.val = ((msig & 0x3) << 1) + 1; \
     573                                        (ostype &)(os | fmt | "" | fmt2); \
     574                                        sepOff( os ); \
     575                                        fmt2.wd = 21; \
     576                                        fmt2.val = lsig & 0x7fffffffffffffff; \
     577                                        (ostype &)(os | fmt2); \
     578                                } else { \
     579                                        if ( f.flags.left ) { \
     580                                                if ( f.wd > 16 ) fmt2.wd = f.wd - 16;   \
     581                                                fmt.wd = 16;                                                    \
    604582                                        } else { \
    605                                                 if ( fmt.wd > 22 ) fmt.wd -= 22; \
    606                                                 else fmt.wd = 1; \
    607                                                 /* 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 ); */ \
    608                                                 (ostype &)(os | fmt | "" | fmt2); \
    609                                                 sepOff( os ); \
    610                                                 fmt2.wd = 21; \
     583                                                if ( f.wd > 16 ) fmt.wd = f.wd - 16;    \
     584                                                fmt2.wd = 16;                                                   \
    611585                                        } /* if */ \
    612                                         fmt2.val = lsig & 0x7fffffffffffffffU; \
    613                                         /* printf( "\nC %llo %d %d '%c' %x\n", fmt2.val, fmt2.wd, fmt2.pc, fmt2.base, fmt2.all ); */ \
    614                                         (ostype &)(os | fmt2); \
    615                                 } else { /* f.base == 'x'  | f.base == 'X' */ \
    616                                         if ( fmt.flags.pc && fmt.pc > 16 ) fmt.pc -= 16; else { fmt.flags.pc = false; fmt.pc = 0; } \
    617                                         if ( fmt.flags.left ) { \
    618                                                 fmt.flags.left = false; \
    619                                                 fmt.wd = 0; \
    620                                                 /* printf( "L %llo %llo %llo %d %d '%c' %x\n", msig, lsig, fmt.val, fmt.wd, fmt.pc, fmt.base, fmt.all ); */ \
    621                                                 fmt2.flags.left = true; \
    622                                                 int msigd = high1( msig ); \
    623                                                 fmt2.wd = f.wd - (fmt.pc > msigd ? fmt.pc : msigd); \
    624                                                 if ( ! fmt.flags.nobsdp ) fmt2.wd -= 2; /* compensate for 0x base specifier */ \
    625                                                 if ( (int)fmt2.wd < 16 ) fmt2.wd = 16; /* cast deals with negative value */ \
    626                                                 fmt2.flags.pc = true; fmt2.pc = 16; \
    627                                         } else { \
    628                                                 if ( fmt.wd > 16 ) fmt.wd -= 16; \
    629                                                 else fmt.wd = 1; \
    630                                                 /* printf( "R %llo %llo %llo %d %d '%c' %x\n", msig, lsig, fmt.val, fmt.wd, fmt.pc, fmt.base, fmt.all ); */ \
    631                                                 fmt2.wd = 16; \
    632                                         } /* if */ \
    633                                         /* printf( "C %llo %d %d '%c' %x\n", fmt2.val, fmt2.wd, fmt2.pc, fmt2.base, fmt2.all ); */ \
    634586                                        (ostype &)(os | fmt | "" | fmt2); \
    635587                                } /* if */ \
    636588                        } /* if */ \
    637589                } else { \
    638                         if ( CODE == 'd' ) { \
    639                                 if ( f.val < 0 )  { fmt( os, "-" ); sepOff( os ); f.val = -f.val; f.flags.sign = false; } \
    640                         } /* if */ \
    641590                        base10_128( os, f ); \
    642591                } /* if */ \
     
    649598IntegralFMTImpl128( unsigned int128, unsigned, 'u', "%    *ll ", "%    *.*ll " )
    650599#endif // __SIZEOF_INT128__
    651 #endif // 0
    652 
    653 #if 1
    654 #if defined( __SIZEOF_INT128__ )
    655 // Default prefix for non-decimal prints is 0b, 0, 0x.
    656 forall( dtype ostype | ostream( ostype ) )
    657 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 ) {
    658         int wd = 1;                                                                                     // f.wd is never 0 because 0 implies left-pad
    659         if ( val > power ) {                                                            // subdivide value into printable 64-bit values
    660                 base_128( os, val / power, power, f, maxdig, bits, cnt + 1 ); // recursive
    661                 f.val = val % power;
    662                 if ( cnt == 1 && f.flags.left ) { wd = f.wd; f.wd = maxdig; } // copy f.wd and reset for printing middle chunk
    663                 // printf( "R val:%#lx(%lu) wd:%u pc:%u base:%c neg:%d pc:%d left:%d nobsdp:%d sign:%d pad0:%d\n",
    664                 //              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 );
    665                 (ostype &)(os | f);
    666                 if ( cnt == 1 ) {
    667                         if ( f.flags.left ) { wd -= maxdig; f.wd = wd < 0 ? 1 : wd; } // update and restore f.wd for printing end chunk
    668                         sepOff( os );                                                           // no seperator between chunks
    669                 } // if
    670         } else {                                                                                        // print start chunk
    671                 f.val = val;
    672                 // f.pc is unsigned => use wd
    673                 if ( f.flags.pc && f.pc > maxdig * cnt ) { wd = f.pc - maxdig * cnt; f.pc = wd < 0 ? 0 : wd; }
    674                 else { f.flags.pc = false; f.pc = 0; }
    675 
    676                 if ( ! f.flags.left ) {                                                 // right justify
    677                         wd = f.wd - maxdig * cnt;
    678                         f.wd = wd < 0 ? 1 : wd;
    679                         wd = maxdig;
    680                 } else {                                                                                // left justify
    681                         if ( cnt != 0 ) {                                                       // value >= 2^64 ?
    682                                 unsigned int dig, bs = 0;
    683                                 // compute size of prefix digits and base
    684                                 if ( f.base == 'd' || f.base == 'u' ) { // no base prefix
    685                                         dig = ceil( log10( f.val ) );           // use floating-point
    686                                         if ( f.base == 'd' && (f.flags.neg || f.flags.sign) ) bs = 1; // sign ?
    687                                 } else {
    688                                         dig = ceiling_div( high1( f.val ), bits );
    689                                         if ( ! f.flags.nobsdp ) {                       // base prefix ?
    690                                                 if ( f.base == 'o' ) {
    691                                                         // 0 prefix for octal is not added for precision with leading zero
    692                                                         if ( f.pc <= dig ) bs = 1;      // 1 character prefix
    693                                                 } else bs = 2;                                  // 2 character prefix
    694                                         } // if
    695                                 } // if
    696                                 wd = f.wd - (f.pc > dig ? f.pc : dig) - bs; // precision > leading digits ?
    697                                 if ( wd < 0 ) wd = 1;
    698                                 f.wd = 1;
    699                         } // if
    700                         // all manipulators handled implicitly for value < 2^64
    701                 } // if
    702                 // prior checks ensure wd not negative
    703 
    704                 if ( f.flags.neg ) f.val = -f.val;
    705                 // printf( "L val:%#lx(%lu) wd:%u pc:%u base:%c neg:%d pc:%d left:%d nobsdp:%d sign:%d pad0:%d\n",
    706                 //              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 );
    707                 (ostype &)(os | f);
    708 
    709                 // remaining middle and end chunks are padded with 0s on the left
    710                 if ( ! f.flags.left ) { f.flags.pad0 = true; f.flags.pc = false; } // left pad with 0s
    711                 else { f.pc = maxdig; f.flags.pc = true; }              // left pad with precision
    712 
    713                 if ( cnt != 0 ) sepOff( os );                                   // no seperator between chunks
    714                 f.wd = wd;                                                                              // reset f.wd for next chunk
    715                 f.flags.sign = false;                                                   // no leading +/- sign
    716                 f.flags.nobsdp = true;                                                  // no leading base prefix
    717         } // if
    718 } // base_128
    719 
    720 #define IntegralFMTImpl128( T ) \
    721 forall( dtype ostype | ostream( ostype ) ) { \
    722         ostype & ?|?( ostype & os, _Ostream_Manip(T) f ) { \
    723                 _Ostream_Manip(uint64_t) fmt; \
    724                 fmt.[wd, pc, base, all] = f.[wd, pc, base, all]; \
    725                 if ( f.base == 'b' | f.base == 'B' ) { \
    726                         base_128( os, f.val, (unsigned int128)1 << 64, fmt, 64, 1 ); \
    727                 } else if ( f.base == 'o' ) { \
    728                         base_128( os, f.val, (unsigned int128)1 << 63, fmt, 21, 3 ); \
    729                 } else if ( f.base == 'd' || f.base == 'u' ) { \
    730                         if ( f.base == 'd' && f.val < 0 ) { f.val = -f.val; fmt.flags.neg = true; } \
    731                         base_128( os, f.val, (unsigned int128)10_000_000_000_000_000_000UL, fmt, 19, 0 ); \
    732                 } else { \
    733                         base_128( os, f.val, (unsigned int128)1 << 64, fmt, 16, 4 ); \
    734                 } /* if */ \
    735                 return os; \
    736         } /* ?|? */ \
    737         void ?|?( ostype & os, _Ostream_Manip(T) f ) { (ostype &)(os | f); ends( os ); } \
    738 } // distribution
    739 
    740 IntegralFMTImpl128( int128 )
    741 IntegralFMTImpl128( unsigned int128 )
    742 #endif // __SIZEOF_INT128__
    743 #endif // 0
    744 
    745 // *********************************** floating point ***********************************
     600
     601//*********************************** floating point ***********************************
    746602
    747603#define PrintWithDP2( os, format, val, ... ) \
     
    799655FloatingPointFMTImpl( long double, "%    *L ", "%    *.*L " )
    800656
    801 // *********************************** character ***********************************
     657//*********************************** character ***********************************
    802658
    803659forall( dtype ostype | ostream( ostype ) ) {
     
    832688} // distribution
    833689
    834 // *********************************** C string ***********************************
     690//*********************************** C string ***********************************
    835691
    836692forall( dtype ostype | ostream( ostype ) ) {
     
    879735
    880736
    881 // *********************************** istream ***********************************
     737//*********************************** istream ***********************************
    882738
    883739
     
    956812        } // ?|?
    957813
    958 #if defined( __SIZEOF_INT128__ )
    959         istype & ?|?( istype & is, int128 & i128 ) {
    960                 return (istype &)(is | (unsigned int128 &)i128);
    961         } // ?|?
    962 
    963         istype & ?|?( istype & is, unsigned int128 & ui128 ) {
    964                 char s[40];
    965                 bool sign = false;
    966 
    967                 if ( fmt( is, " %[-]", s ) == 1 ) sign = true;  // skip whitespace, negative sign ?
    968                 // If the input is too large, the value returned is undefined. If there is no input, no value is returned
    969                 if ( fmt( is, "%39[0-9]%*[0-9]", s ) == 1 ) {   // take first 39 characters, ignore remaining
    970                         ui128 = 0;
    971                         for ( unsigned int i = 0; s[i] != '\0'; i += 1 ) {
    972                                 ui128 = ui128 * 10 + s[i] - '0';
    973                         } // for
    974                         if ( sign ) ui128 = -ui128;
    975                 } else if ( sign ) ungetc( is, '-' );                   // return minus when no digits
    976                 return is;
    977         } // ?|?
    978 #endif // __SIZEOF_INT128__
    979814
    980815        istype & ?|?( istype & is, float & f ) {
     
    1046881} // distribution
    1047882
    1048 // *********************************** manipulators ***********************************
     883//*********************************** manipulators ***********************************
    1049884
    1050885forall( dtype istype | istream( istype ) )
Note: See TracChangeset for help on using the changeset viewer.