Hi all!
I try to set the velocity and the angular momentum to zero at a certain point in my simulation and I want to use the velocity command.
For the velocity it works fine, but for the angular momentum I couldn't get it working.
==============================================================================================================
Error Message
==============================================================================================================
*** Process received signal ***
Signal: Segmentation fault (11)
Signal code: Address not mapped (1)
Failing at address: (nil)
[ 0] /lib/x86_64-linux-gnu/libpthread.so.0(+0x10340) [0x7f15eba89340]
[ 1] /home/ubuntu/LIGGGHTS/LIGGGHTS-PUBLIC/src/lmp_fedora_vtk_fpic() [0xa9ee05]
[ 2] /home/ubuntu/LIGGGHTS/LIGGGHTS-PUBLIC/src/lmp_fedora_vtk_fpic() [0x964598]
[ 3] /home/ubuntu/LIGGGHTS/LIGGGHTS-PUBLIC/src/lmp_fedora_vtk_fpic() [0x9623b8]
[ 4] /home/ubuntu/LIGGGHTS/LIGGGHTS-PUBLIC/src/lmp_fedora_vtk_fpic() [0x962e68]
[ 5] /home/ubuntu/LIGGGHTS/LIGGGHTS-PUBLIC/src/lmp_fedora_vtk_fpic() [0x4187c8]
[ 6] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf5) [0x7f15eb6d5ec5]
[ 7] /home/ubuntu/LIGGGHTS/LIGGGHTS-PUBLIC/src/lmp_fedora_vtk_fpic() [0x4196f7]
*** End of error message ***
==============================================================================================================
==============================================================================================================
The commands I use:
velocity all set 0. 0. 0. sum no units box --> works
velocity all setAngular 0. 0. 0. sum no units box --> does not work
Any ideas or experiences?
Best regards
Harald
Daniel Queteschiner | Thu, 01/28/2016 - 10:07
set command
As a work around you may try to use the set command instead, i.e. something like
set group all angmom 0 0 0aaigner | Fri, 02/12/2016 - 14:21
setAngular
Hi Harald!
velocity setAngular is designed for another atom_style (e.g. ellipsoid) and the corresponding error check was missing.
I added the error check in the dev-version and I think Daniel's workaround should do the job for you.
Before the next release, I will have a closer look.
Bests
Andreas