A bug in ghostFluidHelper?

[furryrabbit]

Some particles at boundary sometimes get large velocity after solving pressure. Then I found that something weird in function ghostFluidHelper function in pressure.cpp
I think this line

if(alpha < gfClamp) return alpha = gfClamp;

should be changed into

if(alpha < gfClamp)  alpha = gfClamp;

In the case of alpha<gfClamp, the Level set value at fluid cell center ( phi_inner ) is nearly 0, ghostFluidHelper() returns (phi_outer/phi_inner) which would be a big value rather than gfClamp(~1e-4)