Hi all,
I am trying to do a CG simulation with only 16 particles. I have the initial particle positions and I am creating a hoomd system from a snapshot. I am adding the topology information manually. The bond, angle, dihedral and pair potentials are table potentials which are set from the respective files.
Currently I am facing the following issues,
1. Even though my system only contains 16 particles, I have to set a larger box size. Eg: hoomd.data.boxdim(L=200)
2. After starting my equilibration, dynamic data will only be updated for 3 steps. These will not be updated after that. Sometimes the output will show ** run complete ** but it gives the following error. Is this something related to the potentials I've used? Please help me to figure this out.
** run complete **
**ERROR**: an illegal memory access was encountered before /hoomd/GPUArray.h:160
terminate called after throwing an instance of 'std::runtime_error'
what(): CUDA Error
Aborted (core dumped)
I have attached the script, table potentials and the data here.
Any help to get my simulation started is highly appreciated!
Thanks,
Geemi