<HTML><BODY>Dear David!<br><br>Thank you very much.<br>I have already made changes in gmx_cluster.cpp code, having solved both problems. An array of boxes and an array of frame numbers are formed similar to the array of frames times at extraction of trajectories. To my regret, commit for gromacs repository is troubles for me so far, but I hope to solve it soon.<br><br>Boris.<br><br><br><blockquote style="border-left:1px solid #0857A6; margin:10px; padding:0 0 0 10px;">
Суббота, 4 августа 2018, 16:16 +03:00 от David van der Spoel <spoel@xray.bmc.uu.se>:<br>
<br>
<div id="">
<div class="js-helper js-readmsg-msg">
<style type="text/css"></style>
<div>
<base target="_self" href="https://e.mail.ru/">
<div id="style_15333886180000000130_BODY">Den 2018-07-26 kl. 13:05, skrev Boris Timofeev:<br>
> Dear developers!<br>
> <br>
> In processing of results of REMD there were inconveniences which, I think,<br>
> it will be easy to eliminate in future versions of Gromacs.<br>
> <br>
This is a useful suggestions, however it is not entirely trivial to <br>
choose which box one would use. Please start a redmine!<br>
<br>
<br>
> 1. Command<br>
> gmx cluster ... -cl somefile.pdb ...<br>
> unlike command<br>
> gmx trjconv ... -o somefile.pdb ...<br>
> doesn't form record about the box sizes in the PDB file:<br>
> CRYST1 boxx boxy boxz ...<br>
> <br>
> As a result, subsequent call<br>
> gmx hbond .. -f somefile.pdb ....<br>
> deals to the error:<br>
> <br>
> Program: gmx hbond, version 2016.2<br>
> Source file: src\gromacs\gmxana\gmx_hbond.cpp (line 947)<br>
> Fatal error:<br>
> Your computational box has shrunk too much.<br>
> gmx hbond can not handle this situation, sorry.<br>
> <br>
> Manual addition the record CRYST1 to pdb-file solves this problem.<br>
> I suggest to add formation of the record CRYST1 to output pdb,<br>
> formed by the gmx cluster command.<br>
> <br>
> 2. The gmx trjconv command can select parts of a trajectory by means of <br>
> the index file: - sub index.ndx:<br>
> <br>
> "Use option -sub. This assumes that the entries in the index file are <br>
> frame numbers and<br>
> dumps each group in the index file to a separate trajectory file."<br>
> <br>
> It would be very convenient if the "gmx cluster" command could form,<br>
> optionally, index files for the found clusters.<br>
> <br>
> Now this problem is solved by copying of times from gmx cluster <br>
> command's log-file<br>
> and transformation of times to indexes of frames manually. This is very <br>
> uncomfortable.<br>
> <br>
> Sincerely yours,<br>
> Boris Timofeev<br>
> <br>
> <br>
> <br>
> <br>
> -- <br>
> Boris Timofeev<br>
> <br>
> <br>
<br>
<br>
-- <br>
David van der Spoel, Ph.D., Professor of Biology<br>
Head of Department, Cell & Molecular Biology, Uppsala University.<br>
Box 596, SE-75124 Uppsala, Sweden. Phone: <span class="js-phone-number">+46184714205</span>.<br>
<a href="http://www.icm.uu.se" target="_blank">http://www.icm.uu.se</a><br>
-- <br>
Gromacs Developers mailing list<br>
<br>
* Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/GMX-developers_List" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/GMX-developers_List</a> before posting!<br>
<br>
* Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br>
<br>
* For (un)subscribe requests visit<br>
<a href="https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-developers" target="_blank">https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-developers</a> or send a mail to gmx-developers-request@gromacs.org.<br>
</div>
<base target="_self" href="https://e.mail.ru/">
</div>
</div>
</div>
</blockquote>
<br>
<br>-- <br>Boris Timofeev<br></BODY></HTML>