mkl_?tppack

Developer Reference for Intel® oneAPI Math Kernel Library for Fortran

Download PDF

ID 766686

Date 12/16/2022

Version

Public

A newer version of this document is available. Customers should click here to go to the newest version.

Visible to Intel only — GUID: GUID-E72E69C2-316E-4124-A0A2-9213DC535A98

View Details

mkl_?tppack

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

Syntax

call mkl_stppack (uplo, trans, n, ap, i, j, rows, cols, a, lda, info )

call mkl_dtppack (uplo, trans, n, ap, i, j, rows, cols, a, lda, info )

call mkl_ctppack (uplo, trans, n, ap, i, j, rows, cols, a, lda, info )

call mkl_ztppack (uplo, trans, n, ap, i, j, rows, cols, a, lda, info )

call mkl_tppack (ap, i, j, rows, cols, a[, uplo] [, trans] [, info])

Include Files

mkl.fi, lapack.f90

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

NOTE:

Any elements of the copied submatrix rectangular outside of the triangular part of the matrix AP are skipped.

Input Parameters

The data types are given for the Fortran interface.

uplo

CHARACTER*1. 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

CHARACTER*1. 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

INTEGER. The order of the matrix AP; n ≥ 0

i, j

INTEGER. 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

INTEGER. Number of rows in the destination submatrix. 0 ≤rows≤n - i + 1.

cols

INTEGER. Number of columns in the destination submatrix. 0 ≤cols≤n - j + 1.

a

REAL for mkl_stppack

DOUBLE PRECISION for mkl_dtppack

COMPLEX for mkl_ctppack

DOUBLE COMPLEX for mkl_ztppack

Pointer to the source submatrix.

Array a(lda, *) 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 must be at least lda*cols for trans = 'N' or lda*rows for trans='T' or trans='C'.

NOTE:

If there are elements outside of the triangular part of AP, they are skipped and are not copied from a.

lda

INTEGER. The leading dimension of the array a.

lda≥ max(1, rows) for trans = 'N' and and lda≥ max(1, cols) for trans='T' or trans='C'.

Output Parameters

REAL for mkl_stppack

DOUBLE PRECISION for mkl_dtppack

COMPLEX for mkl_ctppack

DOUBLE COMPLEX for mkl_ztppack

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.

info

INTEGER. If info=0, the execution is successful. If info = -i, the i-th parameter had an illegal value.

Parent topic: LAPACK Auxiliary Routines

Select Your Language

Using Intel.com Search

Quick Links

Recent Searches

Advanced Search

Only search in

Developer Reference for Intel® oneAPI Math Kernel Library for Fortran

mkl_?tppack