Hi all,
I’m trying to implement the molecular dynamics of rounded cubes in a cylinder. At the moment, I use "fix wall/body/polyhedron k_n c_n c_t zcylinder radius " command, and an error jumped out as follows,
ERROR: Illegal fix wall/body/polyhedron command (…/fix_wall_body_polyhedron.cpp:122)
Last command: fix 3 all wall/body/polyhedron 2000 50 50 zcylinder 0.0 10
I looked through the source code and found that there is no code corresponding to “zcylinder” in fix_wall_body_polyhedron file which cause the error. I wonder it is the bug in the source code.
Could anyone give me some suggestions to implement the simulation of cubes in a cylinder through other approaches?
Best wishes,
Mengjie