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
-- Generated on 19 May 2009 by Doxygen 1.5.5 -- Copyright 2007 by the Trustees of Indiana University.