Archive from user: Teemu Murtola If you use a tool that understands selections, but does not require the -s flag for basic use, it is impossible to provide the selection through an index file (-n) if -s is not provided. This gives an error that the maximum allowed atom index is zero. One example of an affected tool is gmx rdf. This works in release-2018, but not in release-2019 or master. This seems to have got broken in I2f43e62bc2d97f5e654f15c6e474b9b71d7106ec, which changed the semantics of TopologyInformation. All comments in the class still say that it returns nullptr as the topology when no topology is loaded, and the selection code relies on that, but that no longer matches the implementation. Originally the class was a simple proxy to accessing different topology-related information in the tools, but it got a lot of other responsibilities as well with the mentioned change. The desired behavior of the class should be decided on (and/or the different responsibilities split), the documentation fixed to match the behavior, and all uses gone through to check that they actually use it as the decided behavior requires. It is likely that currently different parts of the code rely on conflicting behavior… Not sure if there are also other corner cases that may be broken. *(from redmine: issue id 2847, created on 2019-01-29 by gmxdefault, closed on 2019-02-14)* * Relations: * relates #1862 * Changesets: * Revision 70d0bfa521fb8ece2abf8c413d0f0163469d7a54 by Mark Abraham on 2019-02-14T10:05:34Z: ``` Revert changes to implementation of insert-molecules and solvate This is a partial reversion of 5181f8d1, which made an inappropriate dependency from preparation tools to TopologyInformation, which is intended for analysis tools. This is needed in order to fix some actual analysis tools that were broken by related changes to TopologyInformation. Unfortunately in 5181f8d1 readConformation was replaced with TopologyInformation, and the former is no longer available. Instead, readConfAndTopology is used in the same way that it was originally used to implement readConformation. So the behaviour is reverted even though the code looks different. Refs #2847, #1862 Change-Id: I6e4900eea38b5c2e1b85b0d48427dabe26b09f9d ``` * Revision 413479392e5ca7157d670c699ef9e30d23e5d3ca by Mark Abraham on 2019-02-14T15:26:09Z: ``` Fix use of selections referring to index groups Selections that refer only to groups in an index file do not require a topology file. This (and other behaviours) was inadvertently broken by the wrong implementation of the default constructor, because that constructor is used in setting up the analysis environment. The default constructor of TopologyInformation should establish the intended invariant, namely that until a topology is loaded, topology getters should return nullptr. Fixed that constructor, and fixed the unit test that verifies that behaviour. Since that test was added at the same time as a change to the implementation of TopologyInformation in commit 264458c9, the broken invariant wasn't noticed. Added a test also to gmx rdf that verifies that the selection behaviour is fixed. Fixes #2847 Change-Id: I2149616e0248c6d5cdd5657a8292d2c3346614b9 ```