@@ -28,42 +28,12 @@ void LOOP_ions::output_HS_R(
2828
2929 if (GlobalV::NSPIN==1 ||GlobalV::NSPIN==4 )
3030 {
31- // this->UHM.calculate_STN_R();
32- // this->UHM.GK.cal_vlocal_R(0);
33- // this->UHM.GK.distribute_pvpR_tr();
34- // HS_Matrix::save_HSR_tr(0);
35-
3631 // jingan add 2021-6-4, modify 2021-12-2
3732 this ->UHM .calculate_HSR_sparse (0 , sparse_threshold);
3833 }
3934 // /*
4035 else if (GlobalV::NSPIN==2 )
4136 {
42- // this->UHM.calculate_STN_R();
43- // for(int ik=0; ik<GlobalC::kv.nks; ik++)
44- // {
45- // if(ik==0 || ik==GlobalC::kv.nks/2)
46- // {
47- // if(GlobalV::NSPIN==2)GlobalV::CURRENT_SPIN = GlobalC::kv.isk[ik];
48- // for(int ir=0; ir<GlobalC::pw.nrxx; ir++)
49- // {
50- // GlobalC::pot.vr_eff1[ir] = GlobalC::pot.vr_eff( GlobalV::CURRENT_SPIN, ir);
51- // }
52-
53- // if(!GlobalV::GAMMA_ONLY_LOCAL)
54- // {
55- // if(GlobalV::VL_IN_H)
56- // {
57- // //this->UHM.GK.cal_vlocal_k(GlobalC::pot.vrs1,GridT);
58- // this->UHM.GK.cal_vlocal_k(GlobalC::pot.vr_eff1, GlobalC::GridT, GlobalV::CURRENT_SPIN);
59- // }
60- // }
61- // this->UHM.GK.cal_vlocal_R(GlobalV::CURRENT_SPIN);
62- // this->UHM.GK.distribute_pvpR_tr();
63- // HS_Matrix::save_HSR_tr(GlobalV::CURRENT_SPIN);
64- // }
65- // }
66-
6737 // jingan add 2021-6-4
6838 for (int ik = 0 ; ik < GlobalC::kv.nks ; ik++)
6939 {
0 commit comments