PLASMA
Parallel Linear Algebra Software for Multicore Architectures
tsqrt: QR factorization of a rectangular matrix of two tiles

Functions

 __attribute__ ((weak))
 

Detailed Description

Function Documentation

__attribute__ ( (weak)  )

Computes a QR factorization of a rectangular matrix formed by coupling an n-by-n upper triangular tile A1 on top of an m-by-n tile A2:

| A1 | = Q * R | A2 |

Parameters
[in]mThe number of columns of the tile A2. m >= 0.
[in]nThe number of rows of the tile A1. The number of columns of the tiles A1 and A2. n >= 0.
[in]ibThe inner-blocking size. ib >= 0.
[in,out]A1On entry, the n-by-n tile A1. On exit, the elements on and above the diagonal of the array contain the n-by-n upper trapezoidal tile R; the elements below the diagonal are not referenced.
[in]lda1The leading dimension of the array A1. LDA1 >= max(1,N).
[in,out]A2On entry, the m-by-n tile A2. On exit, all the elements with the array tau, represent the unitary tile Q as a product of elementary reflectors (see Further Details).
[in]lda2The leading dimension of the tile A2. lda2 >= max(1,m).
[out]TThe ib-by-n triangular factor T of the block reflector. T is upper triangular by block (economic storage); The rest of the array is not referenced.
[in]ldtThe leading dimension of the array T. ldt >= ib.
tauAuxiliary workspace array of length n.
workAuxiliary workspace array of length ib*n.
Return values
PlasmaSuccesssuccessful exit
<0 if -i, the i-th argument had an illegal value

Computes a QR factorization of a rectangular matrix formed by coupling an n-by-n upper triangular tile A1 on top of an m-by-n tile A2:

| A1 | = Q * R | A2 |

Parameters
[in]mThe number of columns of the tile A2. m >= 0.
[in]nThe number of rows of the tile A1. The number of columns of the tiles A1 and A2. n >= 0.
[in]ibThe inner-blocking size. ib >= 0.
[in,out]A1On entry, the n-by-n tile A1. On exit, the elements on and above the diagonal of the array contain the n-by-n upper trapezoidal tile R; the elements below the diagonal are not referenced.
[in]lda1The leading dimension of the array A1. LDA1 >= max(1,N).
[in,out]A2On entry, the m-by-n tile A2. On exit, all the elements with the array tau, represent the orthogonal tile Q as a product of elementary reflectors (see Further Details).
[in]lda2The leading dimension of the tile A2. lda2 >= max(1,m).
[out]TThe ib-by-n triangular factor T of the block reflector. T is upper triangular by block (economic storage); The rest of the array is not referenced.
[in]ldtThe leading dimension of the array T. ldt >= ib.
tauAuxiliary workspace array of length n.
workAuxiliary workspace array of length ib*n.
Return values
PlasmaSuccesssuccessful exit
<0 if -i, the i-th argument had an illegal value