Commit 8da5db71 by Francois Gygi

added rescale_velocities member


git-svn-id: http://qboxcode.org/svn/qb/trunk@775 cba15fb0-1239-40c8-b417-11db7ca47a34
parent c5b2ecc1
......@@ -15,7 +15,7 @@
// AtomSet.C
//
////////////////////////////////////////////////////////////////////////////////
// $Id: AtomSet.C,v 1.28 2010-02-20 23:13:02 fgygi Exp $
// $Id: AtomSet.C,v 1.29 2010-04-16 21:40:50 fgygi Exp $
#include "AtomSet.h"
#include "Species.h"
......@@ -397,16 +397,25 @@ void AtomSet::reset_velocities(void)
{
for ( int is = 0; is < atom_list.size(); is++ )
{
int i = 0;
for ( int ia = 0; ia < atom_list[is].size(); ia++ )
{
atom_list[is][ia]->set_velocity(D3vector(0.0, 0.0, 0.0));
i += 3;
}
}
}
////////////////////////////////////////////////////////////////////////////////
void AtomSet::rescale_velocities(double fac)
{
vector<vector<double> > v;
get_velocities(v);
for ( int is = 0; is < v.size(); is++ )
{
for ( int ia = 0; ia < v[is].size(); ia++ )
v[is][ia] *= fac;
}
set_velocities(v);
}
////////////////////////////////////////////////////////////////////////////////
D3vector AtomSet::vcm(void) const
{
D3vector mvsum;
......
......@@ -15,7 +15,7 @@
// AtomSet.h
//
////////////////////////////////////////////////////////////////////////////////
// $Id: AtomSet.h,v 1.24 2009-10-06 06:29:05 fgygi Exp $
// $Id: AtomSet.h,v 1.25 2010-04-16 21:40:50 fgygi Exp $
#ifndef ATOMSET_H
#define ATOMSET_H
......@@ -83,6 +83,7 @@ class AtomSet
}
void sync(void);
void reset_velocities(void);
void rescale_velocities(double fac);
D3vector vcm(void) const;
D3vector dipole(void) const;
void reset_vcm(void);
......
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