(S|D|C|Z)SYR2K¶

Single, double, single complex, and double complex SYR2K.

Description¶

Performs one of the following rank-2k updates on the symmetric matrix C.

$C \assign \alpha A*\trans{B} + \alpha B*\trans{A} + \beta C$

or

$C \assign \alpha \trans{A}*B + \alpha \trans{B}*A + \beta C$

BLAS Interface¶

void ssyr2k(const char *UPLO, const char *TRANS, const qml_long *N,
            const qml_long *K, const float *ALPHA, const float *A,
            const qml_long *LDA, const float *B, const qml_long *LDB,
            const float *BETA, float *C, const qml_long *LDC);

void dsyr2k(const char *UPLO, const char *TRANS, const qml_long *N,
            const qml_long *K, const double *ALPHA, const double *A,
            const qml_long *LDA, const double *B, const qml_long *LDB,
            const double *BETA, double *C, const qml_long *LDC);

void csyr2k(const char *UPLO, const char *TRANS, const qml_long *N,
            const qml_long *K, const qml_single_complex *ALPHA,
            const qml_single_complex *A, const qml_long *LDA,
            const qml_single_complex *B, const qml_long *LDB,
            const qml_single_complex *BETA, qml_single_complex *C,
            const qml_long *LDC);

void zsyr2k(const char *UPLO, const char *TRANS, const qml_long *N,
            const qml_long *K, const qml_double_complex *ALPHA,
            const qml_double_complex *A, const qml_long *LDA,
            const qml_double_complex *B, const qml_long *LDB,
            const qml_double_complex *BETA, qml_double_complex *C,
            const qml_long *LDC);

CBLAS Interface¶

void cblas_ssyr2k(const CBLAS_ORDER ORDER, const CBLAS_UPLO UPLO,
                  const CBLAS_TRANSPOSE TRANS, const qml_long N,
                  const qml_long K, const float ALPHA, const float *A,
                  const qml_long LDA, const float *B, const qml_long LDB,
                  const float BETA, float *C, const qml_long LDC);

void cblas_dsyr2k(const CBLAS_ORDER ORDER, const CBLAS_UPLO UPLO,
                  const CBLAS_TRANSPOSE TRANS, const qml_long N,
                  const qml_long K, const double ALPHA, const double *A,
                  const qml_long LDA, const double *B, const qml_long LDB,
                  const double BETA, double *C, const qml_long LDC);

void cblas_csyr2k(const CBLAS_ORDER ORDER, const CBLAS_UPLO UPLO,
                  const CBLAS_TRANSPOSE TRANS, const qml_long N,
                  const qml_long K, const qml_single_complex *ALPHA,
                  const qml_single_complex *A, const qml_long LDA,
                  const qml_single_complex *B, const qml_long LDB,
                  const qml_single_complex *BETA, qml_single_complex *C,
                  const qml_long LDC);

void cblas_zsyr2k(const CBLAS_ORDER ORDER, const CBLAS_UPLO UPLO,
                  const CBLAS_TRANSPOSE TRANS, const qml_long N,
                  const qml_long K, const qml_double_complex *ALPHA,
                  const qml_double_complex *A, const qml_long LDA,
                  const qml_double_complex *B, const qml_long LDB,
                  const qml_double_complex *BETA, qml_double_complex *C,
                  const qml_long LDC);

Arguments¶

UPLO	Specify whether the upper or lower triangle of matrix C will be used
TRANS	Specifies which rank-k update is performed:
	Non-Transpose: $C \assign \alpha A\trans{B} + \alpha B\trans{A} + \beta C$
	Transpose: $C \assign \alpha \trans{A}B + \alpha \trans{B}A + \beta C$
	Complex Conjugate Transpose: $C \assign \alpha \trans{A}B + \alpha \trans{B}A + \beta C$
N	Order of matrix C
K	For Non-Transpose, K is the number of columns of matrix A and B
	For Transpose and Complex Conjugate Transpose, K is the number of rows of matrix A and B
ALPHA	Scalar multiplied with the matrix-matrix product
A	Input matrix A
LDA	Leading dimension of matrix A
B	Input matrix B
LDB	Leading dimension of matrix B
BETA	Scalar multiplied with matrix C
C	Symmetric matrix C
LDC	Leading dimension of matrix C