 matrix expression
          beta*op( A )*T*op( A )', where T is the triangular matrix
          defined in the Method section.

  LDR     INTEGER
          The leading dimension of array R.  LDR &gt;= MAX(1,M).

  A       (input) DOUBLE PRECISION array, dimension (LDA,k)
          where k is N when TRANS = 'N' and is M when TRANS = 'T' or
          TRANS = 'C'.
          On entry with TRANS = 'N', the leading M-by-N part of this
          array must contain the matrix A.
          On entry with TRANS = 'T'