class_docstr: Class that defines the canonical metric for Stiefel manifolds.
docstr: Compute the radius of the injectivity domain.
  
  This is is the supremum of radii r for which the exponential map is a
  diffeomorphism from the open ball of radius r centered at the base
  point onto its image.
  In this case the exact injectivity radius is not known, and we use here
  a lower bound given by [Rentmeesters2015]_.
  
  Parameters
  ----------
  base_point : array-like, shape=[..., n, p]
      Point on the manifold.
  
  Returns
  -------
  radius : array-like, shape=[...,]
      Injectivity radius.
  
  References
  ----------
  .. [Rentmeesters2015] Rentmeesters, Quentin. “Algorithms for Data
      Fitting on Some Common Homogeneous Spaces.” UCL - Université
      Catholique de Louvain, 2013.
      https://dial.uclouvain.be/pr/boreal/object/boreal:132587.
code: def injectivity_radius(self, base_point):
      """Compute the radius of the injectivity domain.
  
      This is is the supremum of radii r for which the exponential map is a
      diffeomorphism from the open ball of radius r centered at the base
      point onto its image.
      In this case the exact injectivity radius is not known, and we use here
      a lower bound given by [Rentmeesters2015]_.
  
      Parameters
      ----------
      base_point : array-like, shape=[..., n, p]
          Point on the manifold.
  
      Returns
      -------
      radius : array-like, shape=[...,]
          Injectivity radius.
  
      References
      ----------
      .. [Rentmeesters2015] Rentmeesters, Quentin. “Algorithms for Data
          Fitting on Some Common Homogeneous Spaces.” UCL - Université
          Catholique de Louvain, 2013.
          https://dial.uclouvain.be/pr/boreal/object/boreal:132587.
      """
      radius = gs.array(0.89 * gs.pi)
      return repeat_out(self._space.point_ndim, radius, base_point)
  
imports: import warnings
  import geomstats.backend
  import geomstats.errors
  from geomstats import algebra_utils
  from geomstats.geometry.base import LevelSet
  from geomstats.geometry.hermitian_matrices import powermh
  from geomstats.geometry.matrices import Matrices
  from geomstats.geometry.riemannian_metric import RiemannianMetric
  from geomstats.vectorization import repeat_out
label_desc: Use image node_id 8 for calling the StiefelCanonicalMetric obj's underlying member method code with example usage: obj.injectivity_radius(base_point) and returns: repeat_out