@@ -18,14 +18,7 @@ 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-
22-
23- // #include <xc.h>
24- // added by jghan, 2024-10-10
25- // std::map<xc_func_type, double> scaling_factor_xc = { {XC_GGA_X_ITYH, 1.0} };
26- // std::map<xc_func_type, double> scaling_factor_xc;
27- // scaling_factor_xc[XC_GGA_X_B88] = 1.0;
28- std::map<int , double > scaling_factor_xc = { {XC_GGA_X_ITYH, 1.0 } };
21+ std::map<int , double > scaling_factor_xc = { {XC_GGA_X_ITYH, 1.0 } }; // added by jghan, 2024-10-10
2922
3023void XC_Functional::set_hybrid_alpha (const double alpha_in)
3124{
@@ -233,13 +226,7 @@ void XC_Functional::set_xc_type(const std::string xc_func_in)
233226 func_id.push_back (XC_GGA_X_B88); // complete B88_X, id=106
234227 func_id.push_back (XC_GGA_C_LYP); // complete LYP_C, id=131
235228
236- // according to the order of adding functionals above
237- // scaling_factor_xc.push_back(-1.0);
238- // scaling_factor_xc.push_back(-1.0);
239- // scaling_factor_xc.push_back(1.0);
240- // scaling_factor_xc.push_back(1.0);
241-
242- //
229+ // the scaling factor of CWP22-functionals
243230 scaling_factor_xc[XC_GGA_X_ITYH] = -1.0 ;
244231 scaling_factor_xc[XC_GGA_C_LYPR] = -1.0 ;
245232 scaling_factor_xc[XC_GGA_X_B88] = 1.0 ;
0 commit comments