fix UT numerical errors

Author: jieli-matrix

source/module_hsolver/hsolver_pw.h

@@ -36,7 +36,7 @@ class HSolverPW
               const int diag_iter_max_in,
               const double diag_thr_in,
               const bool need_subspace_in,
-              const bool use_k_continuity_in = true)
+              const bool use_k_continuity_in = false)
         : wfc_basis(wfc_basis_in), calculation_type(calculation_type_in), basis_type(basis_type_in), method(method_in),
           use_paw(use_paw_in), use_uspp(use_uspp_in), nspin(nspin_in), scf_iter(scf_iter_in),
           diag_iter_max(diag_iter_max_in), diag_thr(diag_thr_in), need_subspace(need_subspace_in),

source/module_parameter/input_parameter.h

@@ -89,7 +89,7 @@ struct Input_para
     bool diago_smooth_ethr = false; ///< smooth ethr for iter methods
     int pw_diag_ndim = 4;      ///< dimension of workspace for Davidson diagonalization
     int diago_cg_prec = 1;     ///< mohan add 2012-03-31
-    bool use_k_continuity = true;   ///< whether to use k-point continuity for initializing wave functions
+    bool use_k_continuity = false;   ///< whether to use k-point continuity for initializing wave functions
 
     std::string smearing_method = "gauss"; ///< "gauss",
                                            ///< "mp","methfessel-paxton"