- [mathjax enable]
- linalg[meta header]
- function template[meta id-type]
- std::linalg[meta namespace]
- cpp26[meta cpp]
namespace std::linalg {
template<in-object InObj1,
in-object InObj2,
out-object OutObj>
void add(InObj1 x,
InObj2 y,
OutObj z); // (1)
template<class ExecutionPolicy,
in-object InObj1,
in-object InObj2,
out-object OutObj>
void add(ExecutionPolicy&& exec,
InObj1 x,
InObj2 y,
OutObj z); // (2)
}- in-object[link inout-object.md]
- out-object[link inout-object.md]
同じサイズの行列またはベクトルのxとy、zに対して、$x + y$ をzに代入する:
xとy、zの次元が全て等しくなければならない。x.rank() == y.rank() && y.rank() == z.rank()
- (1), (2):
possibly-addable<InObj1, InObj2, OutObj>()がtrue(行列またはベクトルのx、y、zの各次元の静的要素数が同じ) - (2):
is_execution_policy<ExecutionPolicy>::valueがtrue
行列またはベクトルのx、y、zの各次元が同じであること。
addable(x, y, z) == true
zに代入する。
- (1): 逐次実行する。
- (2): 指定された実行ポリシーに応じて実行する。
なし
zをxまたはyとしてもよい。
[注意] 処理系にあるコンパイラで確認していないため、間違っているかもしれません。
#include <cmath>
#include <execution>
#include <iostream>
#include <linalg>
#include <mdspan>
#include <vector>
template <class Vector>
void print(Vector v) {
for (int i = 0; i < v.extent(0) - 1; ++i) {
std::cout << v[i] << ", ";
}
std::cout << v[v.extent(0) - 1] << std::endl;
}
int main()
{
constexpr size_t N = 3;
std::vector<double> a_vec({1, 2, 3});
std::mdspan a(a_vec.data(), N);
std::vector<double> b_vec({4, 5, 6});
std::mdspan b(b_vec.data(), N);
std::vector<double> c_vec(N);
std::mdspan c(c_vec.data(), N);
// (1)
std::linalg::add(a, b, c);
print(c);
// (2)
std::linalg::add(std::execution::par, a, b, c);
print(c);
return 0;
}- std::linalg::add[color ff0000]
5, 7, 9
5, 7, 9
- C++26
- Clang: ??
- GCC: ??
- ICC: ??
- Visual C++: ??