def METHOD_NAME(Nphi, Ntheta, Nr, dtype, dealias, build_basis, func):
    c, d, b, phi, theta, r, x, y, z = build_basis(Nphi, Ntheta, Nr, dtype, dealias)
    f = field.Field(dist=d, bases=(b,), dtype=dtype)
    if func is np.arccosh:
        f['g'] = 1 + r**2
        a = 2
    else:
        f['g'] = r**2
        a = 0.5
    g0 = func(a*f['g'])
    g = func(a*f).evaluate()
    assert np.allclose(g['g'], g0)