fix bug (confuse evolve_psi with evolve_phi)

Author: lyb9812

source/source_esolver/esolver_of_tddft.cpp

@@ -22,12 +22,12 @@ namespace ModuleESolver
 ESolver_OF_TDDFT::ESolver_OF_TDDFT()
 {
     this->classname = "ESolver_OF_TDDFT";
-    this->evolve_psi=new EVOLVE_PSI();
+    this->evolve_phi=new EVOLVE_PHI();
 }
 
 ESolver_OF_TDDFT::~ESolver_OF_TDDFT()
 {
-    delete this->evolve_psi;
+    delete this->evolve_phi;
     for (int is = 0; is < PARAM.inp.nspin; ++is)
     {
         delete[] this->pphi_td[is];
@@ -102,7 +102,7 @@ void ESolver_OF_TDDFT::runner(UnitCell& ucell, const int istep)
     }
     else
     {
-        this->evolve_psi->propagate_psi(this->pelec, this->chr, ucell, this->pphi_td, this->pw_rho);
+        this->evolve_phi->propagate_psi(this->pelec, this->chr, ucell, this->pphi_td, this->pw_rho);
         for (int is = 0; is < PARAM.inp.nspin; ++is)
         {
             for (int ir = 0; ir < pw_rho->nrxx; ++ir)

source/source_esolver/esolver_of_tddft.h

@@ -2,7 +2,7 @@
 #define ESOLVER_OF_TDDFT_H
 
 #include "esolver_of.h"
-#include "source_pw/module_ofdft/evolve_psi.h"
+#include "source_pw/module_ofdft/evolve_phi.h"
 
 namespace ModuleESolver
 {
@@ -16,7 +16,7 @@ class ESolver_OF_TDDFT : public ESolver_OF
 
   protected:
     std::complex<double>** pphi_td = nullptr;                     // pphi[i] = ppsi.get_pointer(i), which will be freed in ~Psi().
-    EVOLVE_PSI* evolve_psi=nullptr;
+    EVOLVE_PHI* evolve_phi=nullptr;
 };
 } // namespace ModuleESolver