Commit d72fe507 by Francois Gygi

Fixed bug in use of mdionic_stepper pointer when atoms_dyn==LOCKED


git-svn-id: http://qboxcode.org/svn/qb/trunk@1632 cba15fb0-1239-40c8-b417-11db7ca47a34
parent d4776225
......@@ -244,8 +244,9 @@ void CPSampleStepper::step(int niter)
cout << " <temp_ion> " << mdionic_stepper->temp() << " </temp_ion>\n";
cout << " <eta_ion> " << mdionic_stepper->eta() << " </eta_ion>\n";
}
double econst = energy + ekin_ion + ekin_e +
mdionic_stepper->ekin_stepper();
double econst = energy + ekin_ion + ekin_e;
if ( mdionic_stepper )
econst += mdionic_stepper->ekin_stepper();
cout << " <econst> " << econst << " </econst>\n";
cout << " <ekin_ec> " << econst + ekin_e << " </ekin_ec>\n";
}
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment