code: def _margdistphase_loglr(self, mf_snr, opt_snr):
          
          logl = numpy.log(special.i0(mf_snr))
          logl_marg = logl/self._dist_array
          opt_snr_marg = opt_snr/self._dist_array**2
          return special.logsumexp(logl_marg - 0.5*opt_snr_marg,
                                   b=self._deltad*self.dist_prior)
docs: Returns the log likelihood ratio marginalized over distance and
          phase.