anchor.diff:                   ! | det_tmp_bis > = a^{\dagger}_{borb,kspin} a_{aorb,kspin} | det_tmp > 
                    call do_mono_excitation(det_tmp_bis,list_act(aorb),list_act(borb),kspin,i_ok)
                    if(i_ok .ne. 1)cycle
  +                 call get_mono_excitation(det_tmp,det_tmp_bis,exc,phase,N_int)
   !                ! < det_tmp | H | det_tmp_bis >  = F_{aorb,borb}
  -                 call i_H_j(det_tmp_bis,det_tmp,N_int,hab)
  +                 hab = fock_operator_local(aorb,borb,kspin) * phase  
  +!                call i_H_j(det_tmp_bis,det_tmp,N_int,hab)
                    ! < jdet | H | det_tmp_bis > = phase * ( (ir|cv) - (iv|cr) )
   !                call i_H_j(det_tmp_bis,psi_det(1,1,idx(jdet)),N_int,hja)
                    call get_double_excitation(det_tmp_bis,psi_det(1,1,idx(jdet)),exc,phase,N_int)
  
positive.code:  double precision :: delta_e_ja(N_states)
   double precision :: hja
   double precision :: contrib_hij 
   accu_contrib = 0.d0
  !matrix_1h2p = 0.d0
   
  
negative.code: !!!!!!!!!!!!!!!!!!!!!!!!!!! determination of the connections between I and the other J determinants mono excited in the CAS
  !!!!!!!!!!!!!!!!!!!!!!!!!!!! the determinants I and J must be connected by the following operator 
  !!!!!!!!!!!!!!!!!!!!!!!!!!!! <Jdet | a^{\dagger}_b a_{a}  | Idet> 
               integer :: i_hole,i_part
         do jdet = 1, idx(0)
           if(idx(jdet).ne.idet)then
  !         print*, degree(jdet)
  