1/*
2 * Copyright (c) 2000-2001,2011,2014 Apple Inc. All Rights Reserved.
3 *
4 * The contents of this file constitute Original Code as defined in and are
5 * subject to the Apple Public Source License Version 1.2 (the 'License').
6 * You may not use this file except in compliance with the License. Please obtain
7 * a copy of the License at http://www.apple.com/publicsource and read it before
8 * using this file.
9 *
10 * This Original Code and all software distributed under the License are
11 * distributed on an 'AS IS' basis, WITHOUT WARRANTY OF ANY KIND, EITHER EXPRESS
12 * OR IMPLIED, AND APPLE HEREBY DISCLAIMS ALL SUCH WARRANTIES, INCLUDING WITHOUT
13 * LIMITATION, ANY WARRANTIES OF MERCHANTABILITY, FITNESS FOR A PARTICULAR
14 * PURPOSE, QUIET ENJOYMENT OR NON-INFRINGEMENT. Please see the License for the
15 * specific language governing rights and limitations under the License.
16 */
17
18
19/* unused */
20
21#include <stdio.h>
22#include <openssl/tmdiff.h>
23#include "bn_lcl.h"
24
25#define SIZE	256
26#define NUM	(8*8*8)
27#define MOD	(8*8*8*8*8)
28
29main(argc,argv)
30int argc;
31char *argv[];
32	{
33	BN_CTX ctx;
34	BIGNUM a,b,c,r,rr,t,l;
35	int j,i,size=SIZE,num=NUM,mod=MOD;
36	char *start,*end;
37	BN_MONT_CTX mont;
38	double d,md;
39
40	BN_MONT_CTX_init(&mont);
41	BN_CTX_init(&ctx);
42	BN_init(&a);
43	BN_init(&b);
44	BN_init(&c);
45	BN_init(&r);
46
47	start=ms_time_new();
48	end=ms_time_new();
49	while (size <= 1024*8)
50		{
51		BN_rand(&a,size,0,0);
52		BN_rand(&b,size,1,0);
53		BN_rand(&c,size,0,1);
54
55		BN_mod(&a,&a,&c,&ctx);
56
57		ms_time_get(start);
58		for (i=0; i<10; i++)
59			BN_MONT_CTX_set(&mont,&c,&ctx);
60		ms_time_get(end);
61		md=ms_time_diff(start,end);
62
63		ms_time_get(start);
64		for (i=0; i<num; i++)
65			{
66			/* bn_mull(&r,&a,&b,&ctx); */
67			/* BN_sqr(&r,&a,&ctx); */
68			BN_mod_exp_mont(&r,&a,&b,&c,&ctx,&mont);
69			}
70		ms_time_get(end);
71		d=ms_time_diff(start,end)/* *50/33 */;
72		printf("%5d bit:%6.2f %6d %6.4f %4d m_set(%5.4f)\n",size,
73			d,num,d/num,(int)((d/num)*mod),md/10.0);
74		num/=8;
75		mod/=8;
76		if (num <= 0) num=1;
77		size*=2;
78		}
79
80	}
81