code: def test_score(self, X_train, y_train, X_test, y_expected, atol):
      self.estimator.fit(X_train, y_train)
  
      score = self.estimator.score(X_test, y_expected)
      self.assertAllClose(score, 1.0, atol=atol)
  
imports: from geomstats.test_cases.learning._base import BaseEstimatorTestCase
label_desc: Use image node_id 5 for calling the RiemannianMinimumDistanceToMeanTestCase obj's underlying member method code with example usage: obj.test_score(X_train, y_train, X_test, y_expected, atol) without return types