def METHOD_NAME(problem):
    var = problem.create_variables(['us'])['us']
    val = nm.arange(var.n_dof, dtype=var.dtype)
    var.set_data(val)
    val1 = problem.evaluate('ev_grad.i.Omega(us)', us=var, mode='el_avg')
    tst.report('ev_grad(el_avg): min, max:', val1.min(), val1.max())
    aux = problem.evaluate('ev_grad.i.Omega(us)', us=var, mode='qp')
    val2 = _integrate(var, aux)
    val2.shape = val1.shape
    tst.report('ev_grad(qp): min, max:', val2.min(), val2.max())
    ok = tst.compare_vectors(val1, val2,
                             label1='de mode',
                             label2='dq mode')
    if not ok:
        tst.report('failed')
    assert ok