[Pd] Toeplitz solver libraries? Is there a new fast routine anywhere?