Hello!
I want to use potential that requests constraint of all bond lengths in a relatively large molecule. Fix rigid is not the thing I need, I really need fix shake which can handle more than 4 atoms in a cluster. As I understand implementation of such fix shake implies writing a function similar to shake2, shake3, shake4, shake3angle. I understand that it is not a problem of greatest interest for you, probably I could write this function by myself, and you could use it in improved fix shake. Is it possible?
Thank you for your attention,
Gleb Ivanovsky