Te propongo este otro planteamiento, que utiliza la recursividad.<br>

Los parámetros de la función son:<br>
**x, y**: posición, dentro de Data, de la esquina superior izquierda de la submatriz a considerar.<br>
**size**: ancho de la submatriz a considerar. Suponemos que es cuadrada.<br>

    int Data[][]; /* Matriz general, principal, o como se quiera decir. */
    
    int equisuma( int x, int y, int size ) {
      int result = 0;
      int sub = 0;
      int subSize = 0;

      if( size == 2 ) {
        /* Devolvemos la suma de los 4 elementos. */
        result = Data[x][y] +
                 Data[x+1][y] +
                 Data[x][y+1] +
                 Data[x+1][y+1];
      } else {
        /* Tamaño de las submatrices. */
        subSize = size / 2;

        result = equisuma( x, y, subSize );
        
        if( result != -1 ) {
          sub = equisuma( x + subSize, y, subSize );
          if( ( sub == -1 ) || ( sub != result ) )
            result = -1;
        }

        if( result != -1 ) {
          sub = equisuma( x, y + subSize, subSize );
          if( ( sub == -1 ) || ( sub != result ) )
            result = -1;
        }

        if( result != -1 ) {
          sub = equisuma( x + subSize, y + subSize, subSize );
          if( ( sub == -1 ) || ( sub != result ) )
            result = -1;
        }
      }

      return result;
    }

    void main( ) {
      equisuma( 0, 0, AnchoDeLaMatriz );
    }

Como verás, utiliza el clásico esquema de 'divide y vencerás', y creo que se ajusta al enunciado. Tendrás que ajustarlo a los datos exactos de tu problema.