| 166 | 
  | 
  | 
| 167 | 
  | 
} | 
| 168 | 
  | 
 | 
| 169 | 
< | 
void Molecule::getCOMvel( double mtot, double COMvel[3] ) { | 
| 169 | 
> | 
double Molecule::getCOMvel( double COMvel[3] ) { | 
| 170 | 
  | 
 | 
| 171 | 
< | 
  double mass; | 
| 171 | 
> | 
  double mass, mtot; | 
| 172 | 
  | 
  int i; | 
| 173 | 
  | 
 | 
| 174 | 
  | 
  COMvel[0] = 0.0; | 
| 192 | 
  | 
  COMvel[1] /= mtot; | 
| 193 | 
  | 
  COMvel[2] /= mtot; | 
| 194 | 
  | 
  | 
| 195 | 
+ | 
  return mtot; | 
| 196 | 
+ | 
 | 
| 197 | 
  | 
} | 
| 198 | 
+ | 
 | 
| 199 | 
+ | 
void Molecule::atomicRollCall(int* molMembership) { | 
| 200 | 
+ | 
  int i, which; | 
| 201 | 
+ | 
 | 
| 202 | 
+ | 
  for (i=0; i < nAtoms; i++) { | 
| 203 | 
+ | 
    if (myAtoms[i] != NULL) { | 
| 204 | 
+ | 
      which = myAtoms[i]->getIndex(); | 
| 205 | 
+ | 
      molMembership[which] = myIndex; | 
| 206 | 
+ | 
    } | 
| 207 | 
+ | 
  } | 
| 208 | 
+ | 
} |