I have an upper triangular Toeplitz matrix whose elements are block Toeplitz matrices as well. I need to factor this matrix as a poduct of two upper triangular block Toeplitz matrices having the same structure as the original matrix. Furthermore, the original matrix is rank defficient. An SVD will destroy the Toeplitz structure that I need. Does anyone know an algorithm that can do this? Any help will be appreciated.