anchor.diff:        
         do i=1,npre
         do j=1,i
  -        !$OMP ATOMIC
  +        !!$OMP ATOMIC
           delta_cas(ipre(i),ipre(j),i_state) += pre(i) * pre(j) * lambda_mrcc(i_state, k)
         end do
         end do
       end do
  -    !$OMP END PARALLEL DO
  +    !!$OMP END PARALLEL DO
       print *, npres
       npre=0
       do i=1,N_det_ref
  
positive.code: &BEGIN_PROVIDER [ integer(bit_kind), det_ref_active, (N_int,2,N_det_ref) ]
  &BEGIN_PROVIDER [ integer(bit_kind), active_sorb, (N_int,2) ]
  &BEGIN_PROVIDER [ integer(bit_kind), det_cepa0, (N_int,2,N_det_non_ref) ]
    use bitmasks
    implicit none
    
    integer(bit_kind) :: det_noactive(N_int, 2, N_det_non_ref), nonactive_sorb(N_int,2), det(N_int, 2)
    integer i, II, j, k, n, ni, idx(N_det_non_ref), shortcut(0:N_det_non_ref+1)
    logical, external :: detEq
    
    active_sorb(:,:) = 0_8
  
negative.code: !         i_I = J_s
  !         J = I_s
  !       end if
  !       
  !     !$OMP PARALLEL DO default(none) schedule(dynamic,1) shared(delta_ij_old, delta_ii_old)       &
  !     !$OMP private(m,m2,kk, i, k, degree, degree2, l, deg, ni, inac, virt)       &
  !     !$OMP private(ok,p1,p2,h1,h2,s1,s2, blok, wall, I_s, J_s)     &
  !     !$OMP private(phase_iI, phase_Ik, phase_Ji, phase_al, diI, hIi, hJi, delta_JI, dkI, HkI, ci_inv, dia_hla) &
  !     !$OMP private(i_state, contrib, exc_iI, exc_Ik, exc_IJ, det_tmp, det_tmp2)    &
  !     !$OMP shared(N_det_non_ref, N_det_ref, N_int, psi_non_ref, psi_non_ref_coef, psi_ref, psi_ref_coef)   &
  !     !$OMP shared(lambda_mrcc, hf_bitmask, active_sorb,cepa0_shortcut,det_cepa0,N_states) &
  !     !$OMP shared(i_I, J, h_,det_cepa0_idx, linked, blokMwen, nlink, iwall, searchance)
  !       do kk = 1 , nlink(i_I)
  !         k = linked(kk,i_I)
  !         blok = blokMwen(kk,i_I)
  