Initial documentation

Still pretty rough.  We can polish it up and go over some details before we submit it.
A little more than tempted to spout the OO nonsense because even though the preconditions may be different the statement that a tri-diagonal matrix is a sparse matrix is still true in the general sense. It may not be in the way that we design and implement it but generally that statement is still true.