Does boost provide any classes to implement a tridiagonal matrix?
http://www.boost.org/doc/libs/1_44_0/libs/numeric/ublas/doc/banded.htm