求助：IB-LBM中的插值问题

xiangchuang

哪位前辈指点下，在IB-LBM中，对受力和速度使用两点插值，可是在一段程序中：
// Identify the lowest fluid lattice node in interpolation range.
                // 'Lowest' means: its x- and y-values are the smallest.
                // The other fluid nodes in range have coordinates
                // (x_int + 1, y_int), (x_int, y_int + 1), and (x_int + 1, y_int + 1).
                int x_int = (int)(particle.node[n].x - 0.5 + Nx) - Nx;
                int y_int = (int)(particle.node[n].y + 0.5);

                // Run over all neighboring fluid nodes.
                // In the case of the two-point interpolation, it is 2x2 fluid nodes.
                for (int X = x_int; X <= x_int + 1; ++X) {
                        for (int Y = y_int; Y <= y_int + 1; ++Y) {
                                // Compute distance between object node and fluid lattice node.
                                const double dist_x = particle.node[n].x - 0.5 - X;
                                const double dist_y = particle.node[n].y + 0.5 - Y;
为什么X方向要-0.5，Y方向要+0.5，不是直接向下取整就行了么？

xiangchuang

particle.node[n]是流体中粒子点，Nx，Ny是格子坐标

xiangchuang

particle.node[n]是流体中粒子点，Nx，Ny是格子坐标

hustzy

因为上下边界用的半反弹，左右边界用的周期