The corresponding header file is p3m.hpp. More...

#include "config/config.hpp"
#include "electrostatics/p3m.hpp"
#include "electrostatics/coulomb.hpp"
#include "electrostatics/elc.hpp"
#include "electrostatics/p3m_gpu_cuda.cuh"
#include "electrostatics/p3m_gpu_error.hpp"
#include "electrostatics/p3m.impl.hpp"
#include "p3m/P3MFFT.hpp"
#include "p3m/TuningAlgorithm.hpp"
#include "p3m/TuningLogger.hpp"
#include "p3m/field_layout_helpers.hpp"
#include "p3m/for_each_3d.hpp"
#include "p3m/influence_function.hpp"
#include "p3m/math.hpp"
#include "BoxGeometry.hpp"
#include "LocalBox.hpp"
#include "Particle.hpp"
#include "ParticlePropertyIterator.hpp"
#include "ParticleRange.hpp"
#include "PropagationMode.hpp"
#include "actor/visitors.hpp"
#include "cell_system/CellStructure.hpp"
#include "cell_system/CellStructureType.hpp"
#include "communication.hpp"
#include "errorhandling.hpp"
#include "integrators/Propagation.hpp"
#include "npt.hpp"
#include "p3m/send_mesh.hpp"
#include "system/GpuParticleData.hpp"
#include "system/System.hpp"
#include "tuning.hpp"
#include <utils/Vector.hpp>
#include <utils/integral_parameter.hpp>
#include <utils/math/int_pow.hpp>
#include <utils/math/sqr.hpp>
#include <utils/serialization/array.hpp>
#include <boost/mpi/collectives/all_reduce.hpp>
#include <boost/mpi/collectives/broadcast.hpp>
#include <boost/mpi/collectives/reduce.hpp>
#include <boost/mpi/communicator.hpp>
#include <boost/range/combine.hpp>
#include <boost/range/numeric.hpp>
#include <algorithm>
#include <array>
#include <cassert>
#include <complex>
#include <cstddef>
#include <functional>
#include <initializer_list>
#include <numbers>
#include <optional>
#include <span>
#include <sstream>
#include <stdexcept>
#include <string>
#include <tuple>
#include <utility>
#include <vector>

Include dependency graph for p3m.cpp:

Go to the source code of this file.

Classes
struct	anonymous_namespace{p3m.cpp}::AssignCharge< cao >

struct	AssignForces< cao >

class	CoulombTuningAlgorithm< FloatType, Architecture >

Namespaces
namespace	anonymous_namespace{p3m.cpp}

Functions
template<typename FloatType >
std::complex< FloatType >	multiply_complex_by_imaginary (std::complex< FloatType > const &z, FloatType k)

template<typename FloatType >
std::complex< FloatType >	multiply_complex_by_real (std::complex< FloatType > const &z, FloatType k)

template<typename FloatType >
FloatType	complex_norm2 (std::complex< FloatType > const &z)

template<Utils::MemoryOrder order_in, Utils::MemoryOrder order_out, typename T >
auto	transpose (std::span< T > const &flat_array, Utils::Vector3i const &shape)

static bool	is_node_grid_compatible_with_mesh (Utils::Vector3i const &node_grid, Utils::Vector3i const &mesh)

static auto	get_size_from_shape (Utils::Vector3i const &shape)

static auto	p3m_tune_aliasing_sums (Utils::Vector3i const &shift, Utils::Vector3i const &mesh, Utils::Vector3d const &mesh_i, int cao, double alpha_L_i)
	Aliasing sum used by p3m_k_space_error.

static double	p3m_real_space_error (double pref, double r_cut_iL, int n_c_part, double sum_q2, double alpha_L, Utils::Vector3d const &box_l)
	Calculate the real space contribution to the rms error in the force (as described by Kolafa and Perram).

