syr2k: Symmetric rank 2k update

\( C = \alpha A B^T + \alpha B A^T + \beta C \) where \( C \) is symmetric More...

Functions
int	plasma_csyr2k (plasma_enum_t uplo, plasma_enum_t trans, int n, int k, plasma_complex32_t alpha, plasma_complex32_t *pA, int lda, plasma_complex32_t *pB, int ldb, plasma_complex32_t beta, plasma_complex32_t *pC, int ldc)
void	plasma_omp_csyr2k (plasma_enum_t uplo, plasma_enum_t trans, plasma_complex32_t alpha, plasma_desc_t A, plasma_desc_t B, plasma_complex32_t beta, plasma_desc_t C, plasma_sequence_t *sequence, plasma_request_t *request)
int	plasma_dsyr2k (plasma_enum_t uplo, plasma_enum_t trans, int n, int k, double alpha, double *pA, int lda, double *pB, int ldb, double beta, double *pC, int ldc)
void	plasma_omp_dsyr2k (plasma_enum_t uplo, plasma_enum_t trans, double alpha, plasma_desc_t A, plasma_desc_t B, double beta, plasma_desc_t C, plasma_sequence_t *sequence, plasma_request_t *request)
int	plasma_ssyr2k (plasma_enum_t uplo, plasma_enum_t trans, int n, int k, float alpha, float *pA, int lda, float *pB, int ldb, float beta, float *pC, int ldc)
void	plasma_omp_ssyr2k (plasma_enum_t uplo, plasma_enum_t trans, float alpha, plasma_desc_t A, plasma_desc_t B, float beta, plasma_desc_t C, plasma_sequence_t *sequence, plasma_request_t *request)
int	plasma_zsyr2k (plasma_enum_t uplo, plasma_enum_t trans, int n, int k, plasma_complex64_t alpha, plasma_complex64_t *pA, int lda, plasma_complex64_t *pB, int ldb, plasma_complex64_t beta, plasma_complex64_t *pC, int ldc)
void	plasma_omp_zsyr2k (plasma_enum_t uplo, plasma_enum_t trans, plasma_complex64_t alpha, plasma_desc_t A, plasma_desc_t B, plasma_complex64_t beta, plasma_desc_t C, plasma_sequence_t *sequence, plasma_request_t *request)

Detailed Description

\( C = \alpha A B^T + \alpha B A^T + \beta C \) where \( C \) is symmetric

Function Documentation

int plasma_csyr2k	(	plasma_enum_t	uplo,
		plasma_enum_t	trans,
		int	n,
		int	k,
		plasma_complex32_t	alpha,
		plasma_complex32_t *	pA,
		int	lda,
		plasma_complex32_t *	pB,
		int	ldb,
		plasma_complex32_t	beta,
		plasma_complex32_t *	pC,
		int	ldc
	)

Performs one of the symmetric rank 2k operations

\[ C = \alpha A \times B^T + \alpha B \times A^T + \beta C, \]

or

\[ C = \alpha A^T \times B + \alpha B^T \times A + \beta C, \]

where alpha and beta are scalars, C is an n-by-n symmetric matrix, and A and B are n-by-k matrices in the first case and k-by-n matrices in the second case.

Parameters

[in]	uplo	PlasmaUpper: Upper triangle of C is stored; PlasmaLower: Lower triangle of C is stored.
[in]	trans	PlasmaNoTrans: \[ C = \alpha A \times B^T + \alpha B \times A^T + \beta C; \] PlasmaTrans: \[ C = \alpha A^T \times B + \alpha B^T \times A + \beta C. \]
[in]	n	The order of the matrix C. n >= zero.
[in]	k	If trans = PlasmaNoTrans, number of columns of the A and B matrices; if trans = PlasmaTrans, number of rows of the A and B matrices.
[in]	alpha	The scalar alpha.
[in]	pA	An lda-by-ka matrix. If trans = PlasmaNoTrans, ka = k; if trans = PlasmaTrans, ka = n.
[in]	lda	The leading dimension of the array A. If trans = PlasmaNoTrans, lda >= max(1, n); if trans = PlasmaTrans, lda >= max(1, k).
[in]	pB	An ldb-by-kb matrix. If trans = PlasmaNoTrans, kb = k; if trans = PlasmaTrans, kb = n.
[in]	ldb	The leading dimension of the array B. If trans = PlasmaNoTrans, ldb >= max(1, n); if trans = PlasmaTrans, ldb >= max(1, k).
[in]	beta	The scalar beta.
[in,out]	pC	An ldc-by-n matrix. On exit, the uplo part of the matrix is overwritten by the uplo part of the updated matrix.
[in]	ldc	The leading dimension of the array C. ldc >= max(1, n).

