53 |
|
#define __C |
54 |
|
#include "UseTheForce/DarkSide/fInteractionMap.h" |
55 |
|
#include "utils/simError.h" |
56 |
+ |
#include "primitives/Bond.hpp" |
57 |
|
#include "primitives/Bend.hpp" |
57 |
– |
#include "primitives/Bend.hpp" |
58 |
|
namespace oopse { |
59 |
|
|
60 |
|
void ForceManager::calcForces(bool needPotential, bool needStress) { |
251 |
|
longRangePotential[i]=0.0; //Initialize array |
252 |
|
} |
253 |
|
|
254 |
< |
doForceLoop( pos, |
255 |
< |
rc, |
256 |
< |
A, |
257 |
< |
electroFrame, |
258 |
< |
frc, |
259 |
< |
trq, |
260 |
< |
tau.getArrayPointer(), |
261 |
< |
longRangePotential, |
262 |
< |
&passedCalcPot, |
263 |
< |
&passedCalcStress, |
264 |
< |
&isError ); |
265 |
< |
|
254 |
> |
doForceLoop(pos, |
255 |
> |
rc, |
256 |
> |
A, |
257 |
> |
electroFrame, |
258 |
> |
frc, |
259 |
> |
trq, |
260 |
> |
tau.getArrayPointer(), |
261 |
> |
longRangePotential, |
262 |
> |
&passedCalcPot, |
263 |
> |
&passedCalcStress, |
264 |
> |
&isError ); |
265 |
> |
|
266 |
|
if( isError ){ |
267 |
|
sprintf( painCave.errMsg, |
268 |
|
"Error returned from the fortran force calculation.\n" ); |