namespace std::linalg {
template<class Real>
setup_givens_rotation_result<Real>
setup_givens_rotation(Real a, Real b) noexcept; // (1)
template<class Real>
setup_givens_rotation_result<complex<Real>>
setup_givens_rotation(complex<Real> a, complex<Real> b) noexcept; // (2)
}
概要
ギブンス回転を計算する。すなわち、以下の式が成り立つような、Real
型の値c
と s
, r
を計算する。
ただし、s
とr
の型はa
とb
の型による。
- (1)
a
とb
の型がReal
の場合、s
とr
の型もReal
。r
は のユークリッドノルム、つまり である。 - (2)
a
とb
の型がcomplex<Real>
の場合、s
とr
の型もcomplex<Real>
。以下で定義される 関数を用いると、r
は である。
適格要件
Real
はcomplex<Real>
が規定できる型であること。
戻り値
c
とs
が入力r
を入力{c, s, r}
である。
例
[注意] 処理系にあるコンパイラで確認していないため、間違っているかもしれません。
#include <cmath>
#include <complex>
#include <iostream>
#include <linalg>
template <class T>
void print(const std::linalg::setup_givens_rotation_result<T>& result) {
std::cout << "c: " << result.c << '\n'
<< "s: " << result.s << '\n'
<< "r: " << result.r << '\n';
}
int main()
{
// (1)
std::cout << "(1)\n";
auto result1 = std::linalg::setup_givens_rotation(1.0, std::sqrt(3.0));
print(result1);
// (2)
std::cout << "(2)\n";
auto result2 = std::linalg::setup_givens_rotation(std::complex<double>(1.0, 0), std::complex<double>(0, std::sqrt(3.0)));
print(result2);
return 0;
}
27
auto result2 = std::linalg::setup_givens_rotation(std::complex<double>(1.0, 0), std::complex<double>(0, std::sqrt(3.0)));
#include <cmath>
#include <complex>
#include <iostream>
#include <linalg>
template <class T>
void print(const std::linalg::setup_givens_rotation_result<T>& result) {
std::cout << "c: " << result.c << '\n'
<< "s: " << result.s << '\n'
<< "r: " << result.r << '\n';
}
int main()
{
// (1)
std::cout << "(1)\n";
auto result1 = std::linalg::setup_givens_rotation(1.0, std::sqrt(3.0));
print(result1);
出力
(1)
c: 0.5
s: -0.866025
r: 2.0
(2)
c: 0.5
s: (0,-0.866025)
r: (2,0)
バージョン
言語
- C++26
処理系
- Clang: ??
- GCC: ??
- ICC: ??
- Visual C++: ??