Dgemm Dsymm Dsyr2k Dsyrk Dtrmm Dtrsm

Dsyrk() Sub Dsyrk ( Uplo As  String, Trans As  String, N As  Long, K As  Long, Alpha As  Double, A() As  Double, Beta As  Double, C() As  Double, Optional Info As  Long  ) Rank k operation: C <- αAAT + βC or C <- αATA + βC (symmetric matrices) (BLAS 3) Purpose This routine performs one of the symmetric rank k operations C <- αAA^T + βC or C <- αA^TA + βC where α and β are scalars, C is an n x n symmetric matrix and A is an n x k matrix in the first case and a k x n matrix in the second case. Parameters [in] Uplo Specifies whether the upper or lower triangular part of the array C() is to be referenced as follows: = "U": Only the upper triangular part of C() is to be referenced. = "L": Only the lower triangular part of C() is to be referenced. [in] Trans Specifies the operation to be performed as follows: = "N": C <- αAA^T + βC. = "T" or "C": C <- αA^TA + βC. [in] N Order of the matrix C. (N >= 0) (If N = 0, returns without computation) [in] K Number of columns of the matrix A when trans = "N", and number of rows of the matrix A when trans = "T" or "C". (K >= 0) [in] Alpha Scalar α. [in] A() Array A(LA1 - 1, LA2 - 1) (LA1 >= N, LA2 >= K (if Trans="N"), LA1 >= K, LA2 >= N (otherwise)) N x K matrix A when Trans = "N", or K x N matrix A otherwise. [in] Beta Scalar β. [in,out] C() Array C(LC1 - 1, LC2 - 1) (LC1 >= N, LC2 >= N) [in] N x N symmetric matrix C. Only the upper or lower triangular part is to be referenced in accordance with Uplo. [out] αAA^T + βC (if Trans="N") or αA^TA + βC (otherwise). Only the upper or lower triangular part is overwritten in accordance with Uplo. [out] Info (Optional) = 0: Successful exit. = -1: The argument Uplo had an illegal value. (Uplo <> "U" nor "L") = -2: The argument Trans had an illegal value. (Trans <> "N", "T" nor "C") = -3: The argument N had an illegal value. (N < 0) = -4: The argument K had an illegal value. (K < 0) = -6: The argument A() is invalid. = -8: The argument C() is invalid. Reference BLAS