SM2算法第二十六篇：openSSL 中BIGNUM的使用记录

出处：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;

}