docs: def _antisymm_12(C):
      """To get rid of NaNs produced by _scalar2array, antisymmetrize the first
      two indices of operators where C_ijkl = -C_jikl"""
      nans = np.isnan(C)
      C[nans] = -np.einsum('jikl', C)[nans]
      return C