Test: add coverage of force/stress for SOC, and delete useless code

Author: dyzheng

source/Makefile.Objects

@@ -571,8 +571,6 @@ OBJS_LCAO=evolve_elec.o\
       stress_tools.o\
       edm.o\
       pulay_force_stress_center2.o\
-	  fvnl_dbeta_gamma.o\
-	  fvnl_dbeta_k.o\
 	  grid_init.o\
       spar_dh.o\
       spar_exx.o\

source/module_hamilt_lcao/hamilt_lcaodft/CMakeLists.txt

@@ -21,8 +21,6 @@ if(ENABLE_LCAO)
         FORCE_k.cpp
         stress_tools.cpp
         edm.cpp
-        fvnl_dbeta_gamma.cpp
-        fvnl_dbeta_k.cpp
         grid_init.cpp
         spar_dh.cpp
         spar_exx.cpp

source/module_hamilt_lcao/hamilt_lcaodft/FORCE.h

@@ -116,17 +116,6 @@ class Force_LCAO
                         ModuleBase::matrix& stvnl_dphi,
                         Record_adj* ra = nullptr);
 
-    void cal_fvnl_dbeta(const elecstate::DensityMatrix<T, double>* dm,
-                        const Parallel_Orbitals& pv,
-                        const UnitCell& ucell,
-                        const LCAO_Orbitals& orb,
-                        const TwoCenterIntegrator& intor_orb_beta,
-                        Grid_Driver& gd,
-                        const bool isforce,
-                        const bool isstress,
-                        ModuleBase::matrix& fvnl_dbeta,
-                        ModuleBase::matrix& svnl_dbeta);
-
     //-------------------------------------------
     // forces related to local pseudopotentials
     //-------------------------------------------

source/module_hamilt_lcao/hamilt_lcaodft/FORCE_gamma.cpp

@@ -219,17 +219,6 @@ void Force_LCAO<double>::ftable(const bool isforce,
     //tvnl_dphi
     PulayForceStress::cal_pulay_fs(ftvnl_dphi, stvnl_dphi, *dm, ucell, pv, dHx, dHxy, isforce, isstress);
 
-    /*this->cal_fvnl_dbeta(dm,
-                         pv,
-                         ucell,
-                         orb,
-                         *(two_center_bundle.overlap_orb_beta),
-                         GlobalC::GridD,
-                         isforce,
-                         isstress,
-                         fvnl_dbeta,
-                         svnl_dbeta);*/
-
     // vl_dphi
     PulayForceStress::cal_pulay_fs(fvl_dphi, svl_dphi, *dm, ucell, pelec->pot, gint, isforce, isstress, false/*reset dm to gint*/);
 

source/module_hamilt_lcao/hamilt_lcaodft/FORCE_k.cpp

@@ -319,17 +319,6 @@ void Force_LCAO<std::complex<double>>::ftable(const bool isforce,
     // vl_dphi
     PulayForceStress::cal_pulay_fs(fvl_dphi, svl_dphi, *dm, ucell, pelec->pot, gint, isforce, isstress, false/*reset dm to gint*/);
 
-    /*this->cal_fvnl_dbeta(dm,
-                         pv,
-                         ucell,
-                         orb,
-                         *(two_center_bundle.overlap_orb_beta),
-                         GlobalC::GridD,
-                         isforce,
-                         isstress,
-                         fvnl_dbeta,
-                         svnl_dbeta);*/
-
 #ifdef __DEEPKS
     if (PARAM.inp.deepks_scf)
     {