Open YiYangddd opened 2 years ago
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.
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.