LG Muliphase possible code error

edited April 2019 in DualSPHysics v4.2
It is noted that in Liquid/Gas multiphase version, the "GetKernelCubicTensil" function, there are below codes:
//££££££££££££££££££££££££££££££££££££££££££££££££££££££££££££££££££££££?/MP
const float tensilp1 = (pressp1 / (rhopp1*rhopp2))*(pressp1 > 0 ? 0.01f : -0.2f);
const float tensilp2 = (pressp2 / (rhopp1*rhopp2))*(pressp2 > 0 ? 0.01f : -0.2f);
//££££££££££££££££££££££££££££££££££££££££££££££££££££££££££££££££££££££?/MP

While in ver 4.2 single phase version, the code is:
const float tensilp1 = (pressp1 / (rhopp1*rhopp1))*(pressp1 > 0 ? 0.01f : -0.2f);
const float tensilp2 = (pressp2 / (rhopp2*rhopp2))*(pressp2 > 0 ? 0.01f : -0.2f);.

the difference is regarding rhopp1*rhopp1 or rhopp1*rhopp2. I can not find any reason for this change in the corresponding paper, so was this made intentionally or is it a code error?
Sign In or Register to comment.