166 |
|
//integrableObjects_.insert(integrableObjects_.end(), rigidBodies_.begin(), rigidBodies_.end()); |
167 |
|
} |
168 |
|
|
169 |
< |
double Molecule::getMass() { |
169 |
> |
RealType Molecule::getMass() { |
170 |
|
StuntDouble* sd; |
171 |
|
std::vector<StuntDouble*>::iterator i; |
172 |
< |
double mass = 0.0; |
172 |
> |
RealType mass = 0.0; |
173 |
|
|
174 |
|
for (sd = beginIntegrableObject(i); sd != NULL; sd = nextIntegrableObject(i)){ |
175 |
|
mass += sd->getMass(); |
183 |
|
StuntDouble* sd; |
184 |
|
std::vector<StuntDouble*>::iterator i; |
185 |
|
Vector3d com; |
186 |
< |
double totalMass = 0; |
187 |
< |
double mass; |
186 |
> |
RealType totalMass = 0; |
187 |
> |
RealType mass; |
188 |
|
|
189 |
|
for (sd = beginIntegrableObject(i); sd != NULL; sd = nextIntegrableObject(i)){ |
190 |
|
mass = sd->getMass(); |
211 |
|
StuntDouble* sd; |
212 |
|
std::vector<StuntDouble*>::iterator i; |
213 |
|
Vector3d velCom; |
214 |
< |
double totalMass = 0; |
215 |
< |
double mass; |
214 |
> |
RealType totalMass = 0; |
215 |
> |
RealType mass; |
216 |
|
|
217 |
|
for (sd = beginIntegrableObject(i); sd != NULL; sd = nextIntegrableObject(i)){ |
218 |
|
mass = sd->getMass(); |
225 |
|
return velCom; |
226 |
|
} |
227 |
|
|
228 |
< |
double Molecule::getPotential() { |
228 |
> |
RealType Molecule::getPotential() { |
229 |
|
|
230 |
|
Bond* bond; |
231 |
|
Bend* bend; |
234 |
|
Molecule::BendIterator bendIter; |
235 |
|
Molecule::TorsionIterator torsionIter; |
236 |
|
|
237 |
< |
double potential = 0.0; |
237 |
> |
RealType potential = 0.0; |
238 |
|
|
239 |
|
for (bond = beginBond(bondIter); bond != NULL; bond = nextBond(bondIter)) { |
240 |
|
potential += bond->getPotential(); |