mkl_?getrinp_compact
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( A
c
) 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.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 thenmAcmatrices. On entry,apcontains 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 theworkarray.
- 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 ormkl_get_format_compactfor 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, Acis overwritten by inv(Ac).appoints 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.