code: def test_contains_SU2():
      """Test the functionality to check an operation is in an existing approximate set"""
      op = qml.prod(qml.T(0), qml.T(0), qml.T(0), qml.T(0), qml.T(0))
      target = _SU2_transform(op.matrix())[0]
  
      approx_ids, _, _, approx_vec = _approximate_set(
          ("T", "T*", "H"), max_length=3
      )
  
      exists, quaternion = _contains_SU2(target, approx_vec)
      assert exists
  
      result = [
          qml.adjoint(qml.T(0)),
          qml.adjoint(qml.T(0)),
          qml.adjoint(qml.T(0)),
      ]
  
      node_points = qml.math.array(approx_vec)
      seq_node = qml.math.array([quaternion])
      _, [index] = sp.spatial.KDTree(node_points).query(
          seq_node, workers=-1
      )
  
      assert approx_ids[index] == result
  
file_path: pennylane/tests/ops/op_math/test_solovay_kitaev.py
imports: import math
  import pytest
  import scipy
  import pennylane
  from pennylane.ops.op_math.decompositions.solovay_kitaev import _SU2_transform, _quaternion_transform, _contains_SU2, _approximate_set, _group_commutator_decompose, sk_decomposition
label_desc: Use image node_id 3 for calling a global function with example usage: test_contains_SU2() without return types