code: def geodesic_discrete_test_data(self):
      expected = gs.array(
          [
              [
                  [-0.27215095, -0.90632948, 0.32326574],
                  [0.13474934, 0.0, 1.39415392],
              ],
              [
                  [-0.13415037, -0.9507542, 0.27941033],
                  [0.17027678, -0.03133438, 1.24635354],
              ],
              [
                  [0.01180225, -0.96796946, 0.25079044],
                  [0.21334916, -0.0301587, 1.08786833],
              ],
              [
                  [0.15891441, -0.95795476, 0.23889097],
                  [0.2547784, 0.0, 0.92187986],
              ],
          ]
      )
      data = [
          dict(
              initial_point=self.pu2,
              end_point=self.pu3,
              expected=expected,
              atol=6e-06,
          )
      ]
      return self.generate_tests(data)
  
imports: import geomstats.backend
  from geomstats.test.data import TestData
label_desc: Use image node_id 4 for calling the SasakiMetricSphereTestData obj's underlying member method code with example usage: obj.geodesic_discrete_test_data() and returns: self