Basic linear algebra algorithms (since C++26)

基础线性代数算法基于密集基础线性代数子程序（ BLAS ），其对应于 BLAS标准 的一个子集。这些访问数组元素的算法通过表示向量或矩阵的 std::mdspan 来访问这些元素。

BLAS算法被划分为三组称为 级别 的操作集，通常对应算法复杂度中多项式的次数：

• BLAS 1 : 所有带有 std::mdspan 参数的算法，其 std::mdspan 数组访问次数和算术运算次数的量级，与任意 std::mdspan 参数各维度大小的最大乘积呈 线性 关系。这些算法包含点积、范数和向量加法等 向量 操作。
• BLAS 2 : 所有算法具有 二次 时间的一般复杂度。这些算法包含矩阵-向量乘法及三角线性方程组求解等 矩阵-向量 操作。
• BLAS 3 : 所有算法具有 三次 时间的一般复杂度。这些算法包含矩阵-矩阵乘法及多重三角线性方程组求解等 矩阵-矩阵 操作。

注释

示例

#include <cassert>
#include <cstddef>
#include <execution>
#include <linalg>
#include <mdspan>
#include <numeric>
#include <vector>
int main()
{
    std::vector<double> x_vec(42);
    std::ranges::iota(x_vec, 0.0);
    std::mdspan x(x_vec.data(), x_vec.size());
    // x[i] *= 2.0，顺序执行
    std::linalg::scale(2.0, x);
    // x[i] *= 3.0，并行执行
    std::linalg::scale(std::execution::par_unseq, 3.0, x);
    for (std::size_t i{}; i != x.size(); ++i)
        assert(x[i] == 6.0 * static_cast<double>(i));
}

外部链接