We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
There was an error while loading. Please reload this page.
1 parent 1c1d006 commit a45138bCopy full SHA for a45138b
source/module_rdmft/rdmft.cpp
@@ -850,6 +850,7 @@ void RDMFT<TK, TR>::cal_Energy(const int cal_type)
850
E_RDMFT[1] = getEnergy(Ehartree_n_k);
851
852
// for Exc
853
+ E_RDMFT[2] = 0.0;
854
#ifdef __EXX
855
if( GlobalC::exx_info.info_global.cal_exx )
856
{
@@ -860,7 +861,6 @@ void RDMFT<TK, TR>::cal_Energy(const int cal_type)
860
861
Parallel_Reduce::reduce_all(E_RDMFT[2]);
862
}
863
#endif
- // the initial value of E_RDMFT[2] is 0.0
864
E_RDMFT[2] += etxc;
865
866
// add up the results obtained by all processors, or we can do reduce_all(wfcHwfc_) before add_wg() used for Etotal to replace it
0 commit comments