Function cmatrixqrunpackq

src/ortfac.cpp:1004–1122

Partial unpacking of matrix Q from QR decomposition of a complex matrix A. Input parameters: A - matrices Q and R in compact form. Output of CMatrixQR subroutine . M - number of rows in matrix A. M>=0. N - number of columns in matrix A. N>=0. Tau - scalar factors which are used to form Q. Output

     Bochkanov Sergey
*************************************************************************/
void cmatrixqrunpackq(const ap::complex_2d_array& a,
     int m,
     int n,
     const ap::complex_1d_array& tau,
     int qcolumns,
     ap::complex_2d_array& q)
{
    ap::complex_1d_array work;
    ap::complex_1d_array t;
    ap::complex_1d_array taubuf;
    int minmn;
    int refcnt;
    ap::complex_2d_array tmpa;
    ap::complex_2d_array tmpt;
    ap::complex_2d_array tmpr;
    int blockstart;
    int blocksize;
    int rowscount;
    int i;
    int j;
//    int k;
    ap::complex v;

    ap::ap_error::make_assertion(qcolumns<=m, "UnpackQFromQR: QColumns>M!");
    if( m<=0||n<=0 )
    {
        return;
    }
    
    //
    // init
    //
    minmn = ap::minint(m, n);
    refcnt = ap::minint(minmn, qcolumns);
    work.setlength(ap::maxint(m, n)+1);
    t.setlength(ap::maxint(m, n)+1);
    taubuf.setlength(minmn);
    tmpa.setlength(m, ablascomplexblocksize(a));
    tmpt.setlength(ablascomplexblocksize(a), ablascomplexblocksize(a));
    tmpr.setlength(2*ablascomplexblocksize(a), qcolumns);
    q.setlength(m, qcolumns);
    for(i = 0; i <= m-1; i++)
    {
        for(j = 0; j <= qcolumns-1; j++)
        {
            if( i==j )
            {
                q(i,j) = 1;
            }
            else
            {
                q(i,j) = 0;
            }
        }
    }
    
    //
    // Blocked code