python，测试矩阵中的所有值是否小于 numpy.finfo(float).eps

Question

我想测试矩阵的所有元素是否都为零。我有一个 3x3 矩阵，其中包含非常小但非零的元素。像这样:

my_matrix = numpy.mat(
[[  1.32439783e-16,  -1.09384837e-16,   5.22327024e-17],
[ -1.47879836e-17,   1.58560316e-17,  -8.18983137e-19],
[  3.40590503e-17,  -1.36768936e-17,   3.33256165e-17]]
)

我通常会使用:

np.array_equal( my_matrix, np.zeros((3,3)) )

但这会返回 False。我的印象是，由于该矩阵中所有元素的绝对值都小于 numpy.finfo(float).eps，其中 1.0 + eps = 1.0，所以这个应该返回 True。

有没有办法测试矩阵中的所有值是否都小于 numpy.finfo(float).eps？

顺便说一下，numpy.finfo(float).eps = 2.22044604925e-16。

Answer 1

您正在比较两个 float 垫，请改用 np.allclose:

In [119]: np.allclose( my_matrix, np.zeros((3,3)) )
Out[119]: True