switched to a bsnr computation which more closely resembles what's shown in the paper (although the same computationally)