dsa_ossl.c 13 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439
  1. /* crypto/dsa/dsa_ossl.c */
  2. /* Copyright (C) 1995-1998 Eric Young ([email protected])
  3. * All rights reserved.
  4. *
  5. * This package is an SSL implementation written
  6. * by Eric Young ([email protected]).
  7. * The implementation was written so as to conform with Netscapes SSL.
  8. *
  9. * This library is free for commercial and non-commercial use as long as
  10. * the following conditions are aheared to. The following conditions
  11. * apply to all code found in this distribution, be it the RC4, RSA,
  12. * lhash, DES, etc., code; not just the SSL code. The SSL documentation
  13. * included with this distribution is covered by the same copyright terms
  14. * except that the holder is Tim Hudson ([email protected]).
  15. *
  16. * Copyright remains Eric Young's, and as such any Copyright notices in
  17. * the code are not to be removed.
  18. * If this package is used in a product, Eric Young should be given attribution
  19. * as the author of the parts of the library used.
  20. * This can be in the form of a textual message at program startup or
  21. * in documentation (online or textual) provided with the package.
  22. *
  23. * Redistribution and use in source and binary forms, with or without
  24. * modification, are permitted provided that the following conditions
  25. * are met:
  26. * 1. Redistributions of source code must retain the copyright
  27. * notice, this list of conditions and the following disclaimer.
  28. * 2. Redistributions in binary form must reproduce the above copyright
  29. * notice, this list of conditions and the following disclaimer in the
  30. * documentation and/or other materials provided with the distribution.
  31. * 3. All advertising materials mentioning features or use of this software
  32. * must display the following acknowledgement:
  33. * "This product includes cryptographic software written by
  34. * Eric Young ([email protected])"
  35. * The word 'cryptographic' can be left out if the rouines from the library
  36. * being used are not cryptographic related :-).
  37. * 4. If you include any Windows specific code (or a derivative thereof) from
  38. * the apps directory (application code) you must include an acknowledgement:
  39. * "This product includes software written by Tim Hudson ([email protected])"
  40. *
  41. * THIS SOFTWARE IS PROVIDED BY ERIC YOUNG ``AS IS'' AND
  42. * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
  43. * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
  44. * ARE DISCLAIMED. IN NO EVENT SHALL THE AUTHOR OR CONTRIBUTORS BE LIABLE
  45. * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
  46. * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
  47. * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
  48. * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
  49. * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
  50. * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
  51. * SUCH DAMAGE.
  52. *
  53. * The licence and distribution terms for any publically available version or
  54. * derivative of this code cannot be changed. i.e. this code cannot simply be
  55. * copied and put under another distribution licence
  56. * [including the GNU Public Licence.]
  57. */
  58. /* Original version from Steven Schoch <[email protected]> */
  59. #include <stdio.h>
  60. #include "cryptlib.h"
  61. #include <openssl/bn.h>
  62. #include <openssl/sha.h>
  63. #include <openssl/dsa.h>
  64. #include <openssl/rand.h>
  65. #include <openssl/asn1.h>
  66. static DSA_SIG *dsa_do_sign(const unsigned char *dgst, int dlen, DSA *dsa);
  67. static int dsa_sign_setup(DSA *dsa, BN_CTX *ctx_in, BIGNUM **kinvp,
  68. BIGNUM **rp);
  69. static int dsa_do_verify(const unsigned char *dgst, int dgst_len,
  70. DSA_SIG *sig, DSA *dsa);
  71. static int dsa_init(DSA *dsa);
  72. static int dsa_finish(DSA *dsa);
  73. static DSA_METHOD openssl_dsa_meth = {
  74. "OpenSSL DSA method",
  75. dsa_do_sign,
  76. dsa_sign_setup,
  77. dsa_do_verify,
  78. NULL, /* dsa_mod_exp, */
  79. NULL, /* dsa_bn_mod_exp, */
  80. dsa_init,
  81. dsa_finish,
  82. 0,
  83. NULL,
  84. NULL,
  85. NULL
  86. };
  87. /*-
  88. * These macro wrappers replace attempts to use the dsa_mod_exp() and
  89. * bn_mod_exp() handlers in the DSA_METHOD structure. We avoid the problem of
  90. * having a the macro work as an expression by bundling an "err_instr". So;
  91. *
  92. * if (!dsa->meth->bn_mod_exp(dsa, r,dsa->g,&k,dsa->p,ctx,
  93. * dsa->method_mont_p)) goto err;
  94. *
  95. * can be replaced by;
  96. *
  97. * DSA_BN_MOD_EXP(goto err, dsa, r, dsa->g, &k, dsa->p, ctx,
  98. * dsa->method_mont_p);
  99. */
  100. #define DSA_MOD_EXP(err_instr,dsa,rr,a1,p1,a2,p2,m,ctx,in_mont) \
  101. do { \
  102. int _tmp_res53; \
  103. if ((dsa)->meth->dsa_mod_exp) \
  104. _tmp_res53 = (dsa)->meth->dsa_mod_exp((dsa), (rr), (a1), (p1), \
  105. (a2), (p2), (m), (ctx), (in_mont)); \
  106. else \
  107. _tmp_res53 = BN_mod_exp2_mont((rr), (a1), (p1), (a2), (p2), \
  108. (m), (ctx), (in_mont)); \
  109. if (!_tmp_res53) err_instr; \
  110. } while(0)
  111. #define DSA_BN_MOD_EXP(err_instr,dsa,r,a,p,m,ctx,m_ctx) \
  112. do { \
  113. int _tmp_res53; \
  114. if ((dsa)->meth->bn_mod_exp) \
  115. _tmp_res53 = (dsa)->meth->bn_mod_exp((dsa), (r), (a), (p), \
  116. (m), (ctx), (m_ctx)); \
  117. else \
  118. _tmp_res53 = BN_mod_exp_mont((r), (a), (p), (m), (ctx), (m_ctx)); \
  119. if (!_tmp_res53) err_instr; \
  120. } while(0)
  121. const DSA_METHOD *DSA_OpenSSL(void)
  122. {
  123. return &openssl_dsa_meth;
  124. }
  125. static DSA_SIG *dsa_do_sign(const unsigned char *dgst, int dlen, DSA *dsa)
  126. {
  127. BIGNUM *kinv = NULL, *r = NULL, *s = NULL;
  128. BIGNUM m;
  129. BIGNUM xr;
  130. BN_CTX *ctx = NULL;
  131. int reason = ERR_R_BN_LIB;
  132. DSA_SIG *ret = NULL;
  133. int noredo = 0;
  134. BN_init(&m);
  135. BN_init(&xr);
  136. if (!dsa->p || !dsa->q || !dsa->g) {
  137. reason = DSA_R_MISSING_PARAMETERS;
  138. goto err;
  139. }
  140. s = BN_new();
  141. if (s == NULL)
  142. goto err;
  143. ctx = BN_CTX_new();
  144. if (ctx == NULL)
  145. goto err;
  146. redo:
  147. if ((dsa->kinv == NULL) || (dsa->r == NULL)) {
  148. if (!DSA_sign_setup(dsa, ctx, &kinv, &r))
  149. goto err;
  150. } else {
  151. kinv = dsa->kinv;
  152. dsa->kinv = NULL;
  153. r = dsa->r;
  154. dsa->r = NULL;
  155. noredo = 1;
  156. }
  157. if (dlen > BN_num_bytes(dsa->q))
  158. /*
  159. * if the digest length is greater than the size of q use the
  160. * BN_num_bits(dsa->q) leftmost bits of the digest, see fips 186-3,
  161. * 4.2
  162. */
  163. dlen = BN_num_bytes(dsa->q);
  164. if (BN_bin2bn(dgst, dlen, &m) == NULL)
  165. goto err;
  166. /* Compute s = inv(k) (m + xr) mod q */
  167. if (!BN_mod_mul(&xr, dsa->priv_key, r, dsa->q, ctx))
  168. goto err; /* s = xr */
  169. if (!BN_add(s, &xr, &m))
  170. goto err; /* s = m + xr */
  171. if (BN_cmp(s, dsa->q) > 0)
  172. if (!BN_sub(s, s, dsa->q))
  173. goto err;
  174. if (!BN_mod_mul(s, s, kinv, dsa->q, ctx))
  175. goto err;
  176. /*
  177. * Redo if r or s is zero as required by FIPS 186-3: this is very
  178. * unlikely.
  179. */
  180. if (BN_is_zero(r) || BN_is_zero(s)) {
  181. if (noredo) {
  182. reason = DSA_R_NEED_NEW_SETUP_VALUES;
  183. goto err;
  184. }
  185. goto redo;
  186. }
  187. ret = DSA_SIG_new();
  188. if (ret == NULL)
  189. goto err;
  190. ret->r = r;
  191. ret->s = s;
  192. err:
  193. if (ret == NULL) {
  194. DSAerr(DSA_F_DSA_DO_SIGN, reason);
  195. BN_free(r);
  196. BN_free(s);
  197. }
  198. if (ctx != NULL)
  199. BN_CTX_free(ctx);
  200. BN_clear_free(&m);
  201. BN_clear_free(&xr);
  202. if (kinv != NULL) /* dsa->kinv is NULL now if we used it */
  203. BN_clear_free(kinv);
  204. return (ret);
  205. }
  206. static int dsa_sign_setup(DSA *dsa, BN_CTX *ctx_in, BIGNUM **kinvp,
  207. BIGNUM **rp)
  208. {
  209. BN_CTX *ctx;
  210. BIGNUM k, kq, *K, *kinv = NULL, *r = NULL;
  211. BIGNUM l, m;
  212. int ret = 0;
  213. int q_bits;
  214. if (!dsa->p || !dsa->q || !dsa->g) {
  215. DSAerr(DSA_F_DSA_SIGN_SETUP, DSA_R_MISSING_PARAMETERS);
  216. return 0;
  217. }
  218. BN_init(&k);
  219. BN_init(&kq);
  220. BN_init(&l);
  221. BN_init(&m);
  222. if (ctx_in == NULL) {
  223. if ((ctx = BN_CTX_new()) == NULL)
  224. goto err;
  225. } else
  226. ctx = ctx_in;
  227. if ((r = BN_new()) == NULL)
  228. goto err;
  229. /* Preallocate space */
  230. q_bits = BN_num_bits(dsa->q);
  231. if (!BN_set_bit(&k, q_bits)
  232. || !BN_set_bit(&l, q_bits)
  233. || !BN_set_bit(&m, q_bits))
  234. goto err;
  235. /* Get random k */
  236. do
  237. if (!BN_rand_range(&k, dsa->q))
  238. goto err;
  239. while (BN_is_zero(&k));
  240. if ((dsa->flags & DSA_FLAG_NO_EXP_CONSTTIME) == 0) {
  241. BN_set_flags(&k, BN_FLG_CONSTTIME);
  242. }
  243. if (dsa->flags & DSA_FLAG_CACHE_MONT_P) {
  244. if (!BN_MONT_CTX_set_locked(&dsa->method_mont_p,
  245. CRYPTO_LOCK_DSA, dsa->p, ctx))
  246. goto err;
  247. }
  248. /* Compute r = (g^k mod p) mod q */
  249. if ((dsa->flags & DSA_FLAG_NO_EXP_CONSTTIME) == 0) {
  250. /*
  251. * We do not want timing information to leak the length of k, so we
  252. * compute G^k using an equivalent scalar of fixed bit-length.
  253. *
  254. * We unconditionally perform both of these additions to prevent a
  255. * small timing information leakage. We then choose the sum that is
  256. * one bit longer than the modulus.
  257. *
  258. * TODO: revisit the BN_copy aiming for a memory access agnostic
  259. * conditional copy.
  260. */
  261. if (!BN_add(&l, &k, dsa->q)
  262. || !BN_add(&m, &l, dsa->q)
  263. || !BN_copy(&kq, BN_num_bits(&l) > q_bits ? &l : &m))
  264. goto err;
  265. BN_set_flags(&kq, BN_FLG_CONSTTIME);
  266. K = &kq;
  267. } else {
  268. K = &k;
  269. }
  270. DSA_BN_MOD_EXP(goto err, dsa, r, dsa->g, K, dsa->p, ctx,
  271. dsa->method_mont_p);
  272. if (!BN_mod(r, r, dsa->q, ctx))
  273. goto err;
  274. /* Compute part of 's = inv(k) (m + xr) mod q' */
  275. if ((kinv = BN_mod_inverse(NULL, &k, dsa->q, ctx)) == NULL)
  276. goto err;
  277. if (*kinvp != NULL)
  278. BN_clear_free(*kinvp);
  279. *kinvp = kinv;
  280. kinv = NULL;
  281. if (*rp != NULL)
  282. BN_clear_free(*rp);
  283. *rp = r;
  284. ret = 1;
  285. err:
  286. if (!ret) {
  287. DSAerr(DSA_F_DSA_SIGN_SETUP, ERR_R_BN_LIB);
  288. if (r != NULL)
  289. BN_clear_free(r);
  290. }
  291. if (ctx_in == NULL)
  292. BN_CTX_free(ctx);
  293. BN_clear_free(&k);
  294. BN_clear_free(&kq);
  295. BN_clear_free(&l);
  296. BN_clear_free(&m);
  297. return ret;
  298. }
  299. static int dsa_do_verify(const unsigned char *dgst, int dgst_len,
  300. DSA_SIG *sig, DSA *dsa)
  301. {
  302. BN_CTX *ctx;
  303. BIGNUM u1, u2, t1;
  304. BN_MONT_CTX *mont = NULL;
  305. int ret = -1, i;
  306. if (!dsa->p || !dsa->q || !dsa->g) {
  307. DSAerr(DSA_F_DSA_DO_VERIFY, DSA_R_MISSING_PARAMETERS);
  308. return -1;
  309. }
  310. i = BN_num_bits(dsa->q);
  311. /* fips 186-3 allows only different sizes for q */
  312. if (i != 160 && i != 224 && i != 256) {
  313. DSAerr(DSA_F_DSA_DO_VERIFY, DSA_R_BAD_Q_VALUE);
  314. return -1;
  315. }
  316. if (BN_num_bits(dsa->p) > OPENSSL_DSA_MAX_MODULUS_BITS) {
  317. DSAerr(DSA_F_DSA_DO_VERIFY, DSA_R_MODULUS_TOO_LARGE);
  318. return -1;
  319. }
  320. BN_init(&u1);
  321. BN_init(&u2);
  322. BN_init(&t1);
  323. if ((ctx = BN_CTX_new()) == NULL)
  324. goto err;
  325. if (BN_is_zero(sig->r) || BN_is_negative(sig->r) ||
  326. BN_ucmp(sig->r, dsa->q) >= 0) {
  327. ret = 0;
  328. goto err;
  329. }
  330. if (BN_is_zero(sig->s) || BN_is_negative(sig->s) ||
  331. BN_ucmp(sig->s, dsa->q) >= 0) {
  332. ret = 0;
  333. goto err;
  334. }
  335. /*
  336. * Calculate W = inv(S) mod Q save W in u2
  337. */
  338. if ((BN_mod_inverse(&u2, sig->s, dsa->q, ctx)) == NULL)
  339. goto err;
  340. /* save M in u1 */
  341. if (dgst_len > (i >> 3))
  342. /*
  343. * if the digest length is greater than the size of q use the
  344. * BN_num_bits(dsa->q) leftmost bits of the digest, see fips 186-3,
  345. * 4.2
  346. */
  347. dgst_len = (i >> 3);
  348. if (BN_bin2bn(dgst, dgst_len, &u1) == NULL)
  349. goto err;
  350. /* u1 = M * w mod q */
  351. if (!BN_mod_mul(&u1, &u1, &u2, dsa->q, ctx))
  352. goto err;
  353. /* u2 = r * w mod q */
  354. if (!BN_mod_mul(&u2, sig->r, &u2, dsa->q, ctx))
  355. goto err;
  356. if (dsa->flags & DSA_FLAG_CACHE_MONT_P) {
  357. mont = BN_MONT_CTX_set_locked(&dsa->method_mont_p,
  358. CRYPTO_LOCK_DSA, dsa->p, ctx);
  359. if (!mont)
  360. goto err;
  361. }
  362. DSA_MOD_EXP(goto err, dsa, &t1, dsa->g, &u1, dsa->pub_key, &u2, dsa->p,
  363. ctx, mont);
  364. /* BN_copy(&u1,&t1); */
  365. /* let u1 = u1 mod q */
  366. if (!BN_mod(&u1, &t1, dsa->q, ctx))
  367. goto err;
  368. /*
  369. * V is now in u1. If the signature is correct, it will be equal to R.
  370. */
  371. ret = (BN_ucmp(&u1, sig->r) == 0);
  372. err:
  373. if (ret < 0)
  374. DSAerr(DSA_F_DSA_DO_VERIFY, ERR_R_BN_LIB);
  375. if (ctx != NULL)
  376. BN_CTX_free(ctx);
  377. BN_free(&u1);
  378. BN_free(&u2);
  379. BN_free(&t1);
  380. return (ret);
  381. }
  382. static int dsa_init(DSA *dsa)
  383. {
  384. dsa->flags |= DSA_FLAG_CACHE_MONT_P;
  385. return (1);
  386. }
  387. static int dsa_finish(DSA *dsa)
  388. {
  389. if (dsa->method_mont_p)
  390. BN_MONT_CTX_free(dsa->method_mont_p);
  391. return (1);
  392. }