Return values

PlasmaSuccess

successful exit

See also

plasma_omp_csyr2k

plasma_csyr2k

plasma_dsyr2k

plasma_ssyr2k

void plasma_omp_csyr2k	(	plasma_enum_t	uplo,
		plasma_enum_t	trans,
		plasma_complex32_t	alpha,
		plasma_desc_t	A,
		plasma_desc_t	B,
		plasma_complex32_t	beta,
		plasma_desc_t	C,
		plasma_sequence_t *	sequence,
		plasma_request_t *	request
	)

Performs rank 2k update. Non-blocking tile version of plasma_csyr2k(). May return before the computation is finished. Operates on matrices stored by tiles. All matrices are passed through descriptors. All dimensions are taken from the descriptors. Allows for pipelining of operations at runtime.

Parameters

[in]	uplo	PlasmaUpper: Upper triangle of C is stored; PlasmaLower: Lower triangle of C is stored.
[in]	trans	PlasmaNoTrans: \[ C = \alpha A \times B^T + \alpha B \times A^T + \beta C; \] PlasmaTrans: \[ C = \alpha A^T \times B + \alpha B^T \times A + \beta C. \]
[in]	alpha	The scalar alpha.
[in]	A	Descriptor of matrix A.
[in]	B	Descriptor of matrix B.
[in]	beta	The scalar beta.
[in,out]	C	Descriptor of matrix C.
[in]	sequence	Identifies the sequence of function calls that this call belongs to (for completion checks and exception handling purposes). Check the sequence->status for errors.
[out]	request	Identifies this function call (for exception handling purposes).

Return values

void	Errors are returned by setting sequence->status and request->status to error values. The sequence->status and request->status should never be set to PlasmaSuccess (the initial values) since another async call may be setting a failure value at the same time.

See also

plasma_csyr2k

plasma_omp_csyr2k

plasma_omp_csyr2k

int plasma_dsyr2k	(	plasma_enum_t	uplo,
		plasma_enum_t	trans,
		int	n,
		int	k,
		double	alpha,
		double *	pA,
		int	lda,
		double *	pB,
		int	ldb,
		double	beta,
		double *	pC,
		int	ldc
	)

Performs one of the symmetric rank 2k operations

\[ C = \alpha A \times B^T + \alpha B \times A^T + \beta C, \]

or

\[ C = \alpha A^T \times B + \alpha B^T \times A + \beta C, \]

where alpha and beta are scalars, C is an n-by-n symmetric matrix, and A and B are n-by-k matrices in the first case and k-by-n matrices in the second case.

Parameters

[in]	uplo	PlasmaUpper: Upper triangle of C is stored; PlasmaLower: Lower triangle of C is stored.
[in]	trans	PlasmaNoTrans: \[ C = \alpha A \times B^T + \alpha B \times A^T + \beta C; \] PlasmaTrans: \[ C = \alpha A^T \times B + \alpha B^T \times A + \beta C. \]
[in]	n	The order of the matrix C. n >= zero.
[in]	k	If trans = PlasmaNoTrans, number of columns of the A and B matrices; if trans = PlasmaTrans, number of rows of the A and B matrices.
[in]	alpha	The scalar alpha.
[in]	pA	An lda-by-ka matrix. If trans = PlasmaNoTrans, ka = k; if trans = PlasmaTrans, ka = n.
[in]	lda	The leading dimension of the array A. If trans = PlasmaNoTrans, lda >= max(1, n); if trans = PlasmaTrans, lda >= max(1, k).
[in]	pB	An ldb-by-kb matrix. If trans = PlasmaNoTrans, kb = k; if trans = PlasmaTrans, kb = n.
[in]	ldb	The leading dimension of the array B. If trans = PlasmaNoTrans, ldb >= max(1, n); if trans = PlasmaTrans, ldb >= max(1, k).
[in]	beta	The scalar beta.
[in,out]	pC	An ldc-by-n matrix. On exit, the uplo part of the matrix is overwritten by the uplo part of the updated matrix.
[in]	ldc	The leading dimension of the array C. ldc >= max(1, n).

Return values

PlasmaSuccess

successful exit

See also

plasma_omp_dsyr2k

plasma_csyr2k

plasma_dsyr2k

plasma_ssyr2k

void plasma_omp_dsyr2k	(	plasma_enum_t	uplo,
		plasma_enum_t	trans,
		double	alpha,
		plasma_desc_t	A,
		plasma_desc_t	B,
		double	beta,
		plasma_desc_t	C,
		plasma_sequence_t *	sequence,
		plasma_request_t *	request
	)

Performs rank 2k update. Non-blocking tile version of plasma_dsyr2k(). May return before the computation is finished. Operates on matrices stored by tiles. All matrices are passed through descriptors. All dimensions are taken from the descriptors. Allows for pipelining of operations at runtime.

