fix for all-zero in both signals

    if np.all(diff_sig_norm == 0):

        return np.asarray([np.inf] * ref_sig_norm.shape[1])

    channels_identical_idx = np.sum(np.abs(diff_sig_norm), axis=0) == 0

    denom_add = 10**-13 * len_seg

    segment_counter = np.zeros(ref_sig.shape[1])

    for ref_seg, diff_seg in zip(

        )

    ssnr[segment_counter == 0] = -np.inf

    # this prevents all-zero channels in both signals to be reported as -inf

    ssnr[channels_identical_idx] = np.inf

    return ssnr