code: def is_solenoidal(field):
      """
      Checks if a field is solenoidal.
  
      Parameters
      ==========
  
      field : Vector
          The field to check for solenoidal property
  
      Examples
      ========
  
      >>> from sympy.physics.vector import ReferenceFrame
      >>> from sympy.physics.vector import is_solenoidal
      >>> R = ReferenceFrame('R')
      >>> is_solenoidal(R[1]*R[2]*R.x + R[0]*R[2]*R.y + R[0]*R[1]*R.z)
      True
      >>> is_solenoidal(R[1] * R.y)
      False
  
      """
  
      # Field is solenoidal irrespective of frame
      # Take the first frame in the result of the separate method in Vector
      if field == Vector(0):
          return True
      frame = list(field.separate())[0]
      return divergence(field, frame).simplify() is S.Zero
  
file_path: sympy/sympy/physics/vector/fieldfunctions.py
imports: from sympy.core.function import diff
  from sympy.core.singleton import S
  from sympy.integrals.integrals import integrate
  from sympy.physics.vector import Vector, express
  from sympy.physics.vector.frame import _check_frame
  from sympy.physics.vector.vector import _check_vector
label_desc: Use image node_id 5 for calling a global function with example usage: is_solenoidal(field) and returns: unknown, True