1 1.1 christos #include <openssl/bn.h> 2 1.1 christos #include <openssl/rand.h> 3 1.1 christos 4 1.1 christos static int Rand(void) 5 1.1 christos { 6 1.1 christos unsigned char x[2]; 7 1.1 christos RAND_bytes(x, 2); 8 1.1 christos return (x[0] + 2 * x[1]); 9 1.1 christos } 10 1.1 christos 11 1.1 christos static void bug(const char *m, BIGNUM *a, BIGNUM *b) 12 1.1 christos { 13 1.1 christos printf("%s!\na=", m); 14 1.1 christos BN_print_fp(stdout, a); 15 1.1 christos printf("\nb="); 16 1.1 christos BN_print_fp(stdout, b); 17 1.1 christos printf("\n"); 18 1.1 christos fflush(stdout); 19 1.1 christos exit(1); 20 1.1 christos } 21 1.1 christos 22 1.1 christos int 23 1.1 christos main(int argc, char *argv[]) 24 1.1 christos { 25 1.1 christos BIGNUM *a = BN_new(), *b = BN_new(), *c = BN_new(), *d = BN_new(), 26 1.1 christos *C = BN_new(), *D = BN_new(); 27 1.1 christos BN_RECP_CTX *recp = BN_RECP_CTX_new(); 28 1.1 christos BN_CTX *ctx = BN_CTX_new(); 29 1.1 christos int i = 0; 30 1.1 christos 31 1.1 christos for (i = 0; i < 10000; i++) { 32 1.2 rillig BN_rand(a, Rand(), 0, 0); 33 1.2 rillig BN_rand(b, Rand(), 0, 0); 34 1.1 christos if (BN_is_zero(b)) 35 1.1 christos continue; 36 1.1 christos 37 1.1 christos BN_RECP_CTX_set(recp, b, ctx); 38 1.1 christos if (BN_div(C, D, a, b, ctx) != 1) 39 1.1 christos bug("BN_div failed", a, b); 40 1.1 christos if (BN_div_recp(c, d, a, recp, ctx) != 1) 41 1.1 christos bug("BN_div_recp failed", a, b); 42 1.1 christos else if (BN_cmp(c, C) != 0 || BN_cmp(c, C) != 0) 43 1.1 christos bug("mismatch", a, b); 44 1.1 christos } 45 1.1 christos exit(0); 46 1.1 christos } 47