polynomial.c

    int R, togo, size;
//    for(i=0; i<r; i++)
    R=r*c;
    for(i=0; i<R; i=i+8)
    {
        togo = R - i;
        size = (togo > 8) ? 8 : togo;

  //     for(j=0; j<c; j++)
         for(j=0; j<size; j++)

        {
         //   subsref(bt,j,i) = subsref(b,i,j);
            asubsref(bt,j+i) = asubsref(b,i+j);

        }
    }