Fix: nspin2 error

Author: dyzheng

source/module_elecstate/module_dm/density_matrix.cpp

@@ -614,10 +614,6 @@ void DensityMatrix<std::complex<double>, double>::cal_DMR_full(hamilt::HContaine
         // get global indexes of whole matrix for each atom in this process
         int row_ap = this->_paraV->atom_begin_row[iat1];
         int col_ap = this->_paraV->atom_begin_col[iat2];
-        if (row_ap == -1 || col_ap == -1)
-        {
-            throw std::string("Atom-pair not belong this process");
-        }
         for (int ir = 0; ir < tmp_ap.get_R_size(); ++ir)
         {
             const ModuleBase::Vector3<int> r_index = tmp_ap.get_R_index(ir);

source/module_hamilt_lcao/hamilt_lcaodft/FORCE_STRESS.cpp

@@ -185,7 +185,7 @@ void Force_Stress_LCAO<T>::getForceStress(const bool isforce,
         const auto* dm_p = dynamic_cast<const elecstate::ElecStateLCAO<T>*>(pelec)->get_DM();
         if(GlobalV::NSPIN==2)
         {
-            const_cast<elecstate::DensityMatrix<T, double>*>(dm_p)->switch_dmr(2);
+            const_cast<elecstate::DensityMatrix<T, double>*>(dm_p)->switch_dmr(1);
         }
         const hamilt::HContainer<double>* dmr = dm_p->get_DMR_pointer(1);
         tmp_nonlocal.cal_force_stress(isforce, isstress, dmr, fvnl_dbeta, svnl_dbeta);