Blame SOURCES/libgcrypt-1.8.5-elgamal-blinding.patch

9b59ab
commit e8b7f10be275bcedb5fc05ed4837a89bfd605c61
9b59ab
Author: NIIBE Yutaka <gniibe@fsij.org>
9b59ab
Date:   Tue Apr 13 10:00:00 2021 +0900
9b59ab
9b59ab
    cipher: Hardening ElGamal by introducing exponent blinding too.
9b59ab
    
9b59ab
    * cipher/elgamal.c (do_encrypt): Also do exponent blinding.
9b59ab
    
9b59ab
    --
9b59ab
    
9b59ab
    Base blinding had been introduced with USE_BLINDING.  This patch add
9b59ab
    exponent blinding as well to mitigate side-channel attack on mpi_powm.
9b59ab
    
9b59ab
    GnuPG-bug-id: 5328
9b59ab
    Signed-off-by: NIIBE Yutaka <gniibe@fsij.org>
9b59ab
9b59ab
diff --git a/cipher/elgamal.c b/cipher/elgamal.c
9b59ab
index 4eb52d62..9835122f 100644
9b59ab
--- a/cipher/elgamal.c
9b59ab
+++ b/cipher/elgamal.c
9b59ab
@@ -522,8 +522,9 @@ do_encrypt(gcry_mpi_t a, gcry_mpi_t b, gcry_mpi_t input, ELG_public_key *pkey )
9b59ab
 static void
9b59ab
 decrypt (gcry_mpi_t output, gcry_mpi_t a, gcry_mpi_t b, ELG_secret_key *skey )
9b59ab
 {
9b59ab
-  gcry_mpi_t t1, t2, r;
9b59ab
+  gcry_mpi_t t1, t2, r, r1, h;
9b59ab
   unsigned int nbits = mpi_get_nbits (skey->p);
9b59ab
+  gcry_mpi_t x_blind;
9b59ab
 
9b59ab
   mpi_normalize (a);
9b59ab
   mpi_normalize (b);
9b59ab
@@ -534,20 +535,33 @@ decrypt (gcry_mpi_t output, gcry_mpi_t a, gcry_mpi_t b, ELG_secret_key *skey )
9b59ab
 
9b59ab
   t2 = mpi_snew (nbits);
9b59ab
   r  = mpi_new (nbits);
9b59ab
+  r1 = mpi_new (nbits);
9b59ab
+  h  = mpi_new (nbits);
9b59ab
+  x_blind = mpi_snew (nbits);
9b59ab
 
9b59ab
   /* We need a random number of about the prime size.  The random
9b59ab
      number merely needs to be unpredictable; thus we use level 0.  */
9b59ab
   _gcry_mpi_randomize (r, nbits, GCRY_WEAK_RANDOM);
9b59ab
 
9b59ab
+  /* Also, exponent blinding: x_blind = x + (p-1)*r1 */
9b59ab
+  _gcry_mpi_randomize (r1, nbits, GCRY_WEAK_RANDOM);
9b59ab
+  mpi_set_highbit (r1, nbits - 1);
9b59ab
+  mpi_sub_ui (h, skey->p, 1);
9b59ab
+  mpi_mul (x_blind, h, r1);
9b59ab
+  mpi_add (x_blind, skey->x, x_blind);
9b59ab
+
9b59ab
   /* t1 = r^x mod p */
9b59ab
-  mpi_powm (t1, r, skey->x, skey->p);
9b59ab
+  mpi_powm (t1, r, x_blind, skey->p);
9b59ab
   /* t2 = (a * r)^-x mod p */
9b59ab
   mpi_mulm (t2, a, r, skey->p);
9b59ab
-  mpi_powm (t2, t2, skey->x, skey->p);
9b59ab
+  mpi_powm (t2, t2, x_blind, skey->p);
9b59ab
   mpi_invm (t2, t2, skey->p);
9b59ab
   /* t1 = (t1 * t2) mod p*/
9b59ab
   mpi_mulm (t1, t1, t2, skey->p);
9b59ab
 
9b59ab
+  mpi_free (x_blind);
9b59ab
+  mpi_free (h);
9b59ab
+  mpi_free (r1);
9b59ab
   mpi_free (r);
9b59ab
   mpi_free (t2);
9b59ab