fix : init matrix in write_electric_pot function

Author: YuLiu98

source/src_io/write_pot.cpp

@@ -268,7 +268,12 @@ void Potential::write_elecstat_pot(const std::string &fn, const std::string &fn_
     //==========================================
     for (int ir = 0;ir < GlobalC::pw.nrxx;ir++)
     {
-        v_elecstat[ir] = Porter[ir].real() + this->vltot[ir] + v_efield(0, ir);
+        v_elecstat[ir] = Porter[ir].real() + this->vltot[ir];
+
+        if (GlobalV::EFIELD && GlobalV::DIPOLE)
+        {
+            v_elecstat[ir] += v_efield(0, ir);
+        }
     }
 
     //-------------------------------------------

source/src_pw/forces.cpp

@@ -76,6 +76,10 @@ void Forces::init(ModuleBase::matrix& force)
     {
         force_e.create(GlobalC::ucell.nat, 3);
         Efield::compute_force(GlobalC::ucell, force_e);
+        if(GlobalV::TEST_FORCE)
+        {
+            Forces::print("EFIELD      FORCE (Ry/Bohr)", force_e);
+        }
     }
 
     //impose total force = 0