def dotestAdjoint(prb, mesh, sigma):
    m = np.random.rand(prb.mapping.nP)
    d = np.random.rand(prb.survey.nD)

    V1 = d.dot(prb.Jvec(sigma, m))
    V2 = m.dot(prb.Jtvec(sigma, d))
    print 'AdjointTest', V1, V2
    return np.abs(V1-V2)/np.abs(V1), 1e-6
dataset/ETHPy150Open simpeg/simpeg/tests/em/tdem/test_TDEM_combos.py/dotestAdjoint/original