code: def test_dist_random(self, n_points, atol):
      point_a = self.data_generator.random_point(n_points)
      point_b = self.data_generator.random_point(n_points)
  
      self.test_dist(point_a, point_b, atol)
  
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 20 for calling the RiemannianMetricComparisonTestCase obj's underlying member method code with example usage: obj.test_dist_random(n_points, atol) without return types