ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/group/trunk/OOPSE/libmdtools/Molecule.cpp
(Generate patch)

Comparing trunk/OOPSE/libmdtools/Molecule.cpp (file contents):
Revision 450 by mmeineke, Thu Apr 3 21:12:51 2003 UTC vs.
Revision 468 by gezelter, Mon Apr 7 16:56:38 2003 UTC

# Line 165 | Line 165 | void Molecule::getCOM( double COM[3] ) {
165    COM[2] /= mtot;
166  
167   }
168 +
169 + void Molecule::getCOMvel( double mtot, double COMvel[3] ) {
170 +
171 +  double mass;
172 +  int i;
173 +
174 +  COMvel[0] = 0.0;
175 +  COMvel[1] = 0.0;
176 +  COMvel[2] = 0.0;
177 +  mtot   = 0.0;
178 +
179 +  for (i=0; i < nAtoms; i++) {
180 +    if (myAtoms[i] != NULL) {
181 +      
182 +      mass = myAtoms[i]->getMass();
183 +      mtot   += mass;
184 +      COMvel[0] += myAtoms[i]->get_vx() * mass;
185 +      COMvel[1] += myAtoms[i]->get_vy() * mass;
186 +      COMvel[2] += myAtoms[i]->get_vz() * mass;
187 +
188 +    }
189 +  }
190 +
191 +  COMvel[0] /= mtot;
192 +  COMvel[1] /= mtot;
193 +  COMvel[2] /= mtot;
194 +
195 + }

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines