code: def sectional_curvature(tangent_vec_a, tangent_vec_b, base_point):
      x, y = base_point[..., 0], base_point[..., 1]
      detg = gs.polygamma(1, x) * gs.polygamma(1, y) - gs.polygamma(
          1, x + y
      ) * (gs.polygamma(1, x) + gs.polygamma(1, y))
      return (
          gs.polygamma(2, x)
          * gs.polygamma(2, y)
          * gs.polygamma(2, x + y)
          * (
              gs.polygamma(1, x) / gs.polygamma(2, x)
              + gs.polygamma(1, y) / gs.polygamma(2, y)
              - gs.polygamma(1, x + y) / gs.polygamma(2, x + y)
          )
          / (4 * detg**2)
      )
  
file_path: geomstats/geomstats/test_cases/information_geometry/beta.py
imports: import pytest
  import geomstats.backend
  from geomstats.geometry.symmetric_matrices import SymmetricMatrices
  from geomstats.test_cases.information_geometry.dirichlet import DirichletDistributionsTestCase, DirichletMetricTestCase
label_desc: Use image node_id 1 for calling a global function with example usage: sectional_curvature(tangent_vec_a, tangent_vec_b, base_point) and returns: unknown