dsa_ossl.c 11 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416
  1. /* crypto/dsa/dsa_ossl.c */
  2. /* Copyright (C) 1995-1998 Eric Young (eay@cryptsoft.com)
  3. * All rights reserved.
  4. *
  5. * This package is an SSL implementation written
  6. * by Eric Young (eay@cryptsoft.com).
  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 (tjh@cryptsoft.com).
  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 (eay@cryptsoft.com)"
  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 (tjh@cryptsoft.com)"
  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 <schoch@sheba.arc.nasa.gov> */
  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, BIGNUM **rp);
  68. static int dsa_do_verify(const unsigned char *dgst, int dgst_len, DSA_SIG *sig,
  69. DSA *dsa);
  70. static int dsa_init(DSA *dsa);
  71. static int dsa_finish(DSA *dsa);
  72. static DSA_METHOD openssl_dsa_meth = {
  73. "OpenSSL DSA method",
  74. dsa_do_sign,
  75. dsa_sign_setup,
  76. dsa_do_verify,
  77. NULL, /* dsa_mod_exp, */
  78. NULL, /* dsa_bn_mod_exp, */
  79. dsa_init,
  80. dsa_finish,
  81. 0,
  82. NULL,
  83. NULL,
  84. NULL
  85. };
  86. /* These macro wrappers replace attempts to use the dsa_mod_exp() and
  87. * bn_mod_exp() handlers in the DSA_METHOD structure. We avoid the problem of
  88. * having a the macro work as an expression by bundling an "err_instr". So;
  89. *
  90. * if (!dsa->meth->bn_mod_exp(dsa, r,dsa->g,&k,dsa->p,ctx,
  91. * dsa->method_mont_p)) goto err;
  92. *
  93. * can be replaced by;
  94. *
  95. * DSA_BN_MOD_EXP(goto err, dsa, r, dsa->g, &k, dsa->p, ctx,
  96. * dsa->method_mont_p);
  97. */
  98. #define DSA_MOD_EXP(err_instr,dsa,rr,a1,p1,a2,p2,m,ctx,in_mont) \
  99. do { \
  100. int _tmp_res53; \
  101. if((dsa)->meth->dsa_mod_exp) \
  102. _tmp_res53 = (dsa)->meth->dsa_mod_exp((dsa), (rr), (a1), (p1), \
  103. (a2), (p2), (m), (ctx), (in_mont)); \
  104. else \
  105. _tmp_res53 = BN_mod_exp2_mont((rr), (a1), (p1), (a2), (p2), \
  106. (m), (ctx), (in_mont)); \
  107. if(!_tmp_res53) err_instr; \
  108. } while(0)
  109. #define DSA_BN_MOD_EXP(err_instr,dsa,r,a,p,m,ctx,m_ctx) \
  110. do { \
  111. int _tmp_res53; \
  112. if((dsa)->meth->bn_mod_exp) \
  113. _tmp_res53 = (dsa)->meth->bn_mod_exp((dsa), (r), (a), (p), \
  114. (m), (ctx), (m_ctx)); \
  115. else \
  116. _tmp_res53 = BN_mod_exp_mont((r), (a), (p), (m), (ctx), (m_ctx)); \
  117. if(!_tmp_res53) err_instr; \
  118. } while(0)
  119. const DSA_METHOD *DSA_OpenSSL(void)
  120. {
  121. return &openssl_dsa_meth;
  122. }
  123. static DSA_SIG *dsa_do_sign(const unsigned char *dgst, int dlen, DSA *dsa)
  124. {
  125. BIGNUM *kinv=NULL,*r=NULL,*s=NULL;
  126. BIGNUM m;
  127. BIGNUM xr;
  128. BN_CTX *ctx=NULL;
  129. int reason=ERR_R_BN_LIB;
  130. DSA_SIG *ret=NULL;
  131. BN_init(&m);
  132. BN_init(&xr);
  133. if (!dsa->p || !dsa->q || !dsa->g)
  134. {
  135. reason=DSA_R_MISSING_PARAMETERS;
  136. goto err;
  137. }
  138. s=BN_new();
  139. if (s == NULL) goto err;
  140. /* reject a excessive digest length (currently at most
  141. * dsa-with-SHA256 is supported) */
  142. if (dlen > SHA256_DIGEST_LENGTH)
  143. {
  144. reason=DSA_R_DATA_TOO_LARGE_FOR_KEY_SIZE;
  145. goto err;
  146. }
  147. ctx=BN_CTX_new();
  148. if (ctx == NULL) goto err;
  149. if ((dsa->kinv == NULL) || (dsa->r == NULL))
  150. {
  151. if (!DSA_sign_setup(dsa,ctx,&kinv,&r)) goto err;
  152. }
  153. else
  154. {
  155. kinv=dsa->kinv;
  156. dsa->kinv=NULL;
  157. r=dsa->r;
  158. dsa->r=NULL;
  159. }
  160. if (dlen > BN_num_bytes(dsa->q))
  161. /* if the digest length is greater than the size of q use the
  162. * BN_num_bits(dsa->q) leftmost bits of the digest, see
  163. * fips 186-3, 4.2 */
  164. dlen = BN_num_bytes(dsa->q);
  165. if (BN_bin2bn(dgst,dlen,&m) == NULL)
  166. goto err;
  167. /* Compute s = inv(k) (m + xr) mod q */
  168. if (!BN_mod_mul(&xr,dsa->priv_key,r,dsa->q,ctx)) goto err;/* s = xr */
  169. if (!BN_add(s, &xr, &m)) goto err; /* s = m + xr */
  170. if (BN_cmp(s,dsa->q) > 0)
  171. BN_sub(s,s,dsa->q);
  172. if (!BN_mod_mul(s,s,kinv,dsa->q,ctx)) goto err;
  173. ret=DSA_SIG_new();
  174. if (ret == NULL) goto err;
  175. ret->r = r;
  176. ret->s = s;
  177. err:
  178. if (!ret)
  179. {
  180. DSAerr(DSA_F_DSA_DO_SIGN,reason);
  181. BN_free(r);
  182. BN_free(s);
  183. }
  184. if (ctx != NULL) BN_CTX_free(ctx);
  185. BN_clear_free(&m);
  186. BN_clear_free(&xr);
  187. if (kinv != NULL) /* dsa->kinv is NULL now if we used it */
  188. BN_clear_free(kinv);
  189. return(ret);
  190. }
  191. static int dsa_sign_setup(DSA *dsa, BN_CTX *ctx_in, BIGNUM **kinvp, BIGNUM **rp)
  192. {
  193. BN_CTX *ctx;
  194. BIGNUM k,kq,*K,*kinv=NULL,*r=NULL;
  195. int ret=0;
  196. if (!dsa->p || !dsa->q || !dsa->g)
  197. {
  198. DSAerr(DSA_F_DSA_SIGN_SETUP,DSA_R_MISSING_PARAMETERS);
  199. return 0;
  200. }
  201. BN_init(&k);
  202. BN_init(&kq);
  203. if (ctx_in == NULL)
  204. {
  205. if ((ctx=BN_CTX_new()) == NULL) goto err;
  206. }
  207. else
  208. ctx=ctx_in;
  209. if ((r=BN_new()) == NULL) goto err;
  210. /* Get random k */
  211. do
  212. if (!BN_rand_range(&k, dsa->q)) goto err;
  213. while (BN_is_zero(&k));
  214. if ((dsa->flags & DSA_FLAG_NO_EXP_CONSTTIME) == 0)
  215. {
  216. BN_set_flags(&k, BN_FLG_CONSTTIME);
  217. }
  218. if (dsa->flags & DSA_FLAG_CACHE_MONT_P)
  219. {
  220. if (!BN_MONT_CTX_set_locked(&dsa->method_mont_p,
  221. CRYPTO_LOCK_DSA,
  222. dsa->p, ctx))
  223. goto err;
  224. }
  225. /* Compute r = (g^k mod p) mod q */
  226. if ((dsa->flags & DSA_FLAG_NO_EXP_CONSTTIME) == 0)
  227. {
  228. if (!BN_copy(&kq, &k)) goto err;
  229. /* We do not want timing information to leak the length of k,
  230. * so we compute g^k using an equivalent exponent of fixed length.
  231. *
  232. * (This is a kludge that we need because the BN_mod_exp_mont()
  233. * does not let us specify the desired timing behaviour.) */
  234. if (!BN_add(&kq, &kq, dsa->q)) goto err;
  235. if (BN_num_bits(&kq) <= BN_num_bits(dsa->q))
  236. {
  237. if (!BN_add(&kq, &kq, dsa->q)) goto err;
  238. }
  239. K = &kq;
  240. }
  241. else
  242. {
  243. K = &k;
  244. }
  245. DSA_BN_MOD_EXP(goto err, dsa, r, dsa->g, K, dsa->p, ctx,
  246. dsa->method_mont_p);
  247. if (!BN_mod(r,r,dsa->q,ctx)) goto err;
  248. /* Compute part of 's = inv(k) (m + xr) mod q' */
  249. if ((kinv=BN_mod_inverse(NULL,&k,dsa->q,ctx)) == NULL) goto err;
  250. if (*kinvp != NULL) BN_clear_free(*kinvp);
  251. *kinvp=kinv;
  252. kinv=NULL;
  253. if (*rp != NULL) BN_clear_free(*rp);
  254. *rp=r;
  255. ret=1;
  256. err:
  257. if (!ret)
  258. {
  259. DSAerr(DSA_F_DSA_SIGN_SETUP,ERR_R_BN_LIB);
  260. if (r != NULL)
  261. BN_clear_free(r);
  262. }
  263. if (ctx_in == NULL) BN_CTX_free(ctx);
  264. BN_clear_free(&k);
  265. BN_clear_free(&kq);
  266. return(ret);
  267. }
  268. static int dsa_do_verify(const unsigned char *dgst, int dgst_len, DSA_SIG *sig,
  269. DSA *dsa)
  270. {
  271. BN_CTX *ctx;
  272. BIGNUM u1,u2,t1;
  273. BN_MONT_CTX *mont=NULL;
  274. int ret = -1, i;
  275. if (!dsa->p || !dsa->q || !dsa->g)
  276. {
  277. DSAerr(DSA_F_DSA_DO_VERIFY,DSA_R_MISSING_PARAMETERS);
  278. return -1;
  279. }
  280. i = BN_num_bits(dsa->q);
  281. /* fips 186-3 allows only different sizes for q */
  282. if (i != 160 && i != 224 && i != 256)
  283. {
  284. DSAerr(DSA_F_DSA_DO_VERIFY,DSA_R_BAD_Q_VALUE);
  285. return -1;
  286. }
  287. if (BN_num_bits(dsa->p) > OPENSSL_DSA_MAX_MODULUS_BITS)
  288. {
  289. DSAerr(DSA_F_DSA_DO_VERIFY,DSA_R_MODULUS_TOO_LARGE);
  290. return -1;
  291. }
  292. /* reject a excessive digest length (currently at most
  293. * dsa-with-SHA256 is supported) */
  294. if (dgst_len > SHA256_DIGEST_LENGTH)
  295. {
  296. DSAerr(DSA_F_DSA_DO_VERIFY,DSA_R_DATA_TOO_LARGE_FOR_KEY_SIZE);
  297. return -1;
  298. }
  299. BN_init(&u1);
  300. BN_init(&u2);
  301. BN_init(&t1);
  302. if ((ctx=BN_CTX_new()) == NULL) goto err;
  303. if (BN_is_zero(sig->r) || BN_is_negative(sig->r) ||
  304. BN_ucmp(sig->r, dsa->q) >= 0)
  305. {
  306. ret = 0;
  307. goto err;
  308. }
  309. if (BN_is_zero(sig->s) || BN_is_negative(sig->s) ||
  310. BN_ucmp(sig->s, dsa->q) >= 0)
  311. {
  312. ret = 0;
  313. goto err;
  314. }
  315. /* Calculate W = inv(S) mod Q
  316. * save W in u2 */
  317. if ((BN_mod_inverse(&u2,sig->s,dsa->q,ctx)) == NULL) goto err;
  318. /* save M in u1 */
  319. if (dgst_len > (i >> 3))
  320. /* if the digest length is greater than the size of q use the
  321. * BN_num_bits(dsa->q) leftmost bits of the digest, see
  322. * fips 186-3, 4.2 */
  323. dgst_len = (i >> 3);
  324. if (BN_bin2bn(dgst,dgst_len,&u1) == NULL) goto err;
  325. /* u1 = M * w mod q */
  326. if (!BN_mod_mul(&u1,&u1,&u2,dsa->q,ctx)) goto err;
  327. /* u2 = r * w mod q */
  328. if (!BN_mod_mul(&u2,sig->r,&u2,dsa->q,ctx)) goto err;
  329. if (dsa->flags & DSA_FLAG_CACHE_MONT_P)
  330. {
  331. mont = BN_MONT_CTX_set_locked(&dsa->method_mont_p,
  332. CRYPTO_LOCK_DSA, dsa->p, ctx);
  333. if (!mont)
  334. goto err;
  335. }
  336. DSA_MOD_EXP(goto err, dsa, &t1, dsa->g, &u1, dsa->pub_key, &u2, dsa->p, ctx, mont);
  337. /* BN_copy(&u1,&t1); */
  338. /* let u1 = u1 mod q */
  339. if (!BN_mod(&u1,&t1,dsa->q,ctx)) goto err;
  340. /* V is now in u1. If the signature is correct, it will be
  341. * equal to R. */
  342. ret=(BN_ucmp(&u1, sig->r) == 0);
  343. err:
  344. /* XXX: surely this is wrong - if ret is 0, it just didn't verify;
  345. there is no error in BN. Test should be ret == -1 (Ben) */
  346. if (ret != 1) DSAerr(DSA_F_DSA_DO_VERIFY,ERR_R_BN_LIB);
  347. if (ctx != NULL) BN_CTX_free(ctx);
  348. BN_free(&u1);
  349. BN_free(&u2);
  350. BN_free(&t1);
  351. return(ret);
  352. }
  353. static int dsa_init(DSA *dsa)
  354. {
  355. dsa->flags|=DSA_FLAG_CACHE_MONT_P;
  356. return(1);
  357. }
  358. static int dsa_finish(DSA *dsa)
  359. {
  360. if(dsa->method_mont_p)
  361. BN_MONT_CTX_free(dsa->method_mont_p);
  362. return(1);
  363. }