bignum.c 2.27 KB
Newer Older
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36
/* bignum.c
 *
 * bignum operations that are missing from gmp.
 */

/* nettle, low-level cryptographics library
 *
 * Copyright (C) 2001 Niels Mller
 *  
 * The nettle library is free software; you can redistribute it and/or modify
 * it under the terms of the GNU Lesser General Public License as published by
 * the Free Software Foundation; either version 2.1 of the License, or (at your
 * option) any later version.
 * 
 * The nettle library is distributed in the hope that it will be useful, but
 * WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
 * or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public
 * License for more details.
 * 
 * You should have received a copy of the GNU Lesser General Public License
 * along with the nettle library; see the file COPYING.LIB.  If not, write to
 * the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
 * MA 02111-1307, USA.
 */

#if HAVE_CONFIG_H
#include "config.h"
#endif

#if HAVE_LIBGMP

#include "bignum.h"

#include <assert.h>
#include <string.h>

37 38 39 40 41 42
unsigned
nettle_mpz_sizeinbase_256(const mpz_t x)
{
  return (mpz_sizeinbase(x, 2) + 7) / 8;
}

43
void
44
nettle_mpz_get_str_256(unsigned length, uint8_t *s, const mpz_t x)
45 46 47 48 49 50 51 52 53 54 55 56 57
{
  uint8_t *dst = s + length - 1;
  unsigned size = mpz_size(x);
  unsigned i;

  if (!length)
    {
      /* x must be zero */
      assert(!mpz_sgn(x));
      return;
    }
  
  assert(mpz_sgn(x) >= 0);
58
  assert(nettle_mpz_sizeinbase_256(x) <= length);
59 60 61 62 63 64 65 66 67 68 69 70 71 72 73

  for (i = 0; i<size; i++)
    {
      mp_limb_t limb = mpz_getlimbn(x, i);
      unsigned j;

      for (j = 0; length && j < sizeof(mp_limb_t); j++)
        {
          *dst-- = limb & 0xff;
          limb >>= 8;
          length--;
	}
    }
  
  if (length)
74
    memset(s, 0, length);
75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103
}

void
nettle_mpz_set_str_256(mpz_t x,
                       unsigned length, const uint8_t *s)
{
  size_t i;
  mpz_t digit;

  mpz_init(digit);
  mpz_set_ui(x, 0);
  for (i = 0; i < length; i++)
    {
      mpz_set_ui(digit, s[i]);
      mpz_mul_2exp(digit, digit, (length - i - 1) * 8);
      mpz_ior(x, x, digit);
    }
  mpz_clear(digit);
}

void
nettle_mpz_init_set_str_256(mpz_t x,
                            unsigned length, const uint8_t *s)
{
  mpz_init(x);
  nettle_mpz_set_str_256(x, length, s);
}

#endif /* HAVE_LIBGMP */