Method mcopyunblock

src/ialglib.cpp:537–570

This subroutine copies matrix from aligned contigous storage to non- aligned non-contigous storage A: * MxN * aligned * contigous * stride is alglib_r_block * may be transformed during copying (as prescribed by op) B: * alglib_r_block*alglib_r_block (only MxN/NxM submatrix is used) * non-aligned, non-contigous Transformation types: * 0 - no transform * 1 - transposition ***********************

* 1 - transposition
********************************************************************/
void ialglib::mcopyunblock(int m, int n, const double *a, int op, double *b, int stride)
{
    int i, j, n2;
    const double *psrc;
    double *pdst;
    if( op==0 )
    {
        n2 = n/2;
        for(i=0,psrc=a; i<m; i++,a+=alglib_r_block,b+=stride,psrc=a)
        {
            for(j=0,pdst=b; j<n2; j++,pdst+=2,psrc+=2)
            {
                pdst[0] = psrc[0];
                pdst[1] = psrc[1];
            }
            if( n%2!=0 )
                pdst[0] = psrc[0];
        }
    }
    else
    {
        n2 = n/2;
        for(i=0,psrc=a; i<m; i++,a++,b+=stride,psrc=a)
        {
            for(j=0,pdst=b; j<n2; j++,pdst+=2,psrc+=alglib_twice_r_block)
            {
                pdst[0] = psrc[0];
                pdst[1] = psrc[alglib_r_block];
            }
            if( n%2!=0 )
                pdst[0] = psrc[0];
        }
    }
}


/********************************************************************