code: def test_solenoidal():
      assert is_solenoidal(0) is True
      assert is_solenoidal(R.x) is True
      assert is_solenoidal(2 * R.x + 3 * R.y + 4 * R.z) is True
      assert (
          is_solenoidal(
              R[1] * R[2] * R.x + R[0] * R[2] * R.y + R[0] * R[1] * R.z
          )
          is True
      )
      assert is_solenoidal(R[1] * R.y) is False
      assert is_solenoidal(grad_field) is False
      assert is_solenoidal(curl_field) is True
      assert is_solenoidal((-2 * R[1] + 3) * R.z) is True
      assert (
          is_solenoidal(cos(q) * R.x + sin(q) * R.y + cos(q) * P.z)
          is True
      )
      assert is_solenoidal(R[2] * P.x + P[0] * R.z) is True
  
file_path: sympy/sympy/physics/vector/tests/test_fieldfunctions.py
imports: from sympy.core.singleton import S
  from sympy.core.symbol import Symbol
  from sympy.functions.elementary.trigonometric import cos, sin
  from sympy.physics.vector import ReferenceFrame, Vector, Point, dynamicsymbols
  from sympy.physics.vector.fieldfunctions import divergence, gradient, curl, is_conservative, is_solenoidal, scalar_potential, scalar_potential_difference
  from sympy.testing.pytest import raises
label_desc: Use image node_id 5 for calling a global function with example usage: test_solenoidal() without return types