|
Packit |
5c3484 |
/* mpz_out_raw -- write an mpz_t in raw format.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
Copyright 2001, 2002 Free Software Foundation, Inc.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
This file is part of the GNU MP Library.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
The GNU MP Library is free software; you can redistribute it and/or modify
|
|
Packit |
5c3484 |
it under the terms of either:
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
* the GNU Lesser General Public License as published by the Free
|
|
Packit |
5c3484 |
Software Foundation; either version 3 of the License, or (at your
|
|
Packit |
5c3484 |
option) any later version.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
or
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
* the GNU General Public License as published by the Free Software
|
|
Packit |
5c3484 |
Foundation; either version 2 of the License, or (at your option) any
|
|
Packit |
5c3484 |
later version.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
or both in parallel, as here.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
The GNU MP Library is distributed in the hope that it will be useful, but
|
|
Packit |
5c3484 |
WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
|
|
Packit |
5c3484 |
or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
|
Packit |
5c3484 |
for more details.
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
You should have received copies of the GNU General Public License and the
|
|
Packit |
5c3484 |
GNU Lesser General Public License along with the GNU MP Library. If not,
|
|
Packit |
5c3484 |
see https://www.gnu.org/licenses/. */
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
#include <stdio.h>
|
|
Packit |
5c3484 |
#include "gmp.h"
|
|
Packit |
5c3484 |
#include "gmp-impl.h"
|
|
Packit |
5c3484 |
#include "longlong.h"
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
/* HTON_LIMB_STORE takes a normal host byte order limb and stores it as
|
|
Packit |
5c3484 |
network byte order (ie. big endian). */
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
#if HAVE_LIMB_BIG_ENDIAN
|
|
Packit |
5c3484 |
#define HTON_LIMB_STORE(dst, limb) do { *(dst) = (limb); } while (0)
|
|
Packit |
5c3484 |
#endif
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
#if HAVE_LIMB_LITTLE_ENDIAN
|
|
Packit |
5c3484 |
#define HTON_LIMB_STORE(dst, limb) BSWAP_LIMB_STORE (dst, limb)
|
|
Packit |
5c3484 |
#endif
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
#ifndef HTON_LIMB_STORE
|
|
Packit |
5c3484 |
#define HTON_LIMB_STORE(dst, limb) \
|
|
Packit |
5c3484 |
do { \
|
|
Packit |
5c3484 |
mp_limb_t __limb = (limb); \
|
|
Packit |
5c3484 |
char *__p = (char *) (dst); \
|
|
Packit |
5c3484 |
int __i; \
|
|
Packit |
5c3484 |
for (__i = 0; __i < GMP_LIMB_BYTES; __i++) \
|
|
Packit |
5c3484 |
__p[__i] = (char) (__limb >> ((GMP_LIMB_BYTES-1 - __i) * 8)); \
|
|
Packit |
5c3484 |
} while (0)
|
|
Packit |
5c3484 |
#endif
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
size_t
|
|
Packit |
5c3484 |
mpz_out_raw (FILE *fp, mpz_srcptr x)
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
mp_size_t xsize, abs_xsize, bytes, i;
|
|
Packit |
5c3484 |
mp_srcptr xp;
|
|
Packit |
5c3484 |
char *tp, *bp;
|
|
Packit |
5c3484 |
mp_limb_t xlimb;
|
|
Packit |
5c3484 |
int zeros;
|
|
Packit |
5c3484 |
size_t tsize, ssize;
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
xsize = SIZ(x);
|
|
Packit |
5c3484 |
abs_xsize = ABS (xsize);
|
|
Packit |
5c3484 |
bytes = (abs_xsize * GMP_NUMB_BITS + 7) / 8;
|
|
Packit |
5c3484 |
tsize = ROUND_UP_MULTIPLE ((unsigned) 4, GMP_LIMB_BYTES) + bytes;
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
tp = __GMP_ALLOCATE_FUNC_TYPE (tsize, char);
|
|
Packit |
5c3484 |
bp = tp + ROUND_UP_MULTIPLE ((unsigned) 4, GMP_LIMB_BYTES);
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
if (bytes != 0)
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
bp += bytes;
|
|
Packit |
5c3484 |
xp = PTR (x);
|
|
Packit |
5c3484 |
i = abs_xsize;
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
if (GMP_NAIL_BITS == 0)
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
/* reverse limb order, and byte swap if necessary */
|
|
Packit |
5c3484 |
#ifdef _CRAY
|
|
Packit |
5c3484 |
_Pragma ("_CRI ivdep");
|
|
Packit |
5c3484 |
#endif
|
|
Packit |
5c3484 |
do
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
bp -= GMP_LIMB_BYTES;
|
|
Packit |
5c3484 |
xlimb = *xp;
|
|
Packit |
5c3484 |
HTON_LIMB_STORE ((mp_ptr) bp, xlimb);
|
|
Packit |
5c3484 |
xp++;
|
|
Packit |
5c3484 |
}
|
|
Packit |
5c3484 |
while (--i > 0);
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
/* strip high zero bytes (without fetching from bp) */
|
|
Packit |
5c3484 |
count_leading_zeros (zeros, xlimb);
|
|
Packit |
5c3484 |
zeros /= 8;
|
|
Packit |
5c3484 |
bp += zeros;
|
|
Packit |
5c3484 |
bytes -= zeros;
|
|
Packit |
5c3484 |
}
|
|
Packit |
5c3484 |
else
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
mp_limb_t new_xlimb;
|
|
Packit |
5c3484 |
int bits;
|
|
Packit |
5c3484 |
ASSERT_CODE (char *bp_orig = bp - bytes);
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
ASSERT_ALWAYS (GMP_NUMB_BITS >= 8);
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
bits = 0;
|
|
Packit |
5c3484 |
xlimb = 0;
|
|
Packit |
5c3484 |
for (;;)
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
while (bits >= 8)
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
ASSERT (bp > bp_orig);
|
|
Packit |
5c3484 |
*--bp = xlimb & 0xFF;
|
|
Packit |
5c3484 |
xlimb >>= 8;
|
|
Packit |
5c3484 |
bits -= 8;
|
|
Packit |
5c3484 |
}
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
if (i == 0)
|
|
Packit |
5c3484 |
break;
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
new_xlimb = *xp++;
|
|
Packit |
5c3484 |
i--;
|
|
Packit |
5c3484 |
ASSERT (bp > bp_orig);
|
|
Packit |
5c3484 |
*--bp = (xlimb | (new_xlimb << bits)) & 0xFF;
|
|
Packit |
5c3484 |
xlimb = new_xlimb >> (8 - bits);
|
|
Packit |
5c3484 |
bits += GMP_NUMB_BITS - 8;
|
|
Packit |
5c3484 |
}
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
if (bits != 0)
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
ASSERT (bp > bp_orig);
|
|
Packit |
5c3484 |
*--bp = xlimb;
|
|
Packit |
5c3484 |
}
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
ASSERT (bp == bp_orig);
|
|
Packit |
5c3484 |
while (*bp == 0)
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
bp++;
|
|
Packit |
5c3484 |
bytes--;
|
|
Packit |
5c3484 |
}
|
|
Packit |
5c3484 |
}
|
|
Packit |
5c3484 |
}
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
/* total bytes to be written */
|
|
Packit |
5c3484 |
ssize = 4 + bytes;
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
/* twos complement negative for the size value */
|
|
Packit |
5c3484 |
bytes = (xsize >= 0 ? bytes : -bytes);
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
/* so we don't rely on sign extension in ">>" */
|
|
Packit |
5c3484 |
ASSERT_ALWAYS (sizeof (bytes) >= 4);
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
bp[-4] = bytes >> 24;
|
|
Packit |
5c3484 |
bp[-3] = bytes >> 16;
|
|
Packit |
5c3484 |
bp[-2] = bytes >> 8;
|
|
Packit |
5c3484 |
bp[-1] = bytes;
|
|
Packit |
5c3484 |
bp -= 4;
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
if (fp == 0)
|
|
Packit |
5c3484 |
fp = stdout;
|
|
Packit |
5c3484 |
if (fwrite (bp, ssize, 1, fp) != 1)
|
|
Packit |
5c3484 |
ssize = 0;
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
(*__gmp_free_func) (tp, tsize);
|
|
Packit |
5c3484 |
return ssize;
|
|
Packit |
5c3484 |
}
|