出处:http://blog.163.com/xiaoting_hu/blog/static/5046477220136215847597/
转载原因:对于我理解(1)BN_CTX_get() (2)BN_set_word() 有很大帮助
功能:求两个向量的内积
第一种方法
int uInnerProduct(int dime, BIGNUM *r,BIGNUM a[],BIGNUM b[])
{
BIGNUM *temp;
BN_CTX *ctx = NULL; //定义一个BN_CTX类型变量
ctx = BN_CTX_new();
BN_CTX_start(ctx); // 如果有大数函数(如BN-mul())中有BN_CTX变量,也需要在加入这三句。使用完用后面的end和free释放。
temp=BN_CTX_get(ctx); //给Temp分配一个空间
BN_set_word(r,0); //给r所指向的变量赋值为0
//计算内积
for(int i=0;i<dime;i++)
{
BN_mul(temp,&a[i],&b[i],ctx);
BN_add(r,r,temp);
}
BN_CTX_end(ctx);
BN_CTX_free(ctx); //配合上面的new和start
return 0;
}
int _tmain(int argc, _TCHAR* argv[])
{
BIGNUM a[4];
BIGNUM b[4];
BIGNUM *ret;
ret=BN_new(); //和BN_free()配对
int dime=4;
for(int i=0;i<dime;i++)
{
BN_init(&a[i]); //如果不是指针定义方法,则必须在使用该变量前对其初始化,类似于BN-new()
BN_init(&b[i]);
BN_rand(&a[i],512,1,0);
BN_rand(&b[i],512,0,1); //但每次产生的都是一样的
BN_print_fp(stdout,&a[i]);
printf("\n");
BN_print_fp(stdout,&b[i]);
printf("\n");
}
uInnerProduct(dime,ret,a,b);
for(int i=0;i<dime;i++)
{
BN_free(&a[i]); //使用完必须释放
BN_free(&b[i]);
}
BN_print_fp(stdout,ret);
BN_free(ret);
system("pause");
return 0;
}