Strassen算法的矩阵划分

Question

假设我有一个NxN矩阵，其中充满了从1到10的随机整数。现在，我想调用PROC(A(1:n/2, 1:n/2)+A(n/2+1:n, n/2+1:n)...，其中n是矩阵的大小。换句话说，我想做一个子矩阵，从A的第一行和第一列开始，一直到A大小的一半，然后把它加到一个子矩阵中，这个子矩阵从A的一半大小加1开始，一直到A的末尾。

我使用的分区函数是这样的：

public Matrix partition(int rowStart, int rowEnd, int colStart, int colEnd) {
    // int r = 0;
    // int c = 0;
    if (this.N%2 != 0) throw new RuntimeException("Illegal matrix dimensions.");
    Matrix C = new Matrix((this.N)/2);
    for (int i=rowStart-1; i<rowEnd; i++) {
        for (int j=colStart-1; j<colEnd; j++) {
            C.data[i][j] = this.data[i][j];
            // C.data[r][c] = this.data[i][j];
            c++;
        }
        r++;
    }
    return C;
}

现在，这适用于查找给定矩阵左上角(Matrix C = m.partition(1, m.size()/2, 1, m.size()/2);)的子矩阵。

 9.00     5.00     0.00     3.00
 0.00     7.00     8.00     3.00
 9.00     3.00    10.00     8.00
 0.00     6.00     2.00     0.00

 9.00     5.00
 0.00     7.00

但是当我尝试得到另一个子矩阵(Matrix D = m.partition(m.size()/2+1, m.size(), m.size()/2+1, m.size());)时，我得到了一个ArrayIndexOutOfBoundsException: 2。我尝试将单独的行和列计数器添加到我的分区函数，但它给出了相同的错误。我如何修改我的分区函数来处理所有的输入，并且仍然给出正确的输出？

Answer 1

C.data[i][j] = this.data[i][j];<-- Culprit

i为rowStart-1，j为rowEnd-1对于C.data，您需要i和j从0开始

for (int i=rowStart-1,p=0; i<rowEnd; i++,p++) {
        for (int j=colStart-1,q=0; j<colEnd; j++,q++) {
            C.data[p][q] = this.data[i][j];
        }
    }