Commit 3276c6e5 by Francois Gygi

Fixed initialization of ekin_ in compute_v0


git-svn-id: http://qboxcode.org/svn/qb/trunk@228 cba15fb0-1239-40c8-b417-11db7ca47a34
parent 5702d7f8
......@@ -3,7 +3,7 @@
// MDIonicStepper.C
//
////////////////////////////////////////////////////////////////////////////////
// $Id: MDIonicStepper.C,v 1.6 2004-05-03 19:14:08 fgygi Exp $
// $Id: MDIonicStepper.C,v 1.7 2004-05-20 00:20:26 fgygi Exp $
#include "MDIonicStepper.h"
......@@ -70,6 +70,7 @@ void MDIonicStepper::compute_v0(const vector<vector< double> >& f0)
// Note: vhalf contains the velocity at t-1/2 since compute_v0 is called
// after a call to update_r()
ekin_ = 0.0;
for ( int is = 0; is < vhalf_.size(); is++ )
{
assert(pmass_[is] > 0.0);
......
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