code: def test_to_grassmannian_belongs_to_grassmannian(
      self, n_points, atol
  ):
      point = self.space.random_point(n_points)
      grass_point = self.space.to_grassmannian(point)
  
      grassmannian = Grassmannian(
          self.space.n, self.space.p, equip=False
      )
      res = grassmannian.belongs(grass_point, atol=atol)
      expected = gs.ones(n_points, dtype=bool)
  
      self.assertAllEqual(res, expected)
  
imports: import pytest
  import geomstats.backend
  from geomstats.geometry.grassmannian import Grassmannian
  from geomstats.test.test_case import TestCase
  from geomstats.test.vectorization import generate_vectorization_data
  from geomstats.test_cases.geometry.base import LevelSetTestCase
  from geomstats.test_cases.geometry.riemannian_metric import RiemannianMetricTestCase
label_desc: Use image node_id 2 for calling the StiefelTestCase obj's underlying member method code with example usage: obj.test_to_grassmannian_belongs_to_grassmannian(n_points, atol) without return types