@@ -693,7 +693,7 @@ void RDMFT<TK, TR>::cal_V_XC()
693693 &Vxc_fromRI_c->Hexxs
694694 );
695695 }
696- V_exx_XC->contributeHR ();
696+ // V_exx_XC->contributeHR();
697697 }
698698}
699699
@@ -888,8 +888,8 @@ void RDMFT<TK, TR>::cal_Energy(const int cal_type)
888888
889889 std::cout << std::fixed << std::setprecision (10 ) << " ******\n E(TV + Hartree + XC) by RDMFT: " << E_RDMFT[3 ] << " \n\n ETV_RDMFT: "
890890 << E_RDMFT[0 ] << " \n Ehartree_RDMFT: " << E_RDMFT[1 ] << " \n Exc_RDMFT: " << E_RDMFT[2 ] << " \n E_Ewald: " << E_Ewald
891- << " \n E_entropy(-TS): " << E_entropy << " \n E_descf: " << E_descf << " \n\n Etotal_RDMFT: " << Etotal << " \n Exc_dft: " << E_xc_KS
892- << " \n E_exx_KS: " << E_exx_KS <<" \n ******\n\n " << std::endl;
891+ << " \n E_entropy(-TS): " << E_entropy << " \n E_descf: " << E_descf << " \n\n Etotal_RDMFT: " << Etotal << " \n Exc_ksdft: " << E_xc_KS
892+ << " \n E_exx_ksdft: " << E_exx_KS <<" \n ******\n\n " << std::endl;
893893
894894 std::cout << " \n E_deband_KS: " << E_deband_KS << " \n E_deband_harris_KS: " << E_deband_harris_KS << " \n\n " << std::endl;
895895 }
@@ -899,8 +899,8 @@ void RDMFT<TK, TR>::cal_Energy(const int cal_type)
899899
900900 std::cout << std::fixed << std::setprecision (10 ) << " ******\n E(TV + Hartree + XC) by RDMFT: " << E_RDMFT[3 ] << " \n\n ETV_RDMFT: "
901901 << E_RDMFT[0 ] << " \n E_hartree_RDMFT: " << E_RDMFT[1 ] << " \n Exc_" << XC_func_rdmft << " _RDMFT: " << E_RDMFT[2 ] << " \n E_Ewald: " << E_Ewald
902- << " \n E_entropy(-TS): " << E_entropy << " \n E_descf: " << E_descf << " \n\n Etotal_RDMFT: " << Etotal << " \n Exc_dft: " << E_xc_KS
903- << " \n E_exx_KS: " << E_exx_KS <<" \n ******\n\n " << std::endl;
902+ << " \n E_entropy(-TS): " << E_entropy << " \n E_descf: " << E_descf << " \n\n Etotal_RDMFT: " << Etotal << " \n Exc_ksdft: " << E_xc_KS
903+ << " \n E_exx_ksdft: " << E_exx_KS <<" \n ******\n\n " << std::endl;
904904
905905 std::cout << " \n etxc: " << etxc << " \n vtxc: " << vtxc << " \n " ;
906906
0 commit comments