def METHOD_NAME(self):
    if self.connectivity_indices.shape[1] != self._connectivity_length:
        raise RuntimeError
    for gi in range(self.connectivity_indices.shape[0]):
        ind = self.connectivity_indices[gi, :] - self._index_offset
        coords = self.connectivity_coords[ind, :]
        for i in range(3):
            assert np.unique(coords[:, i]).size == 2
    mylog.debug("Connectivity is consistent.")