started experimenting with block-tridiagonal solver for implicit schemes