Added function to compute the matrix exponential for higher dimensional data of shape [NE][NS][NM][NO][ND][7][7].
This is done by using numpy.einsum, to make the dot product of the last two axis.
Task #7807 (https://gna.org/task/index.php?7807): Speed-up of dispersion models for Clustered analysis.
Authored by: tlinnet 2014-06-19
Parent: [r24155]
Child: [r24157]