Parameters

[in]	uplo	PlasmaUpper: Upper triangle of C is stored; PlasmaLower: Lower triangle of C is stored.
[in]	trans	PlasmaNoTrans: \[ C = \alpha A \times B^T + \alpha B \times A^T + \beta C; \] PlasmaTrans: \[ C = \alpha A^T \times B + \alpha B^T \times A + \beta C. \]
[in]	alpha	The scalar alpha.
[in]	A	Descriptor of matrix A.
[in]	B	Descriptor of matrix B.
[in]	beta	The scalar beta.
[in,out]	C	Descriptor of matrix C.
[in]	sequence	Identifies the sequence of function calls that this call belongs to (for completion checks and exception handling purposes). Check the sequence->status for errors.
[out]	request	Identifies this function call (for exception handling purposes).

Return values

void	Errors are returned by setting sequence->status and request->status to error values. The sequence->status and request->status should never be set to PlasmaSuccess (the initial values) since another async call may be setting a failure value at the same time.

See also

plasma_dsyr2k

plasma_omp_dsyr2k

plasma_omp_csyr2k

int plasma_ssyr2k	(	plasma_enum_t	uplo,
		plasma_enum_t	trans,
		int	n,
		int	k,
		float	alpha,
		float *	pA,
		int	lda,
		float *	pB,
		int	ldb,
		float	beta,
		float *	pC,
		int	ldc
	)

Performs one of the symmetric rank 2k operations

\[ C = \alpha A \times B^T + \alpha B \times A^T + \beta C, \]

or

\[ C = \alpha A^T \times B + \alpha B^T \times A + \beta C, \]

where alpha and beta are scalars, C is an n-by-n symmetric matrix, and A and B are n-by-k matrices in the first case and k-by-n matrices in the second case.

Parameters

[in]	uplo	PlasmaUpper: Upper triangle of C is stored; PlasmaLower: Lower triangle of C is stored.
[in]	trans	PlasmaNoTrans: \[ C = \alpha A \times B^T + \alpha B \times A^T + \beta C; \] PlasmaTrans: \[ C = \alpha A^T \times B + \alpha B^T \times A + \beta C. \]
[in]	n	The order of the matrix C. n >= zero.
[in]	k	If trans = PlasmaNoTrans, number of columns of the A and B matrices; if trans = PlasmaTrans, number of rows of the A and B matrices.
[in]	alpha	The scalar alpha.
[in]	pA	An lda-by-ka matrix. If trans = PlasmaNoTrans, ka = k; if trans = PlasmaTrans, ka = n.
[in]	lda	The leading dimension of the array A. If trans = PlasmaNoTrans, lda >= max(1, n); if trans = PlasmaTrans, lda >= max(1, k).
[in]	pB	An ldb-by-kb matrix. If trans = PlasmaNoTrans, kb = k; if trans = PlasmaTrans, kb = n.
[in]	ldb	The leading dimension of the array B. If trans = PlasmaNoTrans, ldb >= max(1, n); if trans = PlasmaTrans, ldb >= max(1, k).
[in]	beta	The scalar beta.
[in,out]	pC	An ldc-by-n matrix. On exit, the uplo part of the matrix is overwritten by the uplo part of the updated matrix.
[in]	ldc	The leading dimension of the array C. ldc >= max(1, n).

Return values

PlasmaSuccess

successful exit

See also

plasma_omp_ssyr2k

plasma_csyr2k

plasma_dsyr2k

plasma_ssyr2k

void plasma_omp_ssyr2k	(	plasma_enum_t	uplo,
		plasma_enum_t	trans,
		float	alpha,
		plasma_desc_t	A,
		plasma_desc_t	B,
		float	beta,
		plasma_desc_t	C,
		plasma_sequence_t *	sequence,
		plasma_request_t *	request
	)

Performs rank 2k update. Non-blocking tile version of plasma_ssyr2k(). May return before the computation is finished. Operates on matrices stored by tiles. All matrices are passed through descriptors. All dimensions are taken from the descriptors. Allows for pipelining of operations at runtime.

Parameters

[in]	uplo	PlasmaUpper: Upper triangle of C is stored; PlasmaLower: Lower triangle of C is stored.
[in]	trans	PlasmaNoTrans: \[ C = \alpha A \times B^T + \alpha B \times A^T + \beta C; \] PlasmaTrans: \[ C = \alpha A^T \times B + \alpha B^T \times A + \beta C. \]
[in]	alpha	The scalar alpha.
[in]	A	Descriptor of matrix A.
[in]	B	Descriptor of matrix B.
[in]	beta	The scalar beta.
[in,out]	C	Descriptor of matrix C.
[in]	sequence	Identifies the sequence of function calls that this call belongs to (for completion checks and exception handling purposes). Check the sequence->status for errors.
[out]	request	Identifies this function call (for exception handling purposes).

