Fix unittests of Charge<double>

Author: PeizeLin

source/source_esolver/esolver_of.cpp

@@ -250,7 +250,7 @@ void ESolver_OF::before_opt(const int istep, UnitCell& ucell)
             for (int ibs = 0; ibs < this->pw_rho->nrxx; ++ibs)
             {
                 // Here we initialize rho to be uniform,
-                // because the rho got by pot.init_pot -> Charge<double>::atomic_rho may contain minus elements.
+                // because the rho got by pot.init_pot -> Charge::atomic_rho may contain minus elements.
                 this->chr.rho[is][ibs] = this->nelec_[is] / ucell.omega;
                 this->pphi_[is][ibs] = sqrt(this->chr.rho[is][ibs]);
             }

source/source_estate/test/elecstate_base_test.cpp

@@ -27,10 +27,10 @@ Potential::~Potential()
 {
 }
 } // namespace elecstate
-Charge<double>::Charge<double>()
+template<> Charge<double>::Charge()
 {
 }
-Charge<double>::~Charge<double>()
+template<> Charge<double>::~Charge()
 {
 }
 UnitCell::UnitCell()
@@ -77,24 +77,24 @@ void ModulePW::PW_Basis::initgrids(double, ModuleBase::Matrix3, int, int, int)
 void ModulePW::PW_Basis::distribute_r()
 {
 }
-void Charge<double>::set_rho_core(const UnitCell& ucell, ModuleBase::ComplexMatrix const&, const bool*)
+template<> void Charge<double>::set_rho_core(const UnitCell& ucell, ModuleBase::ComplexMatrix const&, const bool*)
 {
 }
-void Charge<double>::init_rho(const UnitCell&,
+template<> void Charge<double>::init_rho(const UnitCell&,
                       const Parallel_Grid&,
                       ModuleBase::ComplexMatrix const&,
                       ModuleSymmetry::Symmetry& symm,
                       const void*,
                       const void*)
 {
 }
-void Charge<double>::set_rhopw(ModulePW::PW_Basis*)
+template<> void Charge<double>::set_rhopw(ModulePW::PW_Basis*)
 {
 }
-void Charge<double>::renormalize_rho()
+template<> void Charge<double>::renormalize_rho()
 {
 }
-void Charge<double>::check_rho()
+template<> void Charge<double>::check_rho()
 {
 }
 

source/source_estate/test/elecstate_magnetism_test.cpp

@@ -27,10 +27,10 @@
 #define private public
 #include "source_estate/magnetism.h"
 #undef private
-Charge<double>::Charge<double>()
+template<> Charge<double>::Charge()
 {
 }
-Charge<double>::~Charge<double>()
+template<> Charge<double>::~Charge()
 {
 }
 

source/source_estate/test/elecstate_print_test.cpp

@@ -31,10 +31,10 @@ Magnetism::Magnetism(){}
 Magnetism::~Magnetism(){}
 InfoNonlocal::InfoNonlocal(){}
 InfoNonlocal::~InfoNonlocal(){}
-Charge<double>::Charge()
+template<> Charge<double>::Charge()
 {
 }
-Charge<double>::~Charge()
+template<> Charge<double>::~Charge()
 {
 }
 SepPot::SepPot(){}

source/source_estate/test/elecstate_pw_test.cpp

@@ -25,10 +25,10 @@ Potential::~Potential()
 {
 }
 } // namespace elecstate
-Charge<double>::Charge()
+template<> Charge<double>::Charge()
 {
 }
-Charge<double>::~Charge()
+template<> Charge<double>::~Charge()
 {
 }
 UnitCell::UnitCell()
@@ -126,25 +126,24 @@ Fcoef::~Fcoef()
 {
 }
 #include "source_cell/klist.h"
-
-void Charge<double>::set_rho_core(const UnitCell& ucell, ModuleBase::ComplexMatrix const&, const bool*)
+template<> void Charge<double>::set_rho_core(const UnitCell& ucell, ModuleBase::ComplexMatrix const&, const bool*)
 {
 }
-void Charge<double>::init_rho(const UnitCell&,
+template<> void Charge<double>::init_rho(const UnitCell&,
                       const Parallel_Grid&,
                       ModuleBase::ComplexMatrix const&,
                       ModuleSymmetry::Symmetry& symm,
                       const void*,
                       const void*)
 {
 }
-void Charge<double>::set_rhopw(ModulePW::PW_Basis*)
+template<> void Charge<double>::set_rhopw(ModulePW::PW_Basis*)
 {
 }
-void Charge<double>::renormalize_rho()
+template<> void Charge<double>::renormalize_rho()
 {
 }
-void Charge<double>::check_rho()
+template<> void Charge<double>::check_rho()
 {
 }
 

source/source_estate/test/potential_new_test.cpp

@@ -36,10 +36,10 @@ InfoNonlocal::~InfoNonlocal()
 {
 }
 #endif
-Charge<double>::Charge()
+template<> Charge<double>::Charge()
 {
 }
-Charge<double>::~Charge()
+template<> Charge<double>::~Charge()
 {
 }
 surchem::surchem()

source/source_estate/test_mpi/charge_mpi_test.cpp

@@ -8,10 +8,10 @@
 #include "source_io/module_parameter/parameter.h"
 
 bool XC_Functional::ked_flag = false;
-Charge<double>::Charge()
+template<> Charge<double>::Charge()
 {
 }
-Charge<double>::~Charge()
+template<> Charge<double>::~Charge()
 {
     delete[] rec;
     delete[] dis;

source/source_hamilt/module_xc/test/xc3_mock.h

@@ -181,8 +181,8 @@ namespace GlobalC
 UnitCell::UnitCell(){};
 UnitCell::~UnitCell(){};
 
-Charge<double>::Charge(){};
-Charge<double>::~Charge(){};
+template<> Charge<double>::Charge(){};
+template<> Charge<double>::~Charge(){};
 
 Magnetism::Magnetism(){};
 Magnetism::~Magnetism(){};

source/source_hsolver/test/test_hsolver_sdft.cpp

@@ -183,8 +183,8 @@ void Stochastic_Iter<T, Device>::cal_storho(const UnitCell& ucell,
 }
 template class Stochastic_Iter<std::complex<double>, base_device::DEVICE_CPU>;
 
-Charge<double>::Charge(){};
-Charge<double>::~Charge(){};
+template<> Charge<double>::Charge(){};
+template<> Charge<double>::~Charge(){};
 
 // Mock implementations for the template functions causing linking errors
 namespace ModulePW {

source/source_io/test/read_wf2rho_pw_test.cpp

@@ -28,10 +28,10 @@ Parallel_Grid::Parallel_Grid()
 Parallel_Grid::~Parallel_Grid()
 {
 }
-Charge<double>::Charge()
+template<> Charge<double>::Charge()
 {
 }
-Charge<double>::~Charge()
+template<> Charge<double>::~Charge()
 {
 }
 UnitCell::UnitCell()