此代码用于以下参考中介绍的强排名显示QR（Strong RRQR）。
参考： Gu，Ming和Stanley C. Eisenstat。 “用于计算强排名显示 QR 分解的有效算法。” SIAM 科学计算杂志 17.4 (1996): 848-869。
给定目标矩阵 A，A 的强 RRQR 的形式为AP = [Q1, Q2] * [R11, R12; 0 R22] \大约 Q1 [R11, R12] 其中 P 是置换矩阵，[Q1,Q2] 具有正交列，R11 是上三角矩阵。 R11 和 R12 需要满足矩阵 inv(R11)*R12 的每个条目都以不小于 1 的预先指定的常数为界。
该代码包含揭示 QR 的强秩，具有 (1) 固定秩，即 R11 的维数； (2)误差阈值，即R22的每一列的范数都小于阈值。 实现严格遵循上述参考文献中提出的算法。
笔记： 代码中主要耗时的步骤通常是对
