xref: /aosp_15_r20/external/mbedtls/programs/pkey/dh_genprime.c (revision 62c56f9862f102b96d72393aff6076c951fb8148)
1*62c56f98SSadaf Ebrahimi /*
2*62c56f98SSadaf Ebrahimi  *  Diffie-Hellman-Merkle key exchange (prime generation)
3*62c56f98SSadaf Ebrahimi  *
4*62c56f98SSadaf Ebrahimi  *  Copyright The Mbed TLS Contributors
5*62c56f98SSadaf Ebrahimi  *  SPDX-License-Identifier: Apache-2.0 OR GPL-2.0-or-later
6*62c56f98SSadaf Ebrahimi  */
7*62c56f98SSadaf Ebrahimi 
8*62c56f98SSadaf Ebrahimi #include "mbedtls/build_info.h"
9*62c56f98SSadaf Ebrahimi 
10*62c56f98SSadaf Ebrahimi #include "mbedtls/platform.h"
11*62c56f98SSadaf Ebrahimi 
12*62c56f98SSadaf Ebrahimi #if !defined(MBEDTLS_BIGNUM_C) || !defined(MBEDTLS_ENTROPY_C) ||   \
13*62c56f98SSadaf Ebrahimi     !defined(MBEDTLS_FS_IO) || !defined(MBEDTLS_CTR_DRBG_C) ||     \
14*62c56f98SSadaf Ebrahimi     !defined(MBEDTLS_GENPRIME)
main(void)15*62c56f98SSadaf Ebrahimi int main(void)
16*62c56f98SSadaf Ebrahimi {
17*62c56f98SSadaf Ebrahimi     mbedtls_printf("MBEDTLS_BIGNUM_C and/or MBEDTLS_ENTROPY_C and/or "
18*62c56f98SSadaf Ebrahimi                    "MBEDTLS_FS_IO and/or MBEDTLS_CTR_DRBG_C and/or "
19*62c56f98SSadaf Ebrahimi                    "MBEDTLS_GENPRIME not defined.\n");
20*62c56f98SSadaf Ebrahimi     mbedtls_exit(0);
21*62c56f98SSadaf Ebrahimi }
22*62c56f98SSadaf Ebrahimi #else
23*62c56f98SSadaf Ebrahimi 
24*62c56f98SSadaf Ebrahimi #include "mbedtls/bignum.h"
25*62c56f98SSadaf Ebrahimi #include "mbedtls/entropy.h"
26*62c56f98SSadaf Ebrahimi #include "mbedtls/ctr_drbg.h"
27*62c56f98SSadaf Ebrahimi 
28*62c56f98SSadaf Ebrahimi #include <stdio.h>
29*62c56f98SSadaf Ebrahimi #include <string.h>
30*62c56f98SSadaf Ebrahimi 
31*62c56f98SSadaf Ebrahimi #define USAGE \
32*62c56f98SSadaf Ebrahimi     "\n usage: dh_genprime param=<>...\n"                                   \
33*62c56f98SSadaf Ebrahimi     "\n acceptable parameters:\n"                                           \
34*62c56f98SSadaf Ebrahimi     "    bits=%%d           default: 2048\n"
35*62c56f98SSadaf Ebrahimi 
36*62c56f98SSadaf Ebrahimi #define DFL_BITS    2048
37*62c56f98SSadaf Ebrahimi 
38*62c56f98SSadaf Ebrahimi /*
39*62c56f98SSadaf Ebrahimi  * Note: G = 4 is always a quadratic residue mod P,
40*62c56f98SSadaf Ebrahimi  * so it is a generator of order Q (with P = 2*Q+1).
41*62c56f98SSadaf Ebrahimi  */
42*62c56f98SSadaf Ebrahimi #define GENERATOR "4"
43*62c56f98SSadaf Ebrahimi 
44*62c56f98SSadaf Ebrahimi 
main(int argc,char ** argv)45*62c56f98SSadaf Ebrahimi int main(int argc, char **argv)
46*62c56f98SSadaf Ebrahimi {
47*62c56f98SSadaf Ebrahimi     int ret = 1;
48*62c56f98SSadaf Ebrahimi     int exit_code = MBEDTLS_EXIT_FAILURE;
49*62c56f98SSadaf Ebrahimi     mbedtls_mpi G, P, Q;
50*62c56f98SSadaf Ebrahimi     mbedtls_entropy_context entropy;
51*62c56f98SSadaf Ebrahimi     mbedtls_ctr_drbg_context ctr_drbg;
52*62c56f98SSadaf Ebrahimi     const char *pers = "dh_genprime";
53*62c56f98SSadaf Ebrahimi     FILE *fout;
54*62c56f98SSadaf Ebrahimi     int nbits = DFL_BITS;
55*62c56f98SSadaf Ebrahimi     int i;
56*62c56f98SSadaf Ebrahimi     char *p, *q;
57*62c56f98SSadaf Ebrahimi 
58*62c56f98SSadaf Ebrahimi     mbedtls_mpi_init(&G); mbedtls_mpi_init(&P); mbedtls_mpi_init(&Q);
59*62c56f98SSadaf Ebrahimi     mbedtls_ctr_drbg_init(&ctr_drbg);
60*62c56f98SSadaf Ebrahimi     mbedtls_entropy_init(&entropy);
61*62c56f98SSadaf Ebrahimi 
62*62c56f98SSadaf Ebrahimi     if (argc < 2) {
63*62c56f98SSadaf Ebrahimi usage:
64*62c56f98SSadaf Ebrahimi         mbedtls_printf(USAGE);
65*62c56f98SSadaf Ebrahimi         goto exit;
66*62c56f98SSadaf Ebrahimi     }
67*62c56f98SSadaf Ebrahimi 
68*62c56f98SSadaf Ebrahimi     for (i = 1; i < argc; i++) {
69*62c56f98SSadaf Ebrahimi         p = argv[i];
70*62c56f98SSadaf Ebrahimi         if ((q = strchr(p, '=')) == NULL) {
71*62c56f98SSadaf Ebrahimi             goto usage;
72*62c56f98SSadaf Ebrahimi         }
73*62c56f98SSadaf Ebrahimi         *q++ = '\0';
74*62c56f98SSadaf Ebrahimi 
75*62c56f98SSadaf Ebrahimi         if (strcmp(p, "bits") == 0) {
76*62c56f98SSadaf Ebrahimi             nbits = atoi(q);
77*62c56f98SSadaf Ebrahimi             if (nbits < 0 || nbits > MBEDTLS_MPI_MAX_BITS) {
78*62c56f98SSadaf Ebrahimi                 goto usage;
79*62c56f98SSadaf Ebrahimi             }
80*62c56f98SSadaf Ebrahimi         } else {
81*62c56f98SSadaf Ebrahimi             goto usage;
82*62c56f98SSadaf Ebrahimi         }
83*62c56f98SSadaf Ebrahimi     }
84*62c56f98SSadaf Ebrahimi 
85*62c56f98SSadaf Ebrahimi     if ((ret = mbedtls_mpi_read_string(&G, 10, GENERATOR)) != 0) {
86*62c56f98SSadaf Ebrahimi         mbedtls_printf(" failed\n  ! mbedtls_mpi_read_string returned %d\n", ret);
87*62c56f98SSadaf Ebrahimi         goto exit;
88*62c56f98SSadaf Ebrahimi     }
89*62c56f98SSadaf Ebrahimi 
90*62c56f98SSadaf Ebrahimi     mbedtls_printf("  ! Generating large primes may take minutes!\n");
91*62c56f98SSadaf Ebrahimi 
92*62c56f98SSadaf Ebrahimi     mbedtls_printf("\n  . Seeding the random number generator...");
93*62c56f98SSadaf Ebrahimi     fflush(stdout);
94*62c56f98SSadaf Ebrahimi 
95*62c56f98SSadaf Ebrahimi     if ((ret = mbedtls_ctr_drbg_seed(&ctr_drbg, mbedtls_entropy_func, &entropy,
96*62c56f98SSadaf Ebrahimi                                      (const unsigned char *) pers,
97*62c56f98SSadaf Ebrahimi                                      strlen(pers))) != 0) {
98*62c56f98SSadaf Ebrahimi         mbedtls_printf(" failed\n  ! mbedtls_ctr_drbg_seed returned %d\n", ret);
99*62c56f98SSadaf Ebrahimi         goto exit;
100*62c56f98SSadaf Ebrahimi     }
101*62c56f98SSadaf Ebrahimi 
102*62c56f98SSadaf Ebrahimi     mbedtls_printf(" ok\n  . Generating the modulus, please wait...");
103*62c56f98SSadaf Ebrahimi     fflush(stdout);
104*62c56f98SSadaf Ebrahimi 
105*62c56f98SSadaf Ebrahimi     /*
106*62c56f98SSadaf Ebrahimi      * This can take a long time...
107*62c56f98SSadaf Ebrahimi      */
108*62c56f98SSadaf Ebrahimi     if ((ret = mbedtls_mpi_gen_prime(&P, nbits, 1,
109*62c56f98SSadaf Ebrahimi                                      mbedtls_ctr_drbg_random, &ctr_drbg)) != 0) {
110*62c56f98SSadaf Ebrahimi         mbedtls_printf(" failed\n  ! mbedtls_mpi_gen_prime returned %d\n\n", ret);
111*62c56f98SSadaf Ebrahimi         goto exit;
112*62c56f98SSadaf Ebrahimi     }
113*62c56f98SSadaf Ebrahimi 
114*62c56f98SSadaf Ebrahimi     mbedtls_printf(" ok\n  . Verifying that Q = (P-1)/2 is prime...");
115*62c56f98SSadaf Ebrahimi     fflush(stdout);
116*62c56f98SSadaf Ebrahimi 
117*62c56f98SSadaf Ebrahimi     if ((ret = mbedtls_mpi_sub_int(&Q, &P, 1)) != 0) {
118*62c56f98SSadaf Ebrahimi         mbedtls_printf(" failed\n  ! mbedtls_mpi_sub_int returned %d\n\n", ret);
119*62c56f98SSadaf Ebrahimi         goto exit;
120*62c56f98SSadaf Ebrahimi     }
121*62c56f98SSadaf Ebrahimi 
122*62c56f98SSadaf Ebrahimi     if ((ret = mbedtls_mpi_div_int(&Q, NULL, &Q, 2)) != 0) {
123*62c56f98SSadaf Ebrahimi         mbedtls_printf(" failed\n  ! mbedtls_mpi_div_int returned %d\n\n", ret);
124*62c56f98SSadaf Ebrahimi         goto exit;
125*62c56f98SSadaf Ebrahimi     }
126*62c56f98SSadaf Ebrahimi 
127*62c56f98SSadaf Ebrahimi     if ((ret = mbedtls_mpi_is_prime_ext(&Q, 50, mbedtls_ctr_drbg_random, &ctr_drbg)) != 0) {
128*62c56f98SSadaf Ebrahimi         mbedtls_printf(" failed\n  ! mbedtls_mpi_is_prime returned %d\n\n", ret);
129*62c56f98SSadaf Ebrahimi         goto exit;
130*62c56f98SSadaf Ebrahimi     }
131*62c56f98SSadaf Ebrahimi 
132*62c56f98SSadaf Ebrahimi     mbedtls_printf(" ok\n  . Exporting the value in dh_prime.txt...");
133*62c56f98SSadaf Ebrahimi     fflush(stdout);
134*62c56f98SSadaf Ebrahimi 
135*62c56f98SSadaf Ebrahimi     if ((fout = fopen("dh_prime.txt", "wb+")) == NULL) {
136*62c56f98SSadaf Ebrahimi         mbedtls_printf(" failed\n  ! Could not create dh_prime.txt\n\n");
137*62c56f98SSadaf Ebrahimi         goto exit;
138*62c56f98SSadaf Ebrahimi     }
139*62c56f98SSadaf Ebrahimi 
140*62c56f98SSadaf Ebrahimi     if (((ret = mbedtls_mpi_write_file("P = ", &P, 16, fout)) != 0) ||
141*62c56f98SSadaf Ebrahimi         ((ret = mbedtls_mpi_write_file("G = ", &G, 16, fout)) != 0)) {
142*62c56f98SSadaf Ebrahimi         mbedtls_printf(" failed\n  ! mbedtls_mpi_write_file returned %d\n\n", ret);
143*62c56f98SSadaf Ebrahimi         fclose(fout);
144*62c56f98SSadaf Ebrahimi         goto exit;
145*62c56f98SSadaf Ebrahimi     }
146*62c56f98SSadaf Ebrahimi 
147*62c56f98SSadaf Ebrahimi     mbedtls_printf(" ok\n\n");
148*62c56f98SSadaf Ebrahimi     fclose(fout);
149*62c56f98SSadaf Ebrahimi 
150*62c56f98SSadaf Ebrahimi     exit_code = MBEDTLS_EXIT_SUCCESS;
151*62c56f98SSadaf Ebrahimi 
152*62c56f98SSadaf Ebrahimi exit:
153*62c56f98SSadaf Ebrahimi 
154*62c56f98SSadaf Ebrahimi     mbedtls_mpi_free(&G); mbedtls_mpi_free(&P); mbedtls_mpi_free(&Q);
155*62c56f98SSadaf Ebrahimi     mbedtls_ctr_drbg_free(&ctr_drbg);
156*62c56f98SSadaf Ebrahimi     mbedtls_entropy_free(&entropy);
157*62c56f98SSadaf Ebrahimi 
158*62c56f98SSadaf Ebrahimi     mbedtls_exit(exit_code);
159*62c56f98SSadaf Ebrahimi }
160*62c56f98SSadaf Ebrahimi #endif /* MBEDTLS_BIGNUM_C && MBEDTLS_ENTROPY_C && MBEDTLS_FS_IO &&
161*62c56f98SSadaf Ebrahimi           MBEDTLS_CTR_DRBG_C && MBEDTLS_GENPRIME */
162