Bug in uncertainty propagation formulas, within MRF.py. Not using cross-terms of co-variance matrix.
The error propagation formula on lines 654, 655, 658, 660, and 661 in MRF.py is incorrect. Currently it doesn't correctly utilise the cross terms, generally over-estimating the error.
To correct change einsum(str("ij,ii,ij->j"), *args)
for np.einsum(str("ij,ik,kj->j"), *args)
.
The same error is repeated in the class method residuals_error. See lines 677, 678, 681, 683, and 684.