PLASMA
Parallel Linear Algebra Software for Multicore Architectures

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

Functions

int plasma_csyrk (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 beta, plasma_complex32_t *pC, int ldc)
 
void plasma_omp_csyrk (plasma_enum_t uplo, plasma_enum_t trans, plasma_complex32_t alpha, plasma_desc_t A, plasma_complex32_t beta, plasma_desc_t C, plasma_sequence_t *sequence, plasma_request_t *request)
 
int plasma_dsyrk (plasma_enum_t uplo, plasma_enum_t trans, int n, int k, double alpha, double *pA, int lda, double beta, double *pC, int ldc)
 
void plasma_omp_dsyrk (plasma_enum_t uplo, plasma_enum_t trans, double alpha, plasma_desc_t A, double beta, plasma_desc_t C, plasma_sequence_t *sequence, plasma_request_t *request)
 
int plasma_ssyrk (plasma_enum_t uplo, plasma_enum_t trans, int n, int k, float alpha, float *pA, int lda, float beta, float *pC, int ldc)
 
void plasma_omp_ssyrk (plasma_enum_t uplo, plasma_enum_t trans, float alpha, plasma_desc_t A, float beta, plasma_desc_t C, plasma_sequence_t *sequence, plasma_request_t *request)
 
int plasma_zsyrk (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 beta, plasma_complex64_t *pC, int ldc)
 
void plasma_omp_zsyrk (plasma_enum_t uplo, plasma_enum_t trans, plasma_complex64_t alpha, plasma_desc_t A, plasma_complex64_t beta, plasma_desc_t C, plasma_sequence_t *sequence, plasma_request_t *request)
 

Detailed Description

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

Function Documentation

int plasma_csyrk ( 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  beta,
plasma_complex32_t *  pC,
int  ldc 
)

Performs one of the symmetric rank k operations

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

or

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

where alpha and beta are scalars, C is an n-by-n symmetric matrix, and A is an n-by-k matrix in the first case and a k-by-n matrix 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 A^T + \beta C; \]

  • PlasmaTrans:

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

[in]nThe order of the matrix C. n >= 0.
[in]kIf trans = PlasmaNoTrans, number of columns of the A matrix; if trans = PlasmaTrans, number of rows of the A matrix.
[in]alphaThe scalar alpha.
[in]pAA is an lda-by-ka matrix. If trans = PlasmaNoTrans, ka = k; if trans = PlasmaTrans, ka = n.
[in]ldaThe leading dimension of the array A. If trans = PlasmaNoTrans, lda >= max(1, n); if trans = PlasmaTrans, lda >= max(1, k).
[in]betaThe scalar beta.
[in,out]pCC is an ldc-by-n matrix. On exit, the uplo part of the matrix is overwritten by the uplo part of the updated matrix.
[in]ldcThe leading dimension of the array C. ldc >= max(1, n).
Return values
PlasmaSuccesssuccessful exit
See also
plasma_omp_csyrk
plasma_csyrk
plasma_dsyrk
plasma_ssyrk
void plasma_omp_csyrk ( plasma_enum_t  uplo,
plasma_enum_t  trans,
plasma_complex32_t  alpha,
plasma_desc_t  A,
plasma_complex32_t  beta,
plasma_desc_t  C,
plasma_sequence_t *  sequence,
plasma_request_t *  request 
)

Performs rank k update. Non-blocking tile version of plasma_csyrk(). 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 A^T + \beta C; \]

  • PlasmaTrans:

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

[in]alphaThe scalar alpha.
[in]ADescriptor of matrix A.
[in]betaThe scalar beta.
[in,out]CDescriptor of matrix C.
[in]sequenceIdentifies the sequence of function calls that this call belongs to (for completion checks and exception handling purposes). Check the sequence->status for errors.
[out]requestIdentifies this function call (for exception handling purposes).
Return values
voidErrors 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_csyrk
plasma_omp_csyrk
plasma_omp_csyrk
plasma_omp_dsyrk
plasma_omp_ssyrk
int plasma_dsyrk ( plasma_enum_t  uplo,
plasma_enum_t  trans,
int  n,
int  k,
double  alpha,
double *  pA,
int  lda,
double  beta,
double *  pC,
int  ldc 
)

Performs one of the symmetric rank k operations

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

or

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

where alpha and beta are scalars, C is an n-by-n symmetric matrix, and A is an n-by-k matrix in the first case and a k-by-n matrix 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 A^T + \beta C; \]

  • PlasmaTrans:

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

[in]nThe order of the matrix C. n >= 0.
[in]kIf trans = PlasmaNoTrans, number of columns of the A matrix; if trans = PlasmaTrans, number of rows of the A matrix.
[in]alphaThe scalar alpha.
[in]pAA is an lda-by-ka matrix. If trans = PlasmaNoTrans, ka = k; if trans = PlasmaTrans, ka = n.
[in]ldaThe leading dimension of the array A. If trans = PlasmaNoTrans, lda >= max(1, n); if trans = PlasmaTrans, lda >= max(1, k).
[in]betaThe scalar beta.
[in,out]pCC is an ldc-by-n matrix. On exit, the uplo part of the matrix is overwritten by the uplo part of the updated matrix.
[in]ldcThe leading dimension of the array C. ldc >= max(1, n).
Return values
PlasmaSuccesssuccessful exit
See also
plasma_omp_dsyrk
plasma_csyrk
plasma_dsyrk
plasma_ssyrk
void plasma_omp_dsyrk ( plasma_enum_t  uplo,
plasma_enum_t  trans,
double  alpha,
plasma_desc_t  A,
double  beta,
plasma_desc_t  C,
plasma_sequence_t *  sequence,
plasma_request_t *  request 
)

