exp.c

Go to the documentation of this file.

/* unused */

#include <stdio.h>
#include <openssl/tmdiff.h>
#include "bn_lcl.h"

#define SIZE    256
#define NUM (8*8*8)
#define MOD (8*8*8*8*8)

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;
        }

    }