code: def _expected_sphere_inner_product_derivativ_matrix(base_point):
      theta = base_point[0]
  
      # derivative with respect to theta
      expected_1 = gs.array(
          [[0, 0], [0, 2 * gs.cos(theta) * gs.sin(theta)]]
      )
      # derivative with respect to phi
      expected_2 = gs.zeros((2, 2))
  
      return expected_1, expected_2
  
imports: import geomstats.backend
  from geomstats.test.data import TestData
  from .riemannian_metric import RiemannianMetricComparisonTestData, RiemannianMetricTestData
label_desc: Use image node_id 10 for calling a global function with example usage: _expected_sphere_inner_product_derivativ_matrix(base_point) and returns: expected_1,  expected_2