We consider the problem of updating an invariant subspace of a Hermitian, large and struc-tured matrix when the matrix is modiied slightly. The problem can be formulated as that of computing stationary values of a certain function, with orthogonality constraints. The constraint is formulated as the requirement that the solution must be on the Grassmann manifold, and Newton's method on the manif...

Both Sylvester matrix and convolution are defined by two polynomials. If one of them has small degree, then the associated Sylvester linear system can be solved fast by using its PFD interpretation of the convolution equation. This can immediately simplify the refinement of approximate convolution by means of Newton’s iteration, where we also incorporate the PFD refinement techniques or alterna...

Global Krylov subspace methods are the most efficient and robust methods to solve generalized coupled Sylvester matrix equation. In this paper, we propose the nested splitting conjugate gradient process for solving this equation. This method has inner and outer iterations, which employs the generalized conjugate gradient method as an inner iteration to approximate each outer iterate, while each...

Presents a recurrent neural network for solving the Sylvester equation with time-varying coefficient matrices. The recurrent neural network with implicit dynamics is deliberately developed in the way that its trajectory is guaranteed to converge exponentially to the time-varying solution of a given Sylvester equation. Theoretical results of convergence and sensitivity analysis are presented to ...

Both Sylvester matrix and convolution are defined by two polynomials. If one of them has small degree, then the associated Sylvester linear system can be solved fast by using its PFD interpretation of the convolution equation. This can immediately simplify the refinement of approximate convolution by means of Newton’s iteration, where we also incorporate the PFD refinement techniques or alterna...

We consider an analytic perturbation of the Sylvester matrix equation. Mainly we are interested in the singular case, that is, when the null space of the unperturbed Sylvester operator is not trivial, but the perturbed equation has a unique solution. In this case, the solution of the perturbed equation can be given in terms of a Laurent series. Here we provide a necessary and su cient condition...