Return values

void	Errors are returned by setting sequence->status and request->status to error values. The sequence->status and request->status should never be set to PlasmaSuccess (the initial values) since another async call may be setting a failure value at the same time.

See also

plasma_ssyr2k

plasma_omp_ssyr2k

plasma_omp_csyr2k

int plasma_zsyr2k	(	plasma_enum_t	uplo,
		plasma_enum_t	trans,
		int	n,
		int	k,
		plasma_complex64_t	alpha,
		plasma_complex64_t *	pA,
		int	lda,
		plasma_complex64_t *	pB,
		int	ldb,
		plasma_complex64_t	beta,
		plasma_complex64_t *	pC,
		int	ldc
	)

Performs one of the symmetric rank 2k operations

\[ C = \alpha A \times B^T + \alpha B \times A^T + \beta C, \]

or

\[ C = \alpha A^T \times B + \alpha B^T \times A + \beta C, \]

where alpha and beta are scalars, C is an n-by-n symmetric matrix, and A and B are n-by-k matrices in the first case and k-by-n matrices in the second case.

Parameters

[in]	uplo	PlasmaUpper: Upper triangle of C is stored; PlasmaLower: Lower triangle of C is stored.
[in]	trans	PlasmaNoTrans: \[ C = \alpha A \times B^T + \alpha B \times A^T + \beta C; \] PlasmaTrans: \[ C = \alpha A^T \times B + \alpha B^T \times A + \beta C. \]
[in]	n	The order of the matrix C. n >= zero.
[in]	k	If trans = PlasmaNoTrans, number of columns of the A and B matrices; if trans = PlasmaTrans, number of rows of the A and B matrices.
[in]	alpha	The scalar alpha.
[in]	pA	An lda-by-ka matrix. If trans = PlasmaNoTrans, ka = k; if trans = PlasmaTrans, ka = n.
[in]	lda	The leading dimension of the array A. If trans = PlasmaNoTrans, lda >= max(1, n); if trans = PlasmaTrans, lda >= max(1, k).
[in]	pB	An ldb-by-kb matrix. If trans = PlasmaNoTrans, kb = k; if trans = PlasmaTrans, kb = n.
[in]	ldb	The leading dimension of the array B. If trans = PlasmaNoTrans, ldb >= max(1, n); if trans = PlasmaTrans, ldb >= max(1, k).
[in]	beta	The scalar beta.
[in,out]	pC	An ldc-by-n matrix. On exit, the uplo part of the matrix is overwritten by the uplo part of the updated matrix.
[in]	ldc	The leading dimension of the array C. ldc >= max(1, n).

Return values

PlasmaSuccess

successful exit

See also

plasma_omp_zsyr2k

plasma_csyr2k

plasma_dsyr2k

plasma_ssyr2k

void plasma_omp_zsyr2k	(	plasma_enum_t	uplo,
		plasma_enum_t	trans,
		plasma_complex64_t	alpha,
		plasma_desc_t	A,
		plasma_desc_t	B,
		plasma_complex64_t	beta,
		plasma_desc_t	C,
		plasma_sequence_t *	sequence,
		plasma_request_t *	request
	)

Performs rank 2k update. Non-blocking tile version of plasma_zsyr2k(). May return before the computation is finished. Operates on matrices stored by tiles. All matrices are passed through descriptors. All dimensions are taken from the descriptors. Allows for pipelining of operations at runtime.

Parameters

[in]	uplo	PlasmaUpper: Upper triangle of C is stored; PlasmaLower: Lower triangle of C is stored.
[in]	trans	PlasmaNoTrans: \[ C = \alpha A \times B^T + \alpha B \times A^T + \beta C; \] PlasmaTrans: \[ C = \alpha A^T \times B + \alpha B^T \times A + \beta C. \]
[in]	alpha	The scalar alpha.
[in]	A	Descriptor of matrix A.
[in]	B	Descriptor of matrix B.
[in]	beta	The scalar beta.
[in,out]	C	Descriptor of matrix C.
[in]	sequence	Identifies the sequence of function calls that this call belongs to (for completion checks and exception handling purposes). Check the sequence->status for errors.
[out]	request	Identifies this function call (for exception handling purposes).

Return values

void	Errors are returned by setting sequence->status and request->status to error values. The sequence->status and request->status should never be set to PlasmaSuccess (the initial values) since another async call may be setting a failure value at the same time.

See also

plasma_zsyr2k

plasma_omp_zsyr2k

plasma_omp_csyr2k