About the boundary handling with MLS in DFSPH

[patricktian1994]

Hi,
Recently I worked on the boundary handling with your great framework, as I see in your paper, the velocity after the divergence-free solver is update with the stiffness parameter, but in this framework it is update with pressure and density like the Band do in the paper.But when it comes to compute the pressure acceleration, no matter with which kind of boundary handling method, it still only considered about the "p_rho2_i" . Is there some advice for this?