A method for computing hyperbolic matrix functions based on Hermite matrix polynomial expansions is presented.
A new error bound analysis is given.
A theoretical estimate for the optimal value of its parameters is obtained.
An efficient and highly-accurate algorithm (implemented in MATLAB) is presented.
A parallel implementation for large scale problems was developed.