boost compressed_matrix retrieve style is not random

23 篇文章 0 订阅
20 篇文章 0 订阅
typedef boost::numeric::ublas::compressed_matrix<float>::iterator1 it1_t;
typedef boost::numeric::ublas::compressed_matrix<float>::iterator2 it2_t;
<pre name="code" class="cpp">boost::numeric::ublas::compressed_matrix<float> At(dim,dim);
my_time time;
time.tic();
int count = 0;
int index = 0;
for(it1_t it1=diff_x.begin1(); it1 != diff_x.end1(); it1++)
 for(it2_t it2 = it1.begin(); it2 != it1.end(); it2++)
  {
	count++;
        At(it2.index2(),it2.index2()) -= label_diff*alpha*pow(*it2,2);
	At(it2.index2(),it2.index2()) = MAX(0.0,At(it2.index2(),it2.index2()));
  }
cout<<"count is: "<<count<<endl;
time.toc();
 


the code is change the value in diagonal of matrix,  we test two different sparse data, the output show below:

dim is: 29890095
count is: 54
time is: 68229 us

dim is: 180
count is: 79
time is: 18 us

We can get that, the element we retrieve is in a same magnitude, but the time has a difference over several magnitude.


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值