algorithm - 计算实矩阵的实对数的算法？

Question

我想要一种算法来找到实矩阵的实对数。

这里发现的存在定理是：

一个实数矩阵有一个实数对数当且仅当它是可逆的并且属于负特征值的每个 Jordan 块出现偶数次。

我在 MATLAB 中工作，但任何适应性算法都可以。下面的示例矩阵满足这两个条件。

A = [-1 1 0 0; 0 -1 0 0; 0 0 -1 1; 0 0 0 -1];
RealMatrixLogarithm?(A) = X;

编辑：我试图理解本文中给出的算法，但它似乎只对正常矩阵有效。

感谢您的任何帮助。