Abstract
This paper is concerned with a relative perturbation theory and its entrywise relatively accurate numerical solutions of an M-matrix Sylvester equation AX + XB = C by which we mean both A and B have positive diagonal entries and nonpositive off-diagonal entries and $${P=I_m \otimes A+B^{\rm T} \otimes I_n}$$ is a nonsingular M-matrix, and C is entrywise nonnegative. It is proved that small relative perturbations to the entries of A, B, and C introduce small relative errors to the entries of the solution X. Thus the smaller entries of X do not suffer bigger relative errors than its larger entries, unlikely the existing perturbation theory for (general) Sylvester equations. We then discuss some minor but crucial implementation changes to three existing numerical methods so that they can be used to compute X as accurately as the input data deserve.
Sign-in/Register to access full text options 
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a callDisclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.
