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