dyzheng
diff --git a/‎source/source_lcao/FORCE_STRESS.cpp‎
Lines changed: 2 additions & 0 deletions b/‎source/source_lcao/FORCE_STRESS.cpp‎
Lines changed: 2 additions & 0 deletions
@@ -179,6 +179,7 @@ void Force_Stress_LCAO<T>::getForceStress(UnitCell& ucell,
         if (PARAM.inp.nspin == 2)
         {
             dmat.dm->switch_dmr(1);
+            edm.switch_dmr(1);
         }
 
         const hamilt::HContainer<double>* dmR = dmat.dm->get_DMR_pointer(1);
@@ -209,6 +210,7 @@ void Force_Stress_LCAO<T>::getForceStress(UnitCell& ucell,
         if (PARAM.inp.nspin == 2)
         {
             dmat.dm->switch_dmr(0);
+            edm.switch_dmr(0);
         }
 
         // Calculate local potential force/stress (vl_dphi)
Original file line number	Diff line number	Diff line change
`@@ -179,6 +179,7 @@ void Force_Stress_LCAO<T>::getForceStress(UnitCell& ucell,`
`179`	`179`	`if (PARAM.inp.nspin == 2)`
`180`	`180`	`{`
`181`	`181`	`dmat.dm->switch_dmr(1);`
	`182`	`+ edm.switch_dmr(1);`
`182`	`183`	`}`
`183`	`184`
`184`	`185`	`const hamilt::HContainer<double>* dmR = dmat.dm->get_DMR_pointer(1);`
`@@ -209,6 +210,7 @@ void Force_Stress_LCAO<T>::getForceStress(UnitCell& ucell,`
`209`	`210`	`if (PARAM.inp.nspin == 2)`
`210`	`211`	`{`
`211`	`212`	`dmat.dm->switch_dmr(0);`
	`213`	`+ edm.switch_dmr(0);`
`212`	`214`	`}`
`213`	`215`
`214`	`216`	`// Calculate local potential force/stress (vl_dphi)`