code: def test_cometrix_matrix_is_identity(self, n_points, atol):
      base_point = self.data_generator.random_point(n_points)
  
      res = self.space.metric.cometric_matrix(base_point)
  
      batch_shape = (n_points,) if n_points > 1 else ()
      expected = gs.broadcast_to(
          gs.eye(self.space.dim), batch_shape + 2 * (self.space.dim,)
      )
      self.assertAllClose(res, expected, atol=atol)
  
file_path: geomstats/geomstats/test_cases/geometry/euclidean.py
imports: import pytest
  import geomstats.backend
  from geomstats.test_cases.geometry.base import VectorSpaceTestCase
  from geomstats.test_cases.geometry.mixins import GroupExpTestCaseMixins
  from geomstats.test_cases.geometry.riemannian_metric import RiemannianMetricTestCase
label_desc: Use image node_id 1 for calling the EuclideanMetricTestCase obj's underlying member method code with example usage: obj.test_cometrix_matrix_is_identity(n_points, atol) without return types