I'm getting ready to start working on a C/C++ project that will be building and solving large tri-diagonal, block tri-diagonal and triangular matrices. I know there are a lot of libraries available ...

We describe a computational procedure for evaluating the quasi-stationary distributions of a continuous-time Markov chain. Our method, which is an 'iterative version' of Arnoldi's algorithm, is ...