matrix exponential function

concurrent and serial implementation
Does this function exist anywhere? Either way, we don't want it exported.