<div dir="ltr">Hello,<div><br></div><div>I have noticed a bug where COMBINE_COLVAR incorrectly reports COLVAR values.</div><div><br></div><div>For example, one would expect two COLVARs below to be the same:</div><div><br></div><div><div>  &COLVAR</div><div>   &COORDINATION</div><div>    ATOMS_FROM 1</div><div>    KINDS_TO H</div><div>    NN 9</div><div>    ND 14</div><div>    R0 [angstrom] 1.5</div><div>   &END COORDINATION</div><div>  &END COLVAR</div></div><div><br></div><div>and </div><div><br></div><div><div>  &COLVAR</div><div>   &COMBINE_COLVAR</div><div>    &COLVAR</div><div>     &COORDINATION</div><div>      ATOMS_FROM 1</div><div>      KINDS_TO H</div><div>      NN 9</div><div>      ND 14</div><div>      R0 [angstrom] 1.5</div><div>     &END COORDINATION</div><div>    &END COLVAR</div><div><br></div><div>    FUNCTION CV_combined</div><div>    VARIABLES CV_combined</div><div>   &END COMBINE_COLVAR</div><div>  &END COLVAR</div></div><div><br></div><div>However, metadynLog reports correct CV1, while CV2 (combined CV) is listed as 0.</div><div><br></div><div>I think this bug was reported earlier at https://groups.google.com/forum/#!searchin/cp2k/COMBINE_COLVAR%7Csort:relevance/cp2k/_uFtNhs8ptY/la3vSdiX7I4J but apparently it is not fixed yet.</div><div>I would like to confirm that this bug is (at least) specific to COMBINE_COLVAR, where COORDINATION is used together with ATOMS_FROM and KINDS_TO.</div><div><br></div><div>I have attached an input file that would reproduce this result. It should take less than a minute to run on a single core.</div></div>