Bugs found in the code
Floating bodies with Symplectic algorithm:
during corrector step:
JSphCpu.cpp in line 1088: if(Nfloat)RunFloating(DtPre,false);
CudaSphSU.cu in line 444: if(floating)CsRunFloating(dc,cte,dc->dtpre);
during corrector step:
JSphCpu.cpp in line 1088: if(Nfloat)RunFloating(DtPre,false);
CudaSphSU.cu in line 444: if(floating)CsRunFloating(dc,cte,dc->dtpre);
Comments
If use
it will be wrong in initial velocity