Follow-up Comment #25, task #7793 (project relax): It remains still to figure out how to speed up the numerical models: The matrix_exponential() and square_matrix_power() functions are the problem. If an equivalent for an array of matrices could be designed, it could be added as a new function to lib.linear_algebra.matrix_exponential and lib.linear_algebra.matrix_power. This is a difficult problem. Regards, Edward _______________________________________________________ Reply to this item at: <http://gna.org/task/?7793> _______________________________________________ Message sent via/by Gna! http://gna.org/