cvx-linalg

Linear algebra utilities for portfolio optimization, part of the jebel-quant ecosystem.

Installation

pip install cvx-linalg

Usage

from cvx.linalg import (
    a_norm, cholesky,
    inv, inv_a_norm, is_positive_definite, lstsq, pca, rand_cov, solve, valid,
)
from cvx.linalg.ewm_cov import ewm_covariance  # requires polars

Functions

• a_norm(vector, matrix=None) — Euclidean norm or NaN-aware matrix norm
• cholesky(cov, rhs=None) — Upper triangular Cholesky factor R such that R.T @ R = cov; when rhs is given, solves cov @ x = rhs (falls back to LU for non-positive-definite matrices)
• cond(matrix, p=None) — Condition number of a matrix (NaN-aware); accepts the same p norm values as numpy.linalg.cond
• det(matrix, cond_threshold=1e12) — Determinant of a square matrix with NaN-aware submatrix handling; emits IllConditionedMatrixWarning when near-singular
• eigvals(matrix) — Eigenvalues of a general square matrix (supports complex output for non-symmetric matrices)
• eigh(matrix) — Eigenvalues/eigenvectors of the valid symmetric/Hermitian submatrix in ascending eigenvalue order
• eigvalsh(matrix) — Eigenvalues-only convenience wrapper around eigh
• ewm_covariance(data, assets, index_col, window, is_halflife, warmup) — Exponentially weighted covariance matrices from a Polars DataFrame
• inv(matrix, cond_threshold=1e12) — Invert a matrix restricted to valid rows/columns; NaN rows/columns are returned for invalid positions
• inv_a_norm(vector, matrix=None) — Euclidean norm or inverse NaN-aware matrix norm
• lstsq(matrix, rhs, cond_threshold=1e12) — Solve a least-squares system with NaN-aware row filtering; returns (x, residuals, rank, sv) consistent with numpy.linalg.lstsq
• is_positive_definite(matrix) — Return True if the matrix is symmetric positive-definite
• pca(returns, n_components) — Principal Component Analysis via SVD
• qr(matrix) — Reduced QR decomposition, matching np.linalg.qr(mode='reduced')
• rand_cov(n, seed) — Random positive semi-definite covariance matrix
• solve(matrix, rhs, cond_threshold=1e12) — Solve a linear system restricted to valid rows/columns; NaN entries are returned for invalid positions
• svd(matrix) — Raw compact singular value decomposition via np.linalg.svd(full_matrices=False)
• valid(matrix) — Return a boolean mask and valid submatrix by removing rows/columns with non-finite diagonal entries

Exceptions & Warnings

• DimensionMismatchError — Raised when vector length does not match matrix dimension
• IllConditionedMatrixWarning — Emitted when the condition number exceeds a configurable threshold
• NegativeWarmupError — Raised when a negative warmup period is passed to ewm_covariance
• NonSquareMatrixError — Raised when a square matrix is required but the input is not square
• SingularMatrixError — Raised when a matrix is numerically singular
• check_and_warn_condition(matrix, threshold) — Emit IllConditionedMatrixWarning when the condition number exceeds the threshold

Types

• Matrix — Type alias for numpy.ndarray with float64 dtype