ParticleCellHelpers.h

Go to the documentation of this file.

 
#pragma once
#include "autopas/utils/ArrayMath.h"
 
namespace autopas::internal {
template <class CellType>
static bool checkParticleInCellAndUpdateByID(CellType &cell, const typename CellType::ParticleType &particle) {
  for (auto &p : cell) {
    if (p.getID() == particle.getID()) {
      p = particle;
      return true;
    }
  }
  return false;
}
 
template <class CellType>
static bool checkParticleInCellAndUpdateByIDAndPosition(CellType &cell, const typename CellType::ParticleType &particle,
                                                        double absError) {
  using namespace autopas::utils::ArrayMath::literals;
  // This lock is relevant for octree and directSum.
  std::lock_guard<AutoPasLock> cellLock(cell.getCellLock());
  for (auto &p : cell) {
    if (p.getID() == particle.getID()) {
      auto distanceVec = p.getR() - particle.getR();
      auto distanceSqr = autopas::utils::ArrayMath::dot(distanceVec, distanceVec);
      if (distanceSqr < absError * absError) {
        p = particle;
        // found the particle, returning.
        return true;
      }
    }
  }
  return false;
}
}  // namespace autopas::internal