csapp2e 家庭作业 6.46

void transpose(int *dst,int *src,int dim)
{
 int i,j;
 int im,jm;
 for(i=0;i<dim-3;i+=4)
 {
  im = i*dim;
  for(j=0;j<dim-3;j+=4)
  {
   jm = j*dim;

   dst[jm+i]=src[im+j];
   dst[jm+i+1]=src[im+dim+j];
   dst[jm+i+2]=src[im+2*dim+j];
   dst[jm+i+3]=src[im+3*dim+j];
   dst[jm+dim+i]=src[im+j+1];
   dst[jm+dim+i+1]=src[im+dim+j+1];
   dst[jm+dim+i+2] =src[im+2*dim+j+1];
   dst[jm+dim+i+3]=src[im+3*dim+j+1];
   dst[jm+2*dim+i]=src[im+j+2];
   dst[jm+2*dim+i+1]=src[im+dim+j+2];
   dst[jm+2*dim+i+2]=src[im+2*dim+j+2];
   dst[jm+2*dim+i+3]=src[im+3*dim+j+2];
   dst[jm+3*dim+i]=src[im+dim+j+3];
   dst[jm+3*dim+i+1]=src[im+dim+j+3];
   dst[jm+3*dim+i+2]=src[im+2*dim+j+3];
   dst[jm+3*dim+i+3]=src[im+3*dim+j+3];
  }
 }
 for(;i<dim;i++)
  for(;j<dim;j++)
   dst[j*dim+i]=src[i*dim+j];
}

 

 

用10000*10000的数组测试（linux下需要取消stack的大小限制才运行这么大的数组）

书上的程序平均需要 9.301s

优化后评价需要 3.830秒