- Implement fast solvers for the case of full band toepliz matrices
- (maybe) Implement displacement CR for intermediate value of the Toeplitz bandwith
- Implement matrix functions (funm, sqrtm, logm, ...)
- fast multiplication for full band toeplitz matrices (?)