def METHOD_NAME(fdemType, comp, src):
    prb = getFDEMProblem(fdemType, comp, [src], freq)
    # prb.solverOpts = dict(check_accuracy=True)
    print(f"Adjoint {fdemType} formulation - {src} - {comp}")
    m = np.log(np.ones(prb.sigmaMap.nP) * CONDUCTIVITY)
    mu = np.ones(prb.mesh.nC) * MU
    if addrandoms is True:
        m = m + np.random.randn(prb.sigmaMap.nP) * np.log(CONDUCTIVITY) * 1e-1
        mu = mu + np.random.randn(prb.mesh.nC) * MU * 1e-1
    survey = prb.survey
    u = prb.fields(m)
    v = np.random.rand(survey.nD)
    w = np.random.rand(prb.mesh.nC)
    vJw = v.dot(prb.Jvec(m, w, u))
    wJtv = w.dot(prb.Jtvec(m, v, u))
    tol = np.max([TOL * (10 ** int(np.log10(np.abs(vJw)))), FLR])
    print(vJw, wJtv, vJw - wJtv, tol, np.abs(vJw - wJtv) < tol)
    return np.abs(vJw - wJtv) < tol