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!