static double	p3m_k_space_error (double pref, Utils::Vector3i const &mesh, int cao, int n_c_part, double sum_q2, double alpha_L, Utils::Vector3d const &box_l)
	Calculate the analytic expression of the error estimate for the P3M method in [23] (eq.

template<typename combined_ranges >
static auto	calc_dipole_moment (boost::mpi::communicator const &comm, combined_ranges const &p_q_unfolded_pos_range)

Detailed Description

The corresponding header file is p3m.hpp.

Definition in file p3m.cpp.

Function Documentation

◆ calc_dipole_moment()

template<typename combined_ranges >

static auto calc_dipole_moment	(	boost::mpi::communicator const &	comm,
		combined_ranges const &	p_q_unfolded_pos_range
	)

static

Definition at line 435 of file p3m.cpp.

Referenced by CoulombP3MImpl< FloatType, Architecture >::long_range_kernel().

◆ complex_norm2()

template<typename FloatType >

FloatType complex_norm2 ( std::complex< FloatType > const & z )

Definition at line 113 of file p3m.cpp.

References Utils::sqr().

Referenced by CoulombP3MImpl< FloatType, Architecture >::long_range_kernel(), and CoulombP3MImpl< FloatType, Architecture >::long_range_pressure().

◆ get_size_from_shape()

static auto get_size_from_shape ( Utils::Vector3i const & shape )

static

Definition at line 136 of file p3m.cpp.

References Utils::Array< T, N >::cbegin(), and Utils::Array< T, N >::cend().

Referenced by CoulombP3MImpl< FloatType, Architecture >::init_cpu_kernels().

◆ is_node_grid_compatible_with_mesh()

static bool is_node_grid_compatible_with_mesh	(	Utils::Vector3i const &	node_grid,
		Utils::Vector3i const &	mesh
	)

static

Definition at line 130 of file p3m.cpp.

◆ multiply_complex_by_imaginary()

template<typename FloatType >

std::complex< FloatType > multiply_complex_by_imaginary	(	std::complex< FloatType > const &	z,
		FloatType	k
	)

Definition at line 100 of file p3m.cpp.

◆ multiply_complex_by_real()

template<typename FloatType >

std::complex< FloatType > multiply_complex_by_real	(	std::complex< FloatType > const &	z,
		FloatType	k
	)

Definition at line 107 of file p3m.cpp.

◆ p3m_k_space_error()

static double p3m_k_space_error	(	double	pref,
		Utils::Vector3i const &	mesh,
		int	cao,
		int	n_c_part,
		double	sum_q2,
		double	alpha_L,
		Utils::Vector3d const &	box_l
	)

static

Calculate the analytic expression of the error estimate for the P3M method in [23] (eq.

8-23 p. 275) in order to obtain the rms error in the force for a system of N randomly distributed particles in a cubic box (k-space part).

Parameters

pref	Prefactor of Coulomb interaction.
mesh	number of mesh points in one direction.
cao	charge assignment order.
n_c_part	number of charged particles in the system.
sum_q2	sum of square of charges in the system
alpha_L	rescaled Ewald splitting parameter.
box_l	box dimensions.

Returns: reciprocal (k) space error

Definition at line 255 of file p3m.cpp.

References for_each_3d(), math::get_analytic_cotangent_sum_kernel(), p3m_tune_aliasing_sums(), Utils::product(), ROUND_ERROR_PREC, and Utils::sqr().

Referenced by CoulombTuningAlgorithm< FloatType, Architecture >::calculate_accuracy().

◆ p3m_real_space_error()

static double p3m_real_space_error	(	double	pref,
		double	r_cut_iL,
		int	n_c_part,
		double	sum_q2,
		double	alpha_L,
		Utils::Vector3d const &	box_l
	)

static

Calculate the real space contribution to the rms error in the force (as described by Kolafa and Perram).

Parameters

pref	Prefactor of Coulomb interaction.
r_cut_iL	rescaled real space cutoff for p3m method.
n_c_part	number of charged particles in the system.
sum_q2	sum of square of charges in the system
alpha_L	rescaled Ewald splitting parameter.
box_l	box dimensions.

Returns: real space error

Definition at line 234 of file p3m.cpp.

References Utils::product(), and Utils::sqr().

Referenced by CoulombTuningAlgorithm< FloatType, Architecture >::calculate_accuracy().

◆ p3m_tune_aliasing_sums()

static auto p3m_tune_aliasing_sums	(	Utils::Vector3i const &	shift,
		Utils::Vector3i const &	mesh,
		Utils::Vector3d const &	mesh_i,
		int	cao,
		double	alpha_L_i
	)

static

Aliasing sum used by p3m_k_space_error.

Definition at line 190 of file p3m.cpp.

References Utils::Vector< T, N >::broadcast(), for_each_3d(), P3M_BRILLOUIN, Utils::product(), math::sinc(), and Utils::sqr().

Referenced by p3m_k_space_error().

◆ transpose()

template<Utils::MemoryOrder order_in, Utils::MemoryOrder order_out, typename T >

auto transpose	(	std::span< T > const &	flat_array,
		Utils::Vector3i const &	shape
	)

Definition at line 118 of file p3m.cpp.

References Utils::Vector< T, N >::broadcast(), for_each_3d(), and Utils::get_linear_index().

Classes

Namespaces

Functions

Detailed Description

Function Documentation

◆ calc_dipole_moment()

◆ complex_norm2()

◆ get_size_from_shape()

◆ is_node_grid_compatible_with_mesh()

◆ multiply_complex_by_imaginary()

◆ multiply_complex_by_real()

◆ p3m_k_space_error()

◆ p3m_real_space_error()

◆ p3m_tune_aliasing_sums()

◆ transpose()