code: def ordered(self):
          
          active, inactive = self.active_inactive
          order = active + inactive
          return UnitCell(self.matrix[:,order], self.active[order])
docs: An equivalent unit cell with the active cell vectors coming first