Parallel computation of a Krylov matrix for a sparse and structured input

As in our previous work, we reduce parallel computation of a Krylov matrix to solving a parametrized linear system of equations. This time we show that such a method is effective in the cases of banded matrices, sparse and structured matrices and triangular matrices.