标签： python numpy matrix scipy
我有一个非常大的酉矩阵（15625 x 15625），具有相当简单的块对角结构：它将身份矩阵与非平凡矩阵交替。所有这些块都有不同的大小，所以我无法预测非平凡矩阵的确切位置。
我想提取一个数组中的所有非平凡子矩阵，以及另一个整数数组，告诉我身份矩阵之间有多大。
我怀疑由于大小是不可预测的，我必须检查对角线中的1是否对应于其行和列中的所有零（在这种情况下它属于单位矩阵），然后一旦我确定了身份矩阵是，我可以采取剩余的块。但我对Python很新，可能我错过了一个更简单的解决方案。