return (uint32)high << 16;
}
-#ifdef HAVE_FAST_LONGLONG
+#ifdef HAVE_LONGLONG
// High-Word einer 64-Bit-Zahl bestimmen
// high32(wert)
return (uint64)high << 32;
}
-#endif /* HAVE_FAST_LONGLONG */
+#endif /* HAVE_LONGLONG */
// Multipliziert zwei 16-Bit-Zahlen miteinander und liefert eine 32-Bit-Zahl:
);
return highlow32(_hi,_lo);
}
-#elif defined(__sparc__) || defined(__sparc64__)
+#elif (defined(__sparc__) || defined(__sparc64__)) && !defined(NO_ASM)
extern "C" uint32 mulu16_ (uint16 arg1, uint16 arg2);
#define mulu16 mulu16_ // extern in Assembler
#else
var union { double f; uint32 i[2]; uint16 s[4]; } __fi; \
__fi.f = (double)(sint32)(_x)*(double)(sint32)(_y) \
+ (double)(4503599627370496.0L); /* + 2^52, zum Normalisieren */\
- hi_zuweisung __fi.s[1]; /* mittlere 16 Bit herausholen, (benutzt CL_CPU_BIG_ENDIAN_P !) */\
+ unused (hi_zuweisung __fi.s[1]); /* mittlere 16 Bit herausholen, (benutzt CL_CPU_BIG_ENDIAN_P !) */\
lo_zuweisung __fi.i[1]; /* untere 32 Bit herausholen (benutzt CL_CPU_BIG_ENDIAN_P !) */\
}
#else
// mulu32_unchecked(arg1,arg2)
// > arg1, arg2 : zwei 32-Bit-Zahlen
// < ergebnis : eine 32-Bit-Zahl
-// Es wird vorausgesetzt, daß arg1*arg2 < 2^32.
+// Es wird vorausgesetzt, daß arg1*arg2 < 2^32.
#if defined(__GNUC__) && defined(__sparc64__) && !defined(NO_ASM)
inline uint32 mulu32_unchecked (uint32 arg1, uint32 arg2)
{
#elif defined(__sparc__) && !defined(NO_ASM)
extern "C" uint32 mulu32_unchecked (uint32 x, uint32 y); // extern in Assembler
#else
- // Wir können dafür auch die Bibliotheksroutine des C-Compilers nehmen:
+ // Wir können dafür auch die Bibliotheksroutine des C-Compilers nehmen:
inline uint32 mulu32_unchecked (uint32 arg1, uint32 arg2)
{
return arg1 * arg2;
extern "C" uint32 mulu32_high; // -> High-Teil
#if defined(__GNUC__) && defined(__m68k__) && !defined(NO_ASM)
#define mulu32(x,y,hi_zuweisung,lo_zuweisung) \
- ({ var uint32 _x = (x); \
- var uint32 _y = (y); \
- var uint32 _hi; \
- var uint32 _lo; \
+ ({ var uint32 _x = (x); \
+ var uint32 _y = (y); \
+ var uint32 _hi; \
+ var uint32 _lo; \
__asm__("mulul %3,%0:%1" : "=d" (_hi), "=d"(_lo) : "1" (_x), "dm" (_y) ); \
- hi_zuweisung _hi; \
- lo_zuweisung _lo; \
+ unused (hi_zuweisung _hi); \
+ lo_zuweisung _lo; \
})
#elif defined(__GNUC__) && defined(__m68k__)
#define mulu32(x,y,hi_zuweisung,lo_zuweisung) \
_hi += high16(_mid); _mid = highlow32_0(low16(_mid)); \
_lo += _mid; if (_lo < _mid) { _hi += 1; } /* 64-Bit-Addition */\
} \
- hi_zuweisung _hi; \
+ unused (hi_zuweisung _hi); \
lo_zuweisung _lo; \
})
#elif defined(__GNUC__) && defined(__sparc64__) && !defined(NO_ASM)
: "=r" (_prod) \
: "r" ((uint32)(x)), "r" ((uint32)(y)) \
); \
- hi_zuweisung (uint32)(_prod>>32); \
+ unused (hi_zuweisung (uint32)(_prod>>32)); \
lo_zuweisung (uint32)(_prod); \
})
#elif defined(__GNUC__) && defined(__sparc__) && !defined(NO_ASM)
#define mulu32(x,y,hi_zuweisung,lo_zuweisung) \
({ lo_zuweisung mulu32_(x,y); /* extern in Assembler */ \
{var register uint32 _hi __asm__("%g1"); \
- hi_zuweisung _hi; \
+ unused (hi_zuweisung _hi); \
}})
#elif defined(__GNUC__) && defined(__arm__) && 0 // see comment cl_asm_arm.cc
#define mulu32(x,y,hi_zuweisung,lo_zuweisung) \
({ lo_zuweisung mulu32_(x,y); /* extern in Assembler */ \
{var register uint32 _hi __asm__("%r1"/*"%a2"*/); \
- hi_zuweisung _hi; \
+ unused (hi_zuweisung _hi); \
}})
#elif defined(__GNUC__) && (defined(__i386__) || defined(__x86_64__)) && !defined(NO_ASM)
#define mulu32(x,y,hi_zuweisung,lo_zuweisung) \
: "=d" /* %edx */ (_hi), "=a" /* %eax */ (_lo) \
: "g" ((uint32)(x)), "1" /* %eax */ ((uint32)(y)) \
); \
- hi_zuweisung _hi; lo_zuweisung _lo; \
+ unused (hi_zuweisung _hi); lo_zuweisung _lo; \
})
#elif defined(__GNUC__) && defined(__mips__) && !defined(NO_ASM)
#define mulu32(x,y,hi_zuweisung,lo_zuweisung) \
: "=r" (_hi), "=r" (_lo) \
: "r" ((uint32)(x)), "r" ((uint32)(y)) \
); \
- hi_zuweisung _hi; lo_zuweisung _lo; \
+ unused (hi_zuweisung _hi); lo_zuweisung _lo; \
})
#elif defined(__GNUC__) && defined(HAVE_LONGLONG) && !defined(__arm__)
#define mulu32(x,y,hi_zuweisung,lo_zuweisung) \
({ var register uint64 _prod = (uint64)(uint32)(x) * (uint64)(uint32)(y); \
- hi_zuweisung (uint32)(_prod>>32); \
+ unused (hi_zuweisung (uint32)(_prod>>32)); \
lo_zuweisung (uint32)(_prod); \
})
#elif defined(WATCOM) && defined(__i386__) && !defined(NO_ASM)
#define mulu32(x,y,hi_zuweisung,lo_zuweisung) \
- { var register uint32 _hi; \
- var register uint32 _lo; \
- _lo = mulu32_(x,y), _hi = mulu32_high_(); \
- hi_zuweisung _hi; lo_zuweisung _lo; \
+ { var register uint32 _hi; \
+ var register uint32 _lo; \
+ _lo = mulu32_(x,y), _hi = mulu32_high_(); \
+ unused (hi_zuweisung _hi); lo_zuweisung _lo; \
}
extern "C" uint32 mulu32_high_ (void);
#pragma aux mulu32_ = 0xF7 0xE2 /* mull %edx */ parm [eax] [edx] value [eax] modify [eax edx];
#pragma aux mulu32_high_ = /* */ value [edx] modify [];
#else
#define mulu32(x,y,hi_zuweisung,lo_zuweisung) \
- { lo_zuweisung mulu32_(x,y); hi_zuweisung mulu32_high; }
+ { lo_zuweisung mulu32_(x,y); unused (hi_zuweisung mulu32_high); }
#if (defined(__m68k__) || defined(__sparc__) || defined(__sparc64__) || defined(__arm__) || (defined(__i386__) && !defined(WATCOM) && !defined(MICROSOFT)) || defined(__x86_64__) || defined(__mips__) || defined(__hppa__)) && !defined(NO_ASM)
// mulu32_ extern in Assembler
#if defined(__sparc__) || defined(__sparc64__)
extern "C" uint32 _get_g1 (void);
- #define mulu32_high (_get_g1()) // Rückgabe im Register %g1
+ #define mulu32_high (_get_g1()) // Rückgabe im Register %g1
#elif !defined(__hppa__)
#define NEED_VAR_mulu32_high
#endif
#else
#define mulu64(x,y,hi_zuweisung,lo_zuweisung) \
{ lo_zuweisung mulu64_(x,y); hi_zuweisung mulu64_high; }
- #if defined(__sparc64__)
+ #if defined(__sparc64__) && !defined(NO_ASM)
// mulu64_ extern in Assembler
- #if defined(__sparc64__)
- extern "C" uint64 _get_g2 (void);
- #define mulu64_high (_get_g2()) // Rückgabe im Register %g2
- #else
- #define NEED_VAR_mulu64_high
- #endif
+ extern "C" uint64 _get_g2 (void);
+ #define mulu64_high (_get_g2()) // Rückgabe im Register %g2
#else
#define NEED_FUNCTION_mulu64_
#endif
// Dividiert eine 16-Bit-Zahl durch eine 16-Bit-Zahl und
// liefert einen 16-Bit-Quotienten und einen 16-Bit-Rest.
// divu_1616_1616(x,y,q=,r=);
-// > uint16 x: Zähler
+// > uint16 x: Zähler
// > uint16 y: Nenner
// < uint16 q: floor(x/y)
// < uint16 r: x mod y
// Dividiert eine 32-Bit-Zahl durch eine 16-Bit-Zahl und
// liefert einen 16-Bit-Quotienten und einen 16-Bit-Rest.
// divu_3216_1616(x,y,q=,r=);
-// > uint32 x: Zähler
+// > uint32 x: Zähler
// > uint16 y: Nenner
-// > Es sei bekannt, daß 0 <= x < 2^16*y .
+// > Es sei bekannt, daß 0 <= x < 2^16*y .
// < uint16 q: floor(x/y)
// < uint16 r: x mod y
// < x = q*y+r
__asm__ __volatile__ ( \
"wr %%g0,%%g0,%%y\n\t" \
"udiv %2,%3,%0\n\t" \
- "umul %0,%3,%1" \
+ "umul %0,%3,%1\n\t" \
"sub %2,%1,%1" \
: "=&r" (__q), "=&r" (__r) \
: "r" (__x), "r" (__y)); \
// Dividiert eine 32-Bit-Zahl durch eine 16-Bit-Zahl und
// liefert einen 32-Bit-Quotienten und einen 16-Bit-Rest.
// divu_3216_3216(x,y,q=,r=);
-// > uint32 x: Zähler
+// > uint32 x: Zähler
// > uint16 y: Nenner
-// Es sei bekannt, daß y>0.
+// Es sei bekannt, daß y>0.
// < uint32 q: floor(x/y)
// < uint16 r: x mod y
// < x = q*y+r
__asm__ __volatile__ ( \
"wr %%g0,%%g0,%%y\n\t" \
"udiv %2,%3,%0\n\t" \
- "umul %0,%3,%1" \
+ "umul %0,%3,%1\n\t" \
"sub %2,%1,%1" \
: "=&r" (__q), "=&r" (__r) \
: "r" (__x), "r" (__y)); \
// Dividiert eine 32-Bit-Zahl durch eine 32-Bit-Zahl und
// liefert einen 32-Bit-Quotienten und einen 32-Bit-Rest.
// divu_3232_3232(x,y,q=,r=);
-// > uint32 x: Zähler
+// > uint32 x: Zähler
// > uint32 y: Nenner
-// Es sei bekannt, daß y>0.
+// Es sei bekannt, daß y>0.
// < uint32 q: floor(x/y)
// < uint32 r: x mod y
// < x = q*y+r
__asm__ __volatile__ ( \
"wr %%g0,%%g0,%%y\n\t" \
"udiv %2,%3,%0\n\t" \
- "umul %0,%3,%1" \
+ "umul %0,%3,%1\n\t" \
"sub %2,%1,%1" \
: "=&r" (__q), "=&r" (__r) \
: "r" (__x), "r" (__y)); \
q_zuweisung (uint32)__q; \
r_zuweisung (uint32)__r; \
})
+ #define divu_3232_3232_(x,y) divu_6432_3232_(0,x,y)
#elif defined(__sparc__) || defined(__sparc64__) || defined(__i386__) || defined(__x86_64__)
#define divu_3232_3232(x,y,q_zuweisung,r_zuweisung) \
divu_6432_3232(0,x,y,q_zuweisung,r_zuweisung)
// y habe genau n+k Bits (1 <= k <= n), d.h. 2^(n+k-1) <= y < 2^(n+k).
// Schreibe x = 2^k*x1 + x0 mit x1 := floor(x/2^k)
// und y = 2^k*y1 + y0 mit y1 := floor(y/2^k)
- // und bilde den Näherungs-Quotienten floor(x1/y1)
+ // und bilde den Näherungs-Quotienten floor(x1/y1)
// oder (noch besser) floor(x1/(y1+1)).
// Wegen 0 <= x1 < 2^(2n) und 0 < 2^(n-1) <= y1 < 2^n
// und x1/(y1+1) <= x/y < x1/(y1+1) + 2
// gilt floor(x1/(y1+1)) <= floor(x/y) <= floor(x1/(y1+1)) + 2 .
// Man bildet also q:=floor(x1/(y1+1)) (ein Shift um n Bit oder
// eine (2n)-durch-n-Bit-Division, mit Ergebnis q <= floor(x/y) < beta)
- // und x-q*y und muß hiervon noch höchstens 2 mal y abziehen und q
+ // und x-q*y und muß hiervon noch höchstens 2 mal y abziehen und q
// incrementieren, um den Quotienten q = floor(x/y) und den Rest
// x-floor(x/y)*y der Division zu bekommen.
#define divu_3232_3232(x,y,q_zuweisung,r_zuweisung) \
{ var uint32 _x1 = _x; /* x1 := x */ \
var uint32 _y1 = _y; /* y1 := y */ \
var uint16 _q; \
- do { _x1 = floor(_x1,2); _y1 = floor(_y1,2); } /* k erhöhen */\
+ do { _x1 = floor(_x1,2); _y1 = floor(_y1,2); } /* k erhöhen */\
until (_y1 <= (uint32)(bit(16)-1)); /* bis y1 < beta */ \
{ var uint16 _y2 = low16(_y1)+1; /* y1+1 bilden */ \
if (_y2==0) \
{ divu_3216_1616(_x1,_y2,_q=,); } /* Division von x1 durch y1+1 */\
} \
/* _q = q = floor(x1/(y1+1)) */ \
- /* x-q*y bilden (eine 16-mal-32-Bit-Multiplikation ohne Überlauf): */\
+ /* x-q*y bilden (eine 16-mal-32-Bit-Multiplikation ohne Überlauf): */\
_x -= highlow32_0(mulu16(_q,high16(_y))); /* q * high16(y) * beta */\
/* gefahrlos, da q*high16(y) <= q*y/beta <= x/beta < beta */ \
_x -= mulu16(_q,low16(_y)); /* q * low16(y) */ \
/* gefahrlos, da q*high16(y)*beta + q*low16(y) = q*y <= x */ \
- /* Noch höchstens 2 mal y abziehen: */ \
+ /* Noch höchstens 2 mal y abziehen: */ \
if (_x >= _y) \
{ _q += 1; _x -= _y; \
if (_x >= _y) \
// Dividiert eine 64-Bit-Zahl durch eine 32-Bit-Zahl und
// liefert einen 32-Bit-Quotienten und einen 32-Bit-Rest.
// divu_6432_3232(xhi,xlo,y,q=,r=);
-// > uint32 xhi,xlo: x = 2^32*xhi+xlo = Zähler
+// > uint32 xhi,xlo: x = 2^32*xhi+xlo = Zähler
// > uint32 y: Nenner
-// > Es sei bekannt, daß 0 <= x < 2^32*y .
+// > Es sei bekannt, daß 0 <= x < 2^32*y .
// < uint32 q: floor(x/y)
// < uint32 r: x mod y
// < x = q*y+r
__asm__ __volatile__ ( \
"wr %2,%%g0,%%y\n\t" \
"udiv %3,%4,%0\n\t" \
- "umul %0,%4,%1" \
+ "umul %0,%4,%1\n\t" \
"sub %3,%1,%1" \
: "=&r" (__q), "=&r" (__r) \
: "r" (__xhi), "r" (__xlo), "r" (__y)); \
// divu_6432_3232_ extern in Assembler
#if defined(__sparc__) || defined(__sparc64__)
extern "C" uint32 _get_g1 (void);
- #define divu_32_rest (_get_g1()) // Rückgabe im Register %g1
+ #define divu_32_rest (_get_g1()) // Rückgabe im Register %g1
#else
#define NEED_VAR_divu_32_rest
#endif
// Dividiert eine 64-Bit-Zahl durch eine 32-Bit-Zahl und
// liefert einen 32-Bit-Quotienten und einen 32-Bit-Rest.
// divu_6432_3232_w(x,y,q=,r=);
-// > uint64 x: Zähler
+// > uint64 x: Zähler
// > uint32 y: Nenner
-// > Es sei bekannt, daß 0 <= x < 2^32*y .
+// > Es sei bekannt, daß 0 <= x < 2^32*y .
// < uint32 q: floor(x/y)
// < uint32 r: x mod y
// < x = q*y+r
__asm__ __volatile__ ( \
"wr %2,%%g0,%%y\n\t" \
"udiv %3,%4,%0\n\t" \
- "umul %0,%4,%1" \
+ "umul %0,%4,%1\n\t" \
"sub %3,%1,%1" \
: "=&r" (__q), "=&r" (__r) \
: "r" (__xhi), "r" (__xlo), "r" (__y)); \
// Dividiert eine 64-Bit-Zahl durch eine 32-Bit-Zahl und
// liefert einen 64-Bit-Quotienten und einen 32-Bit-Rest.
// divu_6432_6432(x,y,q=,r=);
-// > uint64 x: Zähler
+// > uint64 x: Zähler
// > uint32 y: Nenner
-// > Es sei bekannt, daß y>0.
+// > Es sei bekannt, daß y>0.
// < uint64 q: floor(x/y)
// < uint32 r: x mod y
// < x = q*y+r
// Dividiert eine 64-Bit-Zahl durch eine 64-Bit-Zahl und
// liefert einen 64-Bit-Quotienten und einen 64-Bit-Rest.
// divu_6464_6464(x,y,q=,r=);
-// > uint64 x: Zähler
+// > uint64 x: Zähler
// > uint64 y: Nenner
-// > Es sei bekannt, daß y>0.
+// > Es sei bekannt, daß y>0.
// < uint64 q: floor(x/y)
// < uint64 r: x mod y
// < x = q*y+r
// Dividiert eine 128-Bit-Zahl durch eine 64-Bit-Zahl und
// liefert einen 64-Bit-Quotienten und einen 64-Bit-Rest.
// divu_12864_6464(xhi,xlo,y,q=,r=);
-// > uint64 xhi,xlo: x = 2^64*xhi+xlo = Zähler
+// > uint64 xhi,xlo: x = 2^64*xhi+xlo = Zähler
// > uint64 y: Nenner
-// > Es sei bekannt, daß 0 <= x < 2^64*y .
+// > Es sei bekannt, daß 0 <= x < 2^64*y .
// < uint64 q: floor(x/y)
// < uint64 r: x mod y
// < x = q*y+r
// < boolean sqrtp: /=0, falls x=y^2
// Methode:
// y := 2^16 als Anfangswert,
- // y := floor((y + floor(x/y))/2) als nächster Wert,
+ // y := floor((y + floor(x/y))/2) als nächster Wert,
// solange z := floor(x/y) < y, setze y := floor((y+z)/2).
// y ist fertig; x=y^2 genau dann, wenn z=y und die letzte Division aufging.
// (Beweis:
// 1. Die Folge der y ist streng monoton fallend.
- // 2. Stets gilt y >= floor(sqrt(x)) (denn für alle y>0 ist
+ // 2. Stets gilt y >= floor(sqrt(x)) (denn für alle y>0 ist
// y + x/y >= 2*sqrt(x) und daher floor((y + floor(x/y))/2) =
// floor(y/2 + x/(2*y)) >= floor(sqrt(x)) ).
- // 3. Am Schluß gilt x >= y^2.
+ // 3. Am Schluß gilt x >= y^2.
// )
#define isqrt_32_16(x,y_zuweisung,sqrtp_zuweisung) \
{ var uint32 _x = (x); \
loop \
{ var uint16 _z; \
var uint16 _r; \
- if (_x1 >= _y) /* Division _x/_y ergäbe Überlauf -> _z > _y */\
+ if (_x1 >= _y) /* Division _x/_y ergäbe Überlauf -> _z > _y */\
{ unused (sqrtp_zuweisung FALSE); break; } \
divu_3216_1616(_x,_y, _z=,_r=); /* Dividiere _x/_y */ \
if (_z >= _y) \
{ unused (sqrtp_zuweisung (_z == _y) && (_r == 0)); break; } \
- _y = floor((uint16)(_z+_y),2) | bit(16-1); /* _y muß >= 2^15 bleiben */\
+ _y = floor((uint16)(_z+_y),2) | bit(16-1); /* _y muß >= 2^15 bleiben */\
} \
y_zuweisung _y; \
}
#if defined(__sparc__) || defined(__sparc64__) || defined(__m68k__) || defined(__hppa__)
// Methode:
// y := 2^32 als Anfangswert,
- // y := floor((y + floor(x/y))/2) als nächster Wert,
+ // y := floor((y + floor(x/y))/2) als nächster Wert,
// solange z := floor(x/y) < y, setze y := floor((y+z)/2).
// y ist fertig; x=y^2 genau dann, wenn z=y und die letzte Division aufging.
// (Beweis:
// 1. Die Folge der y ist streng monoton fallend.
- // 2. Stets gilt y >= floor(sqrt(x)) (denn für alle y>0 ist
+ // 2. Stets gilt y >= floor(sqrt(x)) (denn für alle y>0 ist
// y + x/y >= 2*sqrt(x) und daher floor((y + floor(x/y))/2) =
// floor(y/2 + x/(2*y)) >= floor(sqrt(x)) ).
- // 3. Am Schluß gilt x >= y^2.
+ // 3. Am Schluß gilt x >= y^2.
// )
#define isqrt_64_32(xhi,xlo,y_zuweisung,sqrtp_zuweisung) \
{ var uint32 _xhi = (xhi); \
loop \
{ var uint32 _z; \
var uint32 _rest; \
- if (_xhi >= _y) /* Division _x/_y ergäbe Überlauf -> _z > _y */\
+ if (_xhi >= _y) /* Division _x/_y ergäbe Überlauf -> _z > _y */\
{ sqrtp_zuweisung FALSE; break; } \
divu_6432_3232(_xhi,_xlo,_y, _z=,_rest=); /* Dividiere _x/_y */\
if (_z >= _y) \
{ sqrtp_zuweisung (_z == _y) && (_rest == 0); break; } \
- _y = floor(_z+_y,2) | bit(32-1); /* _y muß >= 2^31 bleiben */ \
+ _y = floor(_z+_y,2) | bit(32-1); /* _y muß >= 2^31 bleiben */ \
} \
y_zuweisung _y; \
}
// Methode:
// Wie bei UDS_sqrt mit n=2.
// y = 2^16*yhi + ylo ansetzen.
- // Dann muß
+ // Dann muß
// yhi = floor(y/2^16) = floor(floor(sqrt(x))/2^16)
// = floor(sqrt(x)/2^16) = floor(sqrt(x/2^32)) = isqrt(xhi)
// sein. Es folgt yhi >= 2^15.
- // Danach sucht man das größte ylo >=0 mit
+ // Danach sucht man das größte ylo >=0 mit
// x - 2^32*yhi^2 >= 2*2^16*yhi*ylo + ylo^2.
// Dazu setzen wir xhi*2^32+xlo := x - 2^32*yhi^2
// (also xhi := xhi - yhi^2, das ist >=0, <=2*yhi).
- // Die Schätzung für die zweite Ziffer
+ // Die Schätzung für die zweite Ziffer
// ylo' := min(2^16-1,floor((xhi*2^32+xlo)/(2*2^16*yhi)))
- // erfüllt ylo'-1 <= ylo <= ylo', ist also um höchstens 1 zu groß.
+ // erfüllt ylo'-1 <= ylo <= ylo', ist also um höchstens 1 zu groß.
// (Beweis: Rechte Ungleichung klar, da ylo < 2^16 und
// xhi*2^32+xlo >= 2*2^16*yhi*ylo + ylo^2 >= 2*2^16*yhi*ylo
- // ==> (xhi*2^32+xlo)/(2*2^16*yhi) >= ylo gelten muß.
+ // ==> (xhi*2^32+xlo)/(2*2^16*yhi) >= ylo gelten muß.
// Linke Ungleichung: Falls floor(...)>=2^16, ist
// xhi*2^32+xlo >= 2*2^16*2^16*yhi >= 2*2^16*yhi*(2^16-1) + 2^32
// >= 2*2^16*yhi*(2^16-1) + (2^16-1)^2
isqrt_32_16(_xhi,_yhi=,); /* yhi := isqrt(xhi) */ \
_xhi -= mulu16(_yhi,_yhi); /* jetzt 0 <= xhi <= 2*yhi */ \
/* x = 2^32*yhi^2 + 2^32*xhi + xlo */ \
- /* Schätzung für die zweite Ziffer berechnen: */ \
+ /* Schätzung für die zweite Ziffer berechnen: */ \
/* ylo := min(2^16-1,floor((xhi*2^32+xlo)/(2*2^16*yhi))) bilden: */\
{var uint32 _z = (_xhi << 15) | (_xlo >> 17); /* < 2^15*(2*yhi+1) */\
var uint32 _r = highlow32_0(_yhi); \
/* Versuche vom Rest 2^32*floor(r/2^15) + xlo z zu subtrahieren. */\
/* Falls Rest >= z (d.h. r>=2^15 oder xlo>=z), ist ylo fertig, */ \
/* und es gilt x=y^2 genau dann, wenn r<2^15 und xlo=z. */ \
- /* Sonst (d.h. r<2^15 und xlo<z), muß man ylo erniedrigen. Dazu */\
+ /* Sonst (d.h. r<2^15 und xlo<z), muß man ylo erniedrigen. Dazu */\
/* setzt man ylo := ylo-1, z := z-(2*ylo+1), */ \
/* Rest := Rest + 2^17*yhi = xlo + 2^17*yhi >= 2^32 > z, also x>y^2. */\
if (_r < bit(15)) \
// Methode:
// Wie bei UDS_sqrt mit n=2.
// y = 2^32*yhi + ylo ansetzen.
- // Dann muß
+ // Dann muß
// yhi = floor(y/2^32) = floor(floor(sqrt(x))/2^32)
// = floor(sqrt(x)/2^32) = floor(sqrt(x/2^64)) = isqrt(xhi)
// sein. Es folgt yhi >= 2^31.
- // Danach sucht man das größte ylo >=0 mit
+ // Danach sucht man das größte ylo >=0 mit
// x - 2^64*yhi^2 >= 2*2^32*yhi*ylo + ylo^2.
// Dazu setzen wir xhi*2^64+xlo := x - 2^64*yhi^2
// (also xhi := xhi - yhi^2, das ist >=0, <=2*yhi).
- // Die Schätzung für die zweite Ziffer
+ // Die Schätzung für die zweite Ziffer
// ylo' := min(2^32-1,floor((xhi*2^64+xlo)/(2*2^32*yhi)))
- // erfüllt ylo'-1 <= ylo <= ylo', ist also um höchstens 1 zu groß.
+ // erfüllt ylo'-1 <= ylo <= ylo', ist also um höchstens 1 zu groß.
// (Beweis: Rechte Ungleichung klar, da ylo < 2^32 und
// xhi*2^64+xlo >= 2*2^32*yhi*ylo + ylo^2 >= 2*2^32*yhi*ylo
- // ==> (xhi*2^64+xlo)/(2*2^32*yhi) >= ylo gelten muß.
+ // ==> (xhi*2^64+xlo)/(2*2^32*yhi) >= ylo gelten muß.
// Linke Ungleichung: Falls floor(...)>=2^32, ist
// xhi*2^64+xlo >= 2*2^32*2^32*yhi >= 2*2^32*yhi*(2^32-1) + 2^64
// >= 2*2^32*yhi*(2^32-1) + (2^32-1)^2
isqrt_64_32(high32(xhi),low32(xhi),yhi=,); /* yhi := isqrt(xhi) */\
xhi -= mulu32_w(yhi,yhi); /* jetzt 0 <= xhi <= 2*yhi */ \
/* x = 2^64*yhi^2 + 2^64*xhi + xlo */ \
- /* Schätzung für die zweite Ziffer berechnen: */ \
+ /* Schätzung für die zweite Ziffer berechnen: */ \
/* ylo := min(2^32-1,floor((xhi*2^64+xlo)/(2*2^32*yhi))) bilden: */\
{var uint64 z = (xhi << 31) | (xlo >> 33); /* < 2^31*(2*yhi+1) */ \
var uint64 r = highlow64_0(yhi); \
/* Versuche vom Rest 2^64*floor(r/2^31) + xlo z zu subtrahieren. */\
/* Falls Rest >= z (d.h. r>=2^31 oder xlo>=z), ist ylo fertig, */ \
/* und es gilt x=y^2 genau dann, wenn r<2^31 und xlo=z. */ \
- /* Sonst (d.h. r<2^31 und xlo<z), muß man ylo erniedrigen. Dazu */\
+ /* Sonst (d.h. r<2^31 und xlo<z), muß man ylo erniedrigen. Dazu */\
/* setzt man ylo := ylo-1, z := z-(2*ylo+1), */ \
/* Rest := Rest + 2^33*yhi = xlo + 2^33*yhi >= 2^64 > z, also x>y^2. */\
if (r < bit(31)) \
// < uintL ergebnis : Wurzel, >=0, <2^16
extern uintL isqrt (uintL x);
+#ifdef HAVE_LONGLONG
+// Extracts integer root of a 64-bit number and returns a 32-bit number.
+// isqrt(x)
+// > uintQ x : radicand, >=0, <2^64
+// < uintL result : square root, >=0, <2^32
+ extern uintL isqrt (uintQ x);
+#endif
+
+// Sorry for this. We need an isqrt function taking uintC arguments but we
+// cannot use overloading since this would lead to ambiguities with any of the
+// two signatures above.
+ inline uintL isqrtC (uintC x)
+ {
+#if (intCsize==32)
+ return isqrt((uintL)x);
+#else
+ return isqrt((uintQ)x);
+#endif
+ }
+
+
// Zieht die Ganzzahl-Wurzel aus einer 64-Bit-Zahl und
// liefert eine 32-Bit-Wurzel.
// isqrt(x1,x0)
extern uintL isqrt (uintL x1, uintL x0);
-// Bits einer 8-Bit-Zahl zählen:
+// Bits einer 8-Bit-Zahl zählen:
// integerlength8(digit,size=);
-// setzt size auf die höchste in digit vorkommende Bitnummer.
+// setzt size auf die höchste in digit vorkommende Bitnummer.
// > digit: ein uint8 >0
// < size: >0, <=8, mit 2^(size-1) <= digit < 2^size
#if defined(__GNUC__) && defined(__m68k__) && !defined(NO_ASM)
#define integerlength8(digit,size_zuweisung) \
- { var uintL _zero_counter; /* zählt die führenden Nullbits in digit */\
+ { var uintL _zero_counter; /* zählt die führenden Nullbits in digit */\
__asm__("bfffo %1{#0:#8},%0" : "=d" (_zero_counter) : "dm" ((uint8)(digit)) ); \
size_zuweisung (8-_zero_counter); \
}
#define integerlength8(digit,size_zuweisung) \
{ var uintC _bitsize = 1; \
var uintL _x8 = (uint8)(digit); \
- /* _x8 hat höchstens 8 Bits. */\
+ /* _x8 hat höchstens 8 Bits. */\
if (_x8 >= bit(4)) { _x8 = _x8>>4; _bitsize += 4; } \
- /* _x8 hat höchstens 4 Bits. */\
+ /* _x8 hat höchstens 4 Bits. */\
if (_x8 >= bit(2)) { _x8 = _x8>>2; _bitsize += 2; } \
- /* _x8 hat höchstens 2 Bits. */\
+ /* _x8 hat höchstens 2 Bits. */\
if (_x8 >= bit(1)) { /* _x8 = _x8>>1; */ _bitsize += 1; } \
- /* _x8 hat höchstens 1 Bit. Dieses Bit muß gesetzt sein. */\
+ /* _x8 hat höchstens 1 Bit. Dieses Bit muß gesetzt sein. */\
size_zuweisung _bitsize; \
}
#endif
-// Bits einer 16-Bit-Zahl zählen:
+// Bits einer 16-Bit-Zahl zählen:
// integerlength16(digit,size=);
-// setzt size auf die höchste in digit vorkommende Bitnummer.
+// setzt size auf die höchste in digit vorkommende Bitnummer.
// > digit: ein uint16 >0
// < size: >0, <=16, mit 2^(size-1) <= digit < 2^size
#if defined(__GNUC__) && defined(__m68k__) && !defined(NO_ASM)
#define integerlength16(digit,size_zuweisung) \
- { var uintL _zero_counter; /* zählt die führenden Nullbits in digit */\
+ { var uintL _zero_counter; /* zählt die führenden Nullbits in digit */\
__asm__("bfffo %1{#0:#16},%0" : "=d" (_zero_counter) : "dm" ((uint16)(digit)) ); \
size_zuweisung (16-_zero_counter); \
}
integerlength32((uint32)(digit),size_zuweisung) // siehe unten
#elif defined(__GNUC__) && defined(__i386__) && !defined(NO_ASM)
#define integerlength16(digit,size_zuweisung) \
- { var uintW _one_position; /* Position der führenden 1 */\
+ { var uintW _one_position; /* Position der führenden 1 */\
__asm__("bsrw %1,%0" : "=r" (_one_position) : "r" ((uint16)(digit)) ); \
size_zuweisung (1+_one_position); \
}
// Die weiteren kommen von gcc/longlong.h :
#elif defined(__GNUC__) && defined(__ibm032__) && !defined(NO_ASM) // RT/ROMP
#define integerlength16(digit,size_zuweisung) \
- { var uintL _zero_counter; /* zählt die führenden Nullbits in digit */\
+ { var uintL _zero_counter; /* zählt die führenden Nullbits in digit */\
__asm__("clz %0,%1" : "=r" (_zero_counter) : "r" ((uint32)(digit)) ); \
size_zuweisung (16-_zero_counter); \
}
#define integerlength16(digit,size_zuweisung) \
{ var uintC _bitsize = 1; \
var uintWL _x16 = (uint16)(digit); \
- /* _x16 hat höchstens 16 Bits. */\
+ /* _x16 hat höchstens 16 Bits. */\
if (_x16 >= bit(8)) { _x16 = _x16>>8; _bitsize += 8; } \
- /* _x16 hat höchstens 8 Bits. */\
+ /* _x16 hat höchstens 8 Bits. */\
if (_x16 >= bit(4)) { _x16 = _x16>>4; _bitsize += 4; } \
- /* _x16 hat höchstens 4 Bits. */\
+ /* _x16 hat höchstens 4 Bits. */\
if (_x16 >= bit(2)) { _x16 = _x16>>2; _bitsize += 2; } \
- /* _x16 hat höchstens 2 Bits. */\
+ /* _x16 hat höchstens 2 Bits. */\
if (_x16 >= bit(1)) { /* _x16 = _x16>>1; */ _bitsize += 1; } \
- /* _x16 hat höchstens 1 Bit. Dieses Bit muß gesetzt sein. */\
+ /* _x16 hat höchstens 1 Bit. Dieses Bit muß gesetzt sein. */\
size_zuweisung _bitsize; \
}
#endif
-// Bits einer 32-Bit-Zahl zählen:
+// Bits einer 32-Bit-Zahl zählen:
// integerlength32(digit,size=);
-// setzt size auf die höchste in digit vorkommende Bitnummer.
+// setzt size auf die höchste in digit vorkommende Bitnummer.
// > digit: ein uint32 >0
// < size: >0, <=32, mit 2^(size-1) <= digit < 2^size
#if defined(__GNUC__) && defined(__m68k__) && !defined(NO_ASM)
#define integerlength32(digit,size_zuweisung) \
- { var uintL _zero_counter; /* zählt die führenden Nullbits in digit */\
+ { var uintL _zero_counter; /* zählt die führenden Nullbits in digit */\
__asm__("bfffo %1{#0:#32},%0" : "=d" (_zero_counter) : "dm" ((uint32)(digit)) ); \
size_zuweisung (32-_zero_counter); \
}
}
#elif defined(__GNUC__) && defined(__i386__) && !defined(NO_ASM)
#define integerlength32(digit,size_zuweisung) \
- { var uintL _one_position; /* Position der führenden 1 */\
+ { var uintL _one_position; /* Position der führenden 1 */\
__asm__("bsrl %1,%0" : "=r" (_one_position) : "rm" ((uint32)(digit)) ); \
size_zuweisung (1+_one_position); \
}
// Die weiteren kommen von gcc/longlong.h :
#elif defined(__GNUC__) && (defined(__a29k__) || defined(___AM29K__)) && !defined(NO_ASM)
#define integerlength32(digit,size_zuweisung) \
- { var uintL _zero_counter; /* zählt die führenden Nullbits in digit */\
+ { var uintL _zero_counter; /* zählt die führenden Nullbits in digit */\
__asm__("clz %0,%1" : "=r" (_zero_counter) : "r" ((uint32)(digit)) ); \
size_zuweisung (32-_zero_counter); \
}
#elif defined(__GNUC__) && defined(__gmicro__) && !defined(NO_ASM)
#define integerlength32(digit,size_zuweisung) \
- { var uintL _zero_counter; /* zählt die führenden Nullbits in digit */\
+ { var uintL _zero_counter; /* zählt die führenden Nullbits in digit */\
__asm__("bsch/1 %1,%0" : "=g" (_zero_counter) : "g" ((uint32)(digit)) ); \
size_zuweisung (32-_zero_counter); \
}
#ifdef _AIX
// old assembler syntax
#define integerlength32(digit,size_zuweisung) \
- { var uintL _zero_counter; /* zählt die führenden Nullbits in digit */\
+ { var uintL _zero_counter; /* zählt die führenden Nullbits in digit */\
__asm__("cntlz %0,%1" : "=r" (_zero_counter) : "r" ((uint32)(digit)) ); \
size_zuweisung (32-_zero_counter); \
}
#else
// new assembler syntax
#define integerlength32(digit,size_zuweisung) \
- { var uintL _zero_counter; /* zählt die führenden Nullbits in digit */\
+ { var uintL _zero_counter; /* zählt die führenden Nullbits in digit */\
__asm__("cntlzw %0,%1" : "=r" (_zero_counter) : "r" ((uint32)(digit)) ); \
size_zuweisung (32-_zero_counter); \
}
#endif
#elif defined(__GNUC__) && defined(__m88k__) && !defined(NO_ASM)
#define integerlength32(digit,size_zuweisung) \
- { var uintL _one_position; /* Position der führenden 1 */\
+ { var uintL _one_position; /* Position der führenden 1 */\
__asm__("ff1 %0,%1" : "=r" (_one_position) : "r" ((uint32)(digit)) ); \
size_zuweisung (1+_one_position); \
}
#define integerlength32(digit,size_zuweisung) \
{ var uintC _bitsize = 1; \
var uintL _x32 = (uint32)(digit); \
- /* _x32 hat höchstens 32 Bits. */\
+ /* _x32 hat höchstens 32 Bits. */\
if (_x32 >= bit(16)) { _x32 = _x32>>16; _bitsize += 16; } \
- /* _x32 hat höchstens 16 Bits. */\
+ /* _x32 hat höchstens 16 Bits. */\
if (_x32 >= bit(8)) { _x32 = _x32>>8; _bitsize += 8; } \
- /* _x32 hat höchstens 8 Bits. */\
+ /* _x32 hat höchstens 8 Bits. */\
if (_x32 >= bit(4)) { _x32 = _x32>>4; _bitsize += 4; } \
- /* _x32 hat höchstens 4 Bits. */\
+ /* _x32 hat höchstens 4 Bits. */\
if (_x32 >= bit(2)) { _x32 = _x32>>2; _bitsize += 2; } \
- /* _x32 hat höchstens 2 Bits. */\
+ /* _x32 hat höchstens 2 Bits. */\
if (_x32 >= bit(1)) { /* _x32 = _x32>>1; */ _bitsize += 1; } \
- /* _x32 hat höchstens 1 Bit. Dieses Bit muß gesetzt sein. */\
+ /* _x32 hat höchstens 1 Bit. Dieses Bit muß gesetzt sein. */\
size_zuweisung _bitsize; \
}
#define GENERIC_INTEGERLENGTH32
#endif
-// Bits einer 64-Bit-Zahl zählen:
+// Bits einer 64-Bit-Zahl zählen:
// integerlength64(digit,size=);
-// setzt size auf die höchste in digit vorkommende Bitnummer.
+// setzt size auf die höchste in digit vorkommende Bitnummer.
// > digit: ein uint64 >0
// < size: >0, <=64, mit 2^(size-1) <= digit < 2^size
#ifdef GENERIC_INTEGERLENGTH32
#define integerlength64(digit,size_zuweisung) \
{ var uintC _bitsize = 1; \
var uint64 _x64 = (uint64)(digit); \
- /* _x64 hat höchstens 64 Bits. */\
+ /* _x64 hat höchstens 64 Bits. */\
if (_x64 >= bit(32)) { _x64 = _x64>>32; _bitsize += 32; } \
- /* _x64 hat höchstens 32 Bits. */\
+ /* _x64 hat höchstens 32 Bits. */\
if (_x64 >= bit(16)) { _x64 = _x64>>16; _bitsize += 16; } \
- /* _x64 hat höchstens 16 Bits. */\
+ /* _x64 hat höchstens 16 Bits. */\
if (_x64 >= bit(8)) { _x64 = _x64>>8; _bitsize += 8; } \
- /* _x64 hat höchstens 8 Bits. */\
+ /* _x64 hat höchstens 8 Bits. */\
if (_x64 >= bit(4)) { _x64 = _x64>>4; _bitsize += 4; } \
- /* _x64 hat höchstens 4 Bits. */\
+ /* _x64 hat höchstens 4 Bits. */\
if (_x64 >= bit(2)) { _x64 = _x64>>2; _bitsize += 2; } \
- /* _x64 hat höchstens 2 Bits. */\
+ /* _x64 hat höchstens 2 Bits. */\
if (_x64 >= bit(1)) { /* _x64 = _x64>>1; */ _bitsize += 1; } \
- /* _x64 hat höchstens 1 Bit. Dieses Bit muß gesetzt sein. */\
+ /* _x64 hat höchstens 1 Bit. Dieses Bit muß gesetzt sein. */\
size_zuweisung _bitsize; \
}
#else
}
#endif
-// Hintere Nullbits eines 32-Bit-Wortes zählen:
+// Bits einer uintC-Zahl zählen:
+// integerlengthC(digit,size=);
+// setzt size auf die höchste in digit vorkommende Bitnummer.
+// > digit: ein uintC >0
+// < size: >0, <=intCsize, mit 2^(size-1) <= digit < 2^size
+ #if (intCsize==32)
+ #define integerlengthC integerlength32
+ #endif
+ #if (intCsize==64)
+ #define integerlengthC integerlength64
+ #endif
+
+// Hintere Nullbits eines 32-Bit-Wortes zählen:
// ord2_32(digit,count=);
// setzt size auf die kleinste in digit vorkommende Bitnummer.
// > digit: ein uint32 >0
}
#endif
-// Hintere Nullbits eines 64-Bit-Wortes zählen:
+// Hintere Nullbits eines 64-Bit-Wortes zählen:
// ord2_64(digit,count=);
// setzt size auf die kleinste in digit vorkommende Bitnummer.
// > digit: ein uint64 >0
}
-// Bits eines Wortes zählen.
+// Bits eines Wortes zählen.
// logcount_NN();
// > xNN: ein uintNN
// < xNN: Anzahl der darin gesetzten Bits
- // Bits von x8 zählen: (Input x8, Output x8)
+ // Bits von x8 zählen: (Input x8, Output x8)
#define logcount_8() \
- ( /* x8 besteht aus 8 1-Bit-Zählern (0,1). */\
+ ( /* x8 besteht aus 8 1-Bit-Zählern (0,1). */\
x8 = (x8 & 0x55U) + ((x8 & 0xAAU) >> 1), \
- /* x8 besteht aus 4 2-Bit-Zählern (0,1,2). */\
+ /* x8 besteht aus 4 2-Bit-Zählern (0,1,2). */\
x8 = (x8 & 0x33U) + ((x8 & 0xCCU) >> 2), \
- /* x8 besteht aus 2 4-Bit-Zählern (0,1,2,3,4). */\
+ /* x8 besteht aus 2 4-Bit-Zählern (0,1,2,3,4). */\
x8 = (x8 & 0x0FU) + (x8 >> 4) \
- /* x8 besteht aus 1 8-Bit-Zähler (0,...,8). */\
+ /* x8 besteht aus 1 8-Bit-Zähler (0,...,8). */\
)
- // Bits von x16 zählen: (Input x16, Output x16)
+ // Bits von x16 zählen: (Input x16, Output x16)
#define logcount_16() \
- ( /* x16 besteht aus 16 1-Bit-Zählern (0,1). */\
+ ( /* x16 besteht aus 16 1-Bit-Zählern (0,1). */\
x16 = (x16 & 0x5555U) + ((x16 & 0xAAAAU) >> 1), \
- /* x16 besteht aus 8 2-Bit-Zählern (0,1,2). */\
+ /* x16 besteht aus 8 2-Bit-Zählern (0,1,2). */\
x16 = (x16 & 0x3333U) + ((x16 & 0xCCCCU) >> 2), \
- /* x16 besteht aus 4 4-Bit-Zählern (0,1,2,3,4). */\
+ /* x16 besteht aus 4 4-Bit-Zählern (0,1,2,3,4). */\
x16 = (x16 & 0x0F0FU) + ((x16 & 0xF0F0U) >> 4), \
- /* x16 besteht aus 2 8-Bit-Zählern (0,...,8). */\
+ /* x16 besteht aus 2 8-Bit-Zählern (0,...,8). */\
x16 = (x16 & 0x00FFU) + (x16 >> 8) \
- /* x16 besteht aus 1 16-Bit-Zähler (0,...,16). */\
+ /* x16 besteht aus 1 16-Bit-Zähler (0,...,16). */\
)
- // Bits von x32 zählen: (Input x32, Output x32)
+ // Bits von x32 zählen: (Input x32, Output x32)
#define logcount_32() \
- ( /* x32 besteht aus 32 1-Bit-Zählern (0,1). */\
+ ( /* x32 besteht aus 32 1-Bit-Zählern (0,1). */\
x32 = (x32 & 0x55555555UL) + ((x32 & 0xAAAAAAAAUL) >> 1), \
- /* x32 besteht aus 16 2-Bit-Zählern (0,1,2). */\
+ /* x32 besteht aus 16 2-Bit-Zählern (0,1,2). */\
x32 = (x32 & 0x33333333UL) + ((x32 & 0xCCCCCCCCUL) >> 2), \
- /* x32 besteht aus 8 4-Bit-Zählern (0,1,2,3,4). */\
+ /* x32 besteht aus 8 4-Bit-Zählern (0,1,2,3,4). */\
x32 = high16(x32)+low16(x32), \
- /* x32 besteht aus 4 4-Bit-Zählern (0,...,8). */\
+ /* x32 besteht aus 4 4-Bit-Zählern (0,...,8). */\
x32 = (x32 & 0x0F0FU) + ((x32 & 0xF0F0U) >> 4), \
- /* x32 besteht aus 2 8-Bit-Zählern (0,...,16). */\
+ /* x32 besteht aus 2 8-Bit-Zählern (0,...,16). */\
x32 = (x32 & 0x00FFU) + (x32 >> 8) \
- /* x32 besteht aus 1 16-Bit-Zähler (0,...,32). */\
+ /* x32 besteht aus 1 16-Bit-Zähler (0,...,32). */\
)
- // Bits von x64 zählen: (Input x64, Output x64)
+ // Bits von x64 zählen: (Input x64, Output x64)
#define logcount_64() \
- ( /* x64 besteht aus 64 1-Bit-Zählern (0,1). */\
+ ( /* x64 besteht aus 64 1-Bit-Zählern (0,1). */\
x64 = (x64 & 0x5555555555555555ULL) + ((x64 & 0xAAAAAAAAAAAAAAAAULL) >> 1),\
- /* x64 besteht aus 32 2-Bit-Zählern (0,1,2). */\
+ /* x64 besteht aus 32 2-Bit-Zählern (0,1,2). */\
x64 = (x64 & 0x3333333333333333ULL) + ((x64 & 0xCCCCCCCCCCCCCCCCULL) >> 2),\
- /* x64 besteht aus 16 4-Bit-Zählern (0,1,2,3,4). */\
+ /* x64 besteht aus 16 4-Bit-Zählern (0,1,2,3,4). */\
x64 = (uint32)(x64 + (x64 >> 32)), \
- /* x64 besteht aus 8 4-Bit-Zählern (0,...,8). */\
+ /* x64 besteht aus 8 4-Bit-Zählern (0,...,8). */\
x64 = (x64 & 0x0F0F0F0FUL) + ((x64 & 0xF0F0F0F0UL) >> 4), \
- /* x64 besteht aus 4 8-Bit-Zählern (0,...,16). */\
+ /* x64 besteht aus 4 8-Bit-Zählern (0,...,16). */\
x64 = (x64 & 0x00FF00FFU) + ((x64 & 0xFF00FF00U) >> 8), \
- /* x64 besteht aus 2 16-Bit-Zählern (0,...,32). */\
+ /* x64 besteht aus 2 16-Bit-Zählern (0,...,32). */\
x64 = (x64 & 0x0000FFFFU) + (x64 >> 16) \
- /* x64 besteht aus 1 16-Bit-Zähler (0,...,64). */\
+ /* x64 besteht aus 1 16-Bit-Zähler (0,...,64). */\
)
} // namespace cln
git://www.ginac.de / cln.git / blobdiff