code: def test_score_at_train_is_one(self, n_samples, atol):
      X, y = self.data_generator.random_dataset(n_samples)
  
      self.estimator.fit(X, y)
  
      res = self.estimator.score(X, y)
      self.assertAllClose(res, 1, atol=atol)
  
file_path: geomstats/geomstats/test_cases/learning/wrapped_gaussian_process.py
imports: import pytest
  import geomstats.backend
  from geomstats.test.random import RandomDataGenerator
  from geomstats.test_cases.learning._base import BaseEstimatorTestCase
label_desc: Use image node_id 2 for calling the WrappedGaussianProcessTestCase obj's underlying member method code with example usage: obj.test_score_at_train_is_one(n_samples, atol) without return types