PartCfg.cpp

Go to the documentation of this file.

/*
 * Copyright (C) 2010-2022 The ESPResSo project
 *
 * This file is part of ESPResSo.
 *
 * ESPResSo is free software: you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation, either version 3 of the License, or
 * (at your option) any later version.
 *
 * ESPResSo is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with this program.  If not, see <http://www.gnu.org/licenses/>.
 */
 
#include "PartCfg.hpp"
 
#include "BoxGeometry.hpp"
#include "particle_node.hpp"
#include "system/System.hpp"
 
#include <algorithm>
#include <cstddef>
#include <span>
 
void PartCfg::update() {
  m_parts.clear();
 
  auto const ids = get_particle_ids();
  auto const chunk_size = fetch_cache_max_size();
 
  for (std::size_t offset = 0; offset < ids.size();) {
    auto const this_size = std::clamp(chunk_size, std::size_t{0},
                                      std::size_t{ids.size() - offset});
    auto const chunk_ids = std::span(ids.data() + offset, this_size);
 
    prefetch_particle_data(chunk_ids);
 
    for (auto id : chunk_ids) {
      m_parts.push_back(get_particle_data(id));
 
      auto &p = m_parts.back();
      p.pos() += m_box_geo.image_shift(p.image_box());
      p.image_box() = {};
    }
 
    offset += this_size;
  }
}