code: def _condition_number(self):
          
          ev = np.linalg.eig(np.matmul(self.xwins.swapaxes(1, 2), self.xwins))[0]
          return np.sqrt(ev.max(axis=1) / ev.min(axis=1))
docs: Condition number of x; ratio of largest to smallest eigenvalue.