dh.c

   1 /*
   2  * dh.c
   3  */
   4 /*
   5  * This file is Free Software.  It was originally written for secnet.
   6  *
   7  * Copyright 1995-2003 Stephen Early
   8  * Copyright 2002-2014 Ian Jackson
   9  *
  10  * You may redistribute secnet as a whole and/or modify it under the
  11  * terms of the GNU General Public License as published by the Free
  12  * Software Foundation; either version 3, or (at your option) any
  13  * later version.
  14  *
  15  * You may redistribute this file and/or modify it under the terms of
  16  * the GNU General Public License as published by the Free Software
  17  * Foundation; either version 2, or (at your option) any later
  18  * version.
  19  *
  20  * This software is distributed in the hope that it will be useful,
  21  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  22  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  23  * GNU General Public License for more details.
  24  *
  25  * You should have received a copy of the GNU General Public License
  26  * along with this software; if not, see
  27  * https://www.gnu.org/licenses/gpl.html.
  28  */
  29
  30 #include <stdio.h>
  31 #include <gmp.h>
  32 #include <limits.h>
  33
  34 #include "secnet.h"
  35 #include "util.h"
  36
  37 struct dh {
  38     closure_t cl;
  39     struct dh_if ops;
  40     struct cloc loc;
  41     MP_INT p,g; /* prime modulus and generator */
  42 };
  43
  44 static string_t dh_makepublic(void *sst, uint8_t *secret, int32_t secretlen)
  45 {
  46     struct dh *st=sst;
  47     string_t r;
  48     MP_INT a, b; /* a is secret key, b is public key */
  49
  50     mpz_init(&a);
  51     mpz_init(&b);
  52
  53     read_mpbin(&a, secret, secretlen);
  54
  55     mpz_powm_sec(&b, &st->g, &a, &st->p);
  56
  57     r=write_mpstring(&b);
  58
  59     mpz_clear(&a);
  60     mpz_clear(&b);
  61     return r;
  62 }
  63
  64 static void write_mpbin_anomalous(MP_INT *a, uint8_t *buffer,
  65                                   int32_t buflen)
  66     /* If the BN is smaller than buflen, pads it *at the wrong end* */
  67 {
  68     char *hb = write_mpstring(a);
  69     int32_t len;
  70     hex_decode(buffer, buflen, &len, hb, True);
  71     if (len<buflen)
  72         memset(buffer+len,0,buflen-len);
  73     free(hb);
  74 }
  75
  76 static dh_makeshared_fn dh_makeshared;
  77 static void dh_makeshared(void *sst, uint8_t *secret, int32_t secretlen,
  78                           cstring_t rempublic, uint8_t *sharedsecret,
  79                           int32_t buflen)
  80 {
  81     struct dh *st=sst;
  82     MP_INT a, b, c;
  83
  84     mpz_init(&a);
  85     mpz_init(&b);
  86     mpz_init(&c);
  87
  88     read_mpbin(&a, secret, secretlen);
  89     mpz_set_str(&b, rempublic, 16);
  90
  91     mpz_powm_sec(&c, &b, &a, &st->p);
  92
  93     write_mpbin_anomalous(&c,sharedsecret,buflen);
  94
  95     mpz_clear(&a);
  96     mpz_clear(&b);
  97     mpz_clear(&c);
  98 }
  99
 100 static list_t *dh_apply(closure_t *self, struct cloc loc, dict_t *context,
 101                         list_t *args)
 102 {
 103     struct dh *st;
 104     string_t p,g;
 105     item_t *i;
 106
 107     NEW(st);
 108     st->cl.description="dh";
 109     st->cl.type=CL_DH;
 110     st->cl.apply=NULL;
 111     st->cl.interface=&st->ops;
 112     st->ops.st=st;
 113     st->ops.makepublic=dh_makepublic;
 114     st->ops.makeshared=dh_makeshared;
 115     st->loc=loc;
 116     /* We have two string arguments: the first is the modulus, and the
 117        second is the generator. Both are in hex. */
 118     i=list_elem(args,0);
 119     if (i) {
 120         if (i->type!=t_string) {
 121             cfgfatal(i->loc,"diffie-hellman","first argument must be a "
 122                      "string\n");
 123         }
 124         p=i->data.string;
 125         if (mpz_init_set_str(&st->p,p,16)!=0) {
 126             cfgfatal(i->loc,"diffie-hellman","\"%s\" is not a hex number "
 127                      "string\n",p);
 128         }
 129     } else {
 130         cfgfatal(loc,"diffie-hellman","you must provide a prime modulus\n");
 131     }
 132
 133     i=list_elem(args,1);
 134     if (i) {
 135         if (i->type!=t_string) {
 136             cfgfatal(i->loc,"diffie-hellman","second argument must be a "
 137                      "string\n");
 138         }
 139         g=i->data.string;
 140         if (mpz_init_set_str(&st->g,g,16)!=0) {
 141             cfgfatal(i->loc,"diffie-hellman","\"%s\" is not a hex number "
 142                      "string\n",g);
 143         }
 144     } else {
 145         cfgfatal(loc,"diffie-hellman","you must provide a generator\n");
 146     }
 147
 148     i=list_elem(args,2);
 149     if (i && i->type==t_bool && i->data.bool==False) {
 150         Message(M_INFO,"diffie-hellman (%s:%d): skipping modulus "
 151                 "primality check\n",loc.file,loc.line);
 152     } else {
 153         /* Test that the modulus is really prime */
 154         if (mpz_probab_prime_p(&st->p,5)==0) {
 155             cfgfatal(loc,"diffie-hellman","modulus must be a prime\n");
 156         }
 157     }
 158
 159     size_t sz=mpz_sizeinbase(&st->p,2)/8;
 160     if (sz>INT_MAX) {
 161         cfgfatal(loc,"diffie-hellman","modulus far too large\n");
 162     }
 163     if (mpz_cmp(&st->g,&st->p) >= 0) {
 164         cfgfatal(loc,"diffie-hellman","generator must be less than modulus\n");
 165     }
 166
 167     st->ops.len=sz;
 168
 169     st->ops.ceil_len=(mpz_sizeinbase(&st->p,2)+7)/8;
 170     /* According to the docs, mpz_sizeinbase(,256) is allowed to return
 171      * an answer which is 1 too large.  But mpz_sizeinbase(,2) isn't. */
 172
 173     return new_closure(&st->cl);
 174 }
 175
 176 void dh_module(dict_t *dict)
 177 {
 178     add_closure(dict,"diffie-hellman",dh_apply);
 179 }