Fix: update template class Charge initialization

Author: PeizeLin

source/source_estate/module_charge/charge.cpp

@@ -764,6 +764,5 @@ void Charge<Tr>::init_final_scf()
     return;
 }
 
-
 template class Charge<double>;
-template class Charge<std::complex<double>>;
+template class Charge<std::complex<double>>;

source/source_estate/module_charge/charge_init.cpp

@@ -456,7 +456,5 @@ void Charge<Tr>::non_linear_core_correction
     return;
 }
 
-
-
 template class Charge<double>;
-template class Charge<std::complex<double>>;
+template class Charge<std::complex<double>>;

source/source_estate/module_charge/charge_mpi.cpp

@@ -150,7 +150,7 @@ void Charge<Tr>::rho_mpi()
     return;
 }
 
-
 template class Charge<double>;
 template class Charge<std::complex<double>>;
+
 #endif

source/source_estate/test/charge_extra_test.cpp

@@ -29,8 +29,7 @@ Parallel_Grid::~Parallel_Grid(){};
 
 
 // mock functions for Charge<double>
-template<>
-Charge<double>::Charge()
+template<> Charge<double>::Charge()
 {
     rhopw = new ModulePW::PW_Basis;
     rhopw->nrxx = 8;
@@ -45,21 +44,20 @@ Charge<double>::Charge()
         rho[0][i] = i + 1;
     }
 }
-template<>
-Charge<double>::~Charge()
+template<> Charge<double>::~Charge()
 {
     delete[] rho[0];
     delete[] rho;
     delete rhopw;
 }
-template<>
-void Charge<double>::atomic_rho(const int spin_number_need,
+template<> void Charge<double>::atomic_rho(const int spin_number_need,
                         const double& omega,
                         double** rho_in,
                         const ModuleBase::ComplexMatrix& strucFac,
                         const UnitCell& ucell) const
 {
 }
+template class Charge<double>;
 
 // mock functions for PW_Basis
 namespace ModulePW

source/source_estate/test/charge_mixing_test.cpp

@@ -15,26 +15,23 @@ int XC_Functional::func_type = 1;
 bool XC_Functional::ked_flag = false;
 
 // mock function
-Magnetism::~Magnetism()
+Magnetism::~Magnetism() 
 {
 }
 Magnetism::Magnetism()
 {
 }
-template<>
-Charge<double>::~Charge()
+template<> Charge<double>::~Charge()
 {
 }
-template<>
-Charge<double>::Charge()
+template<> Charge<double>::Charge()
 {
 }
-
-template<>
-void Charge<double>::set_rhopw(ModulePW::PW_Basis* rhopw_in)
+template<> void Charge<double>::set_rhopw(ModulePW::PW_Basis* rhopw_in)
 {
     this->rhopw = rhopw_in;
 }
+template class Charge<double>;
 #ifdef __LCAO
 InfoNonlocal::InfoNonlocal()
 {

source/source_estate/test/elecstate_base_test.cpp

@@ -97,6 +97,7 @@ template<typename Tr> void Charge<Tr>::renormalize_rho()
 template<typename Tr> void Charge<Tr>::check_rho()
 {
 }
+template class Charge<double>;
 
 /************************************************
  *  unit test of elecstate.cpp

source/source_estate/test/elecstate_magnetism_test.cpp

@@ -33,6 +33,7 @@ template<> Charge<double>::Charge()
 template<> Charge<double>::~Charge()
 {
 }
+template class Charge<double>;
 
 
 class MagnetismTest : public ::testing::Test

source/source_estate/test/elecstate_print_test.cpp

@@ -31,12 +31,9 @@ Magnetism::Magnetism(){}
 Magnetism::~Magnetism(){}
 InfoNonlocal::InfoNonlocal(){}
 InfoNonlocal::~InfoNonlocal(){}
-template<> Charge<double>::Charge()
-{
-}
-template<> Charge<double>::~Charge()
-{
-}
+template<> Charge<double>::Charge(){}
+template<> Charge<double>::~Charge(){}
+template class Charge<double>;
 SepPot::SepPot(){}
 SepPot::~SepPot(){}
 Sep_Cell::Sep_Cell() noexcept {}

source/source_estate/test/elecstate_pw_test.cpp

@@ -146,6 +146,7 @@ template<typename Tr> void Charge<Tr>::renormalize_rho()
 template<typename Tr> void Charge<Tr>::check_rho()
 {
 }
+template class Charge<double>;
 
 void Set_GlobalV_Default()
 {

source/source_estate/test/potential_new_test.cpp

@@ -42,6 +42,7 @@ template<> Charge<double>::Charge()
 template<> Charge<double>::~Charge()
 {
 }
+template class Charge<double>;
 surchem::surchem()
 {
 }