michellab / plumed2

Development version of plumed 2
GNU Lesser General Public License v3.0
1 stars 0 forks source link

jedi.cpp #1

Closed carajillu closed 6 years ago

carajillu commented 6 years ago

I think there is a problem between the way atoms are read and the makeWhole() call just before updating the grid. The plumed developers manual says that makeWhole() makes atoms whole (which I am not sure what it means) "assuming they are in the proper order", but the vector that contains the atom numbers to request is built first with the polar apolar atoms, then with the polar atoms (lines 457-467), which means that they are not in order, so that might be altering the way the atoms are "made whole" and that could be affecting the score. Mentioned lines are in the devel branch

jmichel80 commented 6 years ago

Hi @carajillu there may well be a problem that I missed in my tests. Perhaps the issue becomes apparent when a molecule drifts out of the minimum image.

Better to assign the issue to yourself, you are more familiar than me with the code by now and I doubt i can go deep in the code before a while.

Sent from my iPhone

On 20 Jul 2018, at 15:58, carajillu notifications@github.com<mailto:notifications@github.com> wrote:

I think there is a problem between the way atoms are read and the makeWhole() call just before updating the grid. The plumed developers manual says that makeWhole() makes atoms whole (which I am not sure what it means) "assuming they are in the proper order", but the vector that contains the atom numbers to request is built first with the polar apolar atoms, then with the polar atoms (lines 457-467), which means that they are not in order, so that might be altering the way the atoms are "made whole" and that could be affecting the score. Mentioned lines are in the devel branch

— You are receiving this because you were assigned. Reply to this email directly, view it on GitHubhttps://github.com/michellab/plumed2/issues/1, or mute the threadhttps://github.com/notifications/unsubscribe-auth/ALLd5E6WekUZNYAIY99K__yGXIczJf-Vks5uIfAKgaJpZM4VYLN4.

The University of Edinburgh is a charitable body, registered in Scotland, with registration number SC005336.

carajillu commented 6 years ago

Now the vector is sorted and everything is fixed accordingly (check last commit in devel) but it doesn't seem to have an effect in the fluctuations of the score when no clustering is done. It doesn't hurt anyway, so I pushed the modification to devel. Closing issue.