Performs rank k update. Non-blocking tile version of plasma_dsyrk(). 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 A^T + \beta C; \]

  • PlasmaTrans:

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

[in]alphaThe scalar alpha.
[in]ADescriptor of matrix A.
[in]betaThe scalar beta.
[in,out]CDescriptor of matrix C.
[in]sequenceIdentifies the sequence of function calls that this call belongs to (for completion checks and exception handling purposes). Check the sequence->status for errors.
[out]requestIdentifies this function call (for exception handling purposes).
Return values
voidErrors 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_dsyrk
plasma_omp_dsyrk
plasma_omp_csyrk
plasma_omp_dsyrk
plasma_omp_ssyrk
int plasma_ssyrk ( plasma_enum_t  uplo,
plasma_enum_t  trans,
int  n,
int  k,
float  alpha,
float *  pA,
int  lda,
float  beta,
float *  pC,
int  ldc 
)

Performs one of the symmetric rank k operations

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

or

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

where alpha and beta are scalars, C is an n-by-n symmetric matrix, and A is an n-by-k matrix in the first case and a k-by-n matrix 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 A^T + \beta C; \]

  • PlasmaTrans:

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

[in]nThe order of the matrix C. n >= 0.
[in]kIf trans = PlasmaNoTrans, number of columns of the A matrix; if trans = PlasmaTrans, number of rows of the A matrix.
[in]alphaThe scalar alpha.
[in]pAA is an lda-by-ka matrix. If trans = PlasmaNoTrans, ka = k; if trans = PlasmaTrans, ka = n.
[in]ldaThe leading dimension of the array A. If trans = PlasmaNoTrans, lda >= max(1, n); if trans = PlasmaTrans, lda >= max(1, k).
[in]betaThe scalar beta.
[in,out]pCC is an ldc-by-n matrix. On exit, the uplo part of the matrix is overwritten by the uplo part of the updated matrix.
[in]ldcThe leading dimension of the array C. ldc >= max(1, n).
Return values
PlasmaSuccesssuccessful exit
See also
plasma_omp_ssyrk
plasma_csyrk
plasma_dsyrk
plasma_ssyrk
void plasma_omp_ssyrk ( plasma_enum_t  uplo,
plasma_enum_t  trans,
float  alpha,
plasma_desc_t  A,
float  beta,
plasma_desc_t  C,
plasma_sequence_t *  sequence,
plasma_request_t *  request 
)

Performs rank k update. Non-blocking tile version of plasma_ssyrk(). 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 A^T + \beta C; \]

  • PlasmaTrans:

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

[in]alphaThe scalar alpha.
[in]ADescriptor of matrix A.
[in]betaThe scalar beta.
[in,out]CDescriptor of matrix C.
[in]sequenceIdentifies the sequence of function calls that this call belongs to (for completion checks and exception handling purposes). Check the sequence->status for errors.
[out]requestIdentifies this function call (for exception handling purposes).
Return values
voidErrors 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_ssyrk
plasma_omp_ssyrk
plasma_omp_csyrk
plasma_omp_dsyrk
plasma_omp_ssyrk
int plasma_zsyrk ( 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  beta,
plasma_complex64_t *  pC,
int  ldc 
)

Performs one of the symmetric rank k operations

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

or

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

where alpha and beta are scalars, C is an n-by-n symmetric matrix, and A is an n-by-k matrix in the first case and a k-by-n matrix 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 A^T + \beta C; \]

  • PlasmaTrans:

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

[in]nThe order of the matrix C. n >= 0.
[in]kIf trans = PlasmaNoTrans, number of columns of the A matrix; if trans = PlasmaTrans, number of rows of the A matrix.
[in]alphaThe scalar alpha.
[in]pAA is an lda-by-ka matrix. If trans = PlasmaNoTrans, ka = k; if trans = PlasmaTrans, ka = n.
[in]ldaThe leading dimension of the array A. If trans = PlasmaNoTrans, lda >= max(1, n); if trans = PlasmaTrans, lda >= max(1, k).
[in]betaThe scalar beta.
[in,out]pCC is an ldc-by-n matrix. On exit, the uplo part of the matrix is overwritten by the uplo part of the updated matrix.
[in]ldcThe leading dimension of the array C. ldc >= max(1, n).
Return values
PlasmaSuccesssuccessful exit
See also
plasma_omp_zsyrk
plasma_csyrk
plasma_dsyrk
plasma_ssyrk
void plasma_omp_zsyrk ( plasma_enum_t  uplo,
plasma_enum_t  trans,
plasma_complex64_t  alpha,
plasma_desc_t  A,
plasma_complex64_t  beta,
plasma_desc_t  C,
plasma_sequence_t *  sequence,
plasma_request_t *  request 
)

Performs rank k update. Non-blocking tile version of plasma_zsyrk(). 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 A^T + \beta C; \]

  • PlasmaTrans:

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

[in]alphaThe scalar alpha.
[in]ADescriptor of matrix A.
[in]betaThe scalar beta.
[in,out]CDescriptor of matrix C.
[in]sequenceIdentifies the sequence of function calls that this call belongs to (for completion checks and exception handling purposes). Check the sequence->status for errors.
[out]requestIdentifies this function call (for exception handling purposes).
Return values
voidErrors 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_zsyrk
plasma_omp_zsyrk
plasma_omp_csyrk
plasma_omp_dsyrk
plasma_omp_ssyrk