MTL 4: Triangular Solvers
x= upper_trisolve(A, b);
x= unit_upper_trisolve(A, b);
On matrices with non-unit diagonals, the divisions can be circumvented by inverting the diagonal once with invert_diagonal(A) and then using:
x= inverse_upper_trisolve(A, b);
Likewise, the functions for lower triangular matrices are defined:
x= lower_trisolve(A, b); x= unit_lower_trisolve(A, b); x= inverse_lower_trisolve(A, b);
Return to Other Matrix Functions Table of Content Proceed to Introduction Krylov-Subspace Methods
Triangular Solvers -- MTL 4 -- Peter Gottschling and Andrew Lumsdaine
-- Gen. with
rev. 7542
on 7 Apr 2011 by doxygen 1.5.9 -- © 2010 by SimuNova UG.