PLASMA
2.8.0
PLASMA - Parallel Linear Algebra for Scalable Multi-core Architectures
|
int CORE_dparfb | ( | PLASMA_enum | side, |
PLASMA_enum | trans, | ||
PLASMA_enum | direct, | ||
PLASMA_enum | storev, | ||
int | M1, | ||
int | N1, | ||
int | M2, | ||
int | N2, | ||
int | K, | ||
int | L, | ||
double * | A1, | ||
int | LDA1, | ||
double * | A2, | ||
int | LDA2, | ||
const double * | V, | ||
int | LDV, | ||
const double * | T, | ||
int | LDT, | ||
double * | WORK, | ||
int | LDWORK | ||
) |
CORE_dparfb applies a complex upper triangular block reflector H or its transpose H' to a complex rectangular matrix formed by coupling two tiles A1 and A2. Matrix V is:
COLUMNWISE ROWWISE | K | | N2-L | L | __ _____________ __ __ _________________ __ | | | | | \ | | | | | \ L
M2-L | | | K |_______________|_____\ __ | | | M2 | | __ |____| | | | K-L \ | | __ |______________________| __ L \ | | __ |______| __ | N2 |
| L | K-L |
[in] | side |
|
[in] | trans |
|
[in] | direct | Indicates how H is formed from a product of elementary reflectors
|
[in] | storev | Indicates how the vectors which define the elementary reflectors are stored:
|
[in] | M1 | The number of columns of the tile A1. M1 >= 0. |
[in] | N1 | The number of rows of the tile A1. N1 >= 0. |
[in] | M2 | The number of columns of the tile A2. M2 >= 0. |
[in] | N2 | The number of rows of the tile A2. N2 >= 0. |
[in] | K | The order of the matrix T (= the number of elementary reflectors whose product defines the block reflector). |
[in] | L | The size of the triangular part of V |
[in,out] | A1 | On entry, the M1-by-N1 tile A1. On exit, A1 is overwritten by the application of Q. |
[in] | LDA1 | The leading dimension of the array A1. LDA1 >= max(1,N1). |
[in,out] | A2 | On entry, the M2-by-N2 tile A2. On exit, A2 is overwritten by the application of Q. |
[in] | LDA2 | The leading dimension of the tile A2. LDA2 >= max(1,N2). |
[in] | V | (LDV,K) if STOREV = 'C' (LDV,M2) if STOREV = 'R' and SIDE = 'L' (LDV,N2) if STOREV = 'R' and SIDE = 'R' Matrix V. |
[in] | LDV | The leading dimension of the array V. If STOREV = 'C' and SIDE = 'L', LDV >= max(1,M2); if STOREV = 'C' and SIDE = 'R', LDV >= max(1,N2); if STOREV = 'R', LDV >= K. |
[out] | T | The triangular K-by-K matrix T in the representation of the block reflector. T is upper triangular by block (economic storage); The rest of the array is not referenced. |
[in] | LDT | The leading dimension of the array T. LDT >= K. |
[in,out] | WORK | |
[in] | LDWORK | The dimension of the array WORK. |
PLASMA_SUCCESS | successful exit |
<0 | if -i, the i-th argument had an illegal value |