谁能提供求矩阵的秩的算法，高分相送。

矩阵保存在一个double型的二维数组中。
我以前的想法是，进行初等变换，将其化成简化行梯形，然后数其中的主1的个数，但是变换中难免乘除运算，因为在运算中后面的几位被丢掉，如在变换过程中有一个元素为0.0001,很难判断其是正确的结果，还是误差。这对求线性方程组问题不大，但对求秩就有很大影响。不知该怎样解决。