code: def compute_rdp_threshold(log_pr_answered, sigma, orders):
      logq = min(log_pr_answered, _log1mexp(log_pr_answered))
      # The input to the threshold mechanism has stability 1, compared to
      # GNMax, which has stability = 2. Hence the sqrt(2) factor below.
      return rdp_gaussian(logq, 2**0.5 * sigma, orders)
  
imports: from __future__ import absolute_import
  from __future__ import division
  from __future__ import print_function
  import math
  from absl import app
  import numpy
  import scipy.stats
label_desc: Use image node_id 10 for calling a global function with example usage: compute_rdp_threshold(log_pr_answered, sigma, orders) and returns: rdp_gaussian