code: def test_fit(self, estimator, X, expected, atol):
      estimate = estimator.fit(X).estimate_
      self.assertAllClose(estimate, expected, atol=atol)
  
file_path: geomstats/tests/tests_geomstats/test_learning/test_geometric_median.py
imports: import random
  import pytest
  from geomstats.geometry.euclidean import Euclidean
  from geomstats.geometry.hyperboloid import Hyperboloid
  from geomstats.geometry.hypersphere import Hypersphere
  from geomstats.geometry.spd_matrices import SPDAffineMetric, SPDLogEuclideanMetric, SPDMatrices
  from geomstats.learning.geometric_median import GeometricMedian
  from geomstats.test.parametrizers import DataBasedParametrizer
  from geomstats.test.test_case import TestCase
  from geomstats.test_cases.learning._base import BaseEstimatorTestCase, MeanEstimatorMixinsTestCase
  from .data.geometric_median import GeometricMedianFitTestData, GeometricMedianTestData
label_desc: Use image node_id 1 for calling the TestGeometricMedianFit obj's underlying member method code with example usage: obj.test_fit(estimator, X, expected, atol) without return types