int P_dorgqr(
int M,
int N,
int K,
double *A,
int LDA,
double *T,
double *Q,
int LDQ
) {
int info;
#if CHECK_VERSION_BEQ(2,4,5)
info = PLASMA_dorgqr(M, N, K, A, LDA, T, Q, LDQ);
#else
PLASMA_desc *descT;
int NB, IB;
int MT, NT;
/* Get autotuned or set tile size; T matrix allocated with R */
PLASMA_Alloc_Workspace_dgeqrf(1, 1, &descT);
PLASMA_Get(PLASMA_TILE_SIZE, &NB);
PLASMA_Get(PLASMA_INNER_BLOCK_SIZE, &IB);
PLASMA_Dealloc_Handle_Tile(&descT);
MT = (M%NB==0) ? (M/NB) : (M/NB+1);
NT = (N%NB==0) ? (N/NB) : (N/NB+1);
// possibly allocate space for descT in R and keep it in qr object instead
info = PLASMA_Desc_Create(&descT, T, PlasmaComplexDouble,
IB, NB, IB*NB, MT*IB, NT*NB, 0, 0, MT*IB, NT*NB);
info = PLASMA_dorgqr(M, N, K, A, LDA, descT, Q, LDQ);
PLASMA_Desc_Destroy(&descT);
#endif
return(info);
}
void PLASMA_DORGQR(int *M, int *N, int *K, double *A, int *LDA, double **T, double *B, int *LDB, int *INFO)
{ *INFO = PLASMA_dorgqr(*M, *N, *K, A, *LDA, *T, B, *LDB); }
