Merge pull request #447 from dyzheng/ABACUS_2.2.0_beta

Fixed bug for stress can't be calculated without force in k line

Author: dyzheng

source/src_lcao/FORCE_k.cpp

@@ -1272,9 +1272,9 @@ void Force_LCAO_k::cal_fvl_dphi_k(
 		// Grid integration here.
 		//--------------------------------
 		// fvl_dphi can not be set to zero here if Vna is used
-		if(isstress&&isforce) 
+		if(isstress||isforce) 
 		{
-			GlobalC::UHM.GK.svl_k_RealSpace(fvl_dphi,svl_dphi,GlobalC::pot.vr_eff1);
+			GlobalC::UHM.GK.svl_k_RealSpace(isforce, isstress, fvl_dphi,svl_dphi,GlobalC::pot.vr_eff1);
 		}
 		else if(isforce) 
 		{