#include <stdio.h>
#include "gmp.h"
#include "gmp-impl.h"
#if HAVE_LIMB_BIG_ENDIAN
#define HOST_ENDIAN 1
#endif
#if HAVE_LIMB_LITTLE_ENDIAN
#define HOST_ENDIAN (-1)
#endif
#ifndef HOST_ENDIAN
static const mp_limb_t endian_test = (CNST_LIMB(1) << (GMP_LIMB_BITS-7)) - 1;
#define HOST_ENDIAN (* (signed char *) &endian_test)
#endif
void
mpz_import (mpz_ptr z, size_t count, int order,
size_t size, int endian, size_t nail, const void *data)
{
mp_size_t zsize;
mp_ptr zp;
ASSERT (order == 1 || order == -1);
ASSERT (endian == 1 || endian == 0 || endian == -1);
ASSERT (nail <= 8*size);
zsize = BITS_TO_LIMBS (count * (8*size - nail));
zp = MPZ_NEWALLOC (z, zsize);
if (endian == 0)
endian = HOST_ENDIAN;
if (nail == 0 && GMP_NAIL_BITS == 0)
{
unsigned align = ((char *) data - (char *) NULL) % sizeof (mp_limb_t);
if (order == -1
&& size == sizeof (mp_limb_t)
&& endian == HOST_ENDIAN
&& align == 0)
{
MPN_COPY (zp, (mp_srcptr) data, (mp_size_t) count);
goto done;
}
if (order == -1
&& size == sizeof (mp_limb_t)
&& endian == - HOST_ENDIAN
&& align == 0)
{
MPN_BSWAP (zp, (mp_srcptr) data, (mp_size_t) count);
goto done;
}
if (order == 1
&& size == sizeof (mp_limb_t)
&& endian == HOST_ENDIAN
&& align == 0)
{
MPN_REVERSE (zp, (mp_srcptr) data, (mp_size_t) count);
goto done;
}
}
{
mp_limb_t limb, byte, wbitsmask;
size_t i, j, numb, wbytes;
mp_size_t woffset;
unsigned char *dp;
int lbits, wbits;
numb = size * 8 - nail;
wbytes = numb / 8;
wbits = numb % 8;
wbitsmask = (CNST_LIMB(1) << wbits) - 1;
woffset = (numb + 7) / 8;
woffset = (endian >= 0 ? woffset : -woffset)
+ (order < 0 ? size : - (mp_size_t) size);
dp = (unsigned char *) data
+ (order >= 0 ? (count-1)*size : 0) + (endian >= 0 ? size-1 : 0);
#define ACCUMULATE(N) \
do { \
ASSERT (lbits < GMP_NUMB_BITS); \
ASSERT (limb <= (CNST_LIMB(1) << lbits) - 1); \
\
limb |= (mp_limb_t) byte << lbits; \
lbits += (N); \
if (lbits >= GMP_NUMB_BITS) \
{ \
*zp++ = limb & GMP_NUMB_MASK; \
lbits -= GMP_NUMB_BITS; \
ASSERT (lbits < (N)); \
limb = byte >> ((N) - lbits); \
} \
} while (0)
limb = 0;
lbits = 0;
for (i = 0; i < count; i++)
{
for (j = 0; j < wbytes; j++)
{
byte = *dp;
dp -= endian;
ACCUMULATE (8);
}
if (wbits != 0)
{
byte = *dp & wbitsmask;
dp -= endian;
ACCUMULATE (wbits);
}
dp += woffset;
}
if (lbits != 0)
{
ASSERT (lbits <= GMP_NUMB_BITS);
ASSERT_LIMB (limb);
*zp++ = limb;
}
ASSERT (zp == PTR(z) + zsize);
ASSERT (dp == (unsigned char *) data
+ (order < 0 ? count*size : - (mp_size_t) size)
+ (endian >= 0 ? (mp_size_t) size - 1 : 0));
}
done:
zp = PTR(z);
MPN_NORMALIZE (zp, zsize);
SIZ(z) = zsize;
}