We show that resulting operator matrices are banded and sparse due to the compact support of B-splines. Moreover, we use an algorithm based on Krylov subspace to exploit the sparsity of matrices for ...