lu_solveΒΆ

LU solve wrapper.

Computes the solution to

\[\mathbf{Ax} = \mathbf{b}\]

or

\[\mathbf{A}^T\mathbf{x} = \mathbf{b}\]

if trans=1.

It uses the SuperLU.solve() method if the input is a SuperLU or else will revert to the dense methods in scipy.

param lu_A:object or tuple containing the information of the LU factorisation
type lu_A:SuperLU or tuple
param b:Right hand side vector to solve
type b:np.ndarray
param trans:0 or 1 for either solution option.
type trans:int
returns:Solution to the system.
rtype:np.ndarray