]>
git.ipfire.org Git - thirdparty/openssl.git/blob - crypto/bn/bn_x931p.c
3 * Written by Dr Stephen N Henson (steve@openssl.org) for the OpenSSL project
6 /* ====================================================================
7 * Copyright (c) 2005 The OpenSSL Project. All rights reserved.
9 * Redistribution and use in source and binary forms, with or without
10 * modification, are permitted provided that the following conditions
13 * 1. Redistributions of source code must retain the above copyright
14 * notice, this list of conditions and the following disclaimer.
16 * 2. Redistributions in binary form must reproduce the above copyright
17 * notice, this list of conditions and the following disclaimer in
18 * the documentation and/or other materials provided with the
21 * 3. All advertising materials mentioning features or use of this
22 * software must display the following acknowledgment:
23 * "This product includes software developed by the OpenSSL Project
24 * for use in the OpenSSL Toolkit. (http://www.OpenSSL.org/)"
26 * 4. The names "OpenSSL Toolkit" and "OpenSSL Project" must not be used to
27 * endorse or promote products derived from this software without
28 * prior written permission. For written permission, please contact
29 * licensing@OpenSSL.org.
31 * 5. Products derived from this software may not be called "OpenSSL"
32 * nor may "OpenSSL" appear in their names without prior written
33 * permission of the OpenSSL Project.
35 * 6. Redistributions of any form whatsoever must retain the following
37 * "This product includes software developed by the OpenSSL Project
38 * for use in the OpenSSL Toolkit (http://www.OpenSSL.org/)"
40 * THIS SOFTWARE IS PROVIDED BY THE OpenSSL PROJECT ``AS IS'' AND ANY
41 * EXPRESSED OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
42 * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
43 * PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE OpenSSL PROJECT OR
44 * ITS CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
45 * SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT
46 * NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
47 * LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
48 * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT,
49 * STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
50 * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED
51 * OF THE POSSIBILITY OF SUCH DAMAGE.
52 * ====================================================================
54 * This product includes cryptographic software written by Eric Young
55 * (eay@cryptsoft.com). This product includes software written by Tim
56 * Hudson (tjh@cryptsoft.com).
61 #include <openssl/bn.h>
63 /* X9.31 routines for prime derivation */
66 * X9.31 prime derivation. This is used to generate the primes pi (p1, p2,
67 * q1, q2) from a parameter Xpi by checking successive odd integers.
70 static int bn_x931_derive_pi(BIGNUM
*pi
, const BIGNUM
*Xpi
, BN_CTX
*ctx
,
74 if (!BN_copy(pi
, Xpi
))
76 if (!BN_is_odd(pi
) && !BN_add_word(pi
, 1))
80 BN_GENCB_call(cb
, 0, i
);
81 /* NB 27 MR is specificed in X9.31 */
82 if (BN_is_prime_fasttest_ex(pi
, 27, ctx
, 1, cb
))
84 if (!BN_add_word(pi
, 2))
87 BN_GENCB_call(cb
, 2, i
);
92 * This is the main X9.31 prime derivation function. From parameters Xp1, Xp2
93 * and Xp derive the prime p. If the parameters p1 or p2 are not NULL they
94 * will be returned too: this is needed for testing.
97 int BN_X931_derive_prime_ex(BIGNUM
*p
, BIGNUM
*p1
, BIGNUM
*p2
,
98 const BIGNUM
*Xp
, const BIGNUM
*Xp1
,
99 const BIGNUM
*Xp2
, const BIGNUM
*e
, BN_CTX
*ctx
,
104 BIGNUM
*t
, *p1p2
, *pm1
;
106 /* Only even e supported */
112 p1
= BN_CTX_get(ctx
);
115 p2
= BN_CTX_get(ctx
);
119 p1p2
= BN_CTX_get(ctx
);
121 pm1
= BN_CTX_get(ctx
);
123 if (!bn_x931_derive_pi(p1
, Xp1
, ctx
, cb
))
126 if (!bn_x931_derive_pi(p2
, Xp2
, ctx
, cb
))
129 if (!BN_mul(p1p2
, p1
, p2
, ctx
))
132 /* First set p to value of Rp */
134 if (!BN_mod_inverse(p
, p2
, p1
, ctx
))
137 if (!BN_mul(p
, p
, p2
, ctx
))
140 if (!BN_mod_inverse(t
, p1
, p2
, ctx
))
143 if (!BN_mul(t
, t
, p1
, ctx
))
146 if (!BN_sub(p
, p
, t
))
149 if (p
->neg
&& !BN_add(p
, p
, p1p2
))
152 /* p now equals Rp */
154 if (!BN_mod_sub(p
, p
, Xp
, p1p2
, ctx
))
157 if (!BN_add(p
, p
, Xp
))
160 /* p now equals Yp0 */
164 BN_GENCB_call(cb
, 0, i
++);
165 if (!BN_copy(pm1
, p
))
167 if (!BN_sub_word(pm1
, 1))
169 if (!BN_gcd(t
, pm1
, e
, ctx
))
173 * X9.31 specifies 8 MR and 1 Lucas test or any prime test
174 * offering similar or better guarantees 50 MR is considerably
177 && BN_is_prime_fasttest_ex(p
, 50, ctx
, 1, cb
))
179 if (!BN_add(p
, p
, p1p2
))
183 BN_GENCB_call(cb
, 3, 0);
195 * Generate pair of paramters Xp, Xq for X9.31 prime generation. Note: nbits
196 * paramter is sum of number of bits in both.
199 int BN_X931_generate_Xpq(BIGNUM
*Xp
, BIGNUM
*Xq
, int nbits
, BN_CTX
*ctx
)
204 * Number of bits for each prime is of the form 512+128s for s = 0, 1,
207 if ((nbits
< 1024) || (nbits
& 0xff))
211 * The random value Xp must be between sqrt(2) * 2^(nbits-1) and 2^nbits
212 * - 1. By setting the top two bits we ensure that the lower bound is
215 if (!BN_rand(Xp
, nbits
, 1, 0))
221 for (i
= 0; i
< 1000; i
++) {
222 if (!BN_rand(Xq
, nbits
, 1, 0))
224 /* Check that |Xp - Xq| > 2^(nbits - 100) */
226 if (BN_num_bits(t
) > (nbits
- 100))
240 * Generate primes using X9.31 algorithm. Of the values p, p1, p2, Xp1 and
241 * Xp2 only 'p' needs to be non-NULL. If any of the others are not NULL the
242 * relevant parameter will be stored in it. Due to the fact that |Xp - Xq| >
243 * 2^(nbits - 100) must be satisfied Xp and Xq are generated using the
244 * previous function and supplied as input.
247 int BN_X931_generate_prime_ex(BIGNUM
*p
, BIGNUM
*p1
, BIGNUM
*p2
,
248 BIGNUM
*Xp1
, BIGNUM
*Xp2
,
250 const BIGNUM
*e
, BN_CTX
*ctx
, BN_GENCB
*cb
)
256 Xp1
= BN_CTX_get(ctx
);
258 Xp2
= BN_CTX_get(ctx
);
260 if (!BN_rand(Xp1
, 101, 0, 0))
262 if (!BN_rand(Xp2
, 101, 0, 0))
264 if (!BN_X931_derive_prime_ex(p
, p1
, p2
, Xp
, Xp1
, Xp2
, e
, ctx
, cb
))