code: def _coeff_smooth(lam):
      xi = 1 - 96 * lam + 24 * lam * cupy.sqrt(3 + 144 * lam)
      omeg = cupy.arctan2(cupy.sqrt(144 * lam - 1), cupy.sqrt(xi))
      rho = (24 * lam - 1 - cupy.sqrt(xi)) / (24 * lam)
      rho = rho * cupy.sqrt(
          (48 * lam + 24 * lam * cupy.sqrt(3 + 144 * lam)) / xi
      )
      return rho, omeg
  
imports: import cupy
  import cupyx.scipy.ndimage
  from cupyx.scipy.signal._iir_utils import apply_iir_sos
  from cupyx.scipy.signal._splines import _symiirorder1_nd, _symiirorder2_nd
  from cupyx.scipy.interpolate._bspline import BSpline
  import numpy
label_desc: Use image node_id 4 for calling a global function with example usage: _coeff_smooth(lam) and returns: rho,  omeg