用循环平铺从矩阵中提取子矩阵

Question

我有以下矩阵4x4：

1  2  3  4

5  6  7  8

9 10 11 12

13 14 15 16

我想要提取并存储（在一些新闻变量中）以下四个子矩阵2x2：

[1 2

 5 6]

[3 4

 7 8]

[9 10

 13 14]

[11 12

 15 16]

这就像openCV的“Rect”（http://docs.opencv.org/java/org/opencv/core/Rect.html）函数，但我不想使用OpenCV。

我必须使用并行编译器，所以我想用文献中存在的着名循环变换来提取子矩阵：“循环平铺”（也称为“循环阻塞”或“循环展开和阻塞”或“循环剥离和交换”）。 - （http://en.wikipedia.org/wiki/Loop_tiling）

有可能吗？

Answer 1

有可能吗？

当然......

    int n = 4;
    int matrix[4][4] = {1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16};
    int submatrixes[n/2*n/2][2][2];
    int i, j, x, y, z;
    for (z = i = 0; i < n; i += 2)
      for (j = 0; j < n; j += 2, ++z)
        for (x = 0; x < 2; x++)
          for (y = 0; y < 2; y++)
            submatrixes[z][x][y] = matrix[i+x][j+y];