* \param[in] testData Test data structure.
* \param[in] pbc Periodic boundary data.
*/
-void applyLincsCuda(ConstraintsTestData *testData, t_pbc pbc)
+void applyLincsCuda(ConstraintsTestData* testData, t_pbc pbc)
{
- auto lincsCuda = std::make_unique<LincsCuda>(testData->ir_.nLincsIter,
- testData->ir_.nProjOrder,
- nullptr);
+ auto lincsCuda =
+ std::make_unique<LincsCuda>(testData->ir_.nLincsIter, testData->ir_.nProjOrder, nullptr);
bool updateVelocities = true;
int numAtoms = testData->numAtoms_;
lincsCuda->set(testData->idef_, testData->md_);
lincsCuda->setPbc(&pbc);
- allocateDeviceBuffer(&d_x, numAtoms, nullptr);
+ allocateDeviceBuffer(&d_x, numAtoms, nullptr);
allocateDeviceBuffer(&d_xp, numAtoms, nullptr);
- allocateDeviceBuffer(&d_v, numAtoms, nullptr);
+ allocateDeviceBuffer(&d_v, numAtoms, nullptr);
- copyToDeviceBuffer(&d_x, (float3*)(testData->x_.data()), 0, numAtoms, nullptr, GpuApiCallBehavior::Sync, nullptr);
- copyToDeviceBuffer(&d_xp, (float3*)(testData->xPrime_.data()), 0, numAtoms, nullptr, GpuApiCallBehavior::Sync, nullptr);
+ copyToDeviceBuffer(&d_x, (float3*)(testData->x_.data()), 0, numAtoms, nullptr,
+ GpuApiCallBehavior::Sync, nullptr);
+ copyToDeviceBuffer(&d_xp, (float3*)(testData->xPrime_.data()), 0, numAtoms, nullptr,
+ GpuApiCallBehavior::Sync, nullptr);
if (updateVelocities)
{
- copyToDeviceBuffer(&d_v, (float3*)(testData->v_.data()), 0, numAtoms, nullptr, GpuApiCallBehavior::Sync, nullptr);
+ copyToDeviceBuffer(&d_v, (float3*)(testData->v_.data()), 0, numAtoms, nullptr,
+ GpuApiCallBehavior::Sync, nullptr);
}
- lincsCuda->apply(d_x, d_xp,
- updateVelocities, d_v, testData->invdt_,
- testData->computeVirial_, testData->virialScaled_);
+ lincsCuda->apply(d_x, d_xp, updateVelocities, d_v, testData->invdt_, testData->computeVirial_,
+ testData->virialScaled_);
- copyFromDeviceBuffer((float3*)(testData->xPrime_.data()), &d_xp, 0, numAtoms, nullptr, GpuApiCallBehavior::Sync, nullptr);
+ copyFromDeviceBuffer((float3*)(testData->xPrime_.data()), &d_xp, 0, numAtoms, nullptr,
+ GpuApiCallBehavior::Sync, nullptr);
if (updateVelocities)
{
- copyFromDeviceBuffer((float3*)(testData->v_.data()), &d_v, 0, numAtoms, nullptr, GpuApiCallBehavior::Sync, nullptr);
+ copyFromDeviceBuffer((float3*)(testData->v_.data()), &d_v, 0, numAtoms, nullptr,
+ GpuApiCallBehavior::Sync, nullptr);
}
freeDeviceBuffer(&d_x);