1ffc5ad914da39a369337c1606a8ba788dbf9055,geomstats/geometry/grassmannian.py,Grassmannian,random_uniform,#Grassmannian#Any#,55

Before Change

        Returns
        -------
        
        basis_change = SpecialOrthogonal(self.n).random_uniform(n_samples)
        projector = from_vector_to_diagonal_matrix(
            gs.array([1.] * self.k + [0.] * (self.n - self.k)))
        return Matrices.mul(
            basis_change, projector, GeneralLinear.inverse(basis_change))

    def belongs(self, point, tolerance=TOLERANCE):
        Check if the point belongs to the manifold.

After Change

        New York: Springer-Verlag. 2003, 10.1007/978-0-387-21540-2
        
        points = gs.random.normal(size=n_samples * self.k * self. n)
        points = gs.reshape(points, (n_samples, self.n, self.k))
        full_rank = Matrices.mul(Matrices.transpose(points), points)
        projector = Matrices.mul(
            points,
            GeneralLinear.inverse(full_rank),
            Matrices.transpose(points))
        return projector[0] if n_samples == 1 else projector

    def belongs(self, point, tolerance=TOLERANCE):
        Check if the point belongs to the manifold.

Instances