将Boost UBLAS blas-1应用于矩阵

Applying Boost UBLAS blas-1 to Matrices

本文关键字：应用于 blas-1 Boost UBLAS      更新时间：2023-10-16

如何在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>);
相关文章：