Tested the block Thomas tridiagonal solver.

There is a precision issue as the size of the matrix
increases. This will be explored in the future.