code: def test_metric_matrix_is_spd(self, n_points, atol):
      base_point = self.data_generator.random_point(n_points)
  
      metric_matrix = self.space.metric.metric_matrix(base_point)
  
      res = SPDMatrices(n=math.prod(self.space.shape)).belongs(
          metric_matrix, atol=atol
      )
      expected_shape = get_batch_shape(
          self.space.point_ndim, base_point
      )
      expected = gs.ones(expected_shape, dtype=bool)
  
      self.assertAllEqual(res, expected)
  
imports: import math
  import pytest
  import geomstats.backend
  from geomstats.geometry.spd_matrices import SPDMatrices
  from geomstats.test_cases.geometry.connection import ConnectionComparisonTestCase, ConnectionTestCase
  from geomstats.vectorization import get_batch_shape
label_desc: Use image node_id 2 for calling the RiemannianMetricTestCase obj's underlying member method code with example usage: obj.test_metric_matrix_is_spd(n_points, atol) without return types