@@ -79,9 +79,7 @@ void Force_Stress_LCAO<T>::getForceStress(UnitCell& ucell,
7979 ModuleBase::matrix fewalds;
8080 ModuleBase::matrix fcc;
8181 ModuleBase::matrix fscc;
82- #ifdef __MLALGO
8382 ModuleBase::matrix fvnl_dalpha; // deepks
84- #endif
8583
8684 fvl_dphi.create (nat, 3 ); // must do it now, update it later, noted by zhengdy
8785
@@ -95,9 +93,7 @@ void Force_Stress_LCAO<T>::getForceStress(UnitCell& ucell,
9593 fewalds.create (nat, 3 );
9694 fcc.create (nat, 3 );
9795 fscc.create (nat, 3 );
98- #ifdef __MLALGO
9996 fvnl_dalpha.create (nat, 3 ); // deepks
100- #endif
10197
10298 // calculate basic terms in Force, same method with PW base
10399 this ->calForcePwPart (ucell,
@@ -124,9 +120,7 @@ void Force_Stress_LCAO<T>::getForceStress(UnitCell& ucell,
124120 ModuleBase::matrix stvnl_dphi;
125121 ModuleBase::matrix svnl_dbeta;
126122 ModuleBase::matrix svl_dphi;
127- #ifdef __MLALGO
128123 ModuleBase::matrix svnl_dalpha; // deepks
129- #endif
130124
131125 // ! stress
132126 if (isstress)
@@ -142,9 +136,8 @@ void Force_Stress_LCAO<T>::getForceStress(UnitCell& ucell,
142136 stvnl_dphi.create (3 , 3 );
143137 svnl_dbeta.create (3 , 3 );
144138 svl_dphi.create (3 , 3 );
145- #ifdef __MLALGO
146139 svnl_dalpha.create (3 , 3 );
147- # endif
140+
148141 // calculate basic terms in Stress, similar method with PW base
149142 this ->calStressPwPart (ucell,
150143 sigmadvl,
0 commit comments