Hi,
I have introduced a contact model, and it works well when number of core is 1, the results was in coincide with analytical solutions, but when the number of cores larger than 1, the simulation was not in my expect, like particles had large velocity and disappeared. What may cause the problem as I described?
I am not sure if the initialization of values will lead this. In my code, I have defined some unassigned CSimplevectors in the cpp file and h file, just like:
CSimpervector n1, but not CSimpervector n1(0.0,0.0,0.0).
And I also introduced some new custom properties. Are there any experiences about this?