code: def _margtimephase_loglr(self, mf_snr, opt_snr):
          
          return special.logsumexp(numpy.log(special.i0(mf_snr)),
                                   b=self._deltat) - 0.5*opt_snr
docs: Returns the log likelihood ratio marginalized over time and phase.