d=ms_time_diff(start,end) *50/33 ;
d=ms_time_diff(start,end);
main(argc,argv)int argc;char *argv[];	{	BN_CTX ctx;	BIGNUM a,b,c,r,rr,t,l;	int j,i,size=SIZE,num=NUM,mod=MOD;	char *start,*end;	BN_MONT_CTX mont;	double d,md;	BN_MONT_CTX_init(&mont);	BN_CTX_init(&ctx);	BN_init(&a);	BN_init(&b);	BN_init(&c);	BN_init(&r);	start=ms_time_new();	end=ms_time_new();	while (size <= 1024*8)		{		BN_rand(&a,size,0,0);		BN_rand(&b,size,1,0);		BN_rand(&c,size,0,1);		BN_mod(&a,&a,&c,&ctx);		ms_time_get(start);		for (i=0; i<10; i++)			BN_MONT_CTX_set(&mont,&c,&ctx);		ms_time_get(end);		md=ms_time_diff(start,end);		ms_time_get(start);		for (i=0; i<num; i++)			{			//bn_mull(&r,&a,&b,&ctx);			//BN_sqr(&r,&a,&ctx);			BN_mod_exp_mont(&r,&a,&b,&c,&ctx,&mont);			}		ms_time_get(end);		d=ms_time_diff(start,end) *50/33 /**/;		printf("%5d bit:%6.2f %6d %6.4f %4d m_set(%5.4f)\n",size,			d,num,d/num,(int)((d/num)*mod),md/10.0);		num/=8;		mod/=8;		if (num <= 0) num=1;		size*=2;		}	}