Parallel solution of large Sylvester equations