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