File tree Expand file tree Collapse file tree 1 file changed +5
-5
lines changed Expand file tree Collapse file tree 1 file changed +5
-5
lines changed Original file line number Diff line number Diff line change @@ -92,7 +92,7 @@ void MSST::first_half()
9292 vsum = vel_sum ();
9393
9494 // save the velocities
95- for (int i; i<ucell.nat ; ++i)
95+ for (int i= 0 ; i<ucell.nat ; ++i)
9696 {
9797 old_v[i] = vel[i];
9898 }
@@ -103,7 +103,7 @@ void MSST::first_half()
103103 vsum = vel_sum ();
104104
105105 // reset the velocities
106- for (int i; i<ucell.nat ; ++i)
106+ for (int i= 0 ; i<ucell.nat ; ++i)
107107 {
108108 vel[i] = old_v[i];
109109 }
@@ -274,7 +274,7 @@ void MSST::propagate_vel()
274274 }
275275 else
276276 {
277- vel[i][k] += ( const_C[k] + const_D[k] * vel[i][k] ) * dthalf +
277+ vel[i][k] += ( const_C[k] + const_D[k] * vel[i][k] ) * dthalf +
278278 0.5 * (const_D[k] * const_D[k] * vel[i][k] + const_C[k] * const_D[k] ) * dthalf * dthalf;
279279 }
280280 }
@@ -304,7 +304,7 @@ void MSST::propagate_voldot()
304304 }
305305 else
306306 {
307- omega[sd] += (const_A - const_B * omega[sd]) * dthalf +
307+ omega[sd] += (const_A - const_B * omega[sd]) * dthalf +
308308 0.5 * (const_B * const_B * omega[sd] - const_A * const_B) * dthalf * dthalf;
309309 }
310- }
310+ }
You can’t perform that action at this time.
0 commit comments