Added untest blockThomas tridiagonal solver routine