mpi_bcast energies

Author: Qianruipku

source/module_esolver/esolver_ks_pw.cpp

@@ -441,6 +441,7 @@ void ESolver_KS_PW<T, Device>::update_pot(const int istep, const int iter)
         }
         this->pelec->pot->update_from_charge(this->pelec->charge, &GlobalC::ucell);
         this->pelec->f_en.descf = this->pelec->cal_delta_escf();
+        MPI_Bcast(&(this->pelec->f_en.descf), 1, MPI_DOUBLE, 0, PARAPW_WORLD);
     }
     else
     {

source/module_esolver/esolver_sdft_pw.cpp

@@ -231,6 +231,7 @@ void ESolver_SDFT_PW::hamilt2density(int istep, int iter, double ethr)
         }
 #endif
     }
+    MPI_Bcast(&(this->pelec->f_en.deband), 1, MPI_DOUBLE, 0, PARAPW_WORLD);
 }
 
 double ESolver_SDFT_PW::cal_energy()