 rows / columns
          of R.  0 &lt;= P &lt;= N.

  S       (input)  INTEGER
          The number of block rows / columns of R to compute.
          0 &lt;= S &lt;= N-P.

  T       (input/output)  DOUBLE PRECISION array, dimension
          (LDT,(N-P)*K) / (LDT,K)
          On entry, if P = 0, then the leading K-by-N*K / N*K-by-K
          part of this array must contain the first block row /
          column of an s.p.d. block Toeplitz matrix.
          If P &gt; 0, the leading K-by-(N-P)*K / (N-P)*K-by