Developer Reference for Intel® oneAPI Math Kernel Library for C

ID 766684
Date 7/13/2023
Public

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

Document Table of Contents

mkl_?getrinp_compact

Computes the inverse of a set of LU-factorized general matrices, without pivoting, stored in the compact format (see Compact Format for details).

Syntax

void mkl_sgetrinp_compact (MKL_LAYOUT layout, MKL_INT n, float * ap, MKL_INT ldap, float * work, MKL_INT lwork, MKL_INT * info, MKL_COMPACT_PACK format, MKL_INT nm);

void mkl_dgetrinp_compact (MKL_LAYOUT layout, MKL_INT n, double * ap, MKL_INT ldap, double * work, MKL_INT lwork, MKL_INT * info, MKL_COMPACT_PACK format, MKL_INT nm);

void mkl_cgetrinp_compact (MKL_LAYOUT layout, MKL_INT n, float * ap, MKL_INT ldap, float * work, MKL_INT lwork, MKL_INT * info, MKL_COMPACT_PACK format, MKL_INT nm);

void mkl_zgetrinp_compact (MKL_LAYOUT layout, MKL_INT n, double * ap, MKL_INT ldap, double * work, MKL_INT lwork, MKL_INT * info, MKL_COMPACT_PACK format, MKL_INT nm);

Description

This routine computes the inverse inv( Ac) of a set of general, n x n matrices Ac, that have been stored in Compact format. The factorization (output) data will also be stored in Compact format.

NOTE:

Compact routines have some limitations; see Numerical Limitations.

Input Parameters

layout

Specifies whether two-dimensional array storage is row-major (MKL_ROW_MAJOR) or column-major (MKL_COL_MAJOR).

n
The order of Ac; n >= 0.
ap
Points to the beginning of the nm Ac matrices. On entry, ap contains the LU factorizations of Ac, stored in Compact format, as returned by mkl_?getrfnp_compact : Ac=Lc*Uc.

See Compact Format for more details.

ldap
Column stride (column-major layout) or row stride (row-major layout) of Ac.
work
Points to the beginning of the work array.
lwork
The size of the work array. If lwork = -1, a workspace query is assumed; the routine calculates only the optimal size of the work array and returns this value as the first entry of the work array.
format
Specifies the format of the compact matrices. See Compact Format or mkl_get_format_compact for details.
nm
Total number of matrices stored in Compact format.

Application Notes:

Before calling this routine, mkl_?gepack_compact must be called. After calling this routine, mkl_?geunpack_compact should be called, unless another compact routine will be subsequently called on the Compact format matrices.

The total number of floating-point operations is approximately nm* (4/3) n 3 for real flavors and nm* (16/3) n 3 for complex flavors.

Output Parameters

ap
On exit, A c is overwritten by inv(Ac). ap points to the beginning of nm inv(Ac) matrices stored in Compact format.
work[0]
On exit, work[0] contains the minimum value of lwork required for optimum performance. Use this lwork for subsequent runs.
info
The parameter is not currently used in this routine. It is reserved for the future use.