def METHOD_NAME(par):
    """Solve marchenko equations using input Rs in time domain with multiple
    points
    """
    MarchenkoWM = Marchenko(
        R, dt=dt, dr=dr, nfmax=nfmax, wav=wav, toff=toff, nsmooth=nsmooth
    )
    _, _, g_inv_minus, g_inv_plus = MarchenkoWM.apply_multiplepoints(
        trav_multi,
        nfft=2**11,
        rtm=False,
        greens=True,
        dottest=True,
        **dict(iter_lim=par["niter"], show=0)
    )
    ginvsub = (g_inv_minus + g_inv_plus)[:, 1, nt - 1 :].T
    ginvsub_norm = ginvsub / ginvsub.max()
    gsub_norm = gsub / gsub.max()
    assert np.linalg.norm(gsub_norm - ginvsub_norm) / np.linalg.norm(gsub_norm) < 1e-1