Riemannian Optimization on Relaxed Indicator Matrix Manifold