We propose a simple technique that, if combined with algorithms for computing functions of triangular matrices, can make them more efficient. Basically, such a technique consists in a specific scaling similarity transformation that reduces the departure from normality of a triangular matrix, thus decreasing its norm and in general its function condition number. It can easily be extended to non-triangular matrices, provided that it is combined with algorithms involving a prior Schur decomposition. Situations where the technique should be used or not will be discussed in detail. Special attention is devoted to particular algorithms like the inverse scaling and squaring to the matrix logarithm and the scaling and squaring to the matrix exponential. The advantages of our proposal are supported by theoretical results and illustrated with numerical experiments.
翻译:我们建议一种简单的方法,如果与三角矩阵计算函数的算法相结合,可以提高这些算法的效率。基本上,这种技术包括具体的缩放相似性转换,减少三角矩阵偏离正常性的情况,从而降低其规范,并一般地降低其功能条件编号。它很容易推广到非三角矩阵,只要它与先前Schur分解的算法相结合。在使用或不使用该技术的情况下,将进行详细讨论。特别注意特定的算法,例如反缩放和对齐矩阵对数以及缩放和对齐矩阵指数。我们提案的优点得到理论结果的支持,并以数字实验加以说明。