BOOST SPARSE_MATRIX 不同元素访问方式的时间差异

#include <stdio.h>
#include <stdlib.h>
#include <boost/numeric/ublas/matrix.hpp>
#include <boost/numeric/ublas/matrix_sparse.hpp>
#include <boost/numeric/ublas/io.hpp>
#include <fstream>
#include <boost/format.hpp>
#include <boost/tokenizer.hpp>
#include <boost/algorithm/string.hpp>
#include <vector>
#include <boost/lexical_cast.hpp>
#include <math.h>
#include <boost/algorithm/string.hpp>
#include <ctime>
#include <time.h>
#include <sys/time.h>
#include <memory.h>

#include <algorithm>
#include <cstdlib>

typedef boost::numeric::ublas::compressed_matrix<float>::iterator1 it1_t;
typedef boost::numeric::ublas::compressed_matrix<float>::iterator2 it2_t;
using namespace std;
namespace nub = boost::numeric::ublas;


int main()
{

	nub::compressed_matrix<float> s(1,1000000);
	for(int i = 0; i < s.size2(); i++)
	{
		if(i%10000 == 0)
		  s(0,i) = 99;
		else
		  s(0,i) = 0;
	}

	timeval start, end;
	gettimeofday(&start,NULL);
	for(int i = 0; i < s.size2(); i++)
	{
		if(s(0,i) != 0)
		  s(0,i) = pow(s(0,i),8);
	}
	gettimeofday(&end,NULL);
	cout<<"time is: "<<(end.tv_sec-start.tv_sec)*1000000+(end.tv_usec-start.tv_usec)<<" us"<<endl;


	gettimeofday(&start,NULL);
	for(it1_t t1 = s.begin1(); t1 != s.end1(); t1++)
	  for(it2_t t2 = t1.begin(); t2 != t1.end(); t2++)
	  {
		  if(*t2 != 0)
			*t2 = pow(*t2,8);
	  }
	gettimeofday(&end,NULL);
	cout<<"time is: "<<(end.tv_sec-start.tv_sec)*1000000+(end.tv_usec-start.tv_usec)<<" us"<<endl;
}

g++  -o t t.cpp -O3 -g  -lboost_system

输出的时间分别是：

time is: 536655 us
time is: 218887 us

可以看出，在百万维度上访问的普通的访问时间是相差一倍多