fix the wrong HF Coulomb at nspin=2

Author: maki49

source/module_cell/klist.h

@@ -157,8 +157,8 @@ class K_Vectors
 
 private:
     int nks;         // number of symmetry-reduced k points in this pool(processor, up+dw)
-    int nkstot;      /// number of symmetry-reduced k points in full k mesh
-    int nkstot_full; /// number of k points before symmetry reduction in full k mesh
+    int nkstot;      /// number of symmetry-reduced k points in full k mesh (up+dw)
+    int nkstot_full; /// number of k points before symmetry reduction in full k mesh (none spin)
 
     int nspin;
     bool kc_done;

source/module_ri/Exx_LRI.hpp

@@ -78,8 +78,7 @@ void Exx_LRI<Tdata>::init(const MPI_Comm &mpi_comm_in,
 			case Conv_Coulomb_Pot_K::Ccp_Type::Hf:
 			{
 				// 4/3 * pi * Rcut^3 = V_{supercell} = V_{unitcell} * Nk
-				const int nspin0 = (PARAM.inp.nspin==2) ? 2 : 1;
-				const double hf_Rcut = std::pow(0.75 * this->p_kv->get_nkstot_full()/nspin0 * omega / (ModuleBase::PI), 1.0/3.0);
+                const double hf_Rcut = std::pow(0.75 * this->p_kv->get_nkstot_full() * omega / (ModuleBase::PI), 1.0 / 3.0);
 				return {{"hf_Rcut", hf_Rcut}};
 			}
 			case Conv_Coulomb_Pot_K::Ccp_Type::Erfc: