code: def get_perturbed_indices(self, tol=1e-8):
          
          indices = list(zip(*np.where(abs(self - np.eye(3)) > tol)))
          return indices
docs: Gets indices of perturbed elements of the deformation gradient,
          i. e. those that differ from the identity