[SciPy-User] RZ Factorization

Sturla Molden sturla.molden at gmail.com
Thu Sep 4 10:45:40 EDT 2014


Del Citto Francesco <francesco.delcitto at sauber-motorsport.com> wrote:

>> I would also like to change the triangular matrix solver to use level-3
>> BLAS *TRSM instead of LAPACK
>> *TRTRS. *TRSM has the advantage that is can solve AX=B and XA=B, but
>> *TRTRS can only solve AX=B unless
>> we compute the inverse of A explicitely. It often happens that the
>> triangular matrix is on the right-
>> hand side. And *TRTRS is just a shallow wrapper for *TRSM anyway.
> 
> I'm a bit too rusty with LAPACK to give an opinion on this!
> I haven't been used it for too long.

It was my fault for suggesting *TRTRS for the trisngular solver in the
first place. Big mistake :(

So I feel a bit responsible for correcting it.


Sturla




More information about the SciPy-User mailing list