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 int32_t write_mpbin(MP_INT *a, uint8_t *buffer,
  65                                      int32_t buflen)
  66 {
  67     char *hb = write_mpstring(a);
  68     int32_t len;
  69     hex_decode(buffer, buflen, &len, hb, True);
  70     free(hb);
  71     return len;
  72 }
  73
  74 static dh_makeshared_fn dh_makeshared;
  75 static void dh_makeshared(void *sst, uint8_t *secret, int32_t secretlen,
  76                           cstring_t rempublic, uint8_t *sharedsecret,
  77                           int32_t buflen)
  78 {
  79     struct dh *st=sst;
  80     MP_INT a, b, c;
  81
  82     mpz_init(&a);
  83     mpz_init(&b);
  84     mpz_init(&c);
  85
  86     read_mpbin(&a, secret, secretlen);
  87     mpz_set_str(&b, rempublic, 16);
  88
  89     mpz_powm_sec(&c, &b, &a, &st->p);
  90
  91     write_mpbin(&c,sharedsecret,buflen);
  92
  93     mpz_clear(&a);
  94     mpz_clear(&b);
  95     mpz_clear(&c);
  96 }
  97
  98 static list_t *dh_apply(closure_t *self, struct cloc loc, dict_t *context,
  99                         list_t *args)
 100 {
 101     struct dh *st;
 102     string_t p,g;
 103     item_t *i;
 104
 105     NEW(st);
 106     st->cl.description="dh";
 107     st->cl.type=CL_DH;
 108     st->cl.apply=NULL;
 109     st->cl.interface=&st->ops;
 110     st->ops.st=st;
 111     st->ops.makepublic=dh_makepublic;
 112     st->ops.makeshared=dh_makeshared;
 113     st->loc=loc;
 114     /* We have two string arguments: the first is the modulus, and the
 115        second is the generator. Both are in hex. */
 116     i=list_elem(args,0);
 117     if (i) {
 118         if (i->type!=t_string) {
 119             cfgfatal(i->loc,"diffie-hellman","first argument must be a "
 120                      "string\n");
 121         }
 122         p=i->data.string;
 123         if (mpz_init_set_str(&st->p,p,16)!=0) {
 124             cfgfatal(i->loc,"diffie-hellman","\"%s\" is not a hex number "
 125                      "string\n",p);
 126         }
 127     } else {
 128         cfgfatal(loc,"diffie-hellman","you must provide a prime modulus\n");
 129     }
 130
 131     i=list_elem(args,1);
 132     if (i) {
 133         if (i->type!=t_string) {
 134             cfgfatal(i->loc,"diffie-hellman","second argument must be a "
 135                      "string\n");
 136         }
 137         g=i->data.string;
 138         if (mpz_init_set_str(&st->g,g,16)!=0) {
 139             cfgfatal(i->loc,"diffie-hellman","\"%s\" is not a hex number "
 140                      "string\n",g);
 141         }
 142     } else {
 143         cfgfatal(loc,"diffie-hellman","you must provide a generator\n");
 144     }
 145
 146     i=list_elem(args,2);
 147     if (i && i->type==t_bool && i->data.bool==False) {
 148         Message(M_INFO,"diffie-hellman (%s:%d): skipping modulus "
 149                 "primality check\n",loc.file,loc.line);
 150     } else {
 151         /* Test that the modulus is really prime */
 152         if (mpz_probab_prime_p(&st->p,5)==0) {
 153             cfgfatal(loc,"diffie-hellman","modulus must be a prime\n");
 154         }
 155     }
 156
 157     size_t sz=mpz_sizeinbase(&st->p,2)/8;
 158     if (sz>INT_MAX) {
 159         cfgfatal(loc,"diffie-hellman","modulus far too large\n");
 160     }
 161     if (mpz_cmp(&st->g,&st->p) >= 0) {
 162         cfgfatal(loc,"diffie-hellman","generator must be less than modulus\n");
 163     }
 164
 165     st->ops.len=sz;
 166
 167     st->ops.ceil_len=(mpz_sizeinbase(&st->p,2)+7)/8;
 168     /* According to the docs, mpz_sizeinbase(,256) is allowed to return
 169      * an answer which is 1 too large.  But mpz_sizeinbase(,2) isn't. */
 170
 171     return new_closure(&st->cl);
 172 }
 173
 174 void dh_module(dict_t *dict)
 175 {
 176     add_closure(dict,"diffie-hellman",dh_apply);
 177 }