是否有任何方法可以将边界周期条件合并到非结构化数据点的核求和中？

Question

假设有N个粒子分布在一个体积为V=L^3的周期性立方体内部，我想将立方体划分为一个规则的网格，并在每个网格点r_i处计算如下求和，

其中 r_j 是质量为 y_j 的第 j 个粒子的位置向量，K(r) 是一个各向同性核函数。

但是，我真的不知道如何将周期性边界条件合并到这些非结构化数据点的核求和中。有什么方法或算法可以解决这个问题吗？