[SciPy-user] fastest 'solve' for triangular matrix?

[Emanuele Olivetti]

Dear All,

Is there a better/faster way (provided by scipy) to solve a linear system
Lx=b (when L is lower-triangular) then the following?

piv = N.arange(L.shape[0])
x = scipy.linalg.lu_solve((L,piv),b)

The linear system is trivial to solve but using 'solve' is much much slower
and I was no t able to find something ad-hoc for triangular matrices.

Best,

Emanuele