code: def fast_matrix_copy(f, t=None, n_threads=1):
      
      if not f.flags['F_CONTIGUOUS']:
          raise RuntimeError(
              'This will not be fast unless input array f is Fortran-contiguous.')
  
      if t is None:
          t = asmatrix(empty(f.shape, order='F'))
      elif not t.flags['F_CONTIGUOUS']:
          raise RuntimeError(
              'This will not be fast unless input array t is Fortran-contiguous.')
  
      # Figure out how to divide job up between threads.
      dcopy_wrap(ravel(asarray(f.T)), ravel(asarray(t.T)))
      return t
docs: Not any faster than a serial copy so far.