If I have this:
A * f = g;
A: upper triangular matrix (n x n)
f: (n x 1)
g: (n x 1)
Need to solve for f using back substitution algorithm. I would say that it not really that hard to write one myself, but oh well, if there is a library out there, then why not.
Boost uBlas should work. At least if I understand your question correctly, you probably want to start by looking at
lu_substitute()andinplace_solve().