Abstract

We describe a recurrence method for computing primary $p$th roots of a matrix $A$ with a cost, in terms of elementary arithmetic operations and memory, which is logarithmic with respect to $p$. When $A$ is real and the primary root is real as well, the algorithm is based on the real Schur form of $A$ and uses real arithmetic. The numerical experiments confirm the good behavior of the new algorithm in finite arithmetic. The case of arbitrary fractional powers of $A$ is also considered.

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 call

Disclaimer: 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.