[gmx-users] Gromacs 5 and Shake

[luigi]

Dear Gromacs users,

I wanted to perform a simulation with gromacs v.5, using shake algorithm
to constraint bond lengths.
Looking at the mdp options, shake method is still available, however,
the particle decomposition algorithm has been dropped in gromacs 5, and
shake seems not to be compatible with decomposition of the system in
multiple domains.
I was wondering if there is a way to perform simulations with shake on
multiple threads. Every suggestion (related to the use of shake!) is
welcome.

Thank you in advance,
Luigi