ecc-mul-g-eh.c 2.78 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 37 38 39 40 41 42 43 44 45
/* ecc-mul-g-eh.c

   Copyright (C) 2013, 2014 Niels Möller

   This file is part of GNU Nettle.

   GNU Nettle is free software: you can redistribute it and/or
   modify it under the terms of either:

     * the GNU Lesser General Public License as published by the Free
       Software Foundation; either version 3 of the License, or (at your
       option) any later version.

   or

     * the GNU General Public License as published by the Free
       Software Foundation; either version 2 of the License, or (at your
       option) any later version.

   or both in parallel, as here.

   GNU Nettle 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
   General Public License for more details.

   You should have received copies of the GNU General Public License and
   the GNU Lesser General Public License along with this program.  If
   not, see http://www.gnu.org/licenses/.
*/

/* Development of Nettle's ECC support was funded by the .SE Internet Fund. */

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

#include <assert.h>

#include "ecc.h"
#include "ecc-internal.h"

mp_size_t
ecc_mul_g_eh_itch (const struct ecc_curve *ecc)
{
Niels Möller's avatar
Niels Möller committed
46 47
  /* Needs 3*ecc->p.size + scratch for ecc_add_jja. */
  return ECC_MUL_G_EH_ITCH (ecc->p.size);
48 49 50 51 52 53 54
}

void
ecc_mul_g_eh (const struct ecc_curve *ecc, mp_limb_t *r,
	      const mp_limb_t *np, mp_limb_t *scratch)
{
  /* Scratch need determined by the ecc_add_eh call. Current total is
Niels Möller's avatar
Niels Möller committed
55
     9 * ecc->p.size, at most 648 bytes. */
56
#define tp scratch
Niels Möller's avatar
Niels Möller committed
57
#define scratch_out (scratch + 3*ecc->p.size)
58 59 60 61 62 63 64 65

  unsigned k, c;
  unsigned i, j;
  unsigned bit_rows;

  k = ecc->pippenger_k;
  c = ecc->pippenger_c;

Niels Möller's avatar
Niels Möller committed
66
  bit_rows = (ecc->p.bit_size + k - 1) / k;
67 68

  /* x = 0, y = 1, z = 1 */
Niels Möller's avatar
Niels Möller committed
69 70
  mpn_zero (r, 3*ecc->p.size);
  r[ecc->p.size] = r[2*ecc->p.size] = 1;
71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91

  for (i = k; i-- > 0; )
    {
      ecc_dup_eh (ecc, r, r, scratch);
      for (j = 0; j * c < bit_rows; j++)
	{
	  unsigned bits;
	  /* Avoid the mp_bitcnt_t type for compatibility with older GMP
	     versions. */
	  unsigned bit_index;
	  
	  /* Extract c bits from n, stride k, starting at i + kcj,
	     ending at i + k (cj + c - 1)*/
	  for (bits = 0, bit_index = i + k*(c*j+c); bit_index > i + k*c*j; )
	    {
	      mp_size_t limb_index;
	      unsigned shift;
	      
	      bit_index -= k;

	      limb_index = bit_index / GMP_NUMB_BITS;
Niels Möller's avatar
Niels Möller committed
92
	      if (limb_index >= ecc->p.size)
93 94 95 96 97
		continue;

	      shift = bit_index % GMP_NUMB_BITS;
	      bits = (bits << 1) | ((np[limb_index] >> shift) & 1);
	    }
Niels Möller's avatar
Niels Möller committed
98
	  sec_tabselect (tp, 2*ecc->p.size,
99
			 (ecc->pippenger_table
Niels Möller's avatar
Niels Möller committed
100
			  + (2*ecc->p.size * (mp_size_t) j << c)),
101 102 103 104 105 106 107 108
			 1<<c, bits);

	  ecc_add_eh (ecc, r, r, tp, scratch_out);
	}
    }
#undef tp
#undef scratch_out
}