部分 (行) ピボット付きの LU 分解で、x = P L U です。R の Matrix::lu() と同等です。因子は Matrix::expand(Matrix::lu(x)) のもの (LAPACK の dgetrf) です。NaN または無限大の値があると例外を投げます (Matrix::lu は、NaN または無限大の値を含む因子を返します)。
x = P L U
Matrix::lu()
Matrix::expand(Matrix::lu(x))
NaN
Matrix::lu
lu([[1, 2], [3, 4]]);// { P: [[0, 1], [1, 0]], L: [[1, 0], [0.3333333333333333, 1]], U: [[3, 4], [0, 0.6666666666666667]] } コピー
lu([[1, 2], [3, 4]]);// { P: [[0, 1], [1, 0]], L: [[1, 0], [0.3333333333333333, 1]], U: [[3, 4], [0, 0.6666666666666667]] }
部分 (行) ピボット付きの LU 分解で、
x = P L Uです。R のMatrix::lu()と同等です。因子はMatrix::expand(Matrix::lu(x))のもの (LAPACK の dgetrf) です。NaNまたは無限大の値があると例外を投げます (Matrix::luは、NaN または無限大の値を含む因子を返します)。