算法894

算法894:基于sep逆估计的ɕ-函数的块Schur-Parlett算法。fortran95软件用于计算指数积分器中所需的ψ-函数的矩阵值。库中的子例程接受一个完整的、对角的或上拟三角矩阵作为它们的参数,其中包含四种精度之一的实项或复杂项。实现了两种不同的算法,一种是缩放和平方法,另一种是改进的分块Schur—Parlett算法。在后一种算法中,基于sep逆函数的估计,对参数应用递归的三乘三分块。sep逆函数的估计是由Hager-Higham估计器实现的,它是LAPACK中的子程序xLACON实现的。文中给出了我们对block-Schur-Parlett算法的改进,并给出了数值实验结果。

这个软件也是同行评审按日记帐汤姆斯.