Большие числа

long double x = 1.3624742724742527487782123333;
  // У меня в win32 в консоли 1.3624742724742526789, т.е. 15 знаков совпадают, остальные уже нет
  std::cout.precision( 20 );
  std::cout << x << '\n';

struct u128
{
public:
  u128( int v ) : a(v), b(0UL), c(0UL), d(0UL)  {};
  u128( u32 v ) : a(v), b(0UL), c(0UL), d(0UL)  {};
  //u128( const u128 &x ) : a(x.a), b(x.b), c(x.c), d(x.d) {};

public:
  u32 a, b, c, d;
};
#pragma pack( pop )

void inc( u128 *x )
{
  __asm__
  (
    "addl $1, (%%eax)           ;"
    "adcl $0, 4(%%eax)          ;"
    "adcl $0, 8(%%eax)          ;"
    "adcl $0, 12(%%eax)         ;"
    :
    : "a" (x)
  );
};

void dec( u128 *x )
{
  __asm__
  (
    "subl $1, (%%eax)           ;"
    "sbbl $0, 4(%%eax)          ;"
    "sbbl $0, 8(%%eax)          ;"
    "sbbl $0, 12(%%eax)         ;"
    :
    : "a" (x)
  );
};

void add( u128 *x, const u128 *y )
{
  __asm__
  (
    "movl (%%ebx), %%edx        ;"
    "addl %%edx, (%%eax)        ;"
    "movl 4(%%ebx), %%edx       ;"
    "adcl %%edx, 4(%%eax)       ;"
    "movl 8(%%ebx), %%edx       ;"
    "adcl %%edx, 8(%%eax)       ;"
    "movl 12(%%ebx), %%edx      ;"
    "adcl %%edx, 12(%%eax)      ;"
    :
    : "a" (x), "b" (y)
  );
};

void sub( u128 *x, const u128 *y )
{
  __asm__
  (
    "movl (%%ebx), %%edx        ;"
    "subl %%edx, (%%eax)        ;"
    "movl 4(%%ebx), %%edx       ;"
    "sbbl %%edx, 4(%%eax)       ;"
    "movl 8(%%ebx), %%edx       ;"
    "sbbl %%edx, 8(%%eax)       ;"
    "movl 12(%%ebx), %%edx      ;"
    "sbbl %%edx, 12(%%eax)      ;"
    :
    : "a" (x), "b" (y)
  );
};