[BACK]Return to dh.c CVS log [TXT][DIR] Up to [local] / src / usr.bin / ssh

Annotation of src/usr.bin/ssh/dh.c, Revision 1.40

1.40    ! stevesk     1: /* $OpenBSD: dh.c,v 1.39 2006/07/26 02:35:17 stevesk Exp $ */
1.1       provos      2: /*
                      3:  * Copyright (c) 2000 Niels Provos.  All rights reserved.
                      4:  *
                      5:  * Redistribution and use in source and binary forms, with or without
                      6:  * modification, are permitted provided that the following conditions
                      7:  * are met:
                      8:  * 1. Redistributions of source code must retain the above copyright
                      9:  *    notice, this list of conditions and the following disclaimer.
                     10:  * 2. Redistributions in binary form must reproduce the above copyright
                     11:  *    notice, this list of conditions and the following disclaimer in the
                     12:  *    documentation and/or other materials provided with the distribution.
                     13:  *
                     14:  * THIS SOFTWARE IS PROVIDED BY THE AUTHOR ``AS IS'' AND ANY EXPRESS OR
                     15:  * IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES
                     16:  * OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED.
                     17:  * IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR ANY DIRECT, INDIRECT,
                     18:  * INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT
                     19:  * NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
                     20:  * DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
                     21:  * THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
                     22:  * (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF
                     23:  * THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
                     24:  */
                     25:
                     26: #include "includes.h"
1.39      stevesk    27:
                     28: #include <sys/param.h>
1.1       provos     29:
                     30: #include <openssl/bn.h>
                     31: #include <openssl/dh.h>
1.38      stevesk    32:
1.40    ! stevesk    33: #include <stdlib.h>
1.38      stevesk    34: #include <string.h>
1.1       provos     35:
                     36: #include "dh.h"
1.5       markus     37: #include "pathnames.h"
1.6       markus     38: #include "log.h"
                     39: #include "misc.h"
1.1       provos     40:
1.17      itojun     41: static int
1.1       provos     42: parse_prime(int linenum, char *line, struct dhgroup *dhg)
                     43: {
                     44:        char *cp, *arg;
                     45:        char *strsize, *gen, *prime;
1.35      deraadt    46:        const char *errstr = NULL;
1.1       provos     47:
                     48:        cp = line;
1.32      djm        49:        if ((arg = strdelim(&cp)) == NULL)
                     50:                return 0;
1.1       provos     51:        /* Ignore leading whitespace */
                     52:        if (*arg == '\0')
                     53:                arg = strdelim(&cp);
1.22      markus     54:        if (!arg || !*arg || *arg == '#')
1.1       provos     55:                return 0;
                     56:
                     57:        /* time */
                     58:        if (cp == NULL || *arg == '\0')
                     59:                goto fail;
                     60:        arg = strsep(&cp, " "); /* type */
                     61:        if (cp == NULL || *arg == '\0')
                     62:                goto fail;
                     63:        arg = strsep(&cp, " "); /* tests */
                     64:        if (cp == NULL || *arg == '\0')
                     65:                goto fail;
                     66:        arg = strsep(&cp, " "); /* tries */
                     67:        if (cp == NULL || *arg == '\0')
                     68:                goto fail;
                     69:        strsize = strsep(&cp, " "); /* size */
                     70:        if (cp == NULL || *strsize == '\0' ||
1.35      deraadt    71:            (dhg->size = (u_int)strtonum(strsize, 0, 64*1024, &errstr)) == 0 ||
                     72:            errstr)
1.1       provos     73:                goto fail;
1.9       provos     74:        /* The whole group is one bit larger */
                     75:        dhg->size++;
1.1       provos     76:        gen = strsep(&cp, " "); /* gen */
                     77:        if (cp == NULL || *gen == '\0')
                     78:                goto fail;
                     79:        prime = strsep(&cp, " "); /* prime */
                     80:        if (cp != NULL || *prime == '\0')
                     81:                goto fail;
                     82:
1.18      markus     83:        if ((dhg->g = BN_new()) == NULL)
                     84:                fatal("parse_prime: BN_new failed");
                     85:        if ((dhg->p = BN_new()) == NULL)
                     86:                fatal("parse_prime: BN_new failed");
1.14      markus     87:        if (BN_hex2bn(&dhg->g, gen) == 0)
1.10      provos     88:                goto failclean;
                     89:
1.14      markus     90:        if (BN_hex2bn(&dhg->p, prime) == 0)
1.10      provos     91:                goto failclean;
                     92:
                     93:        if (BN_num_bits(dhg->p) != dhg->size)
1.27      dtucker    94:                goto failclean;
                     95:
                     96:        if (BN_is_zero(dhg->g) || BN_is_one(dhg->g))
1.10      provos     97:                goto failclean;
1.1       provos     98:
                     99:        return (1);
1.10      provos    100:
                    101:  failclean:
1.19      markus    102:        BN_clear_free(dhg->g);
                    103:        BN_clear_free(dhg->p);
1.1       provos    104:  fail:
1.7       millert   105:        error("Bad prime description in line %d", linenum);
1.1       provos    106:        return (0);
                    107: }
                    108:
                    109: DH *
