docs: def _unvec(vecA, m=None):
      """inverse of _vec() operator"""
      N = vecA.shape[0]
      if m is None:
          m = np.sqrt(vecA.shape[1] + 0.25).astype(np.int64)
      return vecA.reshape((N, m, -1), order='F')