code: def _squared_dist_grad_point_b(point_a, point_b, metric):
      """Compute gradient of squared_dist wrt point_b.
  
      Compute the Riemannian gradient of the squared geodesic
      distance with respect to the second point point_b.
  
      Parameters
      ----------
      point_a : array-like, shape=[..., dim]
          Point.
      point_b : array-like, shape=[..., dim]
          Point.
      metric : SpecialEuclideanMatrixCanonicalLeftMetric
          Metric defining the distance.
  
      Returns
      -------
      _ : array-like, shape=[..., dim]
          Riemannian gradient, in the form of a tangent
          vector at base point : point_b.
      """
      return -2 * metric.log(point_a, point_b)
  
imports: import geomstats.backend
  import geomstats.errors
  from geomstats.geometry.base import LevelSet
  from geomstats.geometry.general_linear import GeneralLinear
  from geomstats.geometry.matrices import Matrices, MatricesMetric
  from geomstats.geometry.symmetric_matrices import SymmetricMatrices
  from geomstats.vectorization import repeat_out
label_desc: Use image node_id 2 for calling a global function with example usage: _squared_dist_grad_point_b(point_a, point_b, metric) and returns: unknown