Verified computation for the matrix Lambert W function

摘要

Two iterative algorithms are proposed for numerically computing interval matrices containing primary matrix Lambert W functions. The first algorithm is based on a numerical spectral decomposition and involves only cubic complexity per iteration. The second algorithm is based on a numerical Jordan decomposition and applicable even for defective matrices. Numerical results show the effectiveness and robustness of the algorithms.