code: def FermionicSWAP(phi):
      r"""Fermionic SWAP rotation gate.
  
      Args:
          phi (float): rotation angle :math:`\phi`
      Returns:
          array[complex]: unitary 4x4 rotation matrix
      """
      c = math.cos(phi / 2)
      s = math.sin(phi / 2)
      g = math.exp(1j * phi / 2)
      p = math.exp(1j * phi)
      return math.array(
          [
              [1, 0, 0, 0],
              [0, g * c, -1j * g * s, 0],
              [0, -1j * g * s, g * c, 0],
              [0, 0, 0, p],
          ]
      )
  
imports: import pennylane
  from pennylane import numpy
  from pennylane import math
label_desc: Use image node_id 31 for calling a global function with example usage: FermionicSWAP(phi) and returns: math