anchor.diff:    do i=1,dim_1e
         E1e = E1e + 1.0d0*NOEps(i,1)*Hcore(i,i)
       do j=1,dim_1e
  -      E2ec = E2ec + 1.0d0*NOEps(i,1)*NOEps(j,1)*Jmat(i,j,1)
  -      E2ex = E2ex - 0.5d0*NOEps(i,1)*NOEps(j,1)*Kmat(i,j,1)
  +!      E2ec = E2ec + 1.0d0*NOEps(i,1)*NOEps(j,1)*Jmat(i,j,1)
  +!      E2ex = E2ex - 0.5d0*NOEps(i,1)*NOEps(j,1)*Kmat(i,j,1)
  +      E2e = NOEps(i,1)*NOEps(j,1)*(Fock(i,j,1)-Hcore(i,j))      
       enddo
     enddo
   
  -  E2e = E2ex + E2ec
  -  Erdmft = E1e + E2ec + E2ex
  +!  E2e = E2ex + E2ec
  +  Erdmft = E1e + E2e!+ E2ec + E2ex
   
     write(*,*) "         E1e = ", E1e   
     write(*,*) "         E2ex = ", E2ex   
  
positive.code:   Kmat = 0.0d0
  
    call AOMat_to_MO()
  
    Erdmft = 0.0d0
    E1e  = 0.0d0
  
negative.code: !#                 Run 1-RDMFT
  !#############################################
  subroutine run_rdmft()
    use module_data, only      : dim_1e,Hcore,NOEps,NOCoef,Coef,Fock
    use module_wavefun, only   : Fock_to_MO
    implicit none
    integer                   :: i,j
    double precision          :: E1e,E2ec,E2ex,E2e
  