From 72e1c1b3f7e41bca58215a6194a9705f51efc642 Mon Sep 17 00:00:00 2001 From: maki49 <1579492865@qq.com> Date: Tue, 2 Jul 2024 11:44:46 +0800 Subject: [PATCH] minor fixes --- source/module_esolver/esolver.cpp | 4 ++-- source/module_esolver/esolver_ks_lip.cpp | 6 +++--- source/module_esolver/esolver_ks_lip.h | 6 +++--- source/module_esolver/esolver_ks_pw.cpp | 16 ++++++---------- source/module_esolver/esolver_ks_pw.h | 4 ++-- 5 files changed, 16 insertions(+), 20 deletions(-) diff --git a/source/module_esolver/esolver.cpp b/source/module_esolver/esolver.cpp index 6f74986e1c..4774f38baf 100644 --- a/source/module_esolver/esolver.cpp +++ b/source/module_esolver/esolver.cpp @@ -50,7 +50,7 @@ std::string determine_type(void) esolver_type = "ksdft_lip"; } #else - ModuleBase::WARNING_QUIT("ESolver", "LCAO basis type must be compiled with __LCAO"); + ModuleBase::WARNING_QUIT("ESolver", "Calculation involving numerical orbitals must be compiled with __LCAO"); #endif } else if (GlobalV::BASIS_TYPE == "lcao") @@ -65,7 +65,7 @@ std::string determine_type(void) esolver_type = "ksdft_lcao"; } #else - ModuleBase::WARNING_QUIT("ESolver", "LCAO basis type must be compiled with __LCAO"); + ModuleBase::WARNING_QUIT("ESolver", "Calculation involving numerical orbitals must be compiled with __LCAO"); #endif } diff --git a/source/module_esolver/esolver_ks_lip.cpp b/source/module_esolver/esolver_ks_lip.cpp index b47a690a7b..76a6164e00 100644 --- a/source/module_esolver/esolver_ks_lip.cpp +++ b/source/module_esolver/esolver_ks_lip.cpp @@ -51,12 +51,12 @@ namespace ModuleESolver } template - void ESolver_KS_LIP::init_hsolver() + void ESolver_KS_LIP::allocate_hsolver() { this->phsol = new hsolver::HSolverLIP(this->pw_wfc); } template - void ESolver_KS_LIP::delete_hsolver() + void ESolver_KS_LIP::deallocate_hsolver() { if (this->phsol != nullptr) { @@ -113,7 +113,7 @@ namespace ModuleESolver } else { - ModuleBase::WARNING_QUIT("ESolver_KS_PW", "HSolver has not been initialed!"); + ModuleBase::WARNING_QUIT("ESolver_KS_LIP", "HSolver has not been allocated."); } // add exx #ifdef __LCAO diff --git a/source/module_esolver/esolver_ks_lip.h b/source/module_esolver/esolver_ks_lip.h index c6389099f2..793a227ec5 100644 --- a/source/module_esolver/esolver_ks_lip.h +++ b/source/module_esolver/esolver_ks_lip.h @@ -1,6 +1,6 @@ #ifndef ESOLVER_KS_LIP_H #define ESOLVER_KS_LIP_H -#include "./esolver_ks_pw.h" +#include "module_esolver/esolver_ks_pw.h" #include "module_hsolver/hsolver_lip.h" namespace ModuleESolver { @@ -21,8 +21,8 @@ namespace ModuleESolver protected: - virtual void init_hsolver() override; - virtual void delete_hsolver() override; + virtual void allocate_hsolver() override; + virtual void deallocate_hsolver() override; }; } // namespace ModuleESolver diff --git a/source/module_esolver/esolver_ks_pw.cpp b/source/module_esolver/esolver_ks_pw.cpp index f17eea44d2..bb2e36ff15 100644 --- a/source/module_esolver/esolver_ks_pw.cpp +++ b/source/module_esolver/esolver_ks_pw.cpp @@ -68,7 +68,7 @@ template ESolver_KS_PW::~ESolver_KS_PW() { // delete HSolver and ElecState - this->delete_hsolver(); + this->deallocate_hsolver(); if (this->pelec != nullptr) { delete reinterpret_cast*>(this->pelec); @@ -161,7 +161,7 @@ void ESolver_KS_PW::before_all_runners(Input& inp, UnitCell& ucell) // 2) initialize HSolver if (this->phsol == nullptr) { - this->init_hsolver(); + this->allocate_hsolver(); } // 3) initialize ElecState, @@ -221,18 +221,14 @@ void ESolver_KS_PW::before_all_runners(Input& inp, UnitCell& ucell) } } template -void ESolver_KS_PW::init_hsolver() +void ESolver_KS_PW::allocate_hsolver() { this->phsol = new hsolver::HSolverPW(this->pw_wfc, &this->wf); } template -void ESolver_KS_PW::delete_hsolver() +void ESolver_KS_PW::deallocate_hsolver() { - if (this->phsol != nullptr) - { - delete reinterpret_cast*>(this->phsol); - this->phsol = nullptr; - } + delete reinterpret_cast*>(this->phsol); } template void ESolver_KS_PW::init_after_vc(Input& inp, UnitCell& ucell) @@ -266,7 +262,7 @@ void ESolver_KS_PW::init_after_vc(Input& inp, UnitCell& ucell) this->pw_wfc->collect_local_pw(inp.erf_ecut, inp.erf_height, inp.erf_sigma); delete this->phsol; - this->init_hsolver(); + this->allocate_hsolver(); delete this->pelec; this->pelec = new elecstate::ElecStatePW(this->pw_wfc, diff --git a/source/module_esolver/esolver_ks_pw.h b/source/module_esolver/esolver_ks_pw.h index b10b011116..7af684815f 100644 --- a/source/module_esolver/esolver_ks_pw.h +++ b/source/module_esolver/esolver_ks_pw.h @@ -56,8 +56,8 @@ class ESolver_KS_PW : public ESolver_KS // Init Global class void Init_GlobalC(Input& inp, UnitCell& ucell, pseudopot_cell_vnl& ppcell); - virtual void init_hsolver(); - virtual void delete_hsolver(); + virtual void allocate_hsolver(); + virtual void deallocate_hsolver(); //! hide the psi in ESolver_KS for tmp use psi::Psi, base_device::DEVICE_CPU>* psi = nullptr;