mkl_?getrinp_compact

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( A_c) of a set of general, n x n matrices A_c, that have been stored in Compact format. The factorization (output) data will also be stored in Compact format.

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 A_c; n >= 0.

ap

Points to the beginning of the nm A_c matrices. On entry, ap contains the LU factorizations of A_c, stored in Compact format, as returned by mkl_?getrfnp_compact : A_c=L_c*U_c.

See Compact Format for more details.

ldap

Column stride (column-major layout) or row stride (row-major layout) of A_c.

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 ³ for real flavors and nm* (16/3) n ³ for complex flavors.

Output Parameters

ap

On exit, A _c is overwritten by inv(A_c). ap points to the beginning of nm inv(A_c) 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.