code: def axis_aligned_transforms():
      
      I = N.eye(3)
      xy = I[:2]
      xz = N.vstack((I[0],I[2]))
      yz = I[1:]
      return xy,xz,yz
docs: Generates three transformation matrices
      to map three-dimensional data down
      to slices on the xy, xz and yz planes, respectively.
      The typical use case for this is to iterate over the results
      of this method to provide arguments to e.g. the `HyperbolicErrors`
      function.