Linalg module
The module includes a set of advanced linear algebra operations such as singular value decomposition, and etc.
Currently, Linalg module only supports dense matrix. The support for sparse matrices will be provided very soon.
type mat_d = Owl_dense_matrix_d.mattype mat_z = Owl_dense_matrix_z.mattype ('a, 'b) t = ('a, 'b) Owl_dense_matrix_generic.tA general square matrix A has an LU decomposition into upper and lower triangular matrices, P A = L U where P is a permutation matrix, L is unit lower triangular matrix and U is upper triangular matrix. For square matrices this decomposition can be used to convert the linear system A x = b into a pair of triangular systems (L y = P b, U x = y), which can be solved by forward and back-substitution. Note that the LU decomposition is valid for singular matrices.
val det : mat_d -> floatdet x computes the determinant of a matrix x from its LU decomposition.
qr x calculates QR decomposition for an m by n matrix x as x = Q R. Q is an m by m square matrix (where Q^T Q = I) and R is an m by n right-triangular matrix.
svd x calculates the singular value decomposition of x, and returns a tuple (u,s,v). u is an m by n orthogonal matrix, s an n by 1 matrix of singular values, and v is the transpose of an n by n orthogonal square matrix.
val is_posdef : mat_d -> boolis_posdef x checks whether x is a positive semi-definite matrix.