Implement the function safeCompare in getRateThroughTimeMatrix.R to handle numerical precision issues. safeCompare.R now also in BAMMtools.