warp.tile_lower_solve#

warp.tile_lower_solve(
L: Tile[Float, tuple[int, int]],
y: Tile[Float, tuple[int]],
) → Tile[Float, tuple[int]]#

  • Kernel

Solve for z in Lz = y, where L is a lower triangular matrix.

This performs general forward substitution for a lower triangular system.

Note that computing the adjoint is not yet supported.

Supported datatypes are:

  • float32

  • float64

Parameters:

  • L – A square, non-singular, lower triangular matrix

  • y – A 1D or 2D tile with compatible shape

Returns:

A tile of the same shape as y such that Lz = y.