Hi, all. I found something which is not a “bug”, but at least a pitfall which could be addressed through a error message or strong note in the documentation. I have developed a bond style called fene/MC which I use together with regular FENE bonds using bond_style hybrid. I get the following issue with restarting simulations from restart files:
If, in my input script, I put the “read_restart” command before the “bond_style hybrid” command, like
read_restart ramph.res.500000
bond_style hybrid fene fene/MC
the run fails in an avalanche of seg faults. If, however, I just reverse these two commands, like
bond_style hybrid fene fene/MC
read_restart ramph.res.500000
the run restarts with no problem. Now, in hindsight, looking at the documentation, I am not surprised by this, but a simple “BOND_STYLE HYBRID MUST BE DECLARED BEFORE ANY READ_RESTART COMMAND” warning (sorry for shouting) could have saved me a lot of time spent on fruitless “debugging”. Now it may be that I am wrong here, and that this problem in fact means there is something wrong with my modified code, but as mentioned above switching the order of the two commands seems to fix it.
Best,
Rob