fix compile

Author: Qianruipku

source/module_hsolver/hsolver_pw.h

@@ -12,7 +12,7 @@ namespace hsolver
 template <typename T, typename Device = base_device::DEVICE_CPU>
 class HSolverPW
 {
-  private:
+  protected:
     // Note GetTypeReal<T>::type will
     // return T if T is real type(float, double),
     // otherwise return the real type of T(complex<float>, complex<double>)
@@ -84,12 +84,12 @@ class HSolverPW
     const bool need_subspace; // for cg or dav_subspace
     const bool initialed_psi; 
 
-  private:
+  protected:
     Device* ctx = {};
 
     int rank_in_pool = 0;
     int nproc_in_pool = 1;
-
+  private:
     /// @brief calculate the threshold for iterative-diagonalization for each band
     void cal_ethr_band(const double& wk, const double* wg, const double& ethr, std::vector<double>& ethrs);
 

source/module_hsolver/hsolver_pw_sdft.cpp

@@ -44,7 +44,7 @@ void HSolverPW_SDFT<Device>::solve(hamilt::Hamilt<std::complex<double>, Device>*
         {
             this->updatePsiK(pHamilt, psi, ik);
             // template add precondition calculating here
-            update_precondition(precondition, ik, this->wfc_basis->npwk[ik], pes->pot->get_vl_of_0());
+            this->update_precondition(precondition, ik, this->wfc_basis->npwk[ik], pes->pot->get_vl_of_0());
             /// solve eigenvector and eigenvalue for H(k)
             double* p_eigenvalues = &(pes->ekb(ik, 0));
             this->hamiltSolvePsiK(pHamilt, psi, precondition, p_eigenvalues);

source/module_hsolver/test/test_hsolver_sdft.cpp

@@ -138,23 +138,24 @@ class TestHSolverPW_SDFT : public ::testing::Test
     wavefunc wf;
     StoChe<double> stoche;
     hamilt::HamiltSdftPW<std::complex<double>>* p_hamilt_sto = nullptr;
-    hsolver::HSolverPW_SDFT hs_d = hsolver::HSolverPW_SDFT(&kv, 
-                                                           &pwbk, 
-                                                           &wf, 
-                                                           stowf, 
-                                                           stoche, 
-                                                           p_hamilt_sto,
-                                                           "scf",
-                                                           "pw",
-                                                           "cg",
-                                                           false,
-                                                           PARAM.sys.use_uspp,
-                                                           PARAM.input.nspin,
-                     hsolver::DiagoIterAssist<std::complex<double>>::SCF_ITER,
-                     hsolver::DiagoIterAssist<std::complex<double>>::PW_DIAG_NMAX,
-                     hsolver::DiagoIterAssist<std::complex<double>>::PW_DIAG_THR,
-                     hsolver::DiagoIterAssist<std::complex<double>>::need_subspace,
-                     false);
+    hsolver::HSolverPW_SDFT<base_device::DEVICE_CPU> hs_d = hsolver::HSolverPW_SDFT<base_device::DEVICE_CPU>(
+        &kv,
+        &pwbk,
+        &wf,
+        stowf,
+        stoche,
+        p_hamilt_sto,
+        "scf",
+        "pw",
+        "cg",
+        false,
+        PARAM.sys.use_uspp,
+        PARAM.input.nspin,
+        hsolver::DiagoIterAssist<std::complex<double>>::SCF_ITER,
+        hsolver::DiagoIterAssist<std::complex<double>>::PW_DIAG_NMAX,
+        hsolver::DiagoIterAssist<std::complex<double>>::PW_DIAG_THR,
+        hsolver::DiagoIterAssist<std::complex<double>>::need_subspace,
+        false);
 
     hamilt::Hamilt<std::complex<double>> hamilt_test_d;