@@ -18,12 +18,13 @@ std::vector<int> XC_Functional::func_id(1);
1818int XC_Functional::func_type = 0 ;
1919bool XC_Functional::use_libxc = true ;
2020double XC_Functional::hybrid_alpha = 0.25 ;
21- double static XC_Functional::default_dft_threshold=1.0e-10 ;
21+ double XC_Functional::default_dft_threshold=1.0e-10 ;
22+ double XC_Functional::dens_threshold=1.0e-10 ;
23+ double XC_Functional::zeta_threshold=1.0e-10 ;
24+ double XC_Functional::grho_threshold=1.0e-10 ;
25+ double XC_Functional::tau_threshold=1.0e-10 ;
2226
23- void XC_Functional::set_dft_ingred_thrs (const double d_thr=default_dft_threshold,
24- const double z_thr=default_dft_threshold,
25- const double g_thr=default_dft_threshold,
26- const double t_thr=default_dft_threshold)
27+ void XC_Functional::set_dft_ingred_thrs (double d_thr, double z_thr, double g_thr, double t_thr)
2728{ dens_threshold = d_thr;
2829 zeta_threshold = z_thr;
2930 grho_threshold = g_thr;
@@ -50,11 +51,6 @@ double XC_Functional::get_tau_threshold()
5051 return tau_threshold;
5152}
5253
53- void XC_Functional::set_hybrid_alpha (const double alpha_in)
54- {
55- hybrid_alpha = alpha_in;
56- }
57-
5854double XC_Functional::get_hybrid_alpha ()
5955{
6056 return hybrid_alpha;
@@ -75,7 +71,6 @@ void XC_Functional::set_xc_first_loop(const UnitCell& ucell)
7571the first scf iteration only calculate the functional without exact
7672exchange. but in "nscf" calculation, there is no need of "two-level"
7773method. */
78- XC_Functional::set_dft_ingred_thrs ();
7974 if (ucell.atoms [0 ].ncpp .xc_func == " HF"
8075 || ucell.atoms [0 ].ncpp .xc_func == " PBE0"
8176 || ucell.atoms [0 ].ncpp .xc_func == " HSE" ) {
0 commit comments