Developer Reference for Intel® oneAPI Math Kernel Library - C

Developer Reference

Contents

mkl_?tppack

Copies a triangular/symmetric matrix or submatrix from standard full format to standard packed format.

Syntax

lapack_int LAPACKE_mkl_stppack
(
int
matrix_layout
,
char
uplo
,
char
trans
,
lapack_int
n
,
float*
ap
,
lapack_int
i
,
lapack_int
j
,
lapack_int
rows
,
lapack_int
cols
,
const
float*
a
,
lapack_int
lda
);
lapack_int LAPACKE_mkl_dtppack
(
int
matrix_layout
,
char
uplo
,
char
trans
,
lapack_int
n
,
double*
ap
,
lapack_int
i
,
lapack_int
j
,
lapack_int
rows
,
lapack_int
cols
,
const
double*
a
,
lapack_int
lda
);
lapack_int LAPACKE_mkl_ctppack
(
int
matrix_layout
,
char
uplo
,
char
trans
,
lapack_int
n
,
MKL_Complex8*
ap
,
lapack_int
i
,
lapack_int
j
,
lapack_int
rows
,
lapack_int
cols
,
const
MKL_Complex8*
a
,
lapack_int
lda
);
lapack_int LAPACKE_mkl_ztppack
(
int
matrix_layout
,
char
uplo
,
char
trans
,
lapack_int
n
,
MKL_Complex16*
ap
,
lapack_int
i
,
lapack_int
j
,
lapack_int
rows
,
lapack_int
cols
,
const
MKL_Complex16*
a
,
lapack_int
lda
);
Include Files
  • mkl.h
Description
The routine copies a triangular or symmetric matrix or its submatrix from standard full format to packed format 
AP
i:i+rows-1, j:j+cols-1
 := op(A)
Standard packed formats include:
  • TP: triangular packed storage
  • SP: symmetric indefinite packed storage
  • HP: Hermitian indefinite packed storage
  • PP: symmetric or Hermitian positive definite packed storage
Full formats include:
  • GE: general
  • TR: triangular
  • SY: symmetric indefinite
  • HE: Hermitian indefinite
  • PO: symmetric or Hermitian positive definite
Any elements of the copied submatrix rectangular outside of the triangular part of the matrix
AP
 are skipped.
Input Parameters
A
<datatype>
 placeholder, if present, is used for the C interface data types in the C interface section above. See C Interface Conventions for the C interface principal conventions and type definitions.
uplo
Specifies whether the matrix
AP
 is upper or lower triangular.
If
uplo
 = 'U',
AP
 is upper triangular.
If
uplo
 = 'L':
AP
 is lower triangular.
trans
Specifies whether or not the copied block of
A
 is transposed or not.
If
trans
 = 'N', no transpose:
op(
A
) =
A
.
If
trans
 = 'T',transpose:
op(
A
) =
A
T
.
If
trans
 = 'C',conjugate transpose:
op(
A
) =
A
H
. For real data this is the same as
trans
 = 'T'.
n
The order of the matrix
AP
;
n
 ≥ 0
i, j
Coordinates of the left upper corner of the destination submatrix in
AP
.
If
uplo
=’U’, 1
≤
i
≤
j
≤
n
.
If
uplo
=’L’, 1
≤
j
≤
i
≤
n
.
rows
Number of rows in the destination submatrix. 0
≤
rows
≤
n
 -
i
 + 1.
cols
Number of columns in the destination submatrix. 0
≤
cols
≤
n
 -
j
 + 1.
a
Pointer to the source submatrix.
Array
a
 contains the
rows
-by-
cols
 submatrix stored as unpacked rows-by-columns if
trans
 = ’N’, or unpacked columns-by-rows if
trans
 = ’T’ or
trans
 = ’C’.
The size of
a
 is
trans
 = 'N'
trans
='T' or
trans
='C'
matrix_layout
 =
LAPACK_COL_MAJOR
lda
*
cols
lda
*
rows
matrix_layout
 =
LAPACK_ROW_MAJOR
lda
*
rows
lda
*
cols
If there are elements outside of the triangular part of
AP
, they are skipped and are not copied from
a
.
lda
The leading dimension of the array
a
.
trans
 = 'N'
trans
='T' or
trans
='C'
matrix_layout
 =
LAPACK_COL_MAJOR
lda
≥
 max(1,
rows
)
lda
≥
 max(1,
cols
)
matrix_layout
 =
LAPACK_ROW_MAJOR
lda
≥
 max(1,
cols
)
lda
≥
 max(1,
rows
)
Output Parameters
ap
Array of size at least max(1,
n
(
n
+1)/2). The array
ap
 contains either the upper or the lower triangular part of the matrix
AP
 (as specified by
uplo
) in packed storage (see Matrix Storage Schemes). The submatrix of
ap
 from row
i
 to row
i
 +
rows
 - 1 and column
j
 to column
j
 +
cols
 - 1 is overwritten with a copy of the source matrix.
Return Values
This function returns a value
info
. If
info
=0, the execution is successful. If
info
 = -
i
, the
i
-th parameter had an illegal value.

Product and Performance Information

1

Performance varies by use, configuration and other factors. Learn more at www.Intel.com/PerformanceIndex.