fast solver for symmetric positive definite toeplitz system with preconditioned conjugate gradient method.