anchor.diff:    logical :: ok
     logical, external :: detEq
     integer, external :: omp_get_thread_num
  -
  +  double precision :: contrib(N_states)
  +  
  +  contrib = 0d0
     delta_ij_mrcc = 0d0
     delta_ii_mrcc = 0d0
     delta_ij_s2_mrcc = 0d0
  
positive.code:   !$OMP PARALLEL DO default(none)  schedule(dynamic) &
    !$OMP shared(psi_ref, psi_non_ref, hh_exists, pp_exists, N_int, hh_shortcut) &
    !$OMP shared(N_det_generators, coefs,N_det_non_ref, N_det_ref, delta_ii_mrcc_sto, delta_ij_mrcc_sto) &
    !$OMP shared(psi_det_generators, delta_ii_s2_mrcc_sto, delta_ij_s2_mrcc_sto) &
    !$OMP private(i,j,curnorm,myCoef, h, n, mask, omask, buf, ok, iproc) 
    do gen= 1,N_det_generators
      if(coefs(gen) == 0d0) cycle
  
negative.code:   !double precision, external :: get_dij, get_dij_index
    double precision :: Delta_E_inv(N_states)
    double precision, intent(in) :: coef
    double precision, intent(inout) :: contrib(N_states)
    double precision :: sdress, hdress
  
    if (perturbative_triples) then
      PROVIDE one_anhil fock_virt_total fock_core_inactive_total one_creat
  