add const for variable functions

Author: PeizeLin

source/source_lcao/module_deltaspin/spin_constrain.cpp

@@ -34,7 +34,7 @@ double SpinConstrain<TK>::cal_escon()
 }
 
 template <typename TK>
-double SpinConstrain<TK>::get_escon()
+double SpinConstrain<TK>::get_escon() const
 {
     return this->escon_;
 }
@@ -67,7 +67,7 @@ void SpinConstrain<TK>::set_nspin(int nspin_in)
 
 /// get nspin
 template <typename TK>
-int SpinConstrain<TK>::get_nspin()
+int SpinConstrain<TK>::get_nspin() const
 {
     return this->nspin_;
 }
@@ -439,35 +439,35 @@ void SpinConstrain<TK>::set_input_parameters(double sc_thr_in,
 
 /// get sc_thr
 template <typename TK>
-double SpinConstrain<TK>::get_sc_thr()
+double SpinConstrain<TK>::get_sc_thr() const
 {
     return this->sc_thr_;
 }
 
 /// get nsc
 template <typename TK>
-int SpinConstrain<TK>::get_nsc()
+int SpinConstrain<TK>::get_nsc() const
 {
     return this->nsc_;
 }
 
 /// get nsc_min
 template <typename TK>
-int SpinConstrain<TK>::get_nsc_min()
+int SpinConstrain<TK>::get_nsc_min() const
 {
     return this->nsc_min_;
 }
 
 /// get alpha_trial
 template <typename TK>
-double SpinConstrain<TK>::get_alpha_trial()
+double SpinConstrain<TK>::get_alpha_trial() const
 {
     return this->alpha_trial_;
 }
 
 /// get sccut
 template <typename TK>
-double SpinConstrain<TK>::get_sccut()
+double SpinConstrain<TK>::get_sccut() const
 {
     return this->restrict_current_;
 }
@@ -481,7 +481,7 @@ void SpinConstrain<TK>::set_sc_drop_thr(double sc_drop_thr_in)
 
 /// get sc_drop_thr
 template <typename TK>
-double SpinConstrain<TK>::get_sc_drop_thr()
+double SpinConstrain<TK>::get_sc_drop_thr() const
 {
     return this->sc_drop_thr_;
 }

source/source_lcao/module_deltaspin/spin_constrain.h

@@ -66,7 +66,7 @@ class SpinConstrain
    */
   double cal_escon();
 
-  double get_escon();
+  double get_escon() const;
 
   void run_lambda_loop(int outer_step, 
 		  bool rerun = true);
@@ -181,7 +181,7 @@ class SpinConstrain
     /// set nspin
     void set_nspin(int nspin);
     /// get nspin
-    int get_nspin();
+    int get_nspin() const;
     /// zero atomic magnetic moment
     void zero_Mi();
     /// get decay_grad
@@ -202,17 +202,17 @@ class SpinConstrain
                               double sccut_in,
                               double sc_drop_thr_in);
     /// get sc_thr
-    double get_sc_thr();
+    double get_sc_thr() const;
     /// get nsc
-    int get_nsc();
+    int get_nsc() const;
     /// get nsc_min
-    int get_nsc_min();
+    int get_nsc_min() const;
     /// get alpha_trial
-    double get_alpha_trial();
+    double get_alpha_trial() const;
     /// get sccut
-    double get_sccut();
+    double get_sccut() const;
     /// get sc_drop_thr
-    double get_sc_drop_thr();
+    double get_sc_drop_thr() const;
     /// @brief set orbital parallel info
     void set_ParaV(Parallel_Orbitals* ParaV_in);
     /// @brief set parameters for solver

source/source_pw/module_ofdft/kedf_manager.cpp

@@ -166,7 +166,7 @@ void KEDF_Manager::get_potential(
  *
  * @return kinetic energy
  */
-double KEDF_Manager::get_energy()
+double KEDF_Manager::get_energy() const
 {
     double kinetic_energy = 0.0;
 

source/source_pw/module_ofdft/kedf_manager.h

@@ -41,7 +41,7 @@ class KEDF_Manager
         ModuleBase::matrix& rpot
     );
 
-    double get_energy();
+    double get_energy() const;
 
     void get_energy_density(
         const double* const* prho,