ONT></B>
<PRE>
  To compute the incomplete Cholesky (ICC) factor of a symmetric
  positive definite (s.p.d.) block Toeplitz matrix T, defined by
  either its first block row, or its first block column, depending
  on the routine parameter TYPET.

  By subsequent calls of this routine, further rows / columns of
  the Cholesky factor can be added.
  Furthermore, the generator of the Schur complement of the leading
  (P+S)*K-by-(P+S)*K block in T is available, which can be used,
  e.g., for measuring the quali