anchor.diff:  
   END_PROVIDER
   
  -! BEGIN_PROVIDER [ double precision, delta_ij, (N_states,N_det_non_ref,N_det_ref) ]
  -!&BEGIN_PROVIDER [ double precision, delta_ii, (N_states,N_det_ref) ]
  -! implicit none
  -! BEGIN_DOC
  -! ! Dressing matrix in N_det basis
  -! END_DOC
  -! integer :: i,j,m
  -! delta_ij = 0.d0
  -! delta_ii = 0.d0
  -! call H_apply_mrcc(delta_ij,delta_ii,N_states,N_det_non_ref,N_det_ref)
  -!
  -!END_PROVIDER
  -             
   
   BEGIN_PROVIDER [ double precision, h_matrix_dressed, (N_det,N_det,N_states) ]
    implicit none
  
positive.code:         i = active_excitation_to_determinants_idx(j,a_coll)
          if (i==0) exit
          rho_mrcc(i,s) = rho_mrcc(i,s) + active_excitation_to_determinants_val(s,j,a_coll) * dIj_unique(a_col,s)
  !        rho_mrcc(i,s) = rho_mrcc(i,s) + active_excitation_to_determinants_val(s,j,a_coll) * X(a_col)
        enddo
      end do
  
  
negative.code:     end do
      dIj_unique(1:size(X), s) = X(1:size(X))
  
    enddo
  
    do s=1,N_states
  
  