adevaucorbeil/ArealGlide

Simulation time takes too long

Opened this issue · 2 comments

squarehor100.xy0.txt
The simulation seem will stuck when I use this input file, where I design 100 rod shaped particles squarely arranged. if you could take a look? Thanks.

Yi,

What is happening here is that you have points that are on both the right hand side boundary (i.e. x=0) that coincide with points that are on the left hand side boundary (x=1) . Due to the periodic boundary conditions, this creates duplicates and lead to this kinds of problem. Here are two examples: points 2766 and 1997. Note that differences smaller than 1E-12 can be considered as 0.
2766 0.999999999999999 0.331111111111111
1997 3.04E-18 0.331111111111111

Either you make sure that does not happen, or we have to add a function in the code to clean the input file. The later would be efficient, but might create some unwanted results.

Thanks Alban,
I'll try to modify the input with a different algorithm, this should be able to address the issue.