code: def run_riemannian_relaxation(laplacian, initial_guess,
                                intrinsic_dim, relaxation_kwds):
      
      n, s = initial_guess.shape
      relaxation_kwds = initialize_kwds(relaxation_kwds, n, s, intrinsic_dim)
      if relaxation_kwds['save_init']:
          directory = relaxation_kwds['backup_dir']
          np.save(os.path.join(directory, 'Y0.npy'),initial_guess)
          sp.io.mmwrite(os.path.join(directory, 'L_used.mtx'),
                        sp.sparse.csc_matrix(laplacian))
  
      lossf = relaxation_kwds['lossf']
      return RiemannianRelaxation.init(lossf, laplacian, initial_guess,
                                       intrinsic_dim, relaxation_kwds)
docs: Helper function for creating a RiemannianRelaxation class.