Efficient Scaling and Squaring Method for Computing the Matrix Exponential within a Given Tolerance
This work presents a new algorithm that efficiently computes the matrix exponential within a given tolerance by incorporating Taylor, partitioned, and classical Padé methods, and selecting the most suitable scheme based on the matrix norm and the desired tolerance.