LowerTriQSM

LowerTriQSM#

class tinygp.solvers.quasisep.core.LowerTriQSM(diag: DiagQSM, lower: StrictLowerTriQSM)[source]#

Bases: QSM

A lower triangular quasiseparable matrix

Parameters:

  • diag – The diagonal elements.

  • lower – The strictly lower triangular elements.

matmul(x: tinygp.helpers.JAXArray, *, parallel: bool = False) → tinygp.helpers.JAXArray[source]#

The dot product of this matrix with a dense vector or matrix

Parameters:

  • x (n, ...) – A matrix or vector with leading dimension matching this matrix.

  • parallel – If True, use a parallel associative-scan algorithm instead of the default sequential scan.

scale(other: tinygp.helpers.JAXArray) → LowerTriQSM[source]#

The multiplication of this matrix times a scalar, as a QSM

property shape: tuple[int, int]#

The shape of the matrix

solve(y: tinygp.helpers.JAXArray, *, parallel: bool = False) → tinygp.helpers.JAXArray[source]#

Solve a linear system with this matrix

If this matrix is called L, this solves L @ x = y for x given y, using forward substitution.

Parameters:

  • y (n, ...) – A matrix or vector with leading dimension matching this matrix.

  • parallel – If True, use a parallel associative-scan algorithm.

to_dense() → tinygp.helpers.JAXArray#

Render this representation to a dense matrix

This implementation is not optimized and should really only ever be used for testing purposes.

transpose() → UpperTriQSM[source]#

The matrix transpose as a QSM