Default precision is "whatever the reult of the math operation is" and hence can grow quite large. Because
PDL only offers doubles I tried using BigFloat to fill a Math::Matrix which was then exponentiated to a high power.
Without a defined precision the program slowed down for
each progressive iteration. Woohoo!