code: def _compute_rdp(q, sigma, alpha):
      """Compute RDP of the Sampled Gaussian mechanism at order alpha.
  
      Args:
        q: The sampling rate.
        sigma: The std of the additive Gaussian noise.
        alpha: The order at which RDP is computed.
  
      Returns:
        RDP at alpha, can be np.inf.
      """
      if q == 0:
          return 0
  
      if q == 1.0:
          return alpha / (2 * sigma**2)
  
      if np.isinf(alpha):
          return np.inf
  
      return _compute_log_a(q, sigma, alpha) / (alpha - 1)
  
imports: from __future__ import absolute_import
  from __future__ import division
  from __future__ import print_function
  import math
  import sys
  import numpy
  from scipy import special
  import six
label_desc: Use image node_id 14 for calling a global function with example usage: _compute_rdp(q, sigma, alpha) and returns: unknown, int, unknown, np