QR行列分解

詳細説明

\(M\times N\)( \(M>N\)) 行列 \(A = Q*R\) の QR 分解を実行し、行列方程式 \(A*X=B\) を解く。

引数

src1	行優先順で格納された入力行列 \(A\) へのポインタ。処理完了後、src1 には上三角の \(N\times N\) 行列 \(R\) が格納される。src1 の下三角部分には基本反射行列のベクトルが格納される。詳細は [290] および Lapack の DGEQRF のドキュメントを参照すること。
src1_step	行列 \(A\) の連続する2行間のバイト数。
m	行列 \(A\) の行数。
n	行列 \(A\) の列数。
k	\(M\times K\) 行列 \(B\) における右辺ベクトルの数。
src2	連立方程式 \(A*X=B\) の右辺となる \(M\times K\) 行列 \(B\) へのポインタ。\(B\) は行優先順で格納される。src2 が null ポインタの場合は QR 分解のみが実行される。そうでない場合は連立方程式が解かれ、src1 が一時バッファとして使用されるため、処理完了後 src2 には連立方程式 \(A*X=B\) の解 \(X\) が格納される。
src2_step	行列 \(B\) の連続する2つの行の間のバイト数。
dst	基本反射行列のスカラー係数を格納する連続した \(N\times 1\) 配列へのポインタ。詳細は [290] を参照すること。
info	分解の成否を示す。*info が 0 の場合、分解は失敗である。

関数
int	hal_ni_QR32f (float *src1, size_t src1_step, int m, int n, int k, float *src2, size_t src2_step, float *dst, int *info)
int	hal_ni_QR64f (double *src1, size_t src1_step, int m, int n, int k, double *src2, size_t src2_step, double *dst, int *info)

関数詳解

hal_ni_QR32f()

int hal_ni_QR32f ( float * src1, size_t src1_step, int m, int n, int k, float * src2, size_t src2_step, float * dst, int * info )

inline

#include <core/src/hal_replacement.hpp>

hal_ni_QR64f()

int hal_ni_QR64f ( double * src1, size_t src1_step, int m, int n, int k, double * src2, size_t src2_step, double * dst, int * info )

inline

#include <core/src/hal_replacement.hpp>