code: def test_normalize_is_unitary(self, n_points, atol):
      point = self.data_generator.random_point(n_points)
      vec = self.data_generator.random_tangent_vec(point)
  
      normalized_vec = self.space.metric.normalize(vec, point)
      res = self.space.metric.norm(normalized_vec, point)
  
      batch_shape = (n_points,) if n_points > 1 else ()
      expected = gs.ones(batch_shape)
      self.assertAllClose(res, expected, atol=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 13 for calling the RiemannianMetricTestCase obj's underlying member method code with example usage: obj.test_normalize_is_unitary(n_points, atol) without return types