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

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];