1.9       provos    110: choose_dh(int min, int wantbits, int max)
1.1       provos    111: {
                    112:        FILE *f;
1.28      dtucker   113:        char line[4096];
1.1       provos    114:        int best, bestcount, which;
                    115:        int linenum;
                    116:        struct dhgroup dhg;
                    117:
1.15      provos    118:        if ((f = fopen(_PATH_DH_MODULI, "r")) == NULL &&
                    119:            (f = fopen(_PATH_DH_PRIMES, "r")) == NULL) {
1.30      djm       120:                logit("WARNING: %s does not exist, using fixed modulus",
                    121:                    _PATH_DH_MODULI);
                    122:                return (dh_new_group14());
1.1       provos    123:        }
                    124:
                    125:        linenum = 0;
                    126:        best = bestcount = 0;
                    127:        while (fgets(line, sizeof(line), f)) {
                    128:                linenum++;
                    129:                if (!parse_prime(linenum, line, &dhg))
                    130:                        continue;
1.19      markus    131:                BN_clear_free(dhg.g);
                    132:                BN_clear_free(dhg.p);
1.1       provos    133:
1.9       provos    134:                if (dhg.size > max || dhg.size < min)
                    135:                        continue;
                    136:
                    137:                if ((dhg.size > wantbits && dhg.size < best) ||
                    138:                    (dhg.size > best && best < wantbits)) {
1.1       provos    139:                        best = dhg.size;
                    140:                        bestcount = 0;
                    141:                }
                    142:                if (dhg.size == best)
                    143:                        bestcount++;
                    144:        }
1.16      provos    145:        rewind(f);
1.1       provos    146:
                    147:        if (bestcount == 0) {
1.16      provos    148:                fclose(f);
1.24      itojun    149:                logit("WARNING: no suitable primes in %s", _PATH_DH_PRIMES);
1.31      djm       150:                return (dh_new_group14());
1.1       provos    151:        }
                    152:
                    153:        linenum = 0;
                    154:        which = arc4random() % bestcount;
                    155:        while (fgets(line, sizeof(line), f)) {
                    156:                if (!parse_prime(linenum, line, &dhg))
1.9       provos    157:                        continue;
1.13      markus    158:                if ((dhg.size > max || dhg.size < min) ||
                    159:                    dhg.size != best ||
                    160:                    linenum++ != which) {
1.19      markus    161:                        BN_clear_free(dhg.g);
                    162:                        BN_clear_free(dhg.p);
1.1       provos    163:                        continue;
                    164:                }
                    165:                break;
                    166:        }
                    167:        fclose(f);
1.13      markus    168:        if (linenum != which+1)
                    169:                fatal("WARNING: line %d disappeared in %s, giving up",
                    170:                    which, _PATH_DH_PRIMES);
1.1       provos    171:
                    172:        return (dh_new_group(dhg.g, dhg.p));
1.11      markus    173: }
                    174:
1.30      djm       175: /* diffie-hellman-groupN-sha1 */
1.11      markus    176:
                    177: int
                    178: dh_pub_is_valid(DH *dh, BIGNUM *dh_pub)
                    179: {
                    180:        int i;
                    181:        int n = BN_num_bits(dh_pub);
                    182:        int bits_set = 0;
1.36      djm       183:        BIGNUM *tmp;
1.11      markus    184:
                    185:        if (dh_pub->neg) {
1.24      itojun    186:                logit("invalid public DH value: negativ");
1.11      markus    187:                return 0;
                    188:        }
1.36      djm       189:        if (BN_cmp(dh_pub, BN_value_one()) != 1) {      /* pub_exp <= 1 */
                    190:                logit("invalid public DH value: <= 1");
                    191:                return 0;
                    192:        }
                    193:
                    194:        if ((tmp = BN_new()) == NULL)
                    195:                return (-1);
                    196:        if (!BN_sub(tmp, dh->p, BN_value_one()) ||
                    197:            BN_cmp(dh_pub, tmp) != -1) {                /* pub_exp > p-2 */
                    198:                BN_clear_free(tmp);
                    199:                logit("invalid public DH value: >= p-1");
                    200:                return 0;
                    201:        }
                    202:        BN_clear_free(tmp);
                    203:
1.11      markus    204:        for (i = 0; i <= n; i++)
                    205:                if (BN_is_bit_set(dh_pub, i))
                    206:                        bits_set++;
1.23      markus    207:        debug2("bits set: %d/%d", bits_set, BN_num_bits(dh->p));
1.11      markus    208:
                    209:        /* if g==2 and bits_set==1 then computing log_g(dh_pub) is trivial */
1.36      djm       210:        if (bits_set > 1)
1.11      markus    211:                return 1;
1.36      djm       212:
1.24      itojun    213:        logit("invalid public DH value (%d/%d)", bits_set, BN_num_bits(dh->p));
1.11      markus    214:        return 0;
                    215: }
                    216:
                    217: void
                    218: dh_gen_key(DH *dh, int need)
                    219: {
1.29      dtucker   220:        int i, bits_set, tries = 0;
1.11      markus    221:
                    222:        if (dh->p == NULL)
                    223:                fatal("dh_gen_key: dh->p == NULL");
1.25      miod      224:        if (need > INT_MAX / 2 || 2 * need >= BN_num_bits(dh->p))
1.11      markus    225:                fatal("dh_gen_key: group too small: %d (2*need %d)",
                    226:                    BN_num_bits(dh->p), 2*need);
                    227:        do {
                    228:                if (dh->priv_key != NULL)
1.19      markus    229:                        BN_clear_free(dh->priv_key);
1.18      markus    230:                if ((dh->priv_key = BN_new()) == NULL)
1.11      markus    231:                        fatal("dh_gen_key: BN_new failed");
                    232:                /* generate a 2*need bits random private exponent */
1.21      markus    233:                if (!BN_rand(dh->priv_key, 2*need, 0, 0))
                    234:                        fatal("dh_gen_key: BN_rand failed");
1.11      markus    235:                if (DH_generate_key(dh) == 0)
                    236:                        fatal("DH_generate_key");
1.29      dtucker   237:                for (i = 0, bits_set = 0; i <= BN_num_bits(dh->priv_key); i++)
1.11      markus    238:                        if (BN_is_bit_set(dh->priv_key, i))
                    239:                                bits_set++;
1.23      markus    240:                debug2("dh_gen_key: priv key bits set: %d/%d",
1.11      markus    241:                    bits_set, BN_num_bits(dh->priv_key));
                    242:                if (tries++ > 10)
                    243:                        fatal("dh_gen_key: too many bad keys: giving up");
                    244:        } while (!dh_pub_is_valid(dh, dh->pub_key));
                    245: }
                    246:
                    247: DH *
                    248: dh_new_group_asc(const char *gen, const char *modulus)
                    249: {
                    250:        DH *dh;
                    251:
1.18      markus    252:        if ((dh = DH_new()) == NULL)
                    253:                fatal("dh_new_group_asc: DH_new");
1.11      markus    254:
1.14      markus    255:        if (BN_hex2bn(&dh->p, modulus) == 0)
1.11      markus    256:                fatal("BN_hex2bn p");
1.14      markus    257:        if (BN_hex2bn(&dh->g, gen) == 0)
1.11      markus    258:                fatal("BN_hex2bn g");
                    259:
                    260:        return (dh);
                    261: }
                    262:
                    263: /*
                    264:  * This just returns the group, we still need to generate the exchange
                    265:  * value.
                    266:  */
                    267:
                    268: DH *
                    269: dh_new_group(BIGNUM *gen, BIGNUM *modulus)
                    270: {
                    271:        DH *dh;
                    272:
1.18      markus    273:        if ((dh = DH_new()) == NULL)
                    274:                fatal("dh_new_group: DH_new");
1.11      markus    275:        dh->p = modulus;
                    276:        dh->g = gen;
                    277:
                    278:        return (dh);
                    279: }
                    280:
                    281: DH *
                    282: dh_new_group1(void)
                    283: {
                    284:        static char *gen = "2", *group1 =
                    285:            "FFFFFFFF" "FFFFFFFF" "C90FDAA2" "2168C234" "C4C6628B" "80DC1CD1"
                    286:            "29024E08" "8A67CC74" "020BBEA6" "3B139B22" "514A0879" "8E3404DD"
                    287:            "EF9519B3" "CD3A431B" "302B0A6D" "F25F1437" "4FE1356D" "6D51C245"
                    288:            "E485B576" "625E7EC6" "F44C42E9" "A637ED6B" "0BFF5CB6" "F406B7ED"
                    289:            "EE386BFB" "5A899FA5" "AE9F2411" "7C4B1FE6" "49286651" "ECE65381"
                    290:            "FFFFFFFF" "FFFFFFFF";
                    291:
                    292:        return (dh_new_group_asc(gen, group1));
1.30      djm       293: }
                    294:
                    295: DH *
                    296: dh_new_group14(void)
                    297: {
                    298:        static char *gen = "2", *group14 =
                    299:            "FFFFFFFF" "FFFFFFFF" "C90FDAA2" "2168C234" "C4C6628B" "80DC1CD1"
                    300:            "29024E08" "8A67CC74" "020BBEA6" "3B139B22" "514A0879" "8E3404DD"
                    301:            "EF9519B3" "CD3A431B" "302B0A6D" "F25F1437" "4FE1356D" "6D51C245"
                    302:            "E485B576" "625E7EC6" "F44C42E9" "A637ED6B" "0BFF5CB6" "F406B7ED"
                    303:            "EE386BFB" "5A899FA5" "AE9F2411" "7C4B1FE6" "49286651" "ECE45B3D"
                    304:            "C2007CB8" "A163BF05" "98DA4836" "1C55D39A" "69163FA8" "FD24CF5F"
                    305:            "83655D23" "DCA3AD96" "1C62F356" "208552BB" "9ED52907" "7096966D"
                    306:            "670C354E" "4ABC9804" "F1746C08" "CA18217C" "32905E46" "2E36CE3B"
                    307:            "E39E772C" "180E8603" "9B2783A2" "EC07A28F" "B5C55DF0" "6F4C52C9"
                    308:            "DE2BCBF6" "95581718" "3995497C" "EA956AE5" "15D22618" "98FA0510"
                    309:            "15728E5A" "8AACAA68" "FFFFFFFF" "FFFFFFFF";
                    310:
                    311:        return (dh_new_group_asc(gen, group14));
1.12      markus    312: }
                    313:
                    314: /*
                    315:  * Estimates the group order for a Diffie-Hellman group that has an
                    316:  * attack complexity approximately the same as O(2**bits).  Estimate
                    317:  * with:  O(exp(1.9223 * (ln q)^(1/3) (ln ln q)^(2/3)))
                    318:  */
                    319:
                    320: int
                    321: dh_estimate(int bits)
                    322: {
                    323:
1.26      markus    324:        if (bits <= 128)
1.12      markus    325:                return (1024);  /* O(2**86) */
1.26      markus    326:        if (bits <= 192)
1.12      markus    327:                return (2048);  /* O(2**116) */
                    328:        return (4096);          /* O(2**156) */
1.1       provos    329: }