Modify mp

Author: Flying-dragon-boxing

source/module_cell/klist.h

@@ -140,6 +140,12 @@ class K_Vectors
     {
         this->nkstot_full = value;
     }
+
+    bool get_is_mp() const
+    {
+        return is_mp;
+    }
+
     std::vector<int> ik2iktot; ///<[nks] map ik to the global index of k points
 
   private:

source/module_hamilt_pw/hamilt_pwdft/operator_pw/op_exx_pw.cpp

@@ -54,7 +54,7 @@ OperatorEXXPW<T, Device>::OperatorEXXPW(const int* isk_in,
     : isk(isk_in), wfcpw(wfcpw_in), rhopw(rhopw_in), kv(kv_in), ucell(ucell)
 {
     gamma_extrapolation = PARAM.inp.exx_gamma_extrapolation;
-    if (kv_in->nmp[0] == 0 || kv_in->nmp[1] == 0 || kv_in->nmp[2] == 0)
+    if (!kv_in->get_is_mp())
     {
         gamma_extrapolation = false;
     }

source/module_hamilt_pw/hamilt_pwdft/stress_func_exx.cpp

@@ -11,7 +11,7 @@ void Stress_PW<FPTYPE, Device>::stress_exx(ModuleBase::matrix& sigma,
 {
     double nqs_half1 = 0.5 * p_kv->nmp[0], nqs_half2 = 0.5 * p_kv->nmp[1], nqs_half3 = 0.5 * p_kv->nmp[2];
     bool gamma_extrapolation = PARAM.inp.exx_gamma_extrapolation;
-    if (p_kv->nmp[0] == 0 || p_kv->nmp[1] == 0 || p_kv->nmp[2] == 0)
+    if (!p_kv->get_is_mp())
     {
         gamma_extrapolation = false;
     }