anchor.diff:  !               HkI = h_(J,l)
   !               dkI(i_state) = HkI * lambda_mrcc(i_state, l)
   !               contrib = dkI(i_state) * delta_JI
  +!               !$OMP ATOMIC
   !               delta_ij_old(J,det_cepa0_idx(k),i_state) += contrib
   !               if(dabs(psi_ref_coef(i_I,i_state)).ge.5.d-5) then
  +!                 !$OMP ATOMIC
   !                 delta_ii_old(J,i_state) -= contrib * ci_inv(i_state) * psi_non_ref_coef(det_cepa0_idx(k),i_state)
   !               end if
   !             enddo !i_state
  
positive.code: !         i_I = J_s
  !         J = I_s
  !       end if
  ! 
  !       do kk = 1 , nlink(i_I)
  !         k = linked(kk,i_I)
  !         blok = blokMwen(kk,i_I)
  
negative.code: !     enddo !I
  !       
  ! END_PROVIDER
  ! ! 
  ! 
  
  