Fixed computation for different, but broadcastable, shapes of r,z.
*invariants* function now works with arbitrarily shaped arrays.
Bug fix: Only one unkown dimension possible in reshape.
Symmetrize Voigt stress tensor upon conversion. Invariants now uses stress 3x3 to Voigt by default.§
Cosmetics.