Skip to content

Commit

Permalink
The rt-bops regtest was calculating the values of this type of CV inc…
Browse files Browse the repository at this point in the history
…orrectly because the actions have to be a particular order in the input file in that old version. The new reference values that I have put here are consistent with the values output by master when the actions are in the correct order. Furthermore, because of the changes to the mechanism via which the forces are applied the order of operations is no longer important
  • Loading branch information
Gareth Aneurin Tribello authored and Gareth Aneurin Tribello committed Sep 6, 2024
1 parent f46ca43 commit cd81a35
Show file tree
Hide file tree
Showing 3 changed files with 13 additions and 3 deletions.
6 changes: 3 additions & 3 deletions regtest/crystdistrib/rt-bops/colvar.reference
Original file line number Diff line number Diff line change
@@ -1,4 +1,4 @@
#! FIELDS time s
0.000000 10.66
0.025000 922.28
0.050000 993.43
0.000000 925.50
0.025000 1019.81
0.050000 1064.21
8 changes: 8 additions & 0 deletions regtest/crystdistrib/rt-bops/config
Original file line number Diff line number Diff line change
Expand Up @@ -2,3 +2,11 @@ type=driver
plumed_modules="adjmat crystdistrib"
arg="--plumed plumed.dat --mf_xtc em.xtc --mc mcfile --timestep=0.001 --trajectory-stride=25"
extra_files="../../trajectories/em.xtc"

function plumed_regtest_before() {
export PLUMED_FORBID_CHAINS=yes
}

function plumed_regtest_after() {
export PLUMED_FORBID_CHAINS=no
}
2 changes: 2 additions & 0 deletions regtest/crystdistrib/rt-bops/plumed.dat
Original file line number Diff line number Diff line change
Expand Up @@ -260,6 +260,8 @@ ops: CUSTOM ARG=qp.i,qp.j,qp.k,c1.w VAR=i,j,k,d FUNC=((2.26628554234*exp(-(32.27

ones: ONES SIZE=252

cc: MATRIX_VECTOR_PRODUCT ARG=qp.i,ones
PRINT ARG=cc FILE=cc FMT=%8.4f

sdot: MATRIX_VECTOR_PRODUCT ARG=ops,ones
s: SUM ARG=sdot PERIODIC=NO
Expand Down

0 comments on commit cd81a35

Please sign in to comment.