The tridiagonal algorithm is stable for all positive definite matrices. If you have any problems with the tridiagonal solver, you should question the validity of your problem formulation. What is there about your application that seems to demand division by zero?

Stanford Exploration Project

10/31/1997