矩阵上的
如何在boost::numeric::ublas矩阵上应用1级blas ?例如,我想计算最大项或者所有项的和。在矩阵上使用norm_1或norm_inf不会产生编译错误,但会返回(在我看来)任意结果。我使用boost 1.42
norm_inf
计算由底层向量空间上的无穷范数引起的矩阵范数。这恰好是矩阵的最大绝对行和。
如果你看一下hannes的例子,矩阵(i=99, j=0,…99)的最后一行包含:
9900, 9901, 9902, ... , 9999
如果你把这些元素加起来,你得到994950
,这正是norm_inf
产生的
一个简单的例子如下:
#include<iostream>
#include<boost/numeric/ublas/matrix.hpp>
using namespace boost::numeric::ublas;
int main(){
int l = 100;
matrix<double> m(l,l);
for (int i = 0; i < l; ++i) {
for (int j = 0; j < l; ++j) {
m(i,j)=i*l+j;
}
}
std::cout << norm_inf(m)<<std::endl;
return 0;
}
它应该给出99,但结果是994950。
这是一个至少解决给定任务的一行代码:
float infnorm = accumulate(m.data().begin(),m.data().end(),0,std::max<double>);