c++ - 将 Boost UBLAS blas-1 应用于矩阵

Question

如何在 boost::numeric::ublas 矩阵上应用 1 级 blas？例如，我想计算最大条目或所有条目的总和。在矩阵上使用 norm_1 或 norm_inf 不会产生编译器错误，但会返回(在我看来)任意结果。我正在使用 boost 1.42

Answer 1

矩阵上的

norm_inf 计算矩阵范数 induced通过底层 vector 空间的无穷范数。这恰好是矩阵的最大绝对行总和。

如果您查看 hannes 的示例，矩阵的最后一行 (i=99, j=0...99) 包含:

9900, 9901, 9902, ... , 9999

如果对这些条目求和，您将得到 994950，这正是 norm_inf 产生的结果。