A weird thing in judging the position of the particles

edited November 2016 in DualSPHysics_v4.0
Hi all,

I have tried to add some wave dissipation code in Dualsphysics.
The first job is to judge whether the particle is in the end part of the flume.
To see if the judgement of position is right, I forced the accelerates of the particles in aYWDP1&&float(Lposxy.x)0.)Lace.x=float(double(Lace.x)*double(fdamp)*0.);
if(Lace.y*Lvelrhop.y>0.)Lace.y=float(double(Lace.y)*double(fdamp)*0.);
if(Lace.z*Lvelrhop.z>0.)Lace.z=float(double(Lace.z)*double(fdamp)*0.);
ace[p] = Lace;
}
}


Sign In or Register to comment.