code: def is_gaussian_prime(num):
      r"""Test if num is a Gaussian prime number.
  
      References
      ==========
  
      .. [1] https://oeis.org/wiki/Gaussian_primes
      """
  
      num = sympify(num)
      a, b = num.as_real_imag()
      a = as_int(a, strict=False)
      b = as_int(b, strict=False)
      if a == 0:
          b = abs(b)
          return isprime(b) and b % 4 == 3
      elif b == 0:
          a = abs(a)
          return isprime(a) and a % 4 == 3
      return isprime(a**2 + b**2)
  
imports: from itertools import count
  from sympy.core.sympify import sympify
  from sympy.external.gmpy import gmpy, gcd, jacobi, is_square, bit_scan1, is_fermat_prp, is_euler_prp, is_selfridge_prp, is_strong_selfridge_prp, is_strong_bpsw_prp
  from sympy.external.ntheory import _lucas_sequence
  from sympy.utilities.misc import as_int, filldedent
label_desc: Use image node_id 15 for calling a global function with example usage: is_gaussian_prime(num) and returns: isprime, unknown, unknown