Matlab和Python读取二进制文件的方式不同

Question

我正在将相同的二进制文件读入Python和Matlab并将其放入一个矩阵中。当我取这个矩阵的范数，我得到不同的结果。在

我使用各自的smatload函数来加载二进制文件。在

Python：

def smatload(filename):
    #print 'opening: ', filename
    f = open(filename, 'rb')
    m = np.fromfile(f,'q',1)
    n = np.fromfile(f,'q',1)
    nnz = np.fromfile(f,'q',1)
    print 'reading %d x %d with %d non-zeros' % (m,n,nnz)
    S = np.fromfile(f,'d',3*nnz)
    f.close()
    S = S.reshape((nnz,3))
    rows = S[:,0].astype(int) - 1
    cols = S[:,1].astype(int) - 1
    vals = S[:,2]
    return csr_matrix((vals,(rows,cols)),shape=(m,n))

Matlab软件：

^{pr2}$

对于返回矩阵的范数，我得到的结果是

Matlab：标准值（A.“fro”）=0.018317077159881

Python：np.linalg.norm公司（A） =0.018317077159760

我已经确认它们都读取了正确的值数目（6590x7126矩阵，122526个非零），并且我对两者都使用了相同的范数（frobenius）。在

你知道是什么原因造成的吗？在

Answer 1

快速查看Frobenius规范可以发现，它需要将所有值平方并将它们相加。在

由于read命令中有uint64，因此看起来您可能正在填充浮点存储。当你把两个二进制数相乘时，需要两倍的位来存储答案。这意味着您需要128位来存储所有的十进制值。如果Python和MATLAB的做法不同，这就可以解释为什么十进制值不同。在

有关MATLAB和Python如何处理浮点精度的信息，请参见以下两个链接：

Python： https://docs.python.org/2/tutorial/floatingpoint.html

MATLAB软件： http://blogs.mathworks.com/cleve/2014/07/07/floating-point-numbers/