Hi,
I am writing a code for detecting particle aggloemrates based on contacts. I read all the contact pairs in particle contact force main function, put those pairs in a vector and wrote a function to process them into aggloemrate information.
However, I find it difficult to find a place for this process function. I want this process function to run only once at the end of every timestep, or run only once after calculating contact force in a timestep cycle. Currently I put this process function in particle factory main function. However it always runs at the next timestep because particle factory runs before contact force. This will cause some problems since I also want to remove some particles.
Is there a way for solving that problem?
Regards,
Jason