Maybe there's another numerical factor involved that I'm not aware of, but I think there may be a small typo here: https://github.com/AstroJacobLi/mrf/blob/1de806ff6d9c022d31ad9b94d9b3b1688c8f9e5c/mrf/sbcontrast.py#L276