anchor.diff:       
        print *, k, "res : ", res, "norm : ", sqrt(norm)
        
  -     dIj_unique(:size(X), s) = X(:)
  +     !dIj_unique(:size(X), s) = X(:)
        
        norm = 0.d0
        double precision :: f
  
positive.code:      enddo
       ! rho_mrcc now contains the product of the scaling factors and the
       ! normalization constant
  
     end do
     
  END_PROVIDER
  
  
  BEGIN_PROVIDER [ double precision, dij, (N_det_ref, N_det_non_ref, N_states) ]
  
negative.code:   ok = .true.
  end subroutine
  
  
  
  
  