ecc-mul-g.c 2.87 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
/* ecc-mul-g.c

   Copyright (C) 2013 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/.
*/
Niels Möller's avatar
Niels Möller committed
31

32
/* Development of Nettle's ECC support was funded by the .SE Internet Fund. */
Niels Möller's avatar
Niels Möller committed
33 34 35 36 37 38 39 40 41 42 43 44 45 46 47

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

#include <assert.h>

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

void
ecc_mul_g (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_jja call. Current total is
Niels Möller's avatar
Niels Möller committed
48
     9 * ecc->p.size, at most 648 bytes. */
Niels Möller's avatar
Niels Möller committed
49
#define tp scratch
Niels Möller's avatar
Niels Möller committed
50
#define scratch_out (scratch + 3*ecc->p.size)
Niels Möller's avatar
Niels Möller committed
51 52 53 54 55 56 57 58 59 60

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

  int is_zero;

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

Niels Möller's avatar
Niels Möller committed
61
  bit_rows = (ecc->p.bit_size + k - 1) / k;
Niels Möller's avatar
Niels Möller committed
62
  
Niels Möller's avatar
Niels Möller committed
63
  mpn_zero (r, 3*ecc->p.size);
Niels Möller's avatar
Niels Möller committed
64 65 66 67 68 69 70
  
  for (i = k, is_zero = 1; i-- > 0; )
    {
      ecc_dup_jj (ecc, r, r, scratch);
      for (j = 0; j * c < bit_rows; j++)
	{
	  unsigned bits;
71 72 73
	  /* Avoid the mp_bitcnt_t type for compatibility with older GMP
	     versions. */
	  unsigned bit_index;
Niels Möller's avatar
Niels Möller committed
74 75 76 77 78 79 80 81 82 83 84
	  
	  /* 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
85
	      if (limb_index >= ecc->p.size)
Niels Möller's avatar
Niels Möller committed
86 87 88 89 90
		continue;

	      shift = bit_index % GMP_NUMB_BITS;
	      bits = (bits << 1) | ((np[limb_index] >> shift) & 1);
	    }
Niels Möller's avatar
Niels Möller committed
91
	  sec_tabselect (tp, 2*ecc->p.size,
Niels Möller's avatar
Niels Möller committed
92
			 (ecc->pippenger_table
Niels Möller's avatar
Niels Möller committed
93
			  + (2*ecc->p.size * (mp_size_t) j << c)),
Niels Möller's avatar
Niels Möller committed
94
			 1<<c, bits);
Niels Möller's avatar
Niels Möller committed
95 96
	  cnd_copy (is_zero, r, tp, 2*ecc->p.size);
	  cnd_copy (is_zero, r + 2*ecc->p.size, ecc->unit, ecc->p.size);
Niels Möller's avatar
Niels Möller committed
97 98 99 100
	  
	  ecc_add_jja (ecc, tp, r, tp, scratch_out);
	  /* Use the sum when valid. ecc_add_jja produced garbage if
	     is_zero != 0 or bits == 0, . */	  
Niels Möller's avatar
Niels Möller committed
101
	  cnd_copy (bits & (is_zero - 1), r, tp, 3*ecc->p.size);
Niels Möller's avatar
Niels Möller committed
102 103 104 105 106 107
	  is_zero &= (bits == 0);
	}
    }
#undef tp
#undef scratch_out
}