C++で複素数の射影を計算するproj()関数の使い方を徹底解説

C++のproj()関数とは

この記事では、C++において複素数に対する射影（プロジェクション）を計算するための proj() 関数の動作について詳しく解説します。proj() 関数は、与えられた複素数をリーマン球面へ射影した結果を返す標準ライブラリの機能であり、特に無限大を含む複素数を扱う際に重要な役割を果たします。

まずは、proj() メソッドの基本構文を見ていきましょう。

template <class T> complex<T>
proj (const complex<T>& z);

proj()関数の基本構文と動作

proj() 関数は、引数として複素数を表す値を受け取り、その射影結果を complex<T> 型として返します。通常の有限の複素数が渡された場合は、入力された値がそのまま返されます。一方、無限大（INFINITY）を含む複素数が渡された場合には、実部が無限大、虚部が符号付きゼロ（±0）である複素数が返されるという特別な挙動になります。

使用例

以下のコード例では、proj() 関数を使用して、通常の複素数および無限大を含む複素数に対する射影を計算しています。

#include <iostream>
#include <complex>
using namespace std;
int main(){
    std::complex<double> c1(3, 5);
    cout << "Proj" << c1 << " = " << proj(c1) << endl;
    std::complex<double> c2(0, -INFINITY);
    cout << "Proj" << c2 << " = " << proj(c2) << endl;
    std::complex<double> c3(INFINITY, -1);
    cout << "Proj" << c3 << " = " << proj(c3) << endl;
}

proj() 関数の定義を利用するためには、ソースコード内で <complex> ヘッダーをインクルードすることが必須です。これにより、複素数型や射影関数などの実装が利用可能になります。

上記のコードを正常にコンパイルして実行すると、渡されたそれぞれの複素数に対して、以下のような結果が出力されます。

実行結果

Proj(3,5) = (3,5)
Proj(0,-inf) = (inf,-0)
Proj(inf,-1) = (inf,-0)

出力結果のポイント

実行結果から分かるように、有限の複素数 (3, 5) に対しては元の値がそのまま返されています。一方、虚部が負の無限大である (0, -inf) や、実部が正の無限大である (inf, -1) のような無限大を含む複素数に対しては、いずれも (inf, -0) へと射影されています。これは、リーマン球面モデルにおいて、あらゆる方向の無限遠点が同一の「無限大の点」として扱われるためです。proj() 関数はこの性質を利用しており、無限大を含む複素数の計算を一貫性のある形で処理できるようになっています。