ForceCalculator forceCalculator(simulationState, options);
// The listed force calculator is also initialized with the required arguments
- ListedForceCalculator listedForceCalculator(topology.getInteractionData(),
- topology.numParticles(), 4, box);
+ ListedForceCalculator listedForceCalculator(
+ topology.getInteractionData(), topology.numParticles(), 4, box);
// Integrator is initialized with an array of inverse masses (constructed from topology) and
// the bounding box
LeapFrog integrator(topology, box);
// Print some diagnostic info
- printf("initial position of particle 0: x %4f y %4f z %4f\n", simulationState.coordinates()[0][0],
- simulationState.coordinates()[0][1], simulationState.coordinates()[0][2]);
+ printf("initial position of particle 0: x %4f y %4f z %4f\n",
+ simulationState.coordinates()[0][0],
+ simulationState.coordinates()[0][1],
+ simulationState.coordinates()[0][2]);
// MD Loop
int numSteps = 2;
listedForceCalculator.compute(simulationState.coordinates(), simulationState.forces());
// Integrate with a time step of 1 fs, positions, velocities and forces
- integrator.integrate(1.0, simulationState.coordinates(), simulationState.velocities(),
- simulationState.forces());
+ integrator.integrate(
+ 1.0, simulationState.coordinates(), simulationState.velocities(), simulationState.forces());
}
- printf(" final position of particle 9: x %4f y %4f z %4f\n", simulationState.coordinates()[9][0],
- simulationState.coordinates()[9][1], simulationState.coordinates()[9][2]);
+ printf(" final position of particle 9: x %4f y %4f z %4f\n",
+ simulationState.coordinates()[9][0],
+ simulationState.coordinates()[9][1],
+ simulationState.coordinates()[9][2]);
return 0;
} // main