def _add_labels(self):
        coords = self['pore.coords']
        self['pore.front'] = self.tomask(coords[:, 0] <= self._Lc)
        self['pore.left'] = self.tomask(coords[:, 1] <= self._Lc)
        self['pore.bottom'] = self.tomask(coords[:, 2] <= self._Lc)
        self['pore.back'] = self.tomask(coords[:, 0] >= (self._Lc*(self._Nx-1)))
        self['pore.right'] = self.tomask((self._Lc*(self._Ny-1)) >= coords[:, 1])
        self['pore.top'] = self.tomask(coords[:, 2] >= (self._Lc*(self._Nz-1)))
        for item in ['top', 'bottom', 'left', 'right', 'front', 'back']:
            ps = self.pores(item)
            ts = self.find_neighbor_throats(ps)
            ps = self.find_connected_pores(ts)
            ps0 = self['pore.'+item][ps[:, 0]]
            ps1 = self['pore.'+item][ps[:, 1]]
            ts = ts[ps1*ps0]
            self['throat.'+item] = self.tomask(throats=ts)
dataset/ETHPy150Open PMEAL/OpenPNM/OpenPNM/Network/__TestNet__.py/TestNet._add_labels/swapped operands